王 哲, 夏 輝, 袁 浩, 谷書(shū)杰
(河北農(nóng)業(yè)大學(xué) 城鄉(xiāng)建設(shè)學(xué)院, 河北 保定 071000)
農(nóng)業(yè)面源污染是指在農(nóng)業(yè)生產(chǎn)活動(dòng)中,氮素和磷素等營(yíng)養(yǎng)物質(zhì)、農(nóng)藥以及其他有機(jī)或無(wú)機(jī)污染物質(zhì),通過(guò)農(nóng)田的地表徑流和農(nóng)田滲漏,造成的水環(huán)境的污染[1]。《第一次全國(guó)污染源普查公報(bào)》顯示,我國(guó)農(nóng)業(yè)源已成為污染源之首,是氮磷排放的主要來(lái)源,其排放量分別為270.46萬(wàn)t,28.47萬(wàn)t,分別占排放總量的57.2%,67.4%,其中種植業(yè)是主要部分[2-4]。許多農(nóng)民為追求高產(chǎn)大量施用化肥,不僅造成肥料的浪費(fèi),還加大氮磷淋失的風(fēng)險(xiǎn),引發(fā)地下水硝酸鹽污染、水體富營(yíng)養(yǎng)化等一系列負(fù)面環(huán)境問(wèn)題[5-6]。有研究表明我國(guó)主要糧食作物的氮肥和磷肥的平均利用率分別為27.56%,11.6%[7],氮磷施用量需在一定范圍內(nèi)才能促進(jìn)植物對(duì)氮磷養(yǎng)分的吸收[8]。
2015年中央1號(hào)文件強(qiáng)調(diào)要“加強(qiáng)農(nóng)業(yè)生態(tài)治理”和農(nóng)業(yè)生產(chǎn)面源污染的防控。河北省作為我國(guó)的一個(gè)農(nóng)業(yè)大省,農(nóng)業(yè)面源污染問(wèn)題一直存在,公開(kāi)的數(shù)據(jù)顯示,2017年河北省農(nóng)用化肥使用量居于全國(guó)農(nóng)用化肥使用量的第3位。保定市是京津冀地區(qū)中心城市之一,隨著雄安新區(qū)的建立,保定市的戰(zhàn)略地位越來(lái)越重要,該區(qū)域的生態(tài)好壞將深刻影響雄安新區(qū)的未來(lái)。因此,本研究以保定東部平原農(nóng)業(yè)生產(chǎn)區(qū)為研究對(duì)象,探索農(nóng)民傳統(tǒng)田間管理下夏玉米種植農(nóng)田2 m土層范圍內(nèi)土壤氮磷含量變化情況,分析該區(qū)域農(nóng)業(yè)面源污染時(shí)空變化特征,以期為河北省中西部平原區(qū)農(nóng)業(yè)面源污染風(fēng)險(xiǎn)管控和生態(tài)環(huán)境質(zhì)量提升提供科學(xué)依據(jù),同時(shí)為該區(qū)域進(jìn)一步提出合理施用化肥建議提供科學(xué)指導(dǎo)。
1.1.1 研究區(qū)概況 研究區(qū)位于保定東部淀西平原區(qū),北緯38°19′—39°13′,東經(jīng)114°44′—116°15′,按流域劃分屬于白洋淀以西的淀西清南區(qū),東起潴龍河西至京廣鐵路,北起南拒馬河南至沙河,主體位于河北省保定市東部,面積約為3 952 km2,屬暖溫帶大陸性季風(fēng)氣候,四季分明。多年平均氣溫13.4℃;年平均風(fēng)速1.8 m/s;年平均降水量498.9 mm,降水集中在6—8月,7月最多;年平均蒸發(fā)量為1 430.5 mm。
1.1.2 試驗(yàn)場(chǎng)概況 試驗(yàn)場(chǎng)位于河北省保定市安新縣安州鎮(zhèn)白莊村(38°85′N(xiāo),115°47′E),地處保定東部淀西平原區(qū),位于研究區(qū)內(nèi),屬暖溫帶半濕潤(rùn)大陸季風(fēng)氣候,多年平均氣溫為12.1℃,平均風(fēng)速2.5 m/s,年平均降水529.5 mm,年內(nèi)分配不均,主要集中在7—9月份,占全年總降水量的67.9%以上。多年平均蒸發(fā)量1 773.4 mm。土壤類(lèi)型為潮褐土,0—140 cm土壤質(zhì)地為粉質(zhì)黏土,140—200 cm為粉土。0—2 m土壤基本理化性質(zhì)見(jiàn)表1。

