局部秋澆條件下凍融期“秋澆-春灌”農田土壤水鹽歸趨研究

范理權，史海濱，閆建文*，李仙岳，竇 旭，祁 茜，李慧祥

（1.內蒙古農業大學 水利與土木建筑工程學院，呼和浩特 010018；2.內蒙古河套灌區水利發展中心義長分中心 沙河所渠供水所，內蒙古 巴彥淖爾 015100）

0 引 言1

內蒙古河套灌區是我國典型的干旱地區，是亞洲最大的一首制平原引水灌區[1-2]。該灌區每年作物收獲后均需進行秋澆，以實現土壤保墑和壓鹽。隨著節水改造工程的實施，灌區年平均引水量減少20%[3]，秋澆面積大幅度縮減，春灌面積增加[4]。春灌農田受到種植制度和農民種植意愿的影響，表現為插花狀的空間分布特征，且面積差異較大。

灌區明溝排水不健全且塌坡嚴重導致排水和排鹽能力較差[5]。秋澆后地下水位上升顯著，秋澆與非秋澆農田在局部秋澆后存在水位差，為水鹽遷移提供了驅動力。凍結前不進行秋澆的農田在第二年會進行春灌，秋澆后至翌年春灌前土壤會經歷凍融期。凍融過程中，秋澆與非秋澆農田土壤的物理性狀和水鹽運移規律存在差異，由此引發了新的水鹽分布問題，需要進一步分析，以指導當地秋澆和春灌的實施。

【研究進展】虎膽·吐馬爾白等[6]通過凍融期覆蓋探究了土壤水鹽時空動態變化規律。李瑞平等[7]基于SHAW 模型提出不同鹽漬化土壤的秋澆節水灌溉制度。馬睿等[8]對義長灌域土壤凍結過程中的地下水向凍土層的遷移量進行了估算。毛俊等[9]探究了鹽分對季節性凍融土壤蒸發的影響規律。富廣強等[10]研究了季節性凍融期鹽荒地的水鹽運移規律。張殿發等[11]對凍融條件下土壤水鹽遷移機理進行了研究，認為凍融過程中鹽分的遷移受溫度梯度、鹽分濃度梯度及溫度影響。彭振陽等[12]針對內蒙古河套灌區局部秋澆條件下農田土壤水鹽運移的特殊性進行了研究。此外，以往針對土壤凍融過程的水、熱、鹽變化規律也開展了大量研究[13-14]。【切入點】在河套灌區秋澆面積減小、春灌面積增加的趨勢下，局部秋澆灌溉將直接影響秋澆農田土壤的水鹽狀況，進而影響翌年的春灌制度。【擬解決的關鍵問題】鑒于此，本研究開展局部秋澆條件下凍融期“秋澆-春灌”農田土壤水鹽歸趨研究，明確局部秋澆條件下凍融期土壤水鹽遷移規律，為局部秋澆經歷凍融過程后的翌年春灌農田灌溉制度提供理論依據。

1 材料與方法

1.1 試驗區概況

試驗區位于內蒙古河套灌區五原縣屈二紅圪旦（41°07′46.98″N，108°12′38.40″E，海拔1 020 m），屬義長灌域廣勝西支渠控制范圍。試驗區地處季節性凍融區，土壤于每年的11 月中下旬開始凍結，在5月中旬凍土層全部消融，凍融歷時約180 d。試驗區年平均最高氣溫為14.8 ℃，年平均最低氣溫為1.6 ℃，年平均降水量為136.3 mm，年平均蒸發量為1 937.9 mm，日照時間為266 h。風速在2.5～3 m/s 之間變化，地下水埋深在1.6～2.1 m 之間變化。試驗區土壤為中度鹽漬化土壤，主要土質為粉土及粉壤土。試驗區土壤理化性質見表1。

表1 試驗區土壤理化性質Table 1 Test area soil physical properties table

1.2 試驗設計與數據采集

