李 強(qiáng)，戴美玲，向德明，黎 娟，周米良，田 峰，田明慧，張黎明1,

?

湘西喀斯特區(qū)植煙土壤有效磷時空變異及風(fēng)險評估①

李 強(qiáng)1,2，戴美玲2，向德明3，黎 娟1*，周米良3，田 峰3，田明慧3，張黎明1,3

(1湖南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院，長沙 410128；2湖南農(nóng)業(yè)大學(xué)煙草研究院，長沙 410128；3湖南省煙草公司湘西州公司，湖南吉首 416000)

為揭示典型喀斯特?zé)焻^(qū)土壤有效磷的時空變異特征，為喀斯特地貌煙區(qū)合理施磷和控制磷素面源污染提供依據(jù)，以湘西州2000年和2014年耕地土壤有效磷含量為研究對象，利用空間自相關(guān)、半方差函數(shù)及分形維數(shù)等方法，對縣域內(nèi)土壤速效養(yǎng)分時空變異特征和環(huán)境風(fēng)險進(jìn)行了評價。結(jié)果表明，2014年土壤有效磷平均含量增加了30.89 mg/kg，增幅達(dá)425.48%；土壤有效磷Moran’s I 指數(shù)標(biāo)準(zhǔn)化值下降，分形維數(shù)和塊金效應(yīng)增加，說明土壤有效磷空間結(jié)構(gòu)減弱，隨機(jī)變異性增強(qiáng)；土壤有效磷“極高”、“很高”和“高”等級面積顯著增加，分別增加了28.36%、26.89% 和25.94%；而“低”和“極低”等級面積分別下降了57.45% 和27.12%；土壤有效磷含量與海拔、成土母質(zhì)、土地利用類型、輪作模式及有機(jī)質(zhì)含量關(guān)系密切；2014年湘西州有效磷含量超過風(fēng)險臨界值的概率大幅增加，很高風(fēng)險(概率>80%)、高風(fēng)險(概率60% ~ 80%)和中等風(fēng)險(概率40% ~ 60%)區(qū)域分別占研究區(qū)域面積的1.77%、5.40% 和13.96%。其中，龍山縣北部和永順縣東部應(yīng)作為農(nóng)業(yè)面源磷污染的重點(diǎn)控制區(qū)。

喀斯特；植煙土壤；有效磷；時空變異；環(huán)境風(fēng)險

磷是煙草必需的大量元素之一，一般占煙株干重的0.15% ~ 0.60%，在光合作用、呼吸作用、Krebs循環(huán)和氮代謝過程中起著重要的生理作用[1-2]。土壤磷是植物磷素營養(yǎng)的主要來源，其含量高低和有效性是衡量土壤肥力的主要指標(biāo)之一，對植物生長發(fā)育有重要影響。近30多年來，我國磷肥投入大幅增加，據(jù)報道從1980年到2011年，我國磷肥消費(fèi)量由270萬t上升到820萬t，增長了3倍，中國土壤出現(xiàn)磷盈余，磷盈余總量達(dá)9.8×105t，磷盈余強(qiáng)度2.5 kg/hm2，植煙區(qū)土壤磷素的供應(yīng)能力大幅提高[3-4]。由于烤煙特殊的施肥習(xí)慣，一般常采取氮磷鉀比例為1:1 ~ 1.5:3.0，加上烤煙的磷肥利用率較低，導(dǎo)致植煙土壤比一般農(nóng)田的磷素盈余更為嚴(yán)重[5]，甚至造成一些煙區(qū)水體磷污染[6]。此外，當(dāng)土壤磷素過高對煙葉品質(zhì)也有不良影響，如導(dǎo)致煙株早花，葉脈突出、組織粗糙、彈性差、油分少、易破碎。中國煙草種植主要集中在南部和西南部，地形多為山地和丘陵，過多的磷素隨雨水淋失的幾率較大，極易造成環(huán)境污染風(fēng)險。而有關(guān)喀斯特?zé)焻^(qū)土壤磷素環(huán)境風(fēng)險的研究卻鮮見報道。湖南湘西是典型的喀斯特地貌，烤煙種植歷史悠久，過多的磷素極易隨地表徑流淋失。鑒于此，本文擬采用地統(tǒng)計學(xué)、克里格插值法(Kriging)和概率克里格插值法(Probability Kriging)，研究典型喀斯特植煙區(qū)的磷素環(huán)境風(fēng)險，旨在為煙區(qū)烤煙磷素管理和環(huán)境保護(hù)提供理論依據(jù)。

1 材料與方法

1.1 區(qū)域自然概況

