袁成福，馮紹元(.江西水利職業(yè)學(xué)院，南昌 33003；2.揚(yáng)州大學(xué)水利與能源動(dòng)力工程學(xué)院，江蘇 揚(yáng)州 225009；3.中國(guó)農(nóng)業(yè)大學(xué)中國(guó)農(nóng)業(yè)水問題研究中心，北京 00083)

0 引 言

石羊河流域地處西北干旱內(nèi)陸區(qū)，那里地表水資源極其短缺，而農(nóng)業(yè)灌溉占了總水量很大一部分。對(duì)于地表水資源嚴(yán)重短缺的干旱地區(qū)，地下水是農(nóng)業(yè)灌溉的重要水源。然而，近年來由于當(dāng)?shù)厝藗儗?duì)地下水資源的過度開采利用，使當(dāng)?shù)氐叵滤V化度呈現(xiàn)逐年增加的趨勢(shì)，其中石羊河下游民勤縣的湖區(qū)地下水礦化度高達(dá)3～10 g/L，并且每年以0.12 g/L的速度遞增[1,2]。為了彌補(bǔ)緊缺的淡水資源，使干旱地區(qū)的農(nóng)業(yè)生產(chǎn)能夠發(fā)展穩(wěn)定，其中解決水資源緊缺問題的有效方法可以采用咸水資源進(jìn)行灌溉[3,4]。利用咸水資源進(jìn)行農(nóng)田灌溉一方面能夠減少由干旱對(duì)農(nóng)作物產(chǎn)量造成的損失，另一方面也會(huì)把鹽分帶入到土壤，會(huì)使農(nóng)田土壤環(huán)境發(fā)生變化，特別是在土壤耕作層中會(huì)累積大量的鹽分，進(jìn)而會(huì)導(dǎo)致土壤次生鹽堿化現(xiàn)象的發(fā)生[5]。如果長(zhǎng)時(shí)期利用咸水進(jìn)行農(nóng)田灌溉，會(huì)使大量的鹽分富集在作物根系區(qū)附近，受鹽分脅迫作用，會(huì)抑制農(nóng)作物根系吸水，也會(huì)一定程度上影響農(nóng)作物的生長(zhǎng)和產(chǎn)量。因此，在研究區(qū)開展鹽分脅迫條件下制種玉米根系吸水規(guī)律的研究，對(duì)于提高當(dāng)?shù)叵趟Y源的利用效率、制定合理的灌溉施肥方案和改善農(nóng)田水土環(huán)境具有重要的理論意義和實(shí)用價(jià)值[6]。長(zhǎng)時(shí)期利用咸水資源進(jìn)行農(nóng)田灌溉會(huì)導(dǎo)致土壤累積過量的鹽分，且鹽分主要累積在作物的根系區(qū)附近，由于溶質(zhì)勢(shì)阻礙作物根系吸水，從而會(huì)影響作物的正常生長(zhǎng)[7]。國(guó)內(nèi)外研究人員就鹽分脅迫對(duì)作物根系吸水的影響進(jìn)行過研究。Feng等針對(duì)不同灌水礦化度和不同灌水間隔時(shí)間的農(nóng)田土壤鹽分分布采用ENVIRO-GRO模型進(jìn)行數(shù)值模擬，研究表明咸水灌溉下根系深度對(duì)鹽分分布具有重要的影響[8]。羅長(zhǎng)壽等采用數(shù)值迭代反求方法對(duì)苜蓿的根系吸水規(guī)律進(jìn)了模擬和分析，研究表明土壤累積的鹽分會(huì)明顯降低苜蓿的根系吸水速率，當(dāng)土壤溶液的電導(dǎo)率約為5 dS/m左右時(shí)，對(duì)苜蓿的根系吸水具有較大的影響[9]。國(guó)內(nèi)外已做過大量關(guān)于鹽分脅迫條件下作物根系吸水規(guī)律的研究，其中應(yīng)用較普遍的是由Feddes提出的根系吸水模型[10]。

