曹文超 張運(yùn)龍 嚴(yán)正娟 王敬國(guó) 李俊良 陳 清*

（1青島農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，山東青島266109；2中國(guó)農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，北京100193）

我國(guó)設(shè)施蔬菜的播種面積已達(dá)330萬hm2，占我國(guó)農(nóng)作物播種面積的2.1%（中國(guó)農(nóng)業(yè)統(tǒng)計(jì)年鑒，2010），已呈集約化、規(guī)?；l(fā)展趨勢(shì)。山東是我國(guó)的蔬菜種植大省，其設(shè)施蔬菜種植面積占全國(guó)總面積的1/5以上，總產(chǎn)量居全國(guó)第一位（高中強(qiáng)，2010）。設(shè)施蔬菜種植環(huán)境具有常年高溫、高濕、無降水淋洗及高投入、高產(chǎn)出、高強(qiáng)度利用等特點(diǎn)，是一種高度集約化的利用方式（黃錦法 等，2003）。據(jù)調(diào)查，2004年山東設(shè)施蔬菜氮磷鉀投入量分別為2428、2022 kg·hm-2和2033 kg·hm-2（劉兆輝 等，2008），其中壽光地區(qū)氮磷鉀年投入量高達(dá)4088、3656 kg·hm-2和3438 kg·hm-2（余海英 等，2010），而一般蔬菜形成1000 kg產(chǎn)量平均需要吸收氮磷鉀分別為2～4、0.18～1.20、3.5 kg（杜會(huì)英，2007），其養(yǎng)分投入量遠(yuǎn)遠(yuǎn)超過了蔬菜養(yǎng)分吸收量。此外，目前蔬菜生產(chǎn)中主要采用大水漫灌的灌溉方式，據(jù)調(diào)查，山西省鹽湖區(qū)日光溫室蔬菜每季灌溉10～20次，每次47～63 mm，灌溉總量達(dá)470～1200 mm，平均767 mm（王敬國(guó)，2011）；山東壽光地區(qū)全年灌溉水總量為748～1957 mm，平均灌溉量高達(dá)1307 mm（宋效宗，2007）。“一水一肥，大水漫（溝）灌，肥大水勤，肥隨水走”的傳統(tǒng)灌溉施肥管理方式，大大地降低了水肥利用效率，我國(guó)設(shè)施蔬菜的單位面積水資源利用效率僅為以色列的 1/6～1/5，肥料利用率更低（張志斌，2008），這不僅浪費(fèi)了資源，同時(shí)大大增加了環(huán)境風(fēng)險(xiǎn)。與此同時(shí)，設(shè)施菜田長(zhǎng)期高量投入導(dǎo)致磷、鉀等養(yǎng)分在土壤中過量積累，使土壤酸化、鹽漬化（余海英 等，2006）、微生物群落失衡等，土壤退化問題普遍發(fā)生（王敬國(guó)，2011）。

土壤質(zhì)量是蔬菜生產(chǎn)體系可持續(xù)發(fā)展的基礎(chǔ)。為了解不同種植年限設(shè)施菜田土壤質(zhì)量的變化特點(diǎn)，本試驗(yàn)在山東壽光地區(qū)進(jìn)行了取樣調(diào)查和分析，從進(jìn)一步定量化土壤pH和土壤養(yǎng)分對(duì)比入手，探討不同種植年限設(shè)施菜田土壤質(zhì)量狀況，對(duì)指導(dǎo)設(shè)施蔬菜生產(chǎn)體系持續(xù)高效生產(chǎn)有著重要意義。

1 材料與方法

1.1 調(diào)查區(qū)概況

壽光地處魯北濱海平原，位于北緯 36°41′～37°19′，東經(jīng) 118°32′～119°10′之間，總面積2018 km2。氣候?qū)倥瘻貛Ъ撅L(fēng)性大陸氣候，四季分明。平均年降水量約600 mm，主要集中于6～8月，占全年降水總量的1/2以上。主要種植作物為蔬菜、小麥、玉米。