試驗持續時間為2021 年9 月24 日—2022 年5 月7 日。試驗區東西寬50 m，南北長40 m，總面積為0.2 hm2，秋澆和春灌農田各占試驗區面積的50%。試驗區西側為2021 年進行秋澆灌溉的農田（以下簡稱秋澆農田），東側為2022 年進行春灌的農田（以下簡稱春灌農田）（圖1）。2021 年11 月23 日—12 月24日為凍結期Ⅰ，2021 年12 月25 日—2022 年2 月17日為凍結期Ⅱ，2022 年2 月18 日—3 月25 日為消融期Ⅰ，2022 年3 月26 日—4 月20 日為消融期Ⅱ。

圖1 研究區和采樣點分布示意Fig.1 Study area and sampling point distribution map

土壤溫度采用HZR8T 四路土壤溫度記錄儀自動監測，每隔2 h 記錄1 次。監測深度為田面以下0、10、20、40、60、80、100 cm。在5 眼地下水觀測井中放置HOBO 水位觀測計，在秋澆和春灌農田中心觀測井中放置HOBO 水質觀測計，每12 h 記錄1 次地下水的水位和鹽分濃度。

采用土鉆法采集0～10、10～20、20～40、40～60、60～80、80～100 cm 土層的土樣。分別在10 月20 日進行秋澆、11 月8 日進行秋澆后、翌年5 月6 日進行春灌前取樣；此外，在凍融期每隔10～15 d 取1 次土樣，試驗期間共計取樣16 次。每次取樣后利用烘干法測定土壤質量含水率，并將烘干土樣粉碎過篩，取5 g 干土和25 g 蒸餾水按1∶5 的土水比混合后攪拌、離心，將離心后的提取液用電導儀（DDS-307A型）測定土壤EC。

用梯形量水堰測定秋澆水量，并收集灌溉水，用電導儀測定EC。10 月30 日，研究區進行秋澆，秋澆水量約為225 mm。灌溉水平均礦化度為0.55 g/L。

1.3 研究方法

土壤電導率轉換為土壤全鹽量的計算式為[15]：

式中：C為土壤全鹽量（g/kg）；EC1:5為土水比1∶5的土壤浸提液電導率（dS/m）。

地下水EC和礦化度（TDS）的轉換式為[16]：

式中：TDS為地下水礦化度（g/L）；EC為地下水電導率（dS/m）。

土壤積鹽率為0～100 cm土層剖面某一時期與前一時期相比土壤含鹽量的增加率，計算式為[17]：

式中：t為土壤積鹽率（%）；Wi為第i時期土壤含鹽量（kg/hm2）；Wi-1為第i-1時期土壤含鹽量（kg/hm2）。

土壤脫鹽率為0～100 cm 土層剖面在灌溉后某一時期與灌溉前相比土壤含鹽量的減少率，其計算式為[18]：

式中：ω為土壤脫鹽率（%）；St為灌溉后某一時期土壤含鹽量（kg/hm2）；S0為灌溉前土壤含鹽量（kg/hm2）。

采用式（5）計算第j次取樣時第i層土壤單位面積的儲鹽量[19]：

考慮一維條件下的各土層鹽分平衡，土層內儲鹽量的增加量與其流入與流出量之差相等，計算式為[19]：

1.4 數據處理

采用 Microsoft Excel 2016 軟件進行數據分析，利用 Surfer 12.0 繪制土壤水分和鹽分的時空分布圖。

2 結果與分析

2.1 秋澆-非秋澆（翌年春灌）農田土壤溫度變化特征

由圖2 可知，隨著土層深度的增加，在凍結和消融過程中土壤溫度隨時間的變化均表現出滯后效應，滯后時間隨著土層深度的增大而增加。隨著土層深度的增加，凍結過程中，溫度降幅逐漸變小；消融過程中，溫度漲幅逐漸變大。

圖2 秋澆-非秋澆（翌年春灌）農田土壤溫度變化Fig.2 Temperature changes of each layer of autumn watering - non-autumn irrigation (spring irrigation the following year)

