陳 誠(chéng)

(黃河水利職業(yè)技術(shù)學(xué)院，河南 開(kāi)封 475003)

1 項(xiàng)目簡(jiǎn)介

河南省黃泛區(qū)包括全省20個(gè)縣，涉及耕地24×104hm2。在廣闊的黃泛區(qū)，農(nóng)田土壤漬澇障礙是伴生的。由于地上黃河側(cè)滲導(dǎo)致的漬澇、由于黃河改道及引黃淤灌導(dǎo)致的農(nóng)田土壤是導(dǎo)致河南省黃泛區(qū)農(nóng)田生產(chǎn)力低下的重要根源[1]。周振民在《黃河下游引黃灌區(qū)農(nóng)業(yè)生態(tài)環(huán)境脆弱性評(píng)價(jià)及對(duì)策》中指出，為了提高區(qū)域土水資源利用效率及可持續(xù)利用，保障中原經(jīng)濟(jì)區(qū)快速發(fā)展戰(zhàn)略的實(shí)現(xiàn)，必須深入探究農(nóng)田土壤障礙發(fā)展演變規(guī)律，并構(gòu)建農(nóng)田障礙消減與風(fēng)險(xiǎn)防控對(duì)策[2]。本項(xiàng)目實(shí)施區(qū)為河南省封丘縣域。封丘位于黃河沖積扇平原，南東兩面緊鄰黃河，該縣域農(nóng)田土壤質(zhì)地和剖面質(zhì)地層次呈現(xiàn)豐富的多樣性。封丘縣是我國(guó)典型重要農(nóng)業(yè)生產(chǎn)區(qū)，該縣域農(nóng)田土壤漬澇障礙在廣大的黃泛區(qū)具有充分代表性。因此，探究封丘縣域農(nóng)田土壤漬澇障礙的發(fā)生演變機(jī)制，構(gòu)建農(nóng)田土壤障礙加劇風(fēng)險(xiǎn)防控對(duì)策體系，對(duì)于黃泛區(qū)農(nóng)田土壤障礙消減與農(nóng)業(yè)生態(tài)系統(tǒng)質(zhì)量提升，保障國(guó)家糧食安全和加快中原經(jīng)濟(jì)區(qū)發(fā)展戰(zhàn)略，具有重要的理論意義和實(shí)際應(yīng)用價(jià)值。

2 研究現(xiàn)狀

據(jù)國(guó)土資源部提供的資料顯示，1996～2008年13年間我國(guó)耕地面積減少0.083 33×108hm2,平均每年減少0.006 667×108hm2[3]。通過(guò)對(duì)沙化、酸化和漬澇等障礙性農(nóng)田土壤改性治理，實(shí)現(xiàn)耕地資源的提質(zhì)擴(kuò)量，是保證糧食產(chǎn)量的有效途徑[4]。分析農(nóng)田土壤主要障礙并揭示其發(fā)生機(jī)理，對(duì)于消減農(nóng)田土壤障礙、提高農(nóng)田土壤質(zhì)量具有極為重要的意義。農(nóng)田土壤漬澇等障礙導(dǎo)致農(nóng)田土壤質(zhì)量降低的問(wèn)題，已經(jīng)引起國(guó)內(nèi)外學(xué)者的高度關(guān)注[5～12]，一般認(rèn)為研究農(nóng)田土壤障礙發(fā)生演變機(jī)理的關(guān)鍵基礎(chǔ)是對(duì)農(nóng)田土壤障礙信息精確解譯和評(píng)估，其中最為普遍采用的評(píng)估方法和手段是遙感信息解譯[13～15]。如王玉剛等利用遙感數(shù)據(jù)、GIS及地統(tǒng)計(jì)方法，針對(duì)三工河流域中上游阜康綠洲開(kāi)展研究，得到典型綠洲區(qū)農(nóng)田土壤鹽分的空間變異特征[16]。王金哲等通過(guò)實(shí)驗(yàn)室數(shù)據(jù)和整合所搜集資料，查明環(huán)渤海低平原區(qū)的0～40 cm農(nóng)田土壤含鹽量現(xiàn)狀，在此基礎(chǔ)上，以不同鹽分情況下的農(nóng)田土壤安全容鹽值、農(nóng)田土壤化學(xué)類(lèi)型和作物布局結(jié)構(gòu)為參考指標(biāo)，對(duì)研究區(qū)的農(nóng)田土壤安全容鹽潛力進(jìn)行評(píng)價(jià)[17]。周在明等應(yīng)用協(xié)同克立格插值方法，對(duì)環(huán)渤海低平原農(nóng)田土壤鹽分的空間變異特征進(jìn)行研究[18]。李云鵬等分析不同時(shí)期的遙感資料，演算不同時(shí)期歸一化植物指數(shù)NDVI值，基于遙感影像技術(shù)和光譜特征直接提取土地退化指標(biāo)，分析內(nèi)蒙古自治區(qū)植被覆蓋宏觀變化和荒漠化的情況[19]。Bouaziz Moncef等利用遙感影像提取出影響農(nóng)田土壤障礙的光譜指數(shù)，建立評(píng)價(jià)巴西東北部半干旱區(qū)農(nóng)田土壤障礙的LSU模型[20]。除此之外，還有學(xué)者根據(jù)各障礙類(lèi)型的特點(diǎn)，采取相應(yīng)的評(píng)估方法[21～23]。如石進(jìn)校等通過(guò)不同水平的漬澇處理，研究漬澇對(duì)淫羊藿葉片膜脂過(guò)氧化和SOD活性的影響，為耐澇性淫羊藿的栽培提供理論依據(jù)[24]。Visconti F等針對(duì)內(nèi)陸灌區(qū)農(nóng)田土壤質(zhì)量退化問(wèn)題，同時(shí)考慮種植制度、氣象條件、農(nóng)田土壤物理化學(xué)性狀、灌溉水源特征等要素，對(duì)灌區(qū)農(nóng)田土壤漬澇和鹽化障礙進(jìn)行中長(zhǎng)期預(yù)測(cè)[25]。