1.2 調(diào)查方法

2011年3月至4月中旬，隨機(jī)抽取壽光地區(qū)集約化設(shè)施蔬菜生產(chǎn)區(qū)的5個(gè)鎮(zhèn)（文家鎮(zhèn)、古城街道、紀(jì)臺(tái)鎮(zhèn)、稻田鎮(zhèn)、留呂鎮(zhèn)）的38個(gè)日光溫室進(jìn)行了調(diào)查，并進(jìn)行了相應(yīng)土樣的采集，同時(shí)調(diào)查了附近4塊小麥糧田作為對(duì)照。調(diào)查采用農(nóng)戶調(diào)研的形式，調(diào)查的內(nèi)容包括：采樣點(diǎn)基本情況（種植年限、種植面積、種植作物等）；施肥情況（肥料種類、施肥量等）；灌水情況（灌溉方式、灌溉次數(shù)、灌溉量和灌溉時(shí)期）（表1）。

1.3 土樣采集與測(cè)定方法

土樣采集的每個(gè)點(diǎn)均用GPS定位，日光溫室和小麥田的采集為3個(gè)點(diǎn)的混合土樣，采集深度分別為0～20、20～40、40～60 cm。采集的土樣經(jīng)風(fēng)干、研磨、過篩、混勻后進(jìn)行下一步分析。

土壤測(cè)定采用常規(guī)分析方法（鮑士旦，2000）：土壤 pH測(cè)定采用電位法（去離水浸提，水土比5∶1）；有機(jī)質(zhì)含量測(cè)定采用重鉻酸鉀外加熱法；NO3--N和NH4+-N含量測(cè)定采用新鮮土樣過4 mm篩，稱取12 g鮮土加入100 mL 1 mol·L-1KCl溶液浸提，振蕩1 h，過濾后用流動(dòng)注射分析儀（型號(hào)：TRAACS 2000）測(cè)定其含量，同時(shí)采用烘干法測(cè)定土壤含水量；Olsen-P含量測(cè)定采用 0.5 mol·L-1碳酸氫鈉浸提-鉬藍(lán)比色法；速效鉀含量測(cè)定采用 1 mol·L-1NH4OAc浸提-火焰光度法；電導(dǎo)率測(cè)定采用精密電導(dǎo)儀電導(dǎo)法。

1.4 數(shù)據(jù)處理

采用 Excel2010、Sigmaplot11.0、SPSS19.0軟件進(jìn)行數(shù)據(jù)記錄、作圖和統(tǒng)計(jì)分析。

表1 不同種植年限設(shè)施菜田及糧田調(diào)研點(diǎn)基本狀況

2 結(jié)果與分析

2.1 不同種植年限設(shè)施菜田土壤pH值和EC值的變化

設(shè)施菜田土壤pH值隨土層增加而增大，各土層土壤pH值隨種植年限的增加呈下降趨勢(shì)（圖1）。0～20 cm土層土壤在1～3、≥10 a種植年限下平均pH值為8.03、7.75，與相同土層露地小麥土壤相比，分別下降了0.41、0.69個(gè)單位；20～40 cm土層種植10 a后pH值由8.67下降至7.99，降幅為0.68；不同種植年限設(shè)施菜田0～20 cm土層土壤平均pH值為7.90，較露地小麥0～20 cm土層土壤下降了0.54個(gè)單位，這些充分說明了設(shè)施菜田土壤各土層有隨種植年限的增加呈現(xiàn)酸化的趨勢(shì)。