湘西州位于湖南省西北部，地處湘鄂黔渝4省市交界處，是湖南省烤煙主產(chǎn)區(qū)之一，位于109°10′ ~ 110°22.5′E，27°44.5′ ~ 29°38′N 之間，全州東西寬約170 km，南北長約240 km，總面積15 462 km2。湘西州地勢由西北向東南傾斜，平均海拔800 ~ 200 m，地貌形態(tài)的總體輪廓是一個以山原山地為主，兼有丘陵和小平原，并向北西突出的弧形山區(qū)地貌。西南石灰?guī)r分布極廣，巖溶發(fā)育充分，多溶洞、伏流；西北石英砂巖密布，因地殼作用形成小片峰。東西部為低山丘陵區(qū)，平均海撥200 ~ 500 m，溪河縱橫其間，兩岸多沖積平原。湘西州屬亞熱帶季風(fēng)濕潤氣候，具有明顯的大陸性氣候特征。年平均氣溫為16.5 ~ 17.5 ℃，平均降水1 290 ~ 1 600 mm，日照時數(shù)平均在1 406 ~ 1 219 h。土壤類型有水稻土、潮土、紅壤、黃壤、黃棕壤、石灰土和紫色土等。湘西州基本煙田3.07×105hm2，其中水旱輪作煙田1.32×105hm2，旱地?zé)熖锩娣e1.75×105hm2，主要種植烤煙品種為云煙87和K326，全州年均煙葉產(chǎn)量2.25×104t，輪作模式主要為烤煙–綠肥–玉米–綠肥–烤煙和烤煙–綠肥–水稻–綠肥–烤煙，烤煙施肥采用基肥+追肥的管理模式。

1.2 樣品采集和分析

土壤取樣布點(diǎn)按基本煙田面積進(jìn)行，第1次取樣于2000年12月進(jìn)行，每20 ~ 33.33 hm2布置1個樣點(diǎn)，共采集446個土樣(表1)；第2次取樣于2014年11—12月進(jìn)行，每13.33 ~ 26.67 hm2布置1個取樣點(diǎn)，采集土樣1 242個(表1)。采集方法為在土壤冬翻前，選取667 m2以上的田塊進(jìn)行取樣，用手持式GPS定位，記錄田塊中心的地理信息，采用五點(diǎn)取樣法或“S”形取0 ~ 20 cm土樣，每個田塊確保5點(diǎn)以上，充分混勻后用四分法取大約1 kg土樣。兩期土樣的有效磷均采用鉬銻抗比色法測定[7]。

表1 土壤取樣信息

1.3 土壤有效磷評價標(biāo)準(zhǔn)

土壤有效磷分級標(biāo)準(zhǔn)參考全國第二次土壤普查結(jié)果和文獻(xiàn)[8]制定，將湘西植煙土壤有效磷含量分為極低(<5 mg/kg)、低(5 ~ 10 mg/kg)、適宜(10 ~ 20 mg/kg)、高(20 ~ 30 mg/kg)、很高(30 ~ 40 mg/kg)和極高(>40 mg/kg)(環(huán)境風(fēng)險級)6個等級。

1.4 數(shù)據(jù)分析

采用SPSS17.0進(jìn)行描述性統(tǒng)計和K–S檢驗(yàn)(Kolmogorov-Smirnov test)[9]。實(shí)驗(yàn)半方差函數(shù)模型的計算和理論模型擬合采用GS+9.0軟件完成[10]，克里格插值(Kriging)、概率克里格插值(Probability Kriging)、繪圖及面積統(tǒng)計均在ArcGIS 10.2.2軟件實(shí)現(xiàn)[8,11-12]。

2 結(jié)果與分析

2.1 植煙土壤有效磷的基本統(tǒng)計特征及年代變化

2個時期湘西州土壤有效磷含量列于表2。14年間有效磷含量發(fā)生巨大變化，2014年土壤有效磷含量平均值較2000年增加了30.89 mg/kg，增幅達(dá)425.48%，由“低”等級變?yōu)椤昂芨摺钡燃?。有效磷的變異系?shù)由2000年的85.09% 上升到2014年的90.06%，增加了4.97%。極小值變化較小，而極大值則由2000年的52.60 mg/kg上升到2014年的234.00 mg/kg，增加了181.40 mg/kg。極差由2000年的52.30 mg/kg上升到2014年的233.16 mg/kg，增加了180.86 mg/kg。表明不僅有效磷的含量發(fā)生巨大變化，同時有效磷含量的變異也較大。進(jìn)一步對湘西植煙土壤有效磷含量的等級分布情況進(jìn)行分析(表3)，發(fā)現(xiàn)與2000年相比，2014年土壤有效磷“極低”和“低”的樣品比例分別下降了37.60和26.16個百分點(diǎn)，相應(yīng)地“適宜”、“高”、“很高”和“極高”等級的樣品比例分別增加了3.40、15.18、11.97和33.21個百分點(diǎn)。