3 研究方案

圍繞揭示河南黃泛區(qū)農(nóng)田土壤漬澇障礙的發(fā)生演變機(jī)制、構(gòu)建障礙加劇風(fēng)險(xiǎn)防控與減障對(duì)策體系這一核心研究目標(biāo)，本項(xiàng)目集成應(yīng)用多種研究手段與方法、多尺度多時(shí)相數(shù)據(jù)協(xié)同集成融合的整體學(xué)術(shù)思路，以依托于承擔(dān)單位的國(guó)家和省部級(jí)重點(diǎn)實(shí)驗(yàn)室、野外生態(tài)實(shí)驗(yàn)站、研究中心為研究基地，系統(tǒng)開(kāi)展項(xiàng)目研究工作。本項(xiàng)目應(yīng)用研究區(qū)農(nóng)田土壤漬澇障礙季際、年際動(dòng)態(tài)變化的信息調(diào)研采集，獲得區(qū)域農(nóng)田土壤漬澇障礙的現(xiàn)狀特征。通過(guò)對(duì)區(qū)域農(nóng)田土壤障礙要素歷史資料的整理、匯總與詳細(xì)分析，結(jié)合區(qū)域30年來(lái)4個(gè)代表性時(shí)相遙感影像信息反演，利用GIS制圖與空間數(shù)值分析技術(shù)，獲得各研究時(shí)相的區(qū)域農(nóng)田土壤漬澇障礙的特征參數(shù)(含分布面積、障礙程度等)。綜合農(nóng)田土壤漬澇的各時(shí)相特征，解析區(qū)域農(nóng)田土壤漬澇障礙的演變規(guī)律與發(fā)展趨勢(shì)。

布置田間土壤漬澇障礙動(dòng)態(tài)監(jiān)測(cè)，收集氣象數(shù)據(jù)，監(jiān)測(cè)地下水特征信息、農(nóng)田土壤剖面水分等指標(biāo)數(shù)據(jù)，獲得面域原位農(nóng)田土壤漬澇障礙驅(qū)動(dòng)要素的動(dòng)態(tài)監(jiān)測(cè)數(shù)據(jù)；設(shè)置障礙驅(qū)動(dòng)要素的不同組合(不同土體構(gòu)型與黃河水位)條件下的黃河水文水力學(xué)側(cè)滲和引黃淤灌模擬實(shí)驗(yàn)，獲得不同試驗(yàn)條件下農(nóng)田土壤漬澇障礙發(fā)展程度與范圍特征參數(shù)；通過(guò)田間開(kāi)展驅(qū)動(dòng)因子作用試驗(yàn)，研究灌溉、氣象條件等要素對(duì)農(nóng)田土壤漬澇障礙的作用規(guī)律與內(nèi)在機(jī)理；結(jié)合通過(guò)對(duì)歷史漬澇障礙要素資料的整理分析和遙感信息反演，獲得漬澇障礙主導(dǎo)驅(qū)動(dòng)因子演變特征和對(duì)應(yīng)時(shí)相的農(nóng)田土壤漬澇障礙響應(yīng)特征，確定區(qū)域農(nóng)田土壤漬澇障礙的主導(dǎo)驅(qū)動(dòng)因子，并揭示出主導(dǎo)驅(qū)動(dòng)因子的驅(qū)動(dòng)機(jī)理。