設(shè)施菜田各土層土壤電導(dǎo)率（EC值）隨種植年限的增加均有先升高后下降的趨勢(shì)（圖2），大部分設(shè)施菜田土壤已經(jīng)出現(xiàn)次生鹽漬化問題。0～20 cm土層土壤在種植1～3、4～6、7～9、≥10 a后，EC值分別在1.64～6.86、1.91～9.26、1.60～8.36、1.47～4.86 dS·m-1之間，平均值分別為3.66、5.71、4.61、3.50 dS·m-1，與0～20 cm土層露地小麥土壤平均值相比分別增加了1.90、3.95、2.85、1.74 dS·m-1；20～40、40～60 cm 土層土壤平均 EC 值分別為 2.49、2.11 dS·m-1。設(shè)施菜田土壤隨種植年限的增加，鹽分在土壤中的累積量逐漸增大，種植7 a后EC值開始下降，鹽分累積量減少。這種現(xiàn)象的變化可能與設(shè)施菜田的施肥量、施肥種類、灌水量、灌水方式、栽培作物以及鹽基離子在土壤中的移動(dòng)、淋失等有關(guān)。

圖1 種植年限對(duì)不同土層土壤pH值的影響

圖2 種植年限對(duì)不同土層土壤EC值的影響

設(shè)施菜田土壤中EC值的變化不僅表征了鹽分離子的變化，也可以反映鹽分離子對(duì)土壤pH值的影響。由圖 3可知，設(shè)施菜田 0～20 cm土層土壤 pH值與 EC值呈現(xiàn)一定的負(fù)相關(guān)性（y=-0.043x+8.149，r=-0.290，n=38），但并不顯著?；瘜W(xué)肥料尤其是氮肥的大量施用使土壤酸度增加（劉艷軍 等，2006），此外過量陰離子及磷鉀離子的陪伴也可能會(huì)隨土壤水分的蒸發(fā)使鹽分在表層不斷累積。

2.2 不同種植年限設(shè)施菜田土壤養(yǎng)分含量變化

2.2.1 土壤有機(jī)質(zhì) 設(shè)施菜田土壤有機(jī)質(zhì)含量隨種植年限的增加而逐漸上升（圖4）。0～20 cm土層土壤有機(jī)質(zhì)含量平均值為21.20 g·kg-1，種植1～3、4～6、7～9、≥10 a的土壤有機(jī)質(zhì)含量平均值分別為13.96、19.27、21.91、29.65 g·kg-1，與露地小麥相同土層土壤有機(jī)質(zhì)含量（13.16 g·kg-1）相比分別增加了 6.1%、46.4%、66.5%、125.3%。然而調(diào)查表明，壽光設(shè)施菜田全年施用的有機(jī)肥多以含速效養(yǎng)分多且分解較快的雞糞、豬糞、鴨糞為主，盡管施用量較大，達(dá)到每季14.4～100.0 t·hm-2，但肥料中木質(zhì)素和纖維素含量相對(duì)較低，不利于設(shè)施菜田土壤有機(jī)質(zhì)含量較大程度地提高，同時(shí)不同農(nóng)戶對(duì)菜田投入的有機(jī)肥量差異較大。

2.2.2 土壤無機(jī)氮 土壤無機(jī)氮以硝態(tài)氮形式為主，二者在各土層的變化趨勢(shì)均基本一致。設(shè)施菜田0～20 cm土層種植1～3、4～6、7～9、≥10 a的土壤無機(jī)氮含量分別為130.6、141.1、154.2、125.5 mg·kg-1，分別是露地小麥相同土層無機(jī)氮含量的2.83、3.06、3.34、2.72倍；20～40、40～60 cm土層不同種植年限（1～20 a）設(shè)施菜田土壤無機(jī)氮含量平均值分別為74.43、51.59 mg·kg-1，與相同土層露地小麥土壤無機(jī)氮含量相比，分別增加了40.95、17.12 mg·kg-1（圖5）。0～20 cm土層1～3 a設(shè)施菜田土壤硝態(tài)氮含量與相同土層露地小麥土壤差異顯著，且7～9 a設(shè)施菜田土壤硝態(tài)氮含量分別是露地小麥土壤、1～3 a設(shè)施菜田土壤的3.58、1.21倍（圖6）。可見，設(shè)施菜田土壤無機(jī)氮隨種植年限的增加在各土層出現(xiàn)不同程度累積，且設(shè)施菜田土壤在種植7～9 a后，無機(jī)氮和硝態(tài)氮累積量達(dá)到最大值。

