基于1∶5萬土壤數(shù)據(jù)庫的太湖地區(qū)水稻土全氮含量動態(tài)變化研究*

袁 平 張黎明? 喬 婷 謝安乾 于東升 史學正 邢世和 陳翰閱

（1 福建農林大學資源與環(huán)境學院，福州 350002）

（2 土壤生態(tài)系統(tǒng)健康與調控福建省高校重點實驗室，福州 350002）

（3 土壤與農業(yè)可持續(xù)發(fā)展國家重點實驗室（中國科學院南京土壤研究所），南京 210008）

土壤氮庫（N）是陸地生態(tài)系統(tǒng)養(yǎng)分庫的重要組成成分，它與有機碳庫相輔相成，在維持生態(tài)系統(tǒng)的組成、結構和功能中發(fā)揮著重要作用[1]。土壤中氮素由于植物生長、微生物活動、人為干擾及氣候條件改變而發(fā)生轉化，氮素的氣態(tài)損失和淋溶損失也嚴重影響著生態(tài)環(huán)境。有研究表明，土壤氮素在反硝化作用下產生的各種含氮氣體如一氧化二氮（N2O）、一氧化氮（NO）、二氧化氮（NO2）等直接參與溫室效應[2]。此外，氮也是水體富營養(yǎng)化的重要伴隨因子，過量的氮素流入池塘、湖泊、河流和海洋等會造成藻類種群大爆發(fā)[3]。與自然土壤相比，農業(yè)土壤氮庫受到人類活動的影響更大，其可以在較短的時間內發(fā)生轉化和遷移。因此，明確農業(yè)土壤氮庫的動態(tài)變化對研究陸地生態(tài)系統(tǒng)氮循環(huán)，以及提高氮素利用率等問題的解決均具有重要的意義[4]。

當前，基于數(shù)字化土壤圖結合土壤普查資料進行氮素研究的方法越來越廣泛。從全球尺度來看，Batjes[5]基于1∶500萬的世界土壤圖和4 353個土壤剖面數(shù)據(jù)對全球氮儲量進行了估算，結果表明，全球1 m 土層土壤全氮儲量為133～140 Pg（1 Pg= 1015g）。從國家尺度來看，Tian 等[6]基于1∶100 萬全國土壤空間分布圖和2 480 個土壤剖面對我國氮庫儲量進行了研究，結果表明，全國土壤氮儲量為8.29 Pg。Lin 等[7]利用1∶100 萬中國土壤圖和全國第二次土壤普查1 490 個水稻土剖面數(shù)據(jù)對我國水稻土全氮密度和儲量進行了研究，結果表明，水稻土剖面1 m 土層平均氮密度為12.4 t·hm-2，全氮儲量為569 Tg。從區(qū)域尺度來看，焦閃閃等[8]基于福建省第二次土壤普查的2 021 個采樣點和2008年農業(yè)部測土配方施肥29 945 個耕層樣點對全省耕地全氮密度和儲量動態(tài)變化進行了研究，結果表明，1982—2008年該地區(qū)耕地土壤全氮密度和儲量均為上升趨勢。代子俊等[9]基于第二次土壤普查的45 個剖面資料和2015年土壤調查的61 個剖面數(shù)據(jù)研究了湟水流域土壤全氮的空間變異，結果表明，近30年來湟水流域土壤全氮呈現(xiàn)增長趨勢。