表1 試驗(yàn)場(chǎng)土壤基本理化性質(zhì)
為了能夠更加全面了解保定東部平原區(qū)內(nèi)夏玉米種植期農(nóng)業(yè)面源污染情況,本研究分為兩部分,第一部分為定位試驗(yàn)場(chǎng)取土測(cè)定試驗(yàn),在夏玉米生育期進(jìn)行取土,研究土壤氮磷在垂直方向的分布特性,第二部分為研究區(qū)內(nèi)取土測(cè)定試驗(yàn),在夏玉米生育期進(jìn)行取土,探索土壤氮磷在空間上的分布特性。兩部分試驗(yàn)夏玉米播種、灌溉、施肥等田間管理均按農(nóng)民傳統(tǒng)習(xí)慣進(jìn)行,夏玉米播種和收獲時(shí)間為2016年6月14日和10月1日,夏玉米生育期灌2次水,播種后灌蒙頭水70 mm,抽雄期配合施肥灌第二次水60 mm,均采用當(dāng)?shù)剞r(nóng)民傳統(tǒng)地面灌溉方式,夏玉米播種前底施復(fù)合肥,抽雄期追施復(fù)合肥,玉米季施用氮肥155.8 kg/hm2,磷肥62.3 kg/hm2,鉀肥62.3 kg/hm2。
1.3.1 研究區(qū)取土點(diǎn)及取土?xí)r間的選擇 研究區(qū)土壤樣品的采集應(yīng)能夠如實(shí)反映研究區(qū)土壤氮磷的分布情況及變化規(guī)律。因此,取土點(diǎn)的選擇主要考慮以下因素:(1) 具有代表性。包括其土壤質(zhì)地、種植作物及施肥方式等,農(nóng)民田間管理均按傳統(tǒng)習(xí)慣進(jìn)行;(2) 分布均勻。盡可能將取土點(diǎn)較為均勻地分布在研究區(qū)域內(nèi)。基于以上原則,最終確定10個(gè)取土點(diǎn),取土點(diǎn)分布在高陽(yáng)縣、清苑區(qū)、蠡縣、望都縣、定州市、安國(guó)縣、雄縣、容城縣、徐水區(qū)和安新縣10個(gè)縣。研究區(qū)取土試驗(yàn)取土點(diǎn)基本信息見(jiàn)表2。
由于土壤中氮磷含量受降雨、灌溉、施肥、作物吸收、揮發(fā)等因素影響,為使各取樣點(diǎn)的數(shù)據(jù)具有同步性,研究區(qū)內(nèi)土壤采集定在在2016年9月4—6日集中完成。取土?xí)r間的確定主要考慮以下因素:(1) 該時(shí)間段無(wú)降雨;(2) 9月正處于夏玉米生長(zhǎng)期,這個(gè)時(shí)期一般不會(huì)進(jìn)行灌溉;(3) 6—9月為研究區(qū)的主汛期,降雨較為集中,9月后降雨逐漸減少。經(jīng)過(guò)一個(gè)汛期后,土壤中的氮磷物質(zhì)隨著水分下滲將發(fā)生淋洗現(xiàn)象。因此,該時(shí)間段取土更加能夠反映出研究區(qū)面源污染物進(jìn)入土壤甚至淺層地下水的過(guò)程。