本文通過設(shè)置不同灌水礦化度的咸水灌溉試驗(yàn)，引入Feddes提出的根系吸水模型，根據(jù)田間試驗(yàn)實(shí)測(cè)數(shù)據(jù)，推算該模型的各個(gè)參數(shù)，建立鹽分脅迫條件下制種玉米根系吸水模型，為該研究區(qū)研究制種玉米生長(zhǎng)條件下土壤水鹽運(yùn)動(dòng)規(guī)律奠定基礎(chǔ)。

1 材料與方法

1.1 試驗(yàn)區(qū)位置與土壤物理參數(shù)

田間試驗(yàn)于2013年4-10月在甘肅省武威市涼州區(qū)的中國(guó)農(nóng)業(yè)大學(xué)石羊河試驗(yàn)站(102°52′E，37°52′N，海拔1 581 m)內(nèi)進(jìn)行，該試驗(yàn)站位于石羊河流域中游區(qū)域。田間試驗(yàn)在非稱重式蒸滲儀中進(jìn)行，修建測(cè)坑小區(qū)6個(gè)，每個(gè)測(cè)坑面積為6.66 m2(3.33 m×2 m)，深為3 m。測(cè)坑小區(qū)0～60 cm土壤為砂壤土(國(guó)際制)，其粒徑組成為：砂粒(0.05～1.00 mm)含量(體積分?jǐn)?shù)，下同)為58.90%，粉粒(0.005～0.05 mm)含量為29.02%，黏粒(﹤0.005 mm)含量為11.59%(采用馬爾文MS2000激光粒度分析儀測(cè)定)。測(cè)坑小區(qū)0～60 cm土壤容重為1.5 g/cm3，飽和含水率為0.37 cm3/cm3(體積含水率，下同)，田間持水率為0.30 cm3/cm3。

1.2 試驗(yàn)設(shè)計(jì)

田間試驗(yàn)根據(jù)石羊河流域上游、中游、下游典型地區(qū)地下水的礦化度設(shè)置3種鹽分處理即s1、s2和s3處理，三者分別代表為灌水礦化度為0.71 g/L的淡水、3 g/L的微咸水和6 g/L的咸水。由于受試驗(yàn)地條件所限，每個(gè)處理2個(gè)重復(fù)，共有6個(gè)試驗(yàn)小區(qū)，試驗(yàn)采用裂區(qū)排列方式布置。根據(jù)當(dāng)?shù)刂品N玉米灌溉經(jīng)驗(yàn)，制種玉米生育期內(nèi)灌水5次分別為：6月8日(拔節(jié)期)、6月30 日(抽雄期)，7月20日(開花期)、8月10日(灌漿期)和8月29日(乳熟期)，設(shè)置灌溉制度(表1所示)。試驗(yàn)使用的淡水通過水泵抽取當(dāng)?shù)氐叵滤@得，使用的咸水是根據(jù)當(dāng)?shù)氐叵滤瘜W(xué)組成，采用質(zhì)量比為2∶2∶1的NaCl、MgSO4和CaSO4配置的溶液。供試制種玉米的品種為富農(nóng)340號(hào)，2013年4月20日播種，9月13日收獲，全生育期146 d。測(cè)坑小區(qū)內(nèi)的制種玉米按父本與母本1∶7的比例方式種植，種植密度為每小區(qū)56株。其他農(nóng)藝措施均與當(dāng)?shù)厍闆r保持一致。

表1 各處理灌溉制度Tab.1 Scheme of irrigation under different treatments

1.3 測(cè)定項(xiàng)目與方法

選擇制種玉米拔節(jié)期(2013年6月23日-6月29日)為研究根系吸水時(shí)段1，抽雄期(2013年7月12日-7月19日)為研究根系吸水時(shí)段2。分別在時(shí)段1和時(shí)段2初期和末期通過土鉆獲取土樣，每個(gè)處理土層分為3層，分別為[0,20]，(20,40]，(40,60] cm，每次取樣完后回填鉆孔并做標(biāo)記，采用烘干法測(cè)定土壤含水率。采用電導(dǎo)法測(cè)定土壤溶液電導(dǎo)率EC1∶5，并根據(jù)已有資料的相關(guān)公式(S=0.027 5EC1∶5+0.136 6)把EC1∶5轉(zhuǎn)換成土壤全鹽量S[11]。分別在時(shí)段1、2每隔1 d利用鋼卷尺測(cè)量制種玉米葉子的長(zhǎng)(L)與寬(B)，并采用估算公式得到葉面積指數(shù)：