但是從以上研究來看，我國國家和區(qū)域尺度上的氮素研究多以中小比例尺數(shù)據(jù)庫靜態(tài)研究為主，而結合高精度土壤數(shù)據(jù)庫的兩期大樣本實測樣點動態(tài)變化研究相對較少。水稻土作為我國農業(yè)土壤的重要組成部分，耕作歷史超過10 000年，面積達30 M hm2，水稻產量約占我國糧食總產量的1/2 以上[10-11]。為了獲得作物高產，農民往往加大氮肥的投入，這也導致我國水稻田面臨著高產和生態(tài)環(huán)境保護之間矛盾的多重挑戰(zhàn)[3]。太湖地區(qū)作為我國主要水稻產區(qū)，該地區(qū)高產稻田氮肥施用量已經達到360 kg·hm-2，氮肥的大量施用和高溫多雨的特點導致土壤氮素的大量淋失，造成了該流域頻繁的藍藻爆發(fā)與水體富營養(yǎng)化等環(huán)境問題[3-4]?；诖耍狙芯窟x擇太湖地區(qū)37 個縣（市、區(qū)）2.32 M hm2水稻土為研究對象，以1982年第二次土壤普查和2000年土壤質量演變與持續(xù)利用“973”項目采集的水稻土樣點資料所建立1∶5 萬大比例尺土壤數(shù)據(jù)庫為基礎，揭示水稻土全氮含量的時空動態(tài)變化規(guī)律，明確整個地區(qū)和不同土壤類型、土區(qū)及行政區(qū)的全氮含量動態(tài)變化差異的主控因子，研究結果可為合理制定太湖地區(qū)和我國南方水稻田的農業(yè)管理政策提供參考依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

太 湖 地 區(qū)（118°50′～121°54′E， 29°56′～32°16′N）包括江蘇省鎮(zhèn)江、無錫、常州和蘇州，浙江省嘉興、湖州和杭州的一部分，以及整個上海市，共37 個縣（市、區(qū)），流域總面積達3.65 M hm2（圖1）[10]。氣候類型為亞熱帶季風氣候，溫和濕潤，光照充足，生長季節(jié)較長，干濕季明顯，年平均氣溫在16℃左右，年均降水量為1 100～1 400 mm，≥10℃的積溫在5 000℃左右，日照時數(shù)為1 870～2 225 h，全年無霜期230 d 以上[12]。

經過長期的耕作施肥，太湖地區(qū)形成了以水稻土為主的土壤類型，面積達到2.32 M hm2，占整個研究區(qū)土壤總面積的66%[12]。該地區(qū)水稻土可劃分為6 個亞類（潴育型水稻土、淹育型水稻土、漂洗型水稻土、滲育型水稻土、潛育型水稻土和脫潛型水稻土）和4 個土區(qū)（低山丘陵土區(qū)、太湖平原土區(qū)、沖積平原土區(qū)和低洼圩田土區(qū)）（圖1）。成土 母質為黃土、沖積土和湖泊沉積物等為主。作物管理方式以水稻和冬小麥輪作為主，水稻于6月種植，10月收獲，小麥于11月種植，次年5月收割[11]。

圖1 太湖地區(qū)地理位置和土區(qū)分布圖 Fig.1 Geographic location and paddy soil region map of the Tai-Lake Region，China

1.2 數(shù)據(jù)來源

（1）土壤數(shù)據(jù)。1∶5 萬土壤空間數(shù)據(jù)庫所需基本圖鑒來自于全國第二次土壤普查成果中太湖地區(qū)37 個縣（市、區(qū)）1∶5 萬土壤圖，在ArcGIS 支持下采用雙標準緯線等積圓錐投影，對各個縣（市、區(qū)）1∶5 萬紙質土壤圖進行數(shù)字化、編輯和修改建成，共計水稻土圖斑52 031 個。土壤圖的基本分類單元為土種，共計622 個土種，分別歸屬6 個亞類和137 個土屬，所有土壤圖的制圖單元均采用中國土壤發(fā)生分類系統(tǒng)（Genetic Soil Classification of China，GSCC）[13]。

土壤屬性數(shù)據(jù)庫包括水稻土樣點的位置描述、理化性質和肥力因子，如土壤類型名稱、土層深度、剖面地點、全氮、有機質、全磷和全鉀含量等，其中，土壤全氮采用重鉻酸鉀、硫酸硝化-蒸餾法測定（凱氏定氮法）。1982年表層樣點數(shù)據(jù)來自于太湖地區(qū)37 個縣（市、區(qū)）在《中國土種志》、《省級土種志》、《地市級土種志》和《縣級土種志》中所記錄的水稻土剖面資料，共計1 096 個（圖2）。2000年表層樣點數(shù)據(jù)來自于土壤質量演變與持續(xù)利用“973”項目，共計1 393 個（圖2）?？臻g數(shù)據(jù)和屬性數(shù)據(jù)采用 Shi 等[14]提出的“PKB（Pedological Knowledge Based Method）”法[15-17]連接。