局部秋澆后1 周，秋澆和春灌農田各土層平均溫 度和秋澆前1 周相比均發生了不同程度的降低，秋澆農田平均溫度降幅為0.38～1 ℃；春灌農田平均溫度降幅為1.28～2.38 ℃。由于秋澆農田表面存在隔水層，具有一定的保溫作用，因此秋澆農田降幅小于春灌農田。隨著土層深度的增加，秋澆農田土壤溫度降幅逐漸變大，而春灌農田降幅逐漸變小。凍融過程中，秋澆和春灌農田僅0～60 cm 土層存在負溫。在凍結過程中，秋澆農田0～60 cm 土層溫度降至0 ℃歷時50 d，而春灌農田歷時82 d。在消融階段，秋澆農田0～60 cm土層溫度升至0 ℃歷時55 d，而春灌農田歷時41 d。其中，春灌農田40 cm 深度處的土壤溫度降至0 ℃相比秋澆農田遲8 d，60 cm 深度處的土壤溫度降至0 ℃相比秋澆農田遲23 d。受到上層凍土的阻礙，春灌農田在凍結期Ⅱ的土壤含水率增加主要是由灌溉水滲流所導致，所以春灌農田隨著土層深度的增加溫度降至0 ℃的時間要遲于秋澆農田。

2.2 秋澆-非秋澆（翌年春灌）農田地下水埋深及礦化度變化

由圖3（a）可知，局部秋澆后，地下水位達到最高，秋澆、春灌農田平均地下水埋分別為0.47、0.60 m。凍結期Ⅰ末，地下水位開始下降；翌年3月初，水位降至最低，此時秋澆和春灌農田平均地下水埋深分別為2.06、2.17 m。消融期Ⅰ結束，水位再次上升，秋澆和春灌農田平均地下水埋深分別為0.84、0.82 m。由圖3（b）可知，秋澆開始后，秋澆和春灌農田的地下水礦化度開始提高，而后趨于下降。在局部秋澆后至凍結期Ⅰ末，秋澆和春灌農田地下水位之間相互波動變化，凍結期Ⅰ中12月2—4日、12月15—17日和12月21—23日3個時段的春灌農田水位高于秋澆農田。

圖3 秋澆-非秋澆（翌年春灌）農田地下水埋深和礦化度的動態變化Fig.3 Dynamic changes of groundwater depth and mineralization degree of farmland under autumn watering and non-autumn

2.3 秋澆-非秋澆（翌年春灌）農田土壤剖面含水率時空變化特征

由圖4可知，秋澆后，秋澆農田0～20 cm土層含水率增加74.81%，為春灌農田的2.24倍。凍結期Ⅰ，隨著土壤溫度的降低，秋澆農田0～40 cm土層含水率相比秋澆后增加41.76%。在凍結期Ⅱ，秋澆農田20～80 cm土層含水率增幅為7.04%～40.18%。消融期Ⅰ的秋澆農田0～100 cm土層含水率相比凍結期Ⅱ減少10.44%。由于凍土層尚未融通，因此該時期土壤含水率沒有明顯減小。進入消融期Ⅱ，隨著凍土層逐漸融通，秋澆農田0～80 cm土層含水率相比消融期Ⅰ減少16.35%，而80～100 cm土層含水率增加12.46%。

圖4 秋澆-非秋澆（翌年春灌）農田土壤剖面含水率分布Fig.4 Distribution of soil moisture content in autumn watering - non-autumn irrigation (spring irrigation the following year)

由于春灌農田沒有進行秋澆灌溉，因此主要受灌區地下水位整體上升和灌溉水側向滲流的影響。0～100 cm土層含水率相比秋澆前增加19.91%。凍結期Ⅰ，春灌農田0～40 cm土層含水率相比秋澆后增加57.74%。凍結期Ⅰ秋澆和春灌農田之間的水位相互波動，且地下水會補給土壤水，因此該時期春灌農田0～40 cm土層含水率明顯增大。消融期春灌農田含水率變化規律與秋澆農田基本一致。

經歷凍融過程后，翌年春灌前，土壤含水率由表層至深層呈遞增趨勢，此時秋澆農田與年前秋澆前相比含水率增幅為春灌農田的2.96 倍，春灌農田0～100 cm 土層含水率平均增加8.62%，主要受到灌溉水滲流作用的影響。

2.4 秋澆-非秋澆（翌年春灌）農田土壤剖面鹽分時空變化