2.2 植煙土壤有效磷時空變異特征

土壤養(yǎng)分變異來自結(jié)構(gòu)因子和隨機(jī)因子，結(jié)構(gòu)因子包括土壤母質(zhì)、地形、氣候等非人為因素，隨機(jī)因子主要包括施肥、耕作和灌溉等人為因素；前者引起系統(tǒng)變異，并使變量呈空間相關(guān)性，后者造成隨機(jī)變異，表現(xiàn)為削弱變量的空間相關(guān)性[13]。半方差函數(shù)用于描述區(qū)域化變量的空間變異，可避免樣本方差的影響[11,14]，本試驗(yàn)中分別用不同模型對土壤有效磷進(jìn)行擬合，獲取土壤有效磷的最優(yōu)半方差函數(shù)模型及其相關(guān)參數(shù)(表4)，RMSE接近1且MSE接近0，表明選取的模型具有很高的擬合精度，能夠很好地反映土壤有效磷的空間結(jié)構(gòu)特征。2個時期土壤有效磷最佳半方差函數(shù)模型均為指數(shù)模型，2000年和2014年土壤有效磷的塊金效應(yīng)分別為8.90% 和65.66%，14年間有效磷塊金效應(yīng)增幅較大，反映出結(jié)構(gòu)性因素作用和隨機(jī)性因素作用發(fā)生較大變化[10,15]，表明在過去的14年間，植煙土壤有效磷的空間結(jié)構(gòu)性減弱，隨機(jī)變異性增強(qiáng)，這與磷肥、有機(jī)肥投入等隨機(jī)因子對土壤有效磷含量影響程度增大有關(guān)。

表2 不同時期湘西植煙土壤有效磷含量狀況(mg/kg)

表3 不同時期湘西州土壤樣品有效磷等級分布及變化

表4 土壤有效磷半方差函數(shù)模型及相關(guān)參數(shù)比較

注: RMSE為標(biāo)準(zhǔn)化均方根誤差；MSE 為標(biāo)準(zhǔn)化平均誤差。

此外，從2000年到2014年，土壤有效磷的Moran’s I 值由0.388變?yōu)?.097。經(jīng)標(biāo)準(zhǔn)化計算，兩年Moran’s I 的標(biāo)準(zhǔn)化值均大于2.58，說明2個時期耕地土壤有效磷含量在0.05統(tǒng)計水平下顯著空間自相關(guān)，2014年標(biāo)準(zhǔn)化值低于2000年，表明14年來土壤有效磷的空間自相關(guān)性在減弱，隨機(jī)性因素的影響增強(qiáng)，與前文塊金效應(yīng)的分析一致。

除空間結(jié)構(gòu)和空間自相關(guān)分析以外，進(jìn)一步利用分形維數(shù)對兩個時期土壤有效磷的空間結(jié)構(gòu)特點(diǎn)進(jìn)行定量描述(表4)。土壤有效磷的分形維數(shù)由2000年的1.915變?yōu)?014年的1.984，升高了0.069，14年間研究區(qū)土壤有效磷呈現(xiàn)出更多較小尺度上的變異特點(diǎn)，其隨機(jī)變異比例相對增加，這也印證了前文半方差變異函數(shù)和空間自相關(guān)的分析結(jié)果。

2.3 土壤有效磷時空分布變化