圖2 1982年和2000年太湖地區(qū)水稻土樣點分布圖 Fig.2 Distributions of paddy soil sampling sites in the Tai-Lake Region in 1982 and 2000

（2）氣象和施肥數(shù)據(jù)。氣象數(shù)據(jù)來源于太湖地區(qū)13 個國家氣象站數(shù)據(jù)資料，包括1982—2000年逐日降雨量、最高和最低氣溫[18]。農作物管理數(shù)據(jù)來源于以縣（市、區(qū)）為統(tǒng)計單元的農村統(tǒng)計年鑒，包括1982—2000年37 個縣（市、區(qū)）水稻土年均氮肥和農家肥施用量[19]。氣象數(shù)據(jù)和施肥數(shù)據(jù)均以縣（市、區(qū)）為最小統(tǒng)計單元連接到水稻土數(shù)據(jù)庫中。

1.3 數(shù)據(jù)統(tǒng)計分析

不同土壤類型、土區(qū)和行政區(qū)面積加權平均全氮含量（STN，Soil Total Nitrogen Content，g·kg-1）、全氮含量變化量（VSTN，Variation of Soil Total Nitrogen，g·kg-1）和全氮含量變化幅度（VSTNC，Variability of Soil Total Nitrogen Content，%）計算公式如下：

式中，APS 為研究區(qū)水稻土總面積（hm2），APSi和 STNi分別表示土壤數(shù)據(jù)庫中每個圖斑的面積（hm2）和表層水稻土（0 ～15 cm）全氮含量（g·kg-1），n代表圖斑數(shù)。STNf代表水稻土1982年或2000年土壤全氮含量（g·kg-1），f代表年份（f=1982、2000）。

2 結果與討論

2.1 1982—2000年太湖地區(qū)水稻土全氮含量變化

在1∶5 萬土壤數(shù)據(jù)庫下太湖地區(qū)水稻土圖斑為52 031 個，總面積為2.32 M hm2，1982年和2000年表層水稻土（0～15 cm）全氮含量分別為1.58 g·kg-1和1.79 g·kg-1，近20年的增幅為13.3%（表1），這一結果與很多學者研究認為“第二次土壤普查以來我國表層土壤全氮含量呈上升趨勢的結果一致”[8-10，20-21]。 水稻土全氮含量上升主要驅動因素為該地區(qū)氮肥和農家肥施用量的大幅度增加（圖3）[4，8，21-22]。據(jù)統(tǒng)計，1982年太湖地區(qū)氮肥和農家肥平均施用量分別為N 243 kg·hm-2和C 231 kg·hm-2，而2000年為345 kg·hm-2和270 kg·hm-2，增幅分別達37.66%和19.86%。此外，氣候因子也在長期的全氮含量變化中起著重要作用，據(jù)統(tǒng)計，1982年和2000年太湖地區(qū)年均降雨量分別為1 102 mm 和1 175 mm，該地區(qū)緩慢上升的年均降雨量有利于土壤氮素的積累[23]。