由圖5 可知，秋澆后，秋澆農田0～80 cm 土層脫鹽率為18.45%～49.74%，80～100 cm土層積鹽5.91%。秋澆農田在凍結期Ⅰ的20～40 cm 土層積鹽4.05%，凍結期Ⅱ的40～60 cm 土層積鹽10.71%，其他土層在整個凍結期均呈脫鹽狀態。在消融期Ⅰ，秋澆農田積鹽47.32%。由于凍土層融化，水分釋放會攜帶鹽分，并且該時期凍層尚未徹底融通，故凍層和凍層以下土層鹽分量增加。進入消融期Ⅱ，隨著凍層的融通，秋澆農田和消融期Ⅰ相比0～80 cm 土層脫鹽52.92%，而80～100 cm 土層積鹽8.72%。

圖5 秋澆-非秋澆（翌年春灌）農田土壤剖面鹽分分布Fig.5 Soil cross-sectional salt distribution map of autumn watering-non-autumn irrigation (spring irrigation the following year)

春灌農田秋澆后，0～80 cm 土層脫鹽57.44%，80～100 cm 土層積鹽43.24%。由于春灌農田沒有進行秋澆灌溉，造成其脫鹽的原因是秋澆后地下水位上升，秋澆和春灌農田間水分存在水平側向補給，使春灌農田土壤鹽分得到了淋洗。春灌農田在凍結期Ⅰ積鹽42.59%，由于凍結期Ⅰ秋澆和春灌農田之間水位相互波動，在水勢梯度的作用下，秋澆農田水分通過側向徑流將鹽分傳輸至春灌農田，導致其鹽分量增加。凍結期Ⅱ春灌農田脫鹽19.83%，地下水位的下降是造成春灌農田凍結期Ⅱ脫鹽的主要原因。消融階段春灌農田變化規律和秋澆農田基本一致。在消融期Ⅰ，春灌農田積鹽27.12%。

經歷凍融過程后，在翌年春灌前，春灌農田土層含鹽量為秋澆農田的1.60倍。此時的秋澆農田與秋澆前相比1 m土層整體脫鹽58.35%；而春灌農田0～60 cm土層脫鹽17.29%；60～100 cm土層積鹽41.59%。

2.5 秋澆-非秋澆（翌年春灌）農田鹽分通量

由表2和表3可知，秋澆和春灌農田各時期不同深度土壤鹽分運動方向存在較強的一致性，在秋澆期、凍結期Ⅰ和消融期Ⅰ鹽分整體向下運動，其他時期鹽分則整體向上運動。凍結期Ⅰ，受到秋澆灌溉的影響，秋澆農田各土層鹽分通量為春灌農田的1.61～2.32倍。隨著溫度的持續下降，凍土層逐漸變厚，在凍結作用影響下，凍結期Ⅱ鹽分整體向上運動。消融期Ⅰ，秋澆和春灌農田各土層鹽分通量均隨著土層深度的增大而增大，導致消融階段土層積鹽，秋澆農田各土層鹽分通量為春灌農田3～15倍。消融期Ⅱ，隨著土層逐漸融通，蒸發作用增強，春灌農田各層土壤養分整體向上運動，但60～100 cm土層鹽分通量明顯大于0～60 cm土層，導致深層土壤積鹽。秋澆農田0～60 cm土層鹽分通量隨著土層深度的增加而減小，而春灌農田逐漸增加。

表2 秋澆農田各時期各土層鹽分通量Table 2 The salt flux of each soil layer in each period of autumn watering farmland is divided g/m2

表3 春灌農田各時期各土層鹽分通量Table 3 The salt flux of each soil layer in each period of spring irrigation farmland g/m2

由整個試驗期鹽分凈通量可知，秋澆和春灌農田各土層鹽分通量均為負值，但秋澆農田隨著土層深度的增加，鹽分通量逐漸增加，這說明秋澆農田經歷凍融過程后，在翌年春灌前呈脫鹽狀態；春灌農田隨著土層深度的增加先增加后減小，這也說明春灌農田存在深層積鹽。

3 討 論