表2 取土點(diǎn)基本信息
1.3.2 土壤樣品采集與處理 在試驗(yàn)場(chǎng)夏玉米的拔節(jié)期(2016年7月3日)、小喇叭口期(2016年7月15日)、抽雄期(2016年7月31日)、開(kāi)花期(2016年8月14日)、灌漿期(2016年9月10日)以及研究區(qū)夏玉米的灌漿期(2016年9月4日—6日)的10個(gè)取樣點(diǎn)進(jìn)行土壤樣品采集。取樣方式為土鉆打孔取土,分層取樣,距地表1 m范圍每隔10 cm取樣1次,距地表1—2 m范圍每隔20 cm取樣1次,每次取土垂直方向上共取14次,取樣完畢立即密封并冷凍保存,帶回實(shí)驗(yàn)室測(cè)定。對(duì)夏玉米種植、灌溉、施肥等田間活動(dòng)進(jìn)行監(jiān)測(cè)并記錄;試驗(yàn)場(chǎng)設(shè)置小型氣象站記錄日降水量;土壤硝態(tài)氮測(cè)定采用紫外分光光度法;土壤銨態(tài)氮測(cè)定采用靛酚藍(lán)比色法;土壤有效磷測(cè)定采用鉬銻抗比色法。
利用Excel與SPSS 17.0軟件進(jìn)行統(tǒng)計(jì)分析;以緯度為x軸,經(jīng)度為y軸,土壤指標(biāo)含量為z軸,利用Surfer軟件繪制等值線(xiàn)圖。由于取樣在垂直方向上間隔距離不一樣,在分析數(shù)據(jù)時(shí)采用加權(quán)平均值,文中平均值均指加權(quán)平均值。
表3為試驗(yàn)場(chǎng)夏玉米種植期垂直方向上各土層土壤硝態(tài)氮含量狀況,由數(shù)據(jù)可看出夏玉米種植期土壤硝態(tài)氮含量隨土層深度的增加呈現(xiàn)先減小再增加最后減小的趨勢(shì),大致呈“S”型變化,硝態(tài)氮含量在30 cm土層處出現(xiàn)波谷,最小值為2.66 mg/kg,出現(xiàn)在灌漿期,80 cm土層處出現(xiàn)波峰,最大值為40.25 mg/kg,出現(xiàn)在抽雄期,此土層夏玉米各生育期硝態(tài)氮含量均最大,同時(shí)該層含水量也最大,土壤孔隙小,土壤黏粒增加,吸附作用大,對(duì)硝態(tài)氮向深層土壤滲漏有一定阻礙作用,故儲(chǔ)存在該層的硝態(tài)氮較多,易形成累積,此外7月、8月份降雨充沛,遇較大降雨此土層硝態(tài)氮易隨水分向土壤深處遷移淋溶,增加了該區(qū)域地下水污染的風(fēng)險(xiǎn)。7月31日2 m內(nèi)土層的硝態(tài)氮含量基本都高于其他同等深度土壤硝態(tài)氮含量是因?yàn)?月29日施用復(fù)合肥且7月28日與30日降雨造成的,表層土壤硝態(tài)氮隨水分向土壤深處移動(dòng)。80 cm之后土層硝態(tài)氮含量隨土層深度的增加而逐漸減小,最后在2 m土層深度穩(wěn)定在5 mg/kg左右。夏玉米全生育期表土層硝態(tài)氮平均含量為12.15 mg/kg、心土層為24.82 mg/kg、底土層為7.73 mg/kg,底土層與表土層硝態(tài)氮含量相差4.42 mg/kg,兩層硝態(tài)氮含量較為接近,處于同一水平,但心土層硝態(tài)氮含量明顯大于其他兩層,將近表土層和底土層的兩倍多,因此,在夏玉米種植期間心土層硝態(tài)氮含量大,累積多,淋溶比例大。

表3 試驗(yàn)場(chǎng)夏玉米種植期間各土層硝態(tài)氮含量狀況
圖1為研究區(qū)內(nèi)10個(gè)取樣點(diǎn)表土層、心土層和底土層空間上土壤硝態(tài)氮含量變化情況。在0—2 m土壤范圍內(nèi),位于研究區(qū)中部的2號(hào)取樣點(diǎn)土壤硝態(tài)氮含量最高,平均含量為36.96 mg/kg,1號(hào)取樣點(diǎn)、8號(hào)取樣點(diǎn)、4號(hào)取樣點(diǎn)和3號(hào)取樣點(diǎn)土壤硝態(tài)氮含量較低,平均含量均小于5 mg/kg。由圖1可看出,研究區(qū)表土層硝態(tài)氮含量總體趨勢(shì)是東部低西部高,南部低北部高,且變化幅度較大,極差為39.28 mg/kg;心土層和底土層內(nèi)硝態(tài)氮含量在研究區(qū)中部2號(hào)取樣點(diǎn)均出現(xiàn)了局部高點(diǎn),硝態(tài)氮污染風(fēng)險(xiǎn)大,遇較大灌水或降雨都將會(huì)導(dǎo)致硝態(tài)氮向深層土壤淋洗甚至進(jìn)入地下水,底土層硝態(tài)氮含量從東到西逐漸增加。通過(guò)對(duì)10個(gè)取樣點(diǎn)各土層硝態(tài)氮含量進(jìn)行統(tǒng)計(jì)分析可知,研究區(qū)內(nèi)距離地表60—70 cm處土壤硝態(tài)氮含量較為穩(wěn)定,距地表90 cm處土壤硝態(tài)氮含量變幅最大,0—2 m土壤中硝態(tài)氮含量最大值出現(xiàn)在距地表20 cm位置,到70 cm處降到最小值后再有所回升。