LAI=(K×L×B)∕A

式中：K為葉片擬合系數(shù)，玉米取0.75；A為植株葉片所占地面面積。

采用口徑為8 cm的根鉆，利用“十字法”(每株玉米取5鉆，每鉆分3層，每層20 cm，取至60 cm可以把制種玉米大部分根系取盡)取樣，獲取制種玉米根長(zhǎng)資料，并經(jīng)過根系掃描儀掃描和采用WinRHIZO根系分析軟件獲得根長(zhǎng)密度的分布資料。土面蒸發(fā)量利用自制Microlysimeter(MLS)直接測(cè)定，在研究制種玉米根系吸水時(shí)間段內(nèi)埋設(shè)MLS，每個(gè)處理3個(gè)重復(fù)。MLS內(nèi)徑10 cm，高度20 cm，下端采用不封底的形式，呈圓柱狀，采用PVC材料制成。每天早晨8∶00時(shí)刻對(duì)MLS進(jìn)行稱重，相鄰兩天的重量之差即為前一天的土面蒸發(fā)量。同時(shí)收集好試驗(yàn)站制種玉米全生育期內(nèi)的氣象資料。

1.4 根系吸水模型

Feddes提出的根系吸水模型為：

S(x,t)=β(φ0)Smax

(1)

式中：S(x,t)為僅在鹽分脅迫下作物根系吸水速率，cm3/(cm3·d)；β(φ0)為鹽分脅迫修正因子，φ0為土壤水滲透勢(shì)，cm；Smax為最大根系吸水速率，cm3/(cm3·d)，表示在最優(yōu)土壤水分(水分、養(yǎng)分充足，無鹽分脅迫)條件下的根系吸水速率。

1.4.1最大根系吸水速率Smax的計(jì)算方法

前人研究表明可以采用相對(duì)根長(zhǎng)密度Lnrd(x,t)和潛在蒸騰速率Tp(cm/d)來計(jì)算最大根系吸水速率[12]。即：

(3)

上述兩式中：Lr為最大扎根深度，cm；Ld(x,t)為實(shí)測(cè)根長(zhǎng)密度，cm/cm3。

1.4.2鹽分脅迫修正因子β(φ0)的計(jì)算方法

Van Genuchten等采用一種非線性方程，利用S形曲線函數(shù)來表示鹽分脅迫修正因子[13]。

(4)

式中：φ050為根系吸水鹽分脅迫修正因子等于0.5時(shí)所對(duì)應(yīng)的滲透勢(shì)的值，cm；p為擬合參數(shù)。

2 結(jié)果與分析

2.1 根系吸水模型參數(shù)的推算

2.1.1最大根系吸水速率Smax的計(jì)算

選擇根系吸水時(shí)段2，由于s1處理為淡水充分灌溉且灌溉水量為當(dāng)?shù)刂品N玉米的騰發(fā)量，所以可以認(rèn)為s1處理的根系吸水速率為最大根系吸速率。

潛在蒸騰速率Tp由如下公式[14]估算：

Tp=ET0(1-e-kLAI)

(5)

式中：ET0為參考作物蒸散量，cm，由FAO推薦的Penman-monteith公式[15]計(jì)算，計(jì)算結(jié)果為0.412 cm/d；LAI為葉面積指數(shù)，取值為3.023 cm2/cm2；k是植物冠層輻射衰減系數(shù)，無量綱，對(duì)于玉米而言，k常取0.4[16]。計(jì)算結(jié)果得Tp=0.29 cm/d。

實(shí)測(cè)根長(zhǎng)密度與土層深度的關(guān)系如圖1所示。經(jīng)線性插值擬合后得到Ld(x,t)=-0.011 9x+1.029 2(x為實(shí)測(cè)根長(zhǎng)深度)，制種大扎根深度為100 cm。