從全氮含量空間分布來看（圖4），1982年太湖地區(qū)水稻土全氮含量總體由南向北遞減，南北差異明顯。其中，全氮含量較高（＞2.0 g·kg-1）的水稻土主要分布在東部、東南部以及南部地區(qū)，面積為0.40 M hm2，占太湖地區(qū)水稻土總面積的17%，這些區(qū)域氣候溫和，較低的溫度和較高的降雨量有利于土壤全氮的積累（表2）。全氮含量較低（＜1.0 g·kg-1）的水稻土主要分布在太湖地區(qū)西北部地區(qū)，面積為0.12 M hm2，占水稻土的總面積的5.1%，這些區(qū)域整體溫度較高，且降水量偏少，導致土壤水分不足，加速了氮素的礦化分解（表2）[23-24]。2000年太湖地區(qū)全氮含量空間分布格局與1982年基本相同，但全氮含量較高（＞2.0 g·kg-1）的水稻土面積增加到0.61 M hm2；此外，全氮含量較低（＜1.0 g·kg-1）的面積減少至0.06 M hm2。從全氮含量變化量空間分布來看，太湖地區(qū)1982—2000年不同區(qū)域差異也較大，且增減趨勢不一，約68%面積水稻土表層全氮含量有所上升，主要分布在該地區(qū)的西部、西北部和沿江流域，32%面積水稻土全氮含量有所降低，主要分布在太湖地區(qū)的東部、東南部和中部地區(qū)。其中，全氮含量上升范圍在0～0.5 g·kg-1之間的最多，面積達0.92 M hm2，主要分布在太湖地區(qū)北部和西北部，而全氮含量下降主要集中在-0.5～0 g·kg-1范圍內，面積為0.56 M hm2，主要分布在太湖地區(qū)的東北部和東部區(qū)域。太湖地區(qū)水稻土近20年的全氮含量變化量與初始全氮、有機碳含量呈負相關（表2），主要是因為全氮含量較低的區(qū)域，農民會投入更多的肥料來提高作物產量，而全氮含量較高的區(qū)域農民往往會忽略培肥，影響氮素含量[25]。

表1 太湖地區(qū)水稻土全氮含量統(tǒng)計 Table1 Statistics of soil total nitrogen contents in paddy soils of the Tai-Lake Region，China

圖3 1982—2000年太湖地區(qū)氮肥和農家肥施用量的年際變化 Fig.3 Inter-annual variation of N-fertilizer（a））and manure（b））application rates during the period from 1982 to 2000 in the Tai-Lake Region

圖4 水稻土全氮含量及變化量空間分布 Fig.4 Spatial distribution map of paddy STN and its variation in the Tai-Lake Region，China

表2 不同亞類、土區(qū)和行政區(qū)水稻土的初始土壤屬性、氣候和施肥數(shù)據(jù)統(tǒng)計 Table2 Statistics of initial paddy soil properties，climate data and fertilization data relative to soil subgroup，sub-region and administrative area

2.2 太湖地區(qū)不同亞類水稻土全氮含量差異

從圖5可以看出，太湖地區(qū)不同時期各亞類全氮含量及其變化量差異很大。潴育型水稻土是太湖地區(qū)分布區(qū)域最廣的亞類，主要分布在太湖平原和丘陵地區(qū)，面積為1.23 M hm2，占研究區(qū)水稻土總面積的52.8%。該亞類1982—2000年全氮含量由1.59 g·kg-1上升至1.83 g·kg-1，增幅達15.2%（圖5），這主要與其較高的氮肥和農家肥施用量及較高的降雨量有關[25]。據(jù)統(tǒng)計，潴育型水稻土的氮肥和農家肥年均施用量分別達 N 361 kg·hm-2·a-1和C 278 kg·hm-2·a-1，且年均降雨量超過1 200 mm（表2）。

圖5 表層水稻土（0～15 cm）各亞類全氮含量（a）及全氮含量變化幅度（b） Fig.5 STN（a）and variability of STN（b）in the topsoil layer（0～15 cm）relative to subgroup of paddy soil in the Tai-Lake Region，China

脫潛型水稻土起源于潛育型水稻土，經過長期的耕作施肥和治水改土形成該亞類，面積為0.41 M hm2，占研究區(qū)水稻土總面積的17.7%。該亞類全氮含量變化量是各亞類中最高的，為0.32 g·kg-1，這主要與脫潛型水稻土較高的黏粒含量、氮肥和農家肥施用量和年均降雨量有關[23-25]。據(jù)統(tǒng)計，脫潛型水稻土黏粒含量達30%，氮肥和農家肥施用量分別為N 386 kg·hm-2·a-1和C 267 kg·hm-2·a-1，降雨量為1 256 mm（表2）。