根據(jù)獲取的半方差函數(shù)模型，利用普通Kriging最優(yōu)內(nèi)插法，按照分級標(biāo)準(zhǔn)繪制了兩個時期土壤有效磷含量空間分布圖(圖1)，并利用ArcGIS軟件面積統(tǒng)計模塊對繪制的有效磷空間分布圖進(jìn)行面積統(tǒng)計。兩個時期土壤有效磷含量空間分布規(guī)律均不明顯，2014年湘西植煙土壤有效磷含量分級面積比例與2000年相比有較大變化(圖1和表5)。2000年有效磷含量總體較低，“低”和“極低”的植煙面積比例分別達(dá)58.94% 和27.12%，其中有效磷“低”等級在各個植煙縣均有大面積分布，有效磷“極低”等級主要分布在中部和南部；有效磷“適宜”的植煙面積比例為13.37%，主要分布在東部的永順縣和古丈縣；有效磷“高”和“很高”等級僅零星分布在龍山縣北部。2014年植煙土壤有效磷含量較2000年有大幅增加，新增了2000年未出現(xiàn)的“極高”等級面積達(dá)28.36%，“高”和“很高”等級亦由原來的零星分布分別增加至26.46% 和26.94%，“適宜”等級的面積增加至16.75%，“低”等級面積下降至1.49%，“極低”等級則完全消失。其中，“極高”等級在各植煙縣均有分布，在北部和東部分布較為集中；“很高”和“高”等級的分布較為分散，呈小斑塊狀分布在各縣；“適宜”等級以較大斑塊狀散落分布在各植煙縣；“低”等級僅零星分布在古丈、保靖和龍山。從兩個年份的變化來看，湘西植煙土壤有效磷含量變化劇烈?？梢?，2014年土壤有效磷“極高”等級、“很高”和“高”等級面積顯著增加，分別增加了28.36%、26.89% 和25.94%；而“低”和“極低”等級則大幅下降，分別比2000年下降了57.45% 和27.12%。14年來湘西植煙土壤有效磷含量呈大幅增加趨勢，這與磷肥的長期投入和逐年積累有關(guān)。

圖1 湘西植煙土壤磷素含量空間分布

表5 不同時期土壤有效磷各等級面積統(tǒng)計及變化

2.4 土壤有效磷的影響因素

2.4.1 海拔 湘西煙區(qū)從海拔100 m到1 400 m均有烤煙種植，土壤有效磷與海拔的相關(guān)分析表明，土壤有效磷含量與海拔呈極顯著正相關(guān)關(guān)系(= –0.254，= 0.000)。進(jìn)一步采用決策樹模型分析海拔高度對土壤有效磷含量的影響(圖2)，結(jié)果表明海拔對土壤有機(jī)質(zhì)含量影響的拐點(diǎn)為467、839.7和1 023 m，3個拐點(diǎn)將海拔劃成4個分組：≤467 m、467 ~ 839.7 m、839.7 ~ 1 023 m和>1 023 m，4個海拔分組下土壤有效磷含量分別為27.416、36.013、47.930和62.741 mg/kg，土壤有效磷在海拔分組間差異達(dá)到極顯著水平，呈現(xiàn)隨海拔升高有效磷含量快速升高的趨勢。

2.4.2 成土母質(zhì) 龍山植煙土壤主要成土母質(zhì)為坡積物、洪積物、殘積物和沖積物，有效磷含量在4種成土母質(zhì)間存在極顯著差異(表6)，由高到低依次為坡積物、洪積物、殘積物和沖積物，其中坡積物有效磷含量已基本達(dá)到淋失風(fēng)險臨界值，洪積物為很高等級，殘積物和沖積物為高等級。從變異程度來看，坡積物、洪積物和沖積物為中等強(qiáng)度變異，殘積物表現(xiàn)為強(qiáng)變異。

2.4.3 土地利用類型和輪作制度 湘西植煙土壤的利用類型為旱地和水田，土壤有效磷在2種土地利用類型間差異達(dá)極顯著水平，旱地土壤有效磷含量平均為40.76 mg/kg，水田土壤有效磷含量平均為29.34 mg/kg (表7)。湘西烤煙種植的輪作模式主要有煙–空閑、煙–玉米、煙–綠肥、煙–稻和煙–油菜，不同輪作模式土壤有效磷差異極顯著，以煙–空閑–煙、煙–玉米–煙輪作模式下的土壤有效磷含量最高，顯著高于煙–綠肥–煙、煙–稻–煙、煙–油菜–煙(表8)，這與煙–空閑–煙、煙–玉米–煙2種輪作模式的磷素投入較大有關(guān)。

圖2 海拔高度對土壤有效磷含量影響的決策樹分析

表6 湘西不同成土母質(zhì)土壤有效磷含量(mg/kg)

注： 表中同列數(shù)據(jù)小寫字母不同表示差異達(dá)到<0.05顯著水平，大寫字母不同表示差異達(dá)到<0.01顯著水平，下同。

表7 湘西不同土地利用類型土壤有效磷含量(mg/kg)