圖3 0～20 cm土層土壤pH值和EC值的關(guān)系

圖4 種植年限對(duì)不同土層土壤有機(jī)質(zhì)含量的影響

圖5 種植年限對(duì)不同土層土壤無機(jī)氮含量的影響

圖6 種植年限對(duì)不同土層土壤硝態(tài)氮含量的影響

2.2.3 土壤有效磷 不同種植年限的設(shè)施菜田土壤有效磷含量差異較大（圖 7），種植年限與設(shè)施菜田土壤有效磷含量呈現(xiàn)極顯著的相關(guān)性（0～20、20～40 cm和40～60 cm土層相關(guān)系數(shù)r分別為 0.636**、0.528**和 0.493**，n=38）。設(shè)施菜田 0～20、20～40、40～60 cm土層土壤有效磷平均含量分別為露地小麥土壤的14.0、5.57、3.45倍；0～20 cm土層設(shè)施菜田和露地小麥土壤有效磷含量分別在37.36～553.2 mg·kg-1和7.58～23.1 mg·kg-1之間，平均含量分別為201.2 mg·kg-1和 14.4 mg·kg-1，≥10 a設(shè)施菜田土壤有效磷含量為263.7 mg·kg-1，1～3、4～6、7～9、≥10 a 0～20 cm土層土壤有效磷含量分別是露地小麥土壤的8.9、14.8、18.1、18.3倍。設(shè)施菜田各土層土壤有效磷含量隨種植年限的增加而逐漸累積，并且在 20～40 cm和40～60 cm土層累積明顯，這與磷素淋失有很大關(guān)系。

圖7 種植年限對(duì)不同土層土壤有效磷含量的影響

2.3.4 土壤有效鉀 不同種植年限的設(shè)施菜田土壤中有效鉀含量差異也很大（圖 8），種植年限與設(shè)施菜田土壤有效鉀含量呈現(xiàn)極顯著的相關(guān)性（0～20、20～40 cm和40～60 cm土層相關(guān)系數(shù)r分別為0.479**、0.559**和0.513**，n=38）。1～3 a的設(shè)施菜田土壤在20～40 cm、40～60 cm土層中有效鉀含量與≥10a的設(shè)施菜田土壤差異達(dá)到顯著水平。0～20 cm土層設(shè)施菜田土壤有效鉀含量為250.3～1334 mg·kg-1，平均值為628.6 mg·kg-1，露地小麥土壤有效鉀含量平均值為196.3 mg·kg-1，二者平均值相差3.20倍；20～40 cm土層設(shè)施菜田土壤有效鉀含量平均值為325.3 mg·kg-1，為對(duì)照露地小麥相同土層土壤的3.23倍；40～60 cm土層設(shè)施菜田土壤有效鉀含量最大值可達(dá)366.4 mg·kg-1，平均值為138.8 mg·kg-1，為露地小麥相同土層土壤的 1.34倍。與對(duì)照露地小麥土壤相比，1～3、4～6、7～9、≥10 a 0～20 cm 土層設(shè)施菜田土壤有效鉀含量分別是露地小麥土壤的 2.8、3.2、3.5、3.5倍。隨種植年限的增加設(shè)施菜田土壤中有效鉀含量逐漸累積，有效鉀在20～40 cm和40～60 cm土層累積與鉀在設(shè)施菜田土壤中的淋失有關(guān)。

圖8 種植年限對(duì)不同土層土壤有效鉀含量的影響

3 結(jié)論與討論

3.1 設(shè)施菜田養(yǎng)分累積及其環(huán)境風(fēng)險(xiǎn)