4 農(nóng)田土壤漬澇障礙的原位動(dòng)態(tài)監(jiān)測(cè)方法

在各障礙類(lèi)型區(qū)，設(shè)置農(nóng)田土壤漬澇障礙原位動(dòng)態(tài)監(jiān)測(cè)網(wǎng)，在研究區(qū)設(shè)置4個(gè)監(jiān)測(cè)斷面，每個(gè)斷面設(shè)置6個(gè)點(diǎn)位，每個(gè)點(diǎn)位設(shè)置1眼地下水觀測(cè)井，并分4層埋設(shè)農(nóng)田土壤水分負(fù)壓計(jì)(埋設(shè)分別為20、45、80和125 cm)，每10天測(cè)定1次地下水位，系統(tǒng)獲得農(nóng)田土壤漬澇障礙的動(dòng)態(tài)數(shù)據(jù)。

在研究區(qū)，總體采用網(wǎng)格法布點(diǎn)，采樣間距設(shè)計(jì)為3 km×3 km，采樣層次設(shè)計(jì)為：0～30 cm，30～60 cm，60～100 cm，100～150 cm，150～200 cm，200～250 cm，……，采至恒定地下水位(通過(guò)休灌期人工測(cè)定靜態(tài)井水位確定)。實(shí)驗(yàn)室分析測(cè)定農(nóng)田土壤質(zhì)地組成。連續(xù)進(jìn)行第1年秋季、第2年春季和第2年秋季的3次面域農(nóng)田土壤樣品采集與詳細(xì)調(diào)研，系統(tǒng)獲取田間土壤漬澇障礙的季際和年際的動(dòng)態(tài)變化，以及障礙的現(xiàn)狀特征。

考慮研究區(qū)漬澇的主要成因?yàn)榈厣闲忘S河水側(cè)滲補(bǔ)給，又鑒于當(dāng)前尚未有關(guān)于此類(lèi)型漬澇的統(tǒng)一標(biāo)準(zhǔn)，結(jié)合河南黃泛區(qū)的自然地理等綜合條件，本研究擬按照地下水位埋深小于0.5 m的天數(shù)T分別為：T≤5 d，5 d >T≤15 d，15 d>T≤30 d，30 d>T≤45 d，和T>45 d，相應(yīng)分為非漬澇、輕度、中度、重度和嚴(yán)重漬澇計(jì)5個(gè)等級(jí)。

依據(jù)上述農(nóng)田土壤漬澇障礙的分級(jí)標(biāo)準(zhǔn)，對(duì)季際、年際的農(nóng)田土壤漬澇障礙進(jìn)行特征信息解譯與評(píng)估，并繪制相應(yīng)障礙評(píng)估圖集。

5 區(qū)域農(nóng)田土壤漬澇障礙的規(guī)律分析

采取室內(nèi)判讀與野外考察相結(jié)合的方法，對(duì)研究區(qū)歷史農(nóng)田土壤障礙信息資料(如統(tǒng)計(jì)年鑒、全國(guó)第二次農(nóng)田土壤普查資料、地方農(nóng)田土壤、水利志等)進(jìn)行收集匯總與整理分析，分析農(nóng)田土壤漬澇、障礙與遙感影像真彩色組合特征的對(duì)應(yīng)關(guān)系，確立解譯標(biāo)志；應(yīng)用人機(jī)交互解譯方法，對(duì)發(fā)生漬澇、障礙的農(nóng)田土壤開(kāi)展遙感判讀分類(lèi)；利用GIS軟件對(duì)各時(shí)相的農(nóng)田土壤障礙特征圖進(jìn)行空間疊加，確定出各時(shí)相間的研究區(qū)域農(nóng)田土壤障礙的時(shí)空演變規(guī)律；應(yīng)用轉(zhuǎn)移矩陣統(tǒng)計(jì)學(xué)分析方法，解析農(nóng)田土壤發(fā)生漬澇障礙的總體數(shù)量演變特征。

5.1 不同地下水埋深條件下土壤水分運(yùn)動(dòng)特征