2.4.4 有機(jī)質(zhì)含量 湘西植煙土壤有機(jī)質(zhì)含量平均為28.55 g/kg，處于適宜烤煙生長的范圍。土壤有效磷含量與有機(jī)質(zhì)含量呈極顯著正相關(guān)關(guān)系(= 0.206，= 0.000)。進(jìn)一步采用決策樹模型分析土壤有機(jī)質(zhì)含量對主要土壤養(yǎng)分指標(biāo)的影響(圖3)，結(jié)果表明土壤有機(jī)質(zhì)含量對有效磷影響的拐點(diǎn)為16.5、22.50和27.00 g/kg，該拐點(diǎn)將有機(jī)質(zhì)劃成4個分組：≤16.5 g/kg、16.5 ~ 22.50 g/kg、22.50 ~ 27.00 g/kg和>27.00 g/kg，4組土壤有效磷含量分別為15.034、29.359、38.845和45.965 mg/kg，土壤有效磷在這4個分組間差異達(dá)到極顯著水平，呈現(xiàn)隨有機(jī)質(zhì)增加土壤有效磷含量快速升高的趨勢。

表8 湘西不同輪作模式土壤有效磷含量(mg/kg)

圖3 有機(jī)質(zhì)對土壤有效磷含量影響的決策樹分析

2.5 土壤有效磷含量的環(huán)境風(fēng)險評價

以上分析顯示目前湘西植煙土壤磷素盈余較為嚴(yán)重。磷素移動性較小，且在土壤中易被固定，盈余的磷素在耕層土壤中不斷積累，在地表徑流和土壤侵蝕等因素的作用下，其對水體污染的風(fēng)險增大，有研究發(fā)現(xiàn)土壤有效磷達(dá)33 ~ 36 mg/kg后，土壤磷素便極易發(fā)生磷素淋溶[16]。本文參照這一研究結(jié)果，以40 mg/kg為磷素淋溶臨界值，采用概率克里格插值法對湘西植煙土壤有效磷含量的環(huán)境風(fēng)險進(jìn)行評價，結(jié)果見圖4。

圖4 湘西植煙土壤有效磷含量超過環(huán)境臨界值的概率

與圖1比較發(fā)現(xiàn)，有效磷含量高的區(qū)域其超過風(fēng)險臨界值的概率也較高，與空間分布的趨勢基本一致。2000年全市有效磷含量超過風(fēng)險臨界值的概率普遍較低，概率區(qū)間60% ~ 80%(高風(fēng)險)的面積主要分布在花垣縣，占研究區(qū)面積的1.16%，概率區(qū)間40% ~ 60%(中等風(fēng)險)的面積僅零星分布在永順、花垣和保靖，占研究區(qū)面積的4.41%；概率區(qū)間20% ~ 40%(低風(fēng)險)主要分布在永順、花垣和保靖，占研究區(qū)面積的16.37%；概率小于20%(無風(fēng)險)的區(qū)域面積最大，占總面積的78.06%，分布在全州各縣。2014年全市有效磷含量超過風(fēng)險臨界值的概率大幅增加，概率>80%(很高風(fēng)險)的區(qū)域以小斑塊狀分布龍山縣、永順縣和花垣縣，占研究區(qū)域面積的1.77%；概率區(qū)間60% ~ 80%(高風(fēng)險)的區(qū)域分散分布在全州7個植煙縣中的4個，占研究區(qū)域面積的5.40%；概率區(qū)間40% ~ 60%(中等風(fēng)險)的區(qū)域在全州7個植煙縣均有分布，占研究區(qū)域面積的13.96%；低風(fēng)險和無風(fēng)險區(qū)在全州7個植煙縣均有分布，分布占研究區(qū)域面積的37.42% 和41.44%；14年來土壤磷素環(huán)境風(fēng)險大幅增加，不同風(fēng)險級別呈點(diǎn)狀和片狀同時發(fā)生，說明磷素環(huán)境風(fēng)險的發(fā)生普遍存在。

3 討論

2000年到2014年湘西植煙土壤有效磷平均含量由7.26 mg/kg上升到38.15 mg/kg，增加了30.89 mg/kg，增幅高達(dá)425.48%，整體上由“低”等級變?yōu)椤昂芨摺钡燃?；從不同等級的面積變化來看，2014年土壤有效磷“極高”、“很高”和“高”等級面積顯著增加，分別增加了28.36%、26.89% 和25.94%；而“低”和“極低”等級則大幅下降，分別比2000年下降了57.45% 和27.12%；14年間，研究區(qū)域土壤有效磷環(huán)境風(fēng)險概率顯著上升，其中中等風(fēng)險區(qū)面積增加了9.56%，高風(fēng)險和很高風(fēng)險區(qū)面積分別增加了4.25% 和1.77%，無風(fēng)險區(qū)下降了36.62%，低風(fēng)險區(qū)增加了21.05%。與他人關(guān)于我國很多農(nóng)區(qū)土壤有效磷含量大幅增加的研究結(jié)果基本一致[17-20]?？緹煂τ诹姿毓?yīng)的要求是足而不過多，一般植煙土壤有效磷含量10 ~ 20 mg/kg即可滿足烤煙生長。過多的磷素營養(yǎng)對于煙葉品質(zhì)有顯著的不利影響，研究中發(fā)現(xiàn)高磷土壤上種植烤煙，煙葉的木質(zhì)素和細(xì)胞壁物質(zhì)總量呈增加趨勢[21]，而在高磷土壤上過量施用磷肥則會造成化學(xué)成分協(xié)調(diào)性變差，感官質(zhì)量變差[22]。