設(shè)施菜田單次灌溉量高達(dá)100～150 mm，可入滲到2 m以下的土壤深度（王敬國(guó)，2011），宋效宗（2007）調(diào)查壽光全年施氮肥量為618～1647 kg·hm-2，平均1154 kg·hm-2。設(shè)施菜田長(zhǎng)期較高的水肥投入，不利于水肥資源的高效利用和可持續(xù)發(fā)展，也給土壤本身帶來物理、化學(xué)以及生物學(xué)障礙等負(fù)面影響。

硝酸鹽是設(shè)施蔬菜喜好的氮素形態(tài)。本試驗(yàn)結(jié)果表明，設(shè)施菜田土壤無機(jī)氮以硝態(tài)氮形式為主，無論是對(duì)照露地小麥還是設(shè)施菜田，全土壤剖面0～20 cm土層土壤硝態(tài)氮含量最高，且菜田剖面深層有一定的硝酸鹽累積，試驗(yàn)結(jié)果同黨菊香等（2004）的結(jié)論一致。Cao等（2004）在浙江省嘉興地區(qū)調(diào)查發(fā)現(xiàn)，設(shè)施菜田土壤硝態(tài)氮含量是水稻—小麥輪作土壤的 4倍，本試驗(yàn)結(jié)果表明，1～3 a的0～20 cm土層設(shè)施菜田土壤硝態(tài)氮含量是露地小麥土壤的2.96倍。設(shè)施菜田土壤氮素累積以硝態(tài)氮為主，然而硝態(tài)氮本身帶負(fù)電荷，不易被土壤膠體吸附，易隨灌溉水向下淋洗（Zhang et al.，1996）。硝態(tài)氮在設(shè)施栽培土體內(nèi)的積聚最終會(huì)產(chǎn)生環(huán)境負(fù)效應(yīng)。一方面，表聚的硝態(tài)氮通過植物根系吸收進(jìn)入蔬菜可食部位，有可能導(dǎo)致可食部位硝態(tài)氮的積累量超出安全標(biāo)準(zhǔn)（劉偉 等，2011）；另一方面，由于壽光地區(qū)地下水位埋深較淺，設(shè)施菜田剖面積聚的硝態(tài)氮極易隨灌溉水向下淋洗至地下水系，威脅生態(tài)安全。據(jù)宋效宗等（2008）在2003～2005年監(jiān)測(cè)壽光地下水硝態(tài)氮含量的結(jié)果可知，設(shè)施蔬菜種植區(qū)灌溉水中硝態(tài)氮的平均含量是大田區(qū)的4～7倍，二者差異達(dá)到極顯著水平。在8次的取樣中，硝態(tài)氮含量平均在13.89～40.00 mg·L-1之間，均明顯高于 WHO推薦飲用水和我國(guó)現(xiàn)行飲用水硝態(tài)氮含量上限，并且單井硝態(tài)氮含量最高可達(dá)88.03 mg·L-1，可見設(shè)施菜田灌溉施肥的管理方式對(duì)地下水的影響已經(jīng)不容忽視。