滲育型水稻土多分布于地勢較高的地區(qū)，耕作年限較短，面積為0.37 M hm2，占研究區(qū)水稻土總面積的15.9%，該亞類1982—2000年全氮含量由1.23 g·kg-1上升至1.41 g·kg-1，增幅為15.1%。盡管較高的溫度（16.3℃）和較少的降雨量（1 146 mm）導致該亞類的1982年的土壤全氮含量較低[25-26]，但近20年較多的氮肥（N 324 kg·hm-2·a-1）和農家肥施用量（C 240 kg·hm-2·a-1）也使得該亞類全氮含量呈現(xiàn)緩慢增長趨勢（表2）[24-25]。

漂洗型水稻土主要分布于低山丘陵的梯田，強度淋洗作用導致黏粒淋失嚴重，該亞類太湖地區(qū)分布面積為0.20 M hm2，占研究區(qū)水稻土總面積的8.6%[27]。土壤砂粒多（75%）和黏粒少（16%）導致該亞類1982年全氮含量僅為1.18 g·kg-1，但1982—2000年全氮含量上升了0.24 g·kg-1，這與該亞類近20年氮肥（N 303 kg·hm-2·a-1）和農家肥（C 205 kg·hm-2·a-1）施用量相對較高有關[26]。

潛育型水稻土面積為0.10 M hm2，占研究區(qū)水稻土總面積的4.3%，是太湖地區(qū)唯一全氮含量下降的亞類。該亞類長期處于淹水還原狀態(tài)，通氣性不良導致土壤有機氮的礦化作用較弱而易于積累，1982年的土壤全氮含量高達2.33 g·kg-1[24]。但1982—2000年全氮含量下降了0.61 g·kg-1，變化幅度為-26.3%，這主要與該亞類較高的初始有機氮含量、年均溫度和較低的農業(yè)投入有關。據(jù)統(tǒng)計，該亞類有機肥施用量僅為C 225 kg·hm-2·a-1，年均溫度高達16.7℃（表2）[24-26]。

淹育型水稻土在太湖地區(qū)僅有0.01 M hm2，主要分布在丘陵崗地的上部。1982年該亞類全氮含量較低為1.00 g·kg-1，主要與其高溫少雨的氣候條件和土壤性質有關，據(jù)統(tǒng)計，該亞類溫度和降雨量分別為16.1℃和1 080 mm，黏粒含量僅為15%。后期增加氮肥和農家肥的投入，2000年全氮含量提高了0.22 g·kg-1[25]。

2.3 太湖地區(qū)不同土區(qū)水稻土全氮含量差異

太湖地區(qū)是我國糧食高產地區(qū)，在因地制宜、合理布局方面有著豐富的經驗，根據(jù)成土母質、土地利用方式、耕作制度和農業(yè)利用的特點分為4 個土區(qū)：丘陵低山土區(qū)、太湖平原土區(qū)、沖積平原土區(qū)和低洼圩田土區(qū)[12]。

從圖6可以看出，太湖地區(qū)不同時期各土區(qū)水稻土全氮含量及其變化量差異也很大。低洼圩田土區(qū)位于太湖地區(qū)西北部和南部，水稻土面積為0.69 M hm2，主要土壤類型是潴育型水稻土和脫潛型水稻土，約占該土區(qū)水稻土總面積的90%。1982年和2000年土壤全氮含量均很高，分別為1.83 g·kg-1和2.01 g·kg-1，這主要是該地區(qū)較高的黏粒含量（31%）和氮肥、農家肥施用量（分別為N 385 kg·hm-2· a-1和C 249 kg·hm-2·a-1）導致了全氮的富集（表2）[8，24，29]。但1982—2000年間當?shù)剞r民對該土區(qū)的潛育水稻土進行治水改土，使土壤的結構和肥力等發(fā)生了較大的變化，地下水位的降低、微酸性（6.6）和很高的初始全氮含量（1.83 g·kg-1）降低了土壤全氮富集的速度，近20年全氮含量變幅在各土區(qū)中是最低的，僅為9.8%[25，27-28]。