圖1 研究區(qū)夏玉米種植期土壤硝態(tài)氮含量變化情況
表4為試驗(yàn)場(chǎng)夏玉米種植期垂直方向上各土層土壤銨態(tài)氮含量狀況,夏玉米種植期土壤銨態(tài)氮變化與硝態(tài)氮不同,銨態(tài)氮不易隨水分運(yùn)動(dòng)并能夠被作物直接吸收利用,由數(shù)據(jù)可看出土壤銨態(tài)氮含量在垂直方向上變化不大,分布均勻,7月3日拔節(jié)期20 cm處銨態(tài)氮含量達(dá)到2.5 mg/kg,主要是由于冬小麥?zhǔn)斋@夏玉米播種后,土壤波動(dòng)較大,受施入底肥、灌水與土壤翻動(dòng)等各方面的影響,表層土壤土質(zhì)疏松,蒸發(fā)量大,在田間,氨揮發(fā)隨風(fēng)速和光照強(qiáng)度的增大而增多,因此,30 cm土層銨態(tài)氮含量就恢復(fù)了正常值,除20 cm土層外,其他土層銨態(tài)氮含量基本穩(wěn)定在0.6 mg/kg左右,7月31日20 cm土層處銨態(tài)氮含量高達(dá)3.6 mg/kg是由于7月29日小喇叭口期施用復(fù)合肥所造成的。夏玉米種植期表土層銨態(tài)氮平均含量為0.90 mg/kg、心土層為0.59 mg/kg、底土層為0.51 mg/kg,心土層銨態(tài)氮含量與表土層銨態(tài)氮含量相差0.31 mg/kg,降幅34.44%,底土層硝態(tài)氮含量與表土層硝態(tài)氮含量相差0.39 mg/kg,降幅43.33%。
圖2為研究區(qū)內(nèi)10個(gè)取樣點(diǎn)表土層、心土層和底土層空間上土壤銨態(tài)氮含量變化情況。在0—2 m土壤范圍內(nèi),研究區(qū)內(nèi)1號(hào)取樣點(diǎn)土壤銨態(tài)氮含量最高,平均含量為3.47 mg/kg,6號(hào)取樣點(diǎn)、2號(hào)取樣點(diǎn)和5號(hào)取樣點(diǎn)土壤銨態(tài)氮含量較低,平均含量均小于1 mg/kg。由圖2可看出,研究區(qū)心土層與底土層銨態(tài)氮含量均呈現(xiàn)東部高西部低,南部高北部低的趨勢(shì),兩者變化幅度較小,極差分別為2.85,2.40 mg/kg,而表土層銨態(tài)氮含量則呈波浪形,變化較為明顯,極差為3.40 mg/kg。

表4 試驗(yàn)場(chǎng)夏玉米種植期間各土層銨態(tài)氮含量狀況
表5為試驗(yàn)場(chǎng)夏玉米種植期垂直方向上各土層土壤有效磷含量狀況,由數(shù)據(jù)可看出夏玉米生育期土壤有效磷隨土層深度的增加先減小再增加最后減小,但與土壤硝態(tài)氮變化不同,土壤有效磷增加幅度較小。表土層有效磷含量降幅最大,同時(shí)20 cm土層出現(xiàn)夏玉米種植期間有效磷的最大值16.71 mg/kg,處于開(kāi)花期。土壤有效磷含量在60 cm土層附近出現(xiàn)波動(dòng),主要由于此層腐殖質(zhì)含量較高,團(tuán)粒結(jié)構(gòu)豐富,具有較高的保水保肥能力,富集了部分有效磷,有效磷含量出現(xiàn)短暫增高。夏玉米種植期表土層有效磷平均含量為7.02 mg/kg、心土層為3.44 mg/kg、底土層為2.52 mg/kg,心土層與底土層土壤有效磷含量變化幅度不大,土壤有效磷在20 cm土層保持較高水平,主要有兩點(diǎn):一是與農(nóng)民田間管理采取的措施有關(guān)。農(nóng)民一般施肥僅僅局限在土壤表層,有效磷在土壤表層富集;二是有效磷本身性質(zhì)造成的。有效磷的移動(dòng)性較小,在垂向上只能以水溶態(tài)的形式發(fā)生小范圍的遷移。140—200 cm土層有效磷含量變化不穩(wěn)定,主要由于試驗(yàn)場(chǎng)潛水埋藏深度較淺,7月、8月份降雨充沛,地下水位上升,土壤與地下水之間存在物質(zhì)交換,深層土壤有效磷含量受潛水的影響較大。