圖1 s1處理根長(zhǎng)密度分布Fig.1 Root length density distribution of s1 treatment

由式(2)，(3)計(jì)算得：

Lnrd(x,t)=-0.027 4x+2.37

Smax=-0.795×10-4x+0.687×10-2

2.1.2鹽分脅迫修正因子β(φ0)的計(jì)算

選擇根系吸水時(shí)段2，s2處理的灌溉水礦化度為3 g/L，為充分供水鹽分脅迫處理，計(jì)算相應(yīng)的蒸騰強(qiáng)度：

(7)

式中：Ta為s2處理的實(shí)際蒸騰強(qiáng)度，cm/d；Ta=ETp-Ep，ETp為潛在蒸散強(qiáng)度，可由Penman-monteith公式求得；Ep為潛在土面蒸發(fā)強(qiáng)度，cm/d，可由實(shí)測(cè)的土面蒸發(fā)強(qiáng)度表示。計(jì)算得到：ETp=0.40 cm/d；Ep=0.12 cm/d；所以Tp=0.28 cm/d。

理論上，β(φ0)應(yīng)為S與Smax的比值，由于β(φ0)無法實(shí)測(cè)，因此根據(jù)上式將Ta與Tp的比值近似為鹽分脅迫因子β(φ0)的比值。則：

(8)

φ0為該時(shí)段平均土壤水滲透勢(shì)，cm，由以下公式[5]計(jì)算：

cw=6.7×10-4EC1∶5

(9)

(10)

φ0=-5.63×105cs

(11)

式中：cw為土壤浸出液鹽分質(zhì)量濃度，g/cm3；EC1∶5為土壤浸出液電導(dǎo)率，dS/m；cs為土壤溶液鹽分質(zhì)量濃度，g/cm3；v為土壤浸提液的體積，cm3；γ為土壤容重，g/cm3；m為風(fēng)干土的質(zhì)量，g；θ為體積含水率，cm3/cm3。

由式(9)-(11)得到計(jì)算根系吸水階段的φ0=-4 036.52 cm。則由以上各式可以計(jì)算出p=5.49，取φ050=-7 600 cm[17]。因此，鹽分脅迫修正因子為：

2.2 根系吸水模型參數(shù)的驗(yàn)證

分別選擇時(shí)段1中s2、s3處理和時(shí)段2中s3處理對(duì)鹽分脅迫修正因子β(φ0)的參數(shù)p進(jìn)行驗(yàn)證。其驗(yàn)證結(jié)果見表2。

模型參數(shù)p的推算和驗(yàn)證過程時(shí)的吻合程度可以采用均方根誤差(RMSE)和平均相對(duì)誤差(MRE)來評(píng)價(jià)：

(13)

式中：Oi表示模型參數(shù)推算值；Pi表示模型參數(shù)的驗(yàn)證值。

表2 參數(shù)p推算與驗(yàn)證時(shí)的RMSE和MRE值Tab.2 Root mean square error and mean relative error values in calculation and validation of the parameter (p)

由表2可知，時(shí)段1中s2、s3處理與時(shí)段2中s3處理β(φ0)的RMSE<0.10，MRE<10%，在允許的誤差范圍20%之內(nèi)，因此鹽分脅迫修正因子β(φ0)=0.97是合理的，從而鹽分脅迫修正因子的參數(shù)p為5.49。所以鹽分脅迫條件下制種玉米根系吸水模型為：

3 結(jié) 語

在石羊河流域開展制種玉米生長(zhǎng)條件下的咸水灌溉試驗(yàn)，引入Feddes 提出的根系吸水模型。選擇制種玉米抽雄期為研究制種玉米根系吸水時(shí)段1，設(shè)置淡水充分灌溉試驗(yàn)，通過計(jì)算相對(duì)根長(zhǎng)密度和潛在蒸騰速率，推算得到最優(yōu)灌水條件下的最大根系吸水速率；設(shè)置咸水充分灌溉試驗(yàn)，通過計(jì)算蒸騰強(qiáng)度和土壤滲透勢(shì)，推算得到咸水灌溉條件下的鹽分脅迫修正因子。