沖積平原土區(qū)沿長江與錢塘江呈帶狀分布，水稻土覆蓋面積為0.64 M hm2，占研究區(qū)水稻土總面積的27.5%[13]。該土區(qū)內側沿江過渡地區(qū)石灰性強，砂粒多，黏粒少，土壤空隙大，導致漏水漏肥較嚴重，初始全氮含量較低[29]。1982—2000年全氮含量富集較快，變化量為0.26 g·kg-1，變化幅度為17.9%。這主要與該土區(qū)較低的初始全氮含量及較高的氮肥和農家肥施用量有關，據(jù)統(tǒng)計，沖積平原土區(qū)水稻土初始全氮含量為1.45 g·kg-1，氮肥和農家肥施用量分別達N 378 kg·hm-2·a-1和C 330 kg·hm-2·a-1（表2）[25]。

低山丘陵土區(qū)主要分布在太湖地區(qū)的西部與北部，水稻田面積為0.39 M hm2，占研究區(qū)水稻土總面積的17.0%[13]。該土區(qū)受地形的影響導致水分流失與地表側漏，土壤質地較粗，黏粒含量僅為22%，氮素難以明顯富集[31]。1982年土壤全氮含量僅為1.38 g·kg-1，低于其他3 個土區(qū)，而2000年土壤全氮含量提升至1.82 g·kg-1，主要原因與該土區(qū)較低的年均溫度以及較高的降雨量有關，據(jù)統(tǒng)計，1982—2000年該土區(qū)的年均溫度和降雨量分別為16.1℃和1 217 mm（表2）[24-25]；此外，較低的初始全氮也導致后期氮肥和有機肥的大量投入，導致近20年間全氮含量增幅達31.8%。

太湖平原土區(qū)位于太湖地區(qū)的東部和北部，水稻土面積為0.59 M hm2，占研究區(qū)水稻土總面積的25.5%[13]。該土區(qū)土壤類型以潴育型水稻土（0.28 M hm2）和漂洗型水稻土（0.13 M hm2）為主。1982—2000年該土區(qū)年均氣溫和年均降水量分別為16.5℃和1 154 mm，高溫少雨的氣候條件不利于氮素的富集，全氮含量僅上升0.04 g·kg-1[24-26]。

圖6 表層水稻土（0～15 cm）不同土區(qū)全氮含量（a）及全氮含量變化幅度（b） Fig.6 STN（a）and variability of STN（b）in the topsoil layer（0～15 cm）relative to soil sub-regions in the Tai-Lake Region，China

2.4 太湖地區(qū)不同行政區(qū)水稻土全氮含量差異

從太湖地區(qū)各省份表層水稻土全氮含量變化來看（表3），1982—2000年江蘇省、上海市和浙江省全氮含量均呈富集趨勢。其中，地處南部及西南部的浙江省水稻土全氮含量最高，增幅最大，由1982年的1.79 g·kg-1增長至2000年的2.22 g·kg-1，變化幅度達24.0%；而北部及西北部的江蘇省水稻土全氮含量最低、增長幅度較小，由1982年的1.48 g·kg-1增長至2000年的1.59 g·kg-1，變化幅度僅為7.4%；上海市水稻土全氮含量及增長幅度位于浙江省與江蘇省之間，由1982年的1.64 g·kg-1增長至2000年的1.89 g·kg-1，變化幅度為15.2%。浙江省水稻土氮素富集程度明顯高于江蘇省，一方面由于該省份農業(yè)投入高，據(jù)統(tǒng)計，浙江省氮肥和農家肥施入量分別為N 418 kg·hm-2·a-1和C 302 kg·hm-2·a-1，而江蘇省氮肥和農家肥投入量僅分別為N 315 kg·hm-2·a-1和C 200 kg·hm-2·a-1；另一方面是因為二者年均降雨量和溫度差異較大，浙江省年均降雨量和溫度分別為1 347 mm 和16.4℃，而江蘇省分別為1 156 mm 和16.5℃，高溫少雨的氣候條件會導致土壤水分不足，土壤孔隙度變大，有利于氮素的礦化分解（表2）。