Cao等（2004）在浙江省嘉興地區(qū)調(diào)查研究發(fā)現(xiàn)，設(shè)施菜田土壤有效磷含量是水稻—小麥輪作土壤的4～10倍。本試驗(yàn)中設(shè)施菜田0～20 cm土層土壤有效磷含量平均為露地小麥土壤的13.97倍，同劉兆輝等（2008）的研究結(jié)果相一致。0～20 cm土層土壤有效鉀含量平均值為628.6 mg·kg-1，為露地小麥土壤的3.20倍，這與大慶市設(shè)施菜田有效鉀（449.8～827.5 mg·kg-1）為露地有效鉀含量的 2.1～3.8倍的結(jié)果相類似（趙鳳艷 等，2000）。本試驗(yàn)結(jié)果表明，種植年限與設(shè)施菜田 0～20 cm土層土壤有效磷、有效鉀含量呈現(xiàn)極顯著的相關(guān)性（r分別為 0.636**、0.479**，n=38）。可見，設(shè)施菜田土壤有效磷、有效鉀含量存在著明顯的累積效應(yīng)。Heckrath等（1995）研究認(rèn)為，土壤Olsen-P濃度超過一定值時(shí)（約60 mg·kg-1），磷就會(huì)通過亞地表徑流損失。Hesketh和 Brookes（2000）在洛桑試驗(yàn)站的研究結(jié)果表明，土壤中 Olsen-P與灌溉水中的全磷和水溶性磷呈正相關(guān)，并提出土壤Olsen-P大于57 mg·kg-1時(shí)，灌溉水中磷濃度與土壤水溶性磷含量顯著增加，土壤磷淋失風(fēng)險(xiǎn)顯著增大。據(jù)袁麗金等（2010）對(duì)設(shè)施菜田土壤有效養(yǎng)分累積規(guī)律的研究表明，種植年限小于10 a的設(shè)施菜田在100 cm以下土層土壤有效磷累積峰出現(xiàn)在150 cm剖面處，為12.7 mg·kg-1，大于10 a的設(shè)施菜田在250 cm剖面處出現(xiàn)土壤有效磷的最大值，為19.5 mg·kg-1，同時(shí)也發(fā)現(xiàn)種植年限大于10 a的設(shè)施菜田土壤有效鉀在250～350 cm土層間存在累積峰，且峰值在95.7～119.3 mg·kg-1之間波動(dòng)，說明土壤有效磷、有效鉀在傳統(tǒng)水肥條件下淋洗現(xiàn)象明顯。由于設(shè)施菜田大量施用有機(jī)肥和化肥，土壤磷素和鉀素的逐年累積超過了土壤膠體的吸附能力，菜田耕層土壤磷鉀元素大量累積，有可能加快了土壤磷素通過滲漏淋溶流或地表徑流向土壤深層及水體環(huán)境遷移，加快了鉀素隨灌溉水的向下移動(dòng)，不僅大大降低了元素的利用效率，還給土壤和水體環(huán)境安全帶來了挑戰(zhàn)。

3.2 設(shè)施菜田養(yǎng)分累積與土壤pH值的變化

一般情況下，土壤中磷鉀元素移動(dòng)性較差，易被土壤膠體所吸附固定，硝態(tài)氮易隨水淋失。然而，設(shè)施菜田土壤養(yǎng)分累積已成為一個(gè)不爭(zhēng)的事實(shí)，大量的有效養(yǎng)分累積在土體剖面尤其是表層土壤中，極易隨灌溉水向下淋洗。通常伴隨硝態(tài)氮淋洗的是鹽基離子，而不是質(zhì)子（王敬國(guó)，2011），這樣就導(dǎo)致質(zhì)子在菜田表層土壤中的累積而出現(xiàn)土壤酸化的情況。此外，質(zhì)子與膠體吸附的鹽基離子交換后，在設(shè)施菜田頻繁的灌水情況下，鹽基離子被大量淋洗，而使土壤酸化加劇。

在沒有人類干擾的自然條件下，土壤的酸化過程需要幾百年甚至上百萬年的時(shí)間（Guo et al.，2010）。Guo等（2010）調(diào)查發(fā)現(xiàn)自20世紀(jì)80年代以來，中國(guó)經(jīng)濟(jì)作物類（蔬菜、果樹、茶樹）土壤pH值下降了0.3～0.8個(gè)單位，我國(guó)主要農(nóng)業(yè)地區(qū)的pH值都顯著降低，平均下降了約0.5個(gè)單位。本試驗(yàn)調(diào)查顯示不同種植年限各土層設(shè)施土壤pH值均大于7.0，且各土層土壤pH值隨種植年限的增加而下降。同時(shí)，此次調(diào)查不同種植年限設(shè)施菜田耕層土壤（0～20 cm）平均pH值為7.90，與劉兆輝等（2008）根據(jù)壽光25個(gè)土樣得出耕層設(shè)施土壤平均pH值（7.69）相比高出0.21個(gè)單位；與曾希柏等（2010）得到壽光設(shè)施菜田土壤pH值（6.86）的數(shù)據(jù)也有較大差別，說明壽光地區(qū)土壤pH值不再下降，反而有升高的趨勢(shì)。