本文的研究結(jié)果顯示海拔高度、土壤母質(zhì)、土地利用類型、土壤輪作模式均會對土壤有效磷產(chǎn)生影響。據(jù)文獻(xiàn)報道，我國不少煙區(qū)土壤有效磷含量高于30 mg/kg，甚至超過40 mg/kg[23-29]，一方面與我國農(nóng)業(yè)生產(chǎn)普遍磷肥投入較高有關(guān)，此外還與烤煙種植的特有施肥習(xí)慣關(guān)系密切，我國烤煙施肥普遍采取的氮磷鉀比例為1:1 ~ 1.5:3.0，而烤煙磷肥的當(dāng)季利用率卻僅為15% 左右[30]。據(jù)調(diào)查，湘西植煙土壤每667 m2投入磷素達(dá)8 ~ 10 kg/a，外源磷的大量、長期投入是植煙土壤磷素過高的主要原因，文中不同輪作制度下土壤有效磷含量的差異有力支撐這一論點(diǎn)，這與前人的研究一致[31-32]；據(jù)報道，在我國當(dāng)前的烤煙生產(chǎn)中，磷肥的施用量已遠(yuǎn)遠(yuǎn)超過烤煙需要量，導(dǎo)致土壤磷素積累，增加土壤磷素的流失風(fēng)險[33-34]。當(dāng)土壤磷素超過32 mg/kg時，施磷4 kg/667 m2(約為當(dāng)前常規(guī)施磷水平的一半)，煙葉外觀質(zhì)量和內(nèi)在質(zhì)量均明顯變差[35]，筆者在云南羅平煙區(qū)的調(diào)研有力支持了這一觀點(diǎn)，當(dāng)?shù)夭糠种矡煷髴粼谟行Я?5 mg/kg以上的植煙土壤上均不施用磷肥，但產(chǎn)值產(chǎn)量均和施用磷肥的沒有區(qū)別。因此，針對當(dāng)前研究區(qū)土壤磷素現(xiàn)狀，建議適當(dāng)控制磷肥用量，推行測土施肥技術(shù)，有效磷在30 mg/kg以上的煙田土壤可減少磷肥用量50% 以上，有效磷在40 mg/kg以上的煙田土壤可嘗試開展“啟動性施磷”技術(shù)[36]和隔年施磷技術(shù)[37]，以保護(hù)喀斯特區(qū)植煙生態(tài)環(huán)境和保障煙葉生產(chǎn)的可持續(xù)發(fā)展。

[1] 王海英, 許亦農(nóng). 缺磷脅迫下煙草葉片磷脂酰甘油(PG)含量降低的分子機(jī)理[J]. 植物生理與分子生物學(xué)學(xué)報, 2006, 32(6): 643–648

[2] 胡國松, 鄭偉, 王震東, 等. 烤煙營養(yǎng)原理[M]. 北京: 科學(xué)出版社, 2000: 94–118

[3] 曹志洪, 周健民. 中國土壤質(zhì)量[M]. 北京: 科學(xué)出版社, 2008

[4] 中華人民共和國國家統(tǒng)計局. 中國統(tǒng)計年鑒[M]. 北京: 中國統(tǒng)計出版社, 2012

[5] 許自成, 王林, 肖漢乾. 湖南煙區(qū)烤煙磷含量與土壤磷素的分布特點(diǎn)及關(guān)系分析[J]. 浙江大學(xué)學(xué)報(農(nóng)業(yè)與生命科學(xué)版), 2007, 33(3): 290–297