典型易漬澇區(qū)土壤水分運(yùn)動(dòng)存在顯著的季節(jié)性特征，1 m地下水位和3 m地下水位條件下土壤水分的季節(jié)性運(yùn)動(dòng)規(guī)律整體比較一致，存在顯著的周期性，均表現(xiàn)為夏秋季節(jié)土壤剖面水分含量高、冬春季節(jié)則低，且表層土壤水分含量季節(jié)間波動(dòng)幅度較大，靠近地下水面的底層土壤水分含量季節(jié)間波動(dòng)幅度較小?？刂仆寥榔拭嫠謩?dòng)態(tài)季節(jié)性分布規(guī)律的主要原因是降雨、蒸發(fā)等氣象學(xué)因子，其中研究區(qū)季節(jié)間降雨量動(dòng)態(tài)差異是主要影響因子，研究區(qū)夏季降雨最集中、土壤剖面水分含量最高，秋季降雨少但土壤剖面水分含量往往仍保持較高的原因是土壤水分動(dòng)態(tài)的滯后效應(yīng)，冬春季降雨少、土壤剖面水分含量低。

典型易漬澇區(qū)土壤剖面水分運(yùn)動(dòng)存在顯著的垂直分布特征，1 m地下水位和3 m地下水位條件下土壤水分運(yùn)動(dòng)的垂直分布規(guī)律整體比較一致，土壤水分含量由土壤表層至底層呈現(xiàn)逐漸增加趨勢(shì)，即離地下水位距離愈小、土壤水分含量愈高?？刂仆寥榔拭嫠謩?dòng)態(tài)垂直分布規(guī)律的主要因素是水動(dòng)力學(xué)因子，隨著與地下水面距離的增加，非飽水帶土壤剖面水分補(bǔ)給逐漸降低。同時(shí)，土壤水蒸發(fā)強(qiáng)度逐漸增加，加上植被種植下的根系吸水作用，導(dǎo)致土壤剖面表層含水量最低。見(jiàn)圖1、圖2。

進(jìn)一步研究土壤含水量隨不同地下水埋深條件的變化特征(圖3)，結(jié)果表明，隨著地下水位埋深增加，土壤剖面各層土壤平均含水量呈逐漸下降的規(guī)律。對(duì)土壤剖面不同土層水分含量隨地下水位變化而降低的趨勢(shì)進(jìn)行擬合分析(表1)，結(jié)果表明，土壤剖面各層含水量隨地下水埋深增加均呈現(xiàn)顯著的線性下降趨勢(shì)，其中表土含水量隨地下水埋深增加而下降的幅度最大(回歸方程斜率最大)，底層土壤含水量下降幅度則相對(duì)較小。

圖3 不同地下水位條件對(duì)土壤剖面水分含量變化的影響

項(xiàng) 目0～20 cm20～40 cm40～60 cm整個(gè)土體平均回歸方程Y=-0.080x+0.411Y=-0.068x+0.420Y=-0.061x+0.418Y=-0.045x+0.387R20.9610.9420.9550.950

注:Y為多年平均土壤體積含水量(%)，x為地下水埋深(m)，樣本數(shù)均為5。

5.2 不同土體構(gòu)型條件下土壤水分運(yùn)動(dòng)特征

不同土體構(gòu)型條件下，土壤水分的季節(jié)間運(yùn)動(dòng)規(guī)律整體一致，均表現(xiàn)為夏秋季土壤含水量高、冬春季土壤含水量低。但黏土層水分運(yùn)動(dòng)的季節(jié)間波動(dòng)幅度相對(duì)較小(圖4含厚30 cm和厚1 m黏土夾層土體的40～60 cm層是黏土層)，這是因?yàn)轲ね翆虞^差的透水性能使土壤水分運(yùn)動(dòng)具有一定的滯后效應(yīng)，進(jìn)而影響土壤水分波動(dòng)表現(xiàn)得較平穩(wěn)。

圖4 不同土體構(gòu)型條件下的土壤水分動(dòng)態(tài)

不同土體構(gòu)型條件下，土壤剖面水分含量的垂直分布規(guī)律差異顯著(圖5)。全剖面粉砂壤土體剖面含水量由表層至底層呈逐漸增加的規(guī)律，30～60 cm含黏土夾層土體和0～100 cm為黏土層土體的黏土層部位土壤含水量則明顯高于剖面其他粉砂壤土層，表明黏土層具有良好的保水能力。

圖5 不同土體構(gòu)型條件下的土壤剖面水分分布特征