圖2 研究區(qū)夏玉米種植期土壤銨態(tài)氮含量變化情況

表5 試驗(yàn)場(chǎng)夏玉米種植期間各土層有效磷含量狀況
根據(jù)全國(guó)第二次土壤普查養(yǎng)分分級(jí)標(biāo)準(zhǔn),試驗(yàn)場(chǎng)夏玉米種植期間表土層有效磷為4級(jí),心土層有效磷為5級(jí),底土層有效磷為6級(jí);研究區(qū)夏玉米種植期間表土層有效磷為3級(jí),心土層有效磷為5級(jí),底土層有效磷為4級(jí)。圖3為研究區(qū)內(nèi)10個(gè)取樣點(diǎn)表土層、心土層和底土層空間上土壤有效磷含量變化情況。在0—2 m土壤范圍內(nèi),研究區(qū)內(nèi)3號(hào)取樣點(diǎn)土壤有效磷含量最高,平均含量為17.90 mg/kg,研究區(qū)內(nèi)6號(hào)、8號(hào)和9號(hào)取樣點(diǎn)土壤有效磷含量較低,平均含量在5.0 mg/kg左右。由圖3可看出,表土層有效磷含量在整個(gè)研究區(qū)變化較為明顯,極差為32.66 mg/kg,由南到北表土層有效磷含量大致呈增大趨勢(shì),表土層有效磷易隨水分遷移,形成地表徑流進(jìn)入水體造成水體富營(yíng)養(yǎng)化等,心土層有效磷呈東部高西部低,南部低北部高的趨勢(shì),底土層有效磷呈東部高西部低,南部高北部低的趨勢(shì),心土層和底土層有效磷含量變化幅度較小,極差分別為4.92,8.69 mg/kg。