表3 太湖地區(qū)不同縣（市）水稻土面積、土壤全氮含量、變化量及變化幅度分布 Table3 Area of and soil total nitrogen content，variation of soil total nitrogen in content and variability of soil total nitrogen in content in paddy soil relative to administrative area in the Tai-Lake Region，China

續(xù)表

從太湖地區(qū)各縣（市、區(qū)）表層水稻土全氮含量變化來看（表3），1982—2000年間安吉縣、長興縣、川沙縣、閔行區(qū)、余杭縣和崇明縣全氮含量變化量大于0.6 g·kg-1，變化幅度均超過40%。這可能與閔行區(qū)、川沙縣和余杭縣初始土壤全氮含量較低有關，分別僅為1.24、1.41、1.55 g·kg-1，且安吉縣、余杭縣、長興縣、閔行區(qū)和川沙縣年均降雨量超過1 200 mm，崇明縣和長興縣年均氮肥施用量較高，超過N 350 kg·hm-2·a-1。相反，1982—2000年間吳縣、青浦縣、嘉善縣和丹陽縣4 個縣（市、區(qū)）全氮含量有所下降，變化量分別為-0.42、-0.03、-0.01和-0.004 g·kg-1。這可能與這4 個縣較高的初始全氮含量、初始有機碳含量及氣候條件有關。據(jù)統(tǒng)計，吳縣、嘉善縣、青浦縣和丹陽縣初始全氮含量分別為2.18、2.17、1.98 和1.36 g·kg-1，初始土壤有機碳含量分別為23.45、20.77、20.50 和11.59 g·kg-1，且丹陽縣年均降雨量僅為1 068 mm，年均溫度達到16.2℃。有研究表明，土壤全氮變化與土壤有機碳變化趨勢呈正相關關系[30]。其他各縣（市、區(qū)）全氮含量變化量介于0～0.6 g·kg-1之間，變化幅度介于1%～40%之間。

總體來看，太湖地區(qū)各省份和縣（市、區(qū)）全氮含量變化一方面受到土壤屬性及氣候條件的影響[31]；另一方面也受到化肥施用量的影響，一般情況下高的化肥施用量會導致全氮含量快速富集。相關研究表明，我國肥料利用率平均不到30%，低于世界平均水平，而過量的肥料投入會通過地表進入水體，增加水體營養(yǎng)化的風險[32]。因此，在今后太湖地區(qū)水稻田管理過程中，根據(jù)各個區(qū)域的土壤屬性和肥料投入量制定適宜的農業(yè)管理措施十分重要。

3 結 論

準確評估太湖地區(qū)水稻土全氮含量動態(tài)變化及空間分布可為防控該流域水體富營養(yǎng)化等環(huán)境問題提供重要依據(jù)。從本研究來看，1982—2000年太湖地區(qū)表層水稻土全氮含量整體呈上升趨勢，平均增加了0.21 g·kg-1，氮素富集明顯。從不同水稻土亞類來看，脫潛型、潴育型、漂洗性和淹育型水稻土全氮含量增長均超過0.20 g·kg-1，而潛育型水稻土是唯一下降的亞類，減少了0.61 g·kg-1。從不同土區(qū)來看，低山丘陵土區(qū)和沖擊平原土區(qū)水稻土全氮含量上升較多，分別為0.44 和0.26 g·kg-1，而低洼圩田土區(qū)和太湖平原土區(qū)全氮含量增長緩慢，分別為0.18 和0.04 g·kg-1。從不同行政區(qū)來看，江蘇省、上海市和浙江省太湖地區(qū)水稻土全氮含量均呈不同程度的富集。由于施肥、氣候和土壤理化性質等因素的共同作用，1982—2000年太湖地區(qū)土壤全氮含量產生不同程度的變化，不同亞類、土區(qū)和行政區(qū)全氮含量差異顯著。因此，在今后太湖地區(qū)水稻土種植管理中，根據(jù)不同亞類、土區(qū)和行政區(qū)氮素富集程度適當減少氮肥施用量，以防止氮素流失造成的水體富營養(yǎng)化風險是十分必要的。