進(jìn)一步對(duì)不同土體構(gòu)型條件下的不同土層含水量差異進(jìn)行方差分析和多重比較分析(表2)，0～100 cm為黏土層的土柱0～20 cm層和60～80 cm層含水量均顯著高于另二種土體構(gòu)型土柱的0～20 cm和60～80 cm粉砂壤土層；30～60 cm含黏土夾層土柱和0～100 cm為黏土層土柱的20～60 cm層含水量在0.05顯著水平下不存在顯著差異，但均顯著高于全剖面粉砂壤土土柱20～60 cm的粉砂壤土層含水量。隨著黏土層厚度增加，土體平均含水量也呈顯著增加趨勢(shì)。

表2 不同土體構(gòu)型條件下不同土層含水量的差異顯著性分析(0.05顯著水平下)

5.3 不同植被覆蓋條件下土壤水分運(yùn)動(dòng)特征

不同植被覆蓋條件下，土壤水分含量的季節(jié)性動(dòng)態(tài)規(guī)律基本一致(圖6)，整體均表現(xiàn)為夏秋季土壤含水量高、冬春季土壤含水量低。但植被種植條件下，土壤水分含量的季節(jié)間波動(dòng)幅度相對(duì)更大，這是因?yàn)橹脖桓滴饔玫挠绊?，使土壤水分運(yùn)動(dòng)更加活躍。

圖6 不同植被覆蓋條件下的土壤水分動(dòng)態(tài)

不同植被覆蓋條件下，土壤剖面水分含量的垂直分布規(guī)律較一致(圖7)，整體均表現(xiàn)為表土含水量低、底土含水量高。但受植被耗水作用的影響，植被種植條件下土壤剖面水分含量普遍低于裸土土體，其中表土平均含水量差異最大，平均降低15.0%。進(jìn)一步對(duì)不同植被覆蓋條件下土體各層土壤含水量的差異顯著水平進(jìn)行方差分析和多重比較分析(表3)，在0.01顯著水平下，不同植被覆蓋條件下的表土含水量呈極顯著差異。

圖7 不同植被覆蓋條件下的土壤剖面水分分布特征

項(xiàng) 目0～20 cm20～40 cm40～60 cm60～80 cm80～100 cm100～120 cm含水量0.01Sig.含水量0.01Sig.含水量0.01Sig.含水量0.01Sig.含水量0.01Sig.含水量0.01Sig.有植被覆蓋17.1%A22.3%A23.5%A25.4%A27.2%A28.6%A無(wú)植被覆蓋20.1%B23.4%A24.6%A26.3%A27.6%A29.5%A水分降低15.0%4.4%4.6%3.4%1.4%3.0%

不同植被覆蓋條件下，潛水蒸發(fā)量差異反映了不同的潛水利用效率(表4)。植被生長(zhǎng)需要大量吸收和消耗土壤水分，因此不同地下水位下有無(wú)植被覆蓋土柱的多年平均潛水蒸發(fā)量之比普遍大于1，且隨著地下水埋深增加，該比值呈增大趨勢(shì)，表明隨著地下水埋深增加，植被覆蓋條件下土柱的潛水利用率呈增加趨勢(shì)。

表4 不同植被覆蓋條件下的土體多年平均潛水蒸發(fā)量差異

6 結(jié) 論

土壤水分動(dòng)態(tài)的長(zhǎng)期定位監(jiān)測(cè)試驗(yàn)研究表明，典型易漬澇區(qū)土壤水分運(yùn)動(dòng)年際間差異不大，不同地下水位條件、土體構(gòu)型條件和植被覆蓋條件下，土壤水分運(yùn)動(dòng)均具有明顯的垂直分布特征和季節(jié)性動(dòng)態(tài)規(guī)律，土壤剖面水分含量由表層至底層呈逐漸增加規(guī)律，且表土層水分季節(jié)間波動(dòng)幅度大，夏秋季土壤水分含量高、冬春季土壤水分含量低。

隨地下水埋深增加，土壤剖面含水量呈線性降低趨勢(shì)，其中表土層水分降低幅度較大；受植物吸水作用影響，植被覆蓋條件下土柱的潛水利用效率高于無(wú)植被覆蓋土柱，且潛水利用率隨地下水埋深增加而提高，植被覆蓋條件下土壤剖面水分含量整體偏低于無(wú)植被覆蓋土壤剖面，尤其是表層土壤水分差異顯著，植被覆蓋條件下土壤水分運(yùn)動(dòng)更活躍、季節(jié)性波動(dòng)幅度增加；與通體粉砂壤質(zhì)地土壤相比，含黏土層土體黏土層的水分運(yùn)動(dòng)較平穩(wěn)、季節(jié)性波動(dòng)幅度較小，土壤剖面黏土層部位水分含量顯著增加，黏土層具有良好的保水能力，土體平均含水量隨黏土層厚度增加呈顯著增加趨勢(shì)。