選擇制種玉米拔節(jié)期為研究制種玉米根系吸水時(shí)段2，對(duì)鹽分脅迫修正因子的參數(shù)p進(jìn)行驗(yàn)證，參數(shù)p驗(yàn)證和推算過程時(shí)的RMSE<0.10，MRE<10%，在允許的誤差范圍之內(nèi)。建立了鹽分脅迫條件下制種玉米的根系吸水模型，為研究區(qū)研究制種玉米咸水灌溉條件下土壤水鹽運(yùn)動(dòng)規(guī)律奠定基礎(chǔ)。

□

[1] 黃翠華，薛 嫻，彭 飛，等. 不同礦化度地下水灌溉對(duì)民勤土壤環(huán)境的影響[J]. 中國(guó)沙漠，2013，33(2)：590-595.

[2] 袁成福，馮紹元，蔣 靜，等. 咸水非充分灌溉條件下土壤水鹽運(yùn)動(dòng)SWAP模型模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2014，20(30)：72-81.

[3] 王衛(wèi)光，王修貴，沈榮開，等. 河套灌區(qū)咸水灌溉試驗(yàn)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2004，20(5)：92-96.

[4] Ahmed B A O, Inoue M, Moritani S. Effect of saline water irrigation and manure application on the available water content, soil salinity, and growth of wheat[J]. Agricultural Water Management, 2010,97:165-170.

[5] 王利春，石建初，左 強(qiáng)，等. 鹽分脅迫條件下冬小麥根系吸水模型的構(gòu)建于驗(yàn)證[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2011，1(27)：112-117.

[6] 高曉瑜，霍再琳，馮紹元. 水鹽脅迫條件下作物根系吸水模型研究進(jìn)展及展望[J]. 中國(guó)農(nóng)村水利水電，2013，(1)：45-48.

[7] N K, Van Hoorn J W, A H, et al. Comparison of corn yield response to plant water stress caused by salinity and by drought [J]. Agricultural Water Management., 2004,(65):95-101.

[8] Feng G L, Meiri A, J L. Evaluation of a model for irrigation management under saline condition: Ⅱ. Salt distribution and rooting pattern effects[J]. Soilence Society of America Journal, 2003,67(1):77-80.

[9] 羅長(zhǎng)壽，左 強(qiáng)，李保國(guó)，等. 鹽分脅迫條件下苜蓿根系吸水特性的模擬于分析[J]. 土壤通報(bào)，2001，(S1)：81-84.

[10] Feddes R A, Kowalik P J, Zarandy H. Simulation of field water use and crop yield[M]. Wageningen: Pudoc for Centre for Agricultural Publishing and Documentation, 1978.

[11] 蔣 靜. 石羊河流域咸水非充分灌溉農(nóng)田土壤水鹽運(yùn)移試驗(yàn)與模擬研究[D]. 北京：中國(guó)農(nóng)業(yè)大學(xué)，2011.

[12] Wu J Q, Zhang R D and Gui S X. Modeling Soil water movement with water uptake by roots [J]. Plant and soil, 1999, 215(1):7-17.

[13] Van Genuchten, M T, Hoffman G J. Analysis of crop salt tolerance data[J]. Ecological Studies: Analysis and Synthesus, 1984, 51:258-271.

[14] Prasad R. A liner root water uptake model [J]. Journal of Hydrology, 1988,99(3/4):297-306.

[15] Allen R G, Pereira L S, Raes D, et al. Crop evapotranspiration guidelines for computing crop water requirement[M]. FAO Irrigation and Drainage Food and Agricultural Organization of the United Nations, Rome, Italy, 1998:65-88.

[16] 劉群昌，謝森傳. 華北地區(qū)夏玉米田間水分轉(zhuǎn)化規(guī)律研究[J]. 水利學(xué)報(bào)，1998，(1)：62-68.

[17] Wang Qingming, Huo Zailin, Feng Shaoyuan, et al. Comparison of spring maize root water uptake models under and salinity stressvalidated with field experiment data[J]. Irrigation and Drainage, 2015,5(64):669-682.