菜田土壤 pH值的變化受土壤母質(zhì)、化肥種類數(shù)量、灌水定額、栽培作物、土壤溫度、種植年限以及土壤微生物等因素的影響，有機(jī)物質(zhì)的加入也會(huì)影響到土壤pH值的變化（Feng et al.，1996）。設(shè)施菜田土壤pH值的整體變化與前人的研究成果（李俊良 等，2002；王輝 等，2005；曾希柏 等，2010）不同，原因可能是人們飼養(yǎng)雞、鴨等動(dòng)物的過程中，為了除臭和滅菌而加入了大量的石灰作為墊料，同時(shí)在有機(jī)肥（雞糞、鴨糞、豬糞等）加工包裝過程中為了除臭也摻雜了部分石灰，導(dǎo)致有機(jī)肥含有的碳酸鈣含量較高。據(jù)調(diào)查，壽光地區(qū)施用的風(fēng)干雞糞的碳酸鈣含量平均值為76.63 g·kg-1（n=20），且目前風(fēng)干雞糞的pH值在7.80±0.38（任濤，2011）。有機(jī)肥（雞糞、鴨糞、豬糞等）在設(shè)施菜田中連年大量的施用，不但使有機(jī)質(zhì)在設(shè)施土壤尤其是0～20 cm土層土壤中逐年累積，而且也使菜田土壤pH值增加，一定程度上緩解了設(shè)施菜田土壤的酸化現(xiàn)象。此外，設(shè)施菜田土壤pH值的增加還可能與當(dāng)?shù)毓喔人膲A度、區(qū)域地理位置、作物對(duì)P、S等酸性陰離子基團(tuán)的吸收、農(nóng)戶進(jìn)行的“客土改良”、施用堿性肥料以及加入秸稈等有關(guān)。

鮑士旦．2000．土壤農(nóng)化分析．北京：中國(guó)農(nóng)業(yè)出版社.

黨菊香，郭文龍，郭俊煒，呂家瓏，王軍利．2004．不同種植年限蔬菜大棚土壤鹽分累積及硝態(tài)氮遷移規(guī)律．中國(guó)農(nóng)學(xué)通報(bào)，20（6）：189-191.

杜會(huì)英．2007．保護(hù)地蔬菜氮肥利用、土壤養(yǎng)分和鹽分累積特征研究〔博士論文〕．北京：中國(guó)農(nóng)業(yè)科學(xué)院.

高中強(qiáng)．2010．山東省設(shè)施蔬菜發(fā)展現(xiàn)狀、問題及對(duì)策建議．中國(guó)果菜，（2）：12-14．

黃錦法，曹志洪，李艾芬，張蠶生．2003．稻麥輪作改為保護(hù)地菜田土壤肥力質(zhì)量的演變．植物營(yíng)養(yǎng)與肥料學(xué)報(bào)，9（1）：19-25.

李俊良，崔德杰，孟祥霞，李曉林，張福鎖．2002．山東壽光保護(hù)地蔬菜施肥現(xiàn)狀及問題的研究．土壤通報(bào)，33（2）：126-128.

劉偉，武美燕，胡學(xué)玉，艾天成．2011．設(shè)施栽培地土壤富營(yíng)養(yǎng)化及其潛在的環(huán)境影響．環(huán)境科學(xué)與技術(shù)，34（5）：20-23.