圖3 研究區(qū)夏玉米種植期土壤有效磷含量變化情況
氮磷的淋失不僅直接影響農(nóng)業(yè)生產(chǎn)成本,而且所導(dǎo)致的農(nóng)業(yè)面源污染成為一個(gè)世界性的難題[9]。因灌溉不當(dāng)或夏季集中降雨造成的土壤硝態(tài)氮淋失是我國(guó)華北平原氮損失的重要途徑之一[10]。本研究通過(guò)田間取土測(cè)定試驗(yàn),針對(duì)試驗(yàn)場(chǎng)夏玉米全生育期農(nóng)田土壤氮磷的含量和分布規(guī)律進(jìn)行了垂直方向上的研究分析,并對(duì)保定東部平原區(qū)夏玉米種植農(nóng)田土壤氮磷的含量及變化進(jìn)行了空間上的探索。通過(guò)對(duì)比試驗(yàn)場(chǎng)與研究區(qū)夏玉米灌漿期土壤氮磷在表土層、心土層及底土層的分布特性,以期獲得試驗(yàn)場(chǎng)與研究區(qū)夏玉米灌漿期土壤氮磷含量及分布的關(guān)系,從而更全面地了解農(nóng)業(yè)面源污染狀況。研究表明,試驗(yàn)場(chǎng)垂直方向上夏玉米灌漿期土壤硝態(tài)氮呈先增大后減小的趨勢(shì),銨態(tài)氮及有效磷均呈減小的趨勢(shì),而研究區(qū)垂直方向上土壤氮磷都呈現(xiàn)先減小后增加的趨勢(shì),在含量上夏玉米灌漿期研究區(qū)表土層、底土層比試驗(yàn)場(chǎng)分別高8.65,3.77 mg/kg,但是研究區(qū)心土層硝態(tài)氮含量比試驗(yàn)場(chǎng)低8.66 mg/kg,試驗(yàn)場(chǎng)硝態(tài)氮在心土層累積多;研究區(qū)銨態(tài)氮含量比試驗(yàn)場(chǎng)分別高1.15,0.70,1.12 mg/kg;研究區(qū)表土層、心土層及底土層有效磷含量比試驗(yàn)場(chǎng)分別高14.04,0.18,2.86 mg/kg。綜上可知研究區(qū)土壤硝態(tài)氮含量除心土層以外,均比試驗(yàn)場(chǎng)對(duì)應(yīng)土層含量高,研究區(qū)土壤銨態(tài)氮、有效磷含量除表土層有效磷含量以外,各土層含量較試驗(yàn)場(chǎng)對(duì)應(yīng)土層含量變化幅度小,但是研究區(qū)內(nèi)土壤氮磷與試驗(yàn)場(chǎng)內(nèi)變化規(guī)律不盡一致,主要是由于農(nóng)業(yè)面源污染具有多源性、隨機(jī)性、分散性、分布廣、間歇性和難以監(jiān)測(cè)性的特點(diǎn)[3]。饒靜等[11]從宏觀、中觀和微觀3個(gè)層面分析了我國(guó)農(nóng)業(yè)面源污染發(fā)生機(jī)制;丘雯文[12]、楊濱鍵[13]等的研究表明我國(guó)農(nóng)業(yè)面源污染具有明顯的空間差異,東部和中部地區(qū)的排放強(qiáng)度較高,西部和東北地區(qū)則相對(duì)較低,農(nóng)業(yè)面源污染排放總體差異表現(xiàn)為“下降—上升—下降”的波動(dòng)變化趨勢(shì),農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)、生產(chǎn)規(guī)模的不同也使農(nóng)業(yè)面源污染表現(xiàn)出明顯的地域差異性。針對(duì)農(nóng)業(yè)面源污染造成的污染,目前主要治理措施有源頭控制、遷移過(guò)程攔截、末端治理等[14-15]。土壤中氮磷的大量輸出是引起農(nóng)業(yè)面源污染的主要原因,要防治農(nóng)業(yè)面源污染最關(guān)鍵的問(wèn)題是對(duì)氮磷的來(lái)源加以控制,因此,在研究區(qū)內(nèi)根據(jù)各土層的氮磷含量合理施肥將會(huì)減少農(nóng)業(yè)面源污染,對(duì)保護(hù)環(huán)境具有重要意義。但是農(nóng)業(yè)面源污染的機(jī)制、影響因素以及防治措施都不唯一,并且具有不確定性,對(duì)此還有待進(jìn)一步研究。
(1) 土壤硝態(tài)氮在垂直方向上隨土層深度的增加呈“S”型變化,心土層硝態(tài)氮累積多,淋溶比例大,降雨與灌溉將會(huì)導(dǎo)致此層硝態(tài)氮向深層淋洗。空間上研究區(qū)表土層硝態(tài)氮東部低西部高,南部低北部高,變化幅度較大,心土層和底土層硝態(tài)氮含量在研究區(qū)中部出現(xiàn)了局部高點(diǎn),底土層硝態(tài)氮含量從東到西逐漸增加。
(2) 土壤銨態(tài)氮垂直方向上分布均勻,心土層和底土層銨態(tài)氮含量較為接近。空間上研究區(qū)心土層與底土層銨態(tài)氮均呈現(xiàn)東部高西部低,南部高北部低的趨勢(shì),兩者變化幅度小,表土層銨態(tài)氮?jiǎng)t呈波浪形,變化幅度較為明顯。
(3) 土壤有效磷在垂直方向上隨土層深度的增加先減小再增加最后減小,主要受施肥和地下水埋深的影響。空間上表土層有效磷在整個(gè)研究區(qū)變化較為明顯,由南到北表土層有效磷含量大致呈現(xiàn)增大的趨勢(shì),心土層有效磷呈現(xiàn)東部高西部低,南部低北部高的趨勢(shì),底土層有效磷呈現(xiàn)東部高西部低,南部高北部低的趨勢(shì),心土層和底土層有效磷變化幅度較小。