[6] 褚利平, 王克勤, 宋澤芬, 等. 烤煙坡耕地壤中流氮、磷濃度的動態(tài)特征[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2010, 29(7): 1346–1354

[7] 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國農(nóng)業(yè)科技出版社, 2000

[8] 王淑英, 胡克林, 路蘋, 等. 北京平谷區(qū)土壤有效磷的空間變異特征及其環(huán)境風(fēng)險評價[J]. 中國農(nóng)業(yè)科學(xué), 2009, 42(4): 1290–1298

[9] 郝黎仁, 攀元, 郝哲歐. SPSS實(shí)用統(tǒng)計分析[M]. 北京: 中國水利水電出版社, 2002: 110–135

[10] 王政權(quán). 地統(tǒng)計學(xué)及在生態(tài)學(xué)中的應(yīng)用[M]. 北京: 科學(xué)出版社, 1999: 65–132

[11] 路鵬, 彭佩欽, 宋變蘭, 等. 洞庭湖平原區(qū)土壤全磷含量地統(tǒng)計學(xué)和GIS分析[J]. 中國農(nóng)業(yè)科學(xué), 2005, 38(6): 1204–1212

[12] 吳秀芹, 張洪巖, 李瑞改, 等. ArcGIS9地理信息系統(tǒng)應(yīng)用與實(shí)踐(下冊)[M]. 北京: 清華大學(xué)出版社, 2007: 463–517

[13] 張久明, 宿慶瑞, 遲鳳琴, 等. 嫩江縣土壤速效養(yǎng)分差異性分析[J]. 黑龍江農(nóng)業(yè)科學(xué), 2010(3): 38–40

[14] 苑小勇, 黃元仿, 高如泰, 等. 北京市平谷區(qū)農(nóng)用地土壤有機(jī)質(zhì)空間變異特征[J]. 農(nóng)業(yè)工程學(xué)報, 2008, 24(2): 70–76

[15] 鄭海龍, 陳杰, 鄧文靖, 等. 城市邊緣帶土壤重金屬空間變異及其污染評價[J]. 土壤學(xué)報, 2006, 43(1): 39–45

[16] McDowell R W, Sharpley A N. Approximating phosphorus release from soil to surface drainage[J]. Journal of Environmental Quality[J]. 2001, 30: 508–520

[17] 魯如坤, 時正元, 施建平. 我國南方6省農(nóng)田養(yǎng)分平衡現(xiàn)狀評價和動態(tài)變化研究[J]. 中國農(nóng)業(yè)科學(xué), 2000, 33(2): 63–67

[18] 魯如坤. 土壤磷素水平和水體環(huán)境保護(hù)[J]. 磷肥與復(fù)肥, 2003, 18(1): 4–8

[19] 金明清, 彭月月, 王佩. 四川省鹽源縣植煙土壤氮磷鉀空間變異特征及影響因素[J]. 土壤, 2016, 48(5): 984-991

[20] 吳曉妮, 付登高, 劉興祝, 等. 柴河流域典型景觀類型土壤氮磷含量的空間變異特征[J]. 土壤, 2016, 48(6): 1210–1215

[21] 王艷麗. 豫中褐土區(qū)磷肥用量對烤煙生長發(fā)育和產(chǎn)質(zhì)量的影響[D]. 鄭州: 河南農(nóng)業(yè)大學(xué), 2005

[22] 章新, 李明, 楊碩媛, 等. 磷素水平對煙葉化學(xué)成分和感官評吸質(zhì)量的影響[J]. 安徽農(nóng)業(yè)科學(xué), 2010, 38(10): 5091–5093, 5109

[23] 武德傳, 羅紅香, 宋澤民, 等. 黔南山地植煙土壤主要養(yǎng)分空間變異和管理分區(qū)[J]. 應(yīng)用生態(tài)學(xué)報, 2014, 25(6): 1701–1707

[24] 李強(qiáng), 周冀衡, 楊榮生, 等. 曲靖植煙土壤養(yǎng)分空間變異及土壤肥力適宜性評價[J]. 應(yīng)用生態(tài)學(xué)報, 2011, 22(4): 950–956

[25] 高博超, 婁翼來, 金廣遠(yuǎn), 等. 基于GIS和地統(tǒng)計學(xué)的植煙土壤養(yǎng)分空間分析[J]. 中國煙草學(xué)報, 2009, 15(1): 35–38

[26] 黃瑾, 林北森, 周文亮, 等. 廣西百色植煙土壤主要養(yǎng)分特征及施肥策略[J]. 中國煙草科學(xué), 2010, 31(4): 33–38

[27] 楊美仙, 亞平, 何元勝, 等. 臨滄植煙土壤養(yǎng)分豐缺狀況分析[J]. 云南農(nóng)業(yè)大學(xué)學(xué)報, 2014, 29(2): 269–277

[28] 江厚龍, 張保全, 李鈉鉀, 等. 重慶植煙土壤有效養(yǎng)分含量及其影響因素[J]. 中國煙草科學(xué), 2014, 35(5): 67–73

[29] 宋淑芳, 周冀衡, 鄧小華, 等. 大理植煙土壤養(yǎng)分含量及其對煙葉生產(chǎn)的適宜性[J]. 湖南農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2012, 38(1): 16–21

[30] 李洪斌, 張楊珠, 黃運(yùn)湘, 等. 不同施肥條件對煙草施肥效應(yīng)及肥料利用率的影響[J]. 湖南農(nóng)業(yè)科學(xué), 2013, 29(13): 48–52

[31] 劉建玲, 李仁崗, 廖文華, 等. 白菜–辣椒輪作中磷肥的產(chǎn)量效應(yīng)及土壤磷積累研究[J]. 中國農(nóng)業(yè)科學(xué), 2005, 38(8): 1616–1620

[32] 李中陽, 徐明崗, 李菊梅, 等. 長期施用化肥有機(jī)肥下我國典型土壤無機(jī)磷的變化特征[J]. 土壤通報, 2010, 41(6): 1434–1439

[33] 王艷玲, 何園球, 李成亮, 等. 長期施肥對紅壤磷素持續(xù)供應(yīng)能力的影響[J]. 土壤學(xué)報, 2010, 47(3): 503–507

[34] 顏曉, 王德建, 張剛, 等. 長期施磷稻田土壤磷素累積及其潛在環(huán)境風(fēng)險[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2013, 21(4): 393–400

[35] 周冀衡, 朱小平, 王彥亭, 等. 煙草生理與生物化學(xué)[M]. 合肥: 中國科學(xué)技術(shù)大學(xué)出版社, 1996: 192–193

[36] 柏兆海, 萬其宇, 李海港, 等. 縣域農(nóng)田土壤磷素積累及淋失風(fēng)險分析——以北京市平谷區(qū)為例[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2011, 30(9): 1853–1860

[37] 王海波, 胡軍, 李致新, 等. 三門峽市主要烤煙產(chǎn)區(qū)土壤養(yǎng)分含量分析[J]. 江西農(nóng)業(yè)學(xué)報, 2010, 22(1): 33–36

Temporal and Spatial Variability Risk Assessment of Soil Phosphorous in Karst Tobacco-planting Region of Xiangxi, Hunan

LI Qiang1,2,DAI Meilin2, XIANG Deming3, LI Juan1*, ZHOU Miliang3, TIAN Feng3,TIAN Minghui3, ZHANG Liming1,3

(1 College of Agronomy, Hunan Agricultural University, Changsha 410128, China; 2 Tobacco Research Institute, Hunan Agricultural University, Changsha 410128, China; 3 Xiangxi Branch of Hunan Provincial Tobacco Company, Jishou, Hunan 416000, China)

Based on methods of spatial autocorrelation, semi-variance function and the fractal dimension as well as the data of soil available phosphorus in 2000 and 2014, the spatiotemporal variation and environmental risks of soil available phosphorus in Karst tobacco-planting region in Xiangxi of Hunan were studied. The results showed that, during the 14 years, soil available phosphorus content increased by 30.89 mg/kg; thescores of soil available phosphorus index decreased, fractal dimension and the nugget effect increased, indicating a weaker spatial structure and stronger stochastic variability; the areas of ‘extremely high’, ‘very high’ and ‘high’ levels were increased by 28.36%, 26.89% and 25.94%, while the areas of ‘low’ and ‘very low’ levels decreased by 57.45% and 27.12%, respectively; There were significant correlation between soil available phosphorus and elevation, soil parent material, land use type, rotation mode and SOM content; Probability of soil available phosphorus content exceed the risk critical value increased greatly, the total areas of very high risk (probability >80%), high risk (probability 60%–80%) and medium risk (probability 40%–60%) accounted for 1.77%, 5.40% and 13.96% of the total study area, respectively. North of Longshan and east of Yongshun are the key regions to control agricultural non-point pollution of phosphorus.

Karst area; Tobacco-planting soil; Available phosphorus; Spatiotemporal variation; Environmental risk

湖南省教育廳項目(CX2010B309)，湖南省煙草公司項目(2016-02)和湘西烤煙良好農(nóng)業(yè)操作規(guī)范及推廣應(yīng)用研究項目(州煙司2012-37)資助。

(adalee619@163.com)

李強(qiáng)(1982—)，男，江蘇儀征人，博士，講師，主要從事煙草生理生化方面的研究。E-mail: zqiangli@126.com

10.13758/j.cnki.tr.2018.01.024

S572.06

A