劉兆輝，江麗華，張文君，鄭福麗，王梅，林海濤．2008．山東省設(shè)施蔬菜施肥量演變及土壤養(yǎng)分變化規(guī)律．土壤學(xué)報(bào)，45（2）：296-303．

劉艷軍，李琪，梁文舉．2006．設(shè)施栽培對(duì)土壤環(huán)境質(zhì)量的影響．遼寧工程技術(shù)大學(xué)學(xué)報(bào)，25（3）：465-470.

任濤．2011．不同氮肥及有機(jī)肥投入對(duì)設(shè)施番茄土壤碳氮去向的影響〔博士論文〕．北京：中國(guó)農(nóng)業(yè)大學(xué)．

宋效宗．2007．保護(hù)地生產(chǎn)中硝酸鹽的淋洗及其對(duì)地下水的影響〔博士論文〕．北京：中國(guó)農(nóng)業(yè)大學(xué)．

宋效宗，趙長(zhǎng)星，李季，王小蘭，吳鋼．2008．兩種種植體系下地下水硝態(tài)氮含量變化．生態(tài)學(xué)報(bào)，28（11）：5513-5520.

王輝，董元華，李德成，安瓊．2005．不同種植年限大棚蔬菜地土壤養(yǎng)分狀況研究．土壤，37（4）：460-462.

王敬國(guó)．2011．設(shè)施菜田退化土壤修復(fù)與資源高效利用．北京：中國(guó)農(nóng)業(yè)大學(xué)出版社.

余海英，李廷軒，周健民．2006．典型設(shè)施栽培土壤鹽分變化規(guī)律及潛在的環(huán)境效應(yīng)研究．土壤學(xué)報(bào)，43（4）：571-576．

余海英，李廷軒，張錫洲．2010．溫室栽培系統(tǒng)的養(yǎng)分平衡及土壤養(yǎng)分變化特征．中國(guó)農(nóng)業(yè)科學(xué)，43（3）：514-522．

袁麗金，巨曉棠，張麗娟，王玨，楊志新．2010．設(shè)施蔬菜土壤剖面氮磷鉀積累及對(duì)地下水的影響．中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，18（1）：14-19．曾希柏，白玲玉，蘇世鳴，李蓮芳．2010．山東壽光不同種植年限設(shè)施土壤的酸化與鹽漬化．生態(tài)學(xué)報(bào)，30（7）：1853-1859．

張志斌．2008．我國(guó)設(shè)施蔬菜存在的問題及發(fā)展重點(diǎn)．中國(guó)蔬菜，（5）：1-3．

趙鳳艷，吳鳳芝，劉德，姚禾雨．2000．大棚菜地土壤理化特性的研究．土壤肥料，（2）：11-13．

中國(guó)農(nóng)業(yè)統(tǒng)計(jì)年鑒．2010．北京：中國(guó)農(nóng)業(yè)出版社．

Cao Z H，Huang J F，Zhang C S．2004．Soil quality evolution after land use change from paddy soil to vegetable land．Environmental Geochemistry and Health，26（2）：97-103．

Feng Y，Schubert S，Mengel K．1996．Soil pH increase due to biological decarboxylation of organic antions．Soil Biology & Biochemistry，28：617-624．

Guo J H，Liu X J，Zhang Y，Shen J L，Han W X，Zhang W F，Christie P，Goulding K W T，Vitousek P M，Zhang F S．2010．Significant acidification in major Chinese croplands．Science，327：1008-1010．

Heckrath G，Brookes P C，Poulton P R，Goulding K W T．1995．Phosphorus leaching from soils containing different phosphorus concentrations in the Broadbalk experiment．Journal of Environmental Quality，24：904-910．

Hesketh N，Brookes P C．2000．Development of an indicator for risk of phosphorus leaching．Journal of Environmental Quality，29（1）：105-110．

Zhang W L，Tian Z X，Zhang N，Li X Q．1996．Nitrate pollution of groundwater in northern China．Agriculture，Ecosystems and Environment，59：223-231．