由于節水改造工程的實施，河套灌區引水量逐年減少，多年的水鹽平衡體系被打破[18]，灌區內部鹽分普遍重新分配。隨著引水量的減少，秋澆面積減少，春灌面積相應增加[4]。局部秋澆后，灌溉和非灌溉農田是凍融期耕地鹽分遷移發生的主要區域。

溫度是土壤水分遷移的主要驅動力。局部秋澆后1 周與秋澆前1 周相比，秋澆農田0～100 cm 土層溫度降幅為0.38～1.00 ℃；春灌農田溫度降幅為1.28～2.38 ℃。秋澆農田灌溉后形成隔水層，具有一定的保溫作用，導致秋澆農田溫度降幅小于春灌農田[20]。本研究中，局部秋澆后1 周，隨著土層深度的增加，秋澆農田溫度降幅變大，春灌農田溫度降幅變小。凍結過程中，秋澆農田0～60 cm 土層溫度降至0 ℃相比春灌農田早30 d，消融過程中，秋澆農田0～60 cm 土層溫度升至0 ℃相比春灌農田遲10 d。秋澆農田各土層最低溫度均小于春灌農田，且隨著土層深度的增加，溫差變大。由于春灌農田沒有進行秋澆灌溉，土壤含水率低于秋澆農田。

研究表明，水分在溫度梯度下的運動決定了鹽分的運動趨勢[11]。在凍結期Ⅰ，秋澆農田0～40 cm 土層含水率相比秋澆后增加了41.76%；在凍結期Ⅱ，20～80 cm 土壤含水率增幅為7.04%～40.18%。凍結時，溫度降低導致水分在土水勢作用下向凍層遷移，非凍結層的水分向上層凍結區移動導致凍結層含水率增大，這與以往研究結果一致[11]。春灌農田在凍結期Ⅰ積鹽42.59%，與以往研究結果一致[21-22]。目前，關于凍結過程中鹽分的遷移規律尚不統一。研究表明，水分聚集的凍結層中未發生積鹽[10]，反而呈脫鹽狀態[13,23]。本研究中，秋澆后和凍結期是秋澆農田脫鹽的關鍵時期，并且秋澆農田在凍結期Ⅱ的脫鹽率大于在凍結期Ⅰ的脫鹽率，而春灌農田在凍結期Ⅱ脫鹽19.83%。本研究中，凍結期Ⅱ地下水位開始下降，因此導致春灌農田在凍結期Ⅱ脫鹽。研究表明，凍結前土壤溶液當中的鹽分濃度梯度為正時，凍結層易積鹽；當凍結前土壤溶液鹽分濃度梯度為負時，凍結層易脫鹽[24]。

由凍融期內各土層鹽分凈通量可知，秋澆和春灌農田各土層鹽分凈通量均為負值，但秋澆農田鹽分通量隨著土層深度的增加而增加，春灌農田隨著土層深度的增加先增加后減小。土壤凍結過程中鹽分遷移受土壤類型、初始含水率、土壤溶液濃度、鹽分組成以及凍結溫度等因素共同影響[11]。春灌農田與秋澆農田的相鄰區域，在局部秋澆后經歷凍融過程中，土壤水分和鹽分受到了一定程度影響，秋澆不進行灌溉的農田第二年會進行春灌，因此，對于合理的春灌制度有待進一步研究。

4 結 論

1）局部秋澆后，在凍結過程中，秋澆農田0～60 cm土層溫度降至0 ℃相比春灌農田早30 d，在消融過程中，秋澆農田0～60 cm土層溫度升至0 ℃相比春灌農田遲10 d。

2）局部秋澆后至凍結期Ⅰ末，秋澆和春灌農田水位之間相互波動變化，形成水頭差，成為水分遷移的主要驅動力。

3）秋澆后和凍結期是秋澆農田脫鹽的關鍵時期，秋澆農田從秋澆后至凍結期脫鹽47.38%；而春灌農田在凍結、消融階段分別積鹽35.68%、16.87%。

4）在整個凍融期內，秋澆和春灌農田各層鹽分凈通量均為負值，但秋澆農田鹽分通量隨著土層深度的增加而增加；而春灌農田隨著土層深度的增加先增加后減小。

（作者聲明本文無實際或潛在的利益沖突）