基于双孔隙模型研究刺槐林地和苹果地水量平衡要素的差异性.pdf

第3 0卷第5期2 0 2 3年1 0月水土保持研究R e s e a r c ho fS o i l a n dW a t e rC o n s e r v a t i o nV o l.3 0,N o.5O c t.,2 0 2 3 收稿日期:2 0 2 2-0 7-1 5 修回日期:2 0 2 2-0 8-0 2 资助项目:中国科学院战略性先导专项(B类)课题(X D B 4 0 0 0 0 0 0 0);中国科学院战略性先导科技专项课题(X D B 2 0 0 2 0 2 0 2)第一作者:郭天崎(1 9 9 8),女,山西忻州人,硕士研究生,主要从事土壤物理研究。E-m a i l:g u o t q 2 7 91 6 3.c o m 通信作者:黄明斌(1 9 6 8),男,湖北荆门人,研究员,主要从事生态水文和土壤物理方面研究。E-m a i l:h m b d n w s u a f.e d u.c nh t t p:s t b c y j.p a p e r o n c e.o r gD O I:1 0.1 3 8 6 9/j.c n k i.r s w c.2 0 2 3.0 5.0 1 7.郭天崎,黄明斌,朱谧远,等.基于双孔隙模型研究刺槐林地和苹果地水量平衡要素的差异性J.水土保持研究,2 0 2 3,3 0(5):1 3 8-1 4 5.GUOT i a n q i,HUAN G M i n g b i n,Z HU M i y u a n,e ta l.S t u d yo nt h eD i f f e r e n c eo fW a t e rB a l a n c eB e t w e e nR o b i n i aP s e u d o a c a c i aF o r e s t l a n da n dA p p l eO r c h a r dB a s e do nD u a l-P o r o s i t yM o d e lJ.R e s e a r c ho fS o i l a n dW a t e rC o n s e r v a t i o n,2 0 2 3,3 0(5):1 3 8-1 4 5.基于双孔隙模型研究刺槐林地和苹果地水量平衡要素的差异性郭天崎1,黄明斌2,朱谧远1,张 锐1(1.西北农林科技大学 资源环境学院,陕西 杨凌7 1 2 1 0 0;2.西北农林科技大学水土保持研究所黄土高原土壤侵蚀与旱地农业国家重点实验室,陕西 杨凌7 1 2 1 0 0)摘 要:目的 明晰刺槐林地和苹果园水量平衡要素的动态变化和差异,可为生态可持续发展提供理论指导。方法 利用H y d r u s-1 D中的双孔隙模型和标定的参数模拟了2 0 1 22 0 2 1年刺槐林地和苹果地的水量平衡各要素的动态变化过程。结果 优化期和验证期H y d r u s-1 D中双孔模型的R2分别为0.6 1,0.7 6,RM S E分别为0.0 2 6,0.0 4 8c m3/c m3,模拟期水量平衡误差分别为0.1 8%和0.0 4%;刺槐林地与苹果地的蒸腾量分别为(3 7 27 0)mm和(3 2 03 4)mm,蒸发量为(2 5 6 4 4)mm和(2 8 8 3 6)mm,土壤储水量变化量为(-0.6 1 9 6)mm和(-3.2 68 7)mm,深层渗漏量为(3 44 5)mm和(5 94 8)mm。结论 该模型可以较好地模拟刺槐林地和苹果地的水量平衡要素的动态变化;相比于与苹果地,刺槐林地的蒸腾量偏大,蒸发量和深层渗漏量偏小,S的波动更为明显。关键词:土地利用方式;双孔模型;再分布过程;土壤剖面水分含量;水量平衡中图分类号:S 1 5 2.5;S 1 6 1.4 文献标识码:A 文章编号:1 0 0 5-3 4 0 9(2 0 2 3)0 5-0 1 3 8-0 8S t u d yo nt h eD i f f e r e n c eo fW a t e rB a l a n c eB e t w e e nR o b i n i aP s e u d o a c a c i aF o r e s t l a n da n dA p p l eO r c h a r dB a s e do nD u a l-P o r o s i t yM o d e lGUOT i a n q i1,HUANG M i n g b i n2,Z HU M i y u a n1,Z HANGR u i1(1.C o l l e g eo fN a t u r a lR e s o u r c e sa n dE n v i r o n m e n t,N o r t h w e s tA&FU n i v e r s i t y,Y a n g l i n g,S h a a n x i7 1 2 1 0 0,C h i n a;2.S t a t eK e yL a b o r a t o r yo fS o i lE r o s i o na n dD r y l a n dF a r m i n go nt h eL o e s sP l a t e a u,I n s t i t u t eo fS o i la n dW a t e rC o n s e r v a t i o n,N o r t h w e s tA&FU n i v e r s i t y,Y a n g l i n g,S h a a n x i7 1 2 1 0 0,C h i n a)A b s t r a c t:O b j e c t i v eT h es t u d ya i m st oc l a r i f yt h ed y n a m i cc h a n g e sa n dd i f f e r e n c e si n w a t e rb a l a n c ee l e m e n t sb e t w e e nR o b i n i ap s e u d o a c a c i aw o o d l a n da n da p p l eo r c h a r d,a n d t op r o v i d e t h e o r e t i c a l g u i d a n c e f o re c o l o g i c a l s u s t a i n a b l ed e v e l o p m e n t.M e t h o d sT h i s s t u d yu s e d t h e d u a l p o r em o d e l a n d c a l i b r a t e dp a r a m e t e r s i nH y-d r a s-1 Dt os i m u l a t e t h ed y n a m i cc h a n g e so fv a r i o u se l e m e n t so fw a t e rb a l a n c e i nR o b i n i ap s e u d o a c a c i aw o o d-l a n da n da p p l eo r c h a r df r o m2 0 1 2t o2 0 2 1.R e s u l t sT h eR2a n dRM S Eo f t h e t w oh o l em o d e l i n t h eo p t i m i-z a t i o na n dv a l i d a t i o np e r i o d so fH y d r a s-1 Dw e r e0.6 1a n d0.7 6,r e s p e c t i v e l y,a n d0.0 2 6a n d0.0 4 8c m3/c m3,r e s p e c t i v e l y.T h ew a t e rb a l a n c e e r r o r s i n t h e s i m u l a t i o np e r i o dw e r e 0.1 8%a n d 0.0 4%,r e s p e c t i v e l y;T h e t r a n s p i-r a t i o no fR o b i n i ap s e u d o a c a c i af o r e s t l a n da n da p p l eo r c h a r dw e r e(3 7 27 0)mma n d(3 2 03 4)mm,r e s p e c t i v e l y,a n d t h e e v a p o r a t i o nw a s(2 5 64 4)mma n d(2 8 83 6)mm.V a r i a t i o no f s o i lw a t e r s t o r a g ew a s(-0.6 1 9 6)mma n d(-3.2 6 8 7)mm,a n dd e e pl e a k a g ew a s(3 44 5)mma n d(5 94 8)mm.C o n c l u s i o nT h em o d e l c a nb e t t e r s i m u l a t e t h ed y n a m i c c h a n g e so fw a t e rb a l a n c ee l e m e n t s i nR o b i n i ap s e u d o a c a c i aw o o d-l a n da n da p p l eo r c h a r d.C o m p a r e dw i t ha p p l e f i e l d s,R o b i n i ap s e u d o a c a c i af o r e s th a sa l a r g e r t r a n s p i r a t i o n,as m a l l e re v a p o r a t i o na n dd e e p e r l e a k a g e,a n dam o r es i g n i f i c a n t f l u c t u a t i o n i ns o i lw a t e rs t o r a g e.K e y w o r d s:l a n du s e t y p e;d u a l-p o r o s i t ym o d e l;r e d i s t r i b u t i o np r o c e s s;s o i lw a t e rc o n t e n t;w a t e rb a l a n c e 干旱和半干旱的黄土高原地区植被覆盖度低,水土固持能力差1。为了保持水土、改善生态环境,2 0世纪9 0年代实施了“退耕还林(草)”2的治理措施,导致人工刺槐林和苹果园种植面积大幅度增加3-4。然而,二者均为深根系、高耗水植物5-6,强烈的蒸腾作用会过量消耗深层土壤储水量,导致土壤干燥化,恶化土壤水分生态环境,造成土壤环境恶化,最终引起植被退化7-1 2。因此,明晰人工刺槐林和苹果园耗水强度和水量平衡各要素的动态变化对生态可持续发展具有重要意义。大孔隙通常是指直径位于7 550 0 0m的孔隙1 3,主要是由植物根系、土壤动物、土壤团聚体和土壤裂隙等形成的1 4。大孔隙使土壤导水率显著增加,入渗和再分布过程加快,出现优势流,对土壤剖面水分再分布和水平衡要素产生显著影响。研究表明,人工林增加和生态环境恢复有利于土壤大孔隙形成1 5,因此,研究人工刺槐林和苹果树水分消耗和水分平衡要素动态变化不能忽略大孔隙的影响。例如,M o n i k a等1 6应用H y d r u s-1 D模拟农田灌溉结果表明,考虑优先流的双孔隙模型能更加准确地模拟土壤水分动态过程;L u c a s等1 7人利用H y d r u s-1 D中的单孔和双孔隙模型对巴西半干旱地区水文过程进行模拟,得到双孔隙模型模拟结果更可靠。但是,目前考虑优势流的影响,利用双孔隙模型研究人工刺槐林和苹果地水量平衡差异性的报道还很少。为此,本文以陕西省中国科学院长武黄土高原农业生态试验站附近的刺槐林和苹果地为研究对象,选择H y d r u s-1 D中的双孔水流模型研究优势流存在条件下刺槐林和苹果园水量平衡要素的动态变化和差异,其中,双孔模型中的水力参数用再分布过程和水分交换过程中测定的含水量动态来标定和检验,以期为刺槐林和苹果园的可持续发展提供相应的理论依据。1 试验材料与方法:1.1 研究区概况试验区位于中国科学院长武黄土高原农业生态试验站附近,隶属于陕西省咸阳市长武县(1 0 7.8 E,3 5.2 N),平均海拔为11 8 6.7m,属于暖温带半湿润大陆性季风气候。年平均降水量5 8 4mm,年平均气温9.1,无霜期1 7 1 D,地下水位约5 08 0m。主要土壤类型为黑垆土,母质是马兰黄土。1.2 田间试验1.2.1 刺槐林和苹果地样地布设及土壤剖面水分含量监测 2 0 2 1年5月在试验站附近选择具有代表性的刺槐林地和苹果地。在2块样地中分别选取代表性的位置于安装3根5m和1根3m的中子管,其中3根5m中子管用于土壤水分动态监测,1根3m中子管用于后期土壤水分入渗和再分布试验。同时用铝盒按照01 0 0c m隔1 0c m,1 0 03 0 0c m隔2 0c m取土样,用烘干法测定土壤质量含水量;除此之外,在每块样地中选择代表性区域挖掘5 0c m深的土壤剖面,按照1 0c m间隔划分,共5层,利用环刀取各层原状土,每层3个重复,共1 5个。用于测定土壤容重1 8、饱和含水量1 9,计算出总孔隙度2 0(表1)。假设5 03 0 0c m土壤质地均匀2 1。51 1月用中子仪进行土壤剖面水分含量动态监测,时间间隔为1 01 5d,测定的土层间隔为01 0 0c m隔1 0c m,1 0 03 0 0c m隔2 0c m。1.2.2 入渗与再分布试验 2 0 2 1年7月,利用双环入渗仪(J Z-I N 1 2)进行入渗与再分布试验。双环入渗仪由双环入渗桶和马氏瓶组成,其内环直径为3 0c m,外环直径为6 0c m,与之匹配的马氏瓶容量分别为3L,1 0L。入渗前,根据测定的土壤容重计算出总孔隙度。为确保入渗结束后5 0c m剖面的土壤含水量接近饱和含水量,因此利用05 0c m剖面的总孔隙度计算出刺槐林地和苹果地的总灌溉量分别为1 2 5L和1 0 0L。然后以3m长的中子管为圆心放置双环入渗桶,外环桶壁用土壤填充,确保无水分外流。然后,连接好对应的马氏瓶,开始注水于两环中,打开马氏瓶,内外环水位相同均控制在5c m。灌水结束后,待地表无水层,即入渗过程结束,记录确切的入渗水量和入渗时间,开始水分再分布过程。按照入渗结束后1h,3h,6h,1 2h,2 4h,4 8h,9 6h和1 9 2h用中子仪测定土壤剖面水分含量,其中01 0 0c m隔1 0c m,1 0 03 0 0c m隔2 0c m。1.2.3 叶面积指数的测定及计算 利用L A I-2 2 0 0植物冠层分析仪,每隔7d分别测定刺槐林地和苹果地的叶面积指数(L e a fa r e ai n d e x,L A I)。将最大L A I和刺槐林、苹果树从发芽至落叶天数定为1,以相对L A I(01)和相对时间(01)为参数进行L A I931第5期 郭天崎等:基于双孔隙模型研究刺槐林地和苹果地水量平衡要素的差异性动态模拟。采用C u r v eE x p e r t 1.4曲线拟合软件2 2最终筛选,建立适合刺槐林地(R2=0.9 9)和苹果地(R2=0.9 4)的相对化L A I动态模型:Lf o r e s t=4.0 1t3-9.8 3t2+5.8 4t-0.0 2(1)La p p l e=5.3 3t3-1 1.0 2t2+5.7 5t-0.0 1(2)L A I=LRL A Im a x(3)式中:Lf o r e s t和La p p l e为刺槐林地和苹果地的相对L A I值(LR);t为相对生育时间;L A Im a x为最大叶面积指数。查阅历史资料和实际测定得到2 0 1 02 0 2 1年长武试验站附近刺槐林地和苹果地的L A Im a x,建立L A Im a x的年际变化回归方程,其中刺槐林地(R2为0.9 8):L A Im a x=0.3 0T-6 0 9.5 8(4)苹果地(R2为0.9 8):L A Im a x=0.0 2 3 5T2-9 4.6 4 4T+9 5.3 6 1(5)式中:T为年份,将(4)、(5)式代入(3)式得到每日的叶面积指数。1.3 模型评价与运用1.3.1 双孔水流模型1 7 H y d r u s-1 D模型中的双孔隙模型是基于R i c h a r d s方程和简单的质量平衡方程的组合来描述大孔隙内的水分运动及基质中水分的动态交换过程:m t=aK(h)(h z+1)-Sm-w(6)m t=-Si m+w(7)式中:m,i m分别表示大孔隙域和基质域;=m+mi m为体积含水量(c m3/c m3);Si m,Sm为两个区域的根系吸水速率(1/d);K(h)为非饱和导水率(c m/d);w为大孔隙与基质之间水交换的交换速率(1/d),用如下公式计算:w=i m t=Sme-Si me(8)式中:i m为基质含水量(c m3/c m3);为一级交换速率系数(1/天);Sme,Si me分别为大孔隙域和基质域相对含水量。土壤水分特征曲线和非饱和导水率(K)与土壤含水量的关系由v a nG e n u c h t e n-M u l a e m方程拟合2 3:(h)=r+s-r1+hn()m h0s h0(9)K(h)=KsSe1-1-S1/me()h刺槐林地(1.2 6g/c m3),3 03 0 0c m则苹果地(1.3 5g/c m3)与刺槐林地(1.3 4g/c m3)近似。同一土地利用方式下不同土层的总孔隙度也存在差异,其中刺槐林地表层(03 0c m)的总孔隙度大于深层(3 03 0 0c m),而苹果地表层和深层土壤的总孔隙度相近。测定的饱和含水量的剖面变化规律与总孔隙度一致。表1 刺槐林地和苹果地的土壤物理性质土地利用类型土层土层深度/c m容重/(gc m-3)总孔隙度f饱和含水量/(c m3c m-3)刺槐林地103 01.2 60.0 0 90.5 20.0 0 30.5 00.0 1 123 03 0 01.3 40.0 3 60.4 90.0 1 30.4 50.0 1 8苹果地103 01.3 20.0 0 60.5 00.0 0 20.4 90.0 1 423 03 0 01.3 50.0 1 90.4 90.0 0 70.4 70.0 5 2注:数据组成为平均值标准误差。2.2 模型参数优化基于水分再分布过程土壤剖面含水量监测数据,采用H y d r u s-1 D对土壤再分布过程进行分层模拟,以优化刺槐林地和苹果地的土壤水力参数,结果见表2。从表2可知,刺槐林地和苹果地表层(03 0c m)土壤的大孔隙体积占比分别为2 4%和2 1%,深层(3 03 0 0c m)为7%和1 7%,即随土层深度的增加大孔隙逐步减少,这与根系剖面分布一致;较苹果地而言,刺槐林地表层(03 0c m)和深层(3 03 0 0c m)的大孔隙分布差异更大,表层大孔隙几乎为深层大 孔隙的4倍,且刺槐林地表层的水质交换速率(0.3 4 4)约是深层(0.0 5)的7倍,但苹果地整个土壤的水质交换速率无明显差异,均为0.0 1;此外,饱和导水率(Ks)均表现为随土层深度增加而降低,表层(03 0c m)约为深层(3 03 0 0c m)的2.3倍,二者相比,苹果地的Ks更大,约为刺槐林地的2.3倍。刺槐林地和苹果地的再分布过程模拟结果如图1所示。刺槐林地和苹果地均随时间产生相似的土壤水分动态变化。在3 0c m处由于土壤质地的变化土壤水分发生突然变化。1 9 2h的模拟值明显大于141第5期 郭天崎等:基于双孔隙模型研究刺槐林地和苹果地水量平衡要素的差异性实测值,其他模拟时间段内03 0c m过高的估计的土壤水分含量,可能是中子仪测定表层土壤水分时,中子扩散不稳定造成的3 6,3 03 0 0c m的模拟值与实测值较为一致。1 9 2h实测值与模拟值差异可能来源于入渗水分沿中子管的管壁快速渗漏,导致1 9 2h模拟值始终大于实测值。图2为模拟值与实测值的相关性分析,刺槐林和苹果地的实测值与模拟值的R2,RM S E分别为0.7 6,0.0 4 8c m3/c m3,0.7 6,0.0 3 0c m3/c m3。高的决定系数和低的RM S E值表明该优化参数较合理。表2 刺槐林地和苹果地的土壤水力参数土地利用类型土层土层深度/c mrsnKs/(c mh-1)r i ms i m刺槐林地103 000.1 1 90.0 0 3 62.1 6 31 4.0 00.0 7 80.3 8 10.3 4 423 03 0 000.0 3 00.0 0 4 92.1 5 06.6 30.0 7 80.4 2 00.0 5 0苹果地103 000.1 0 10.0 2 3 71.4 1 94 1.4 00.0 7 80.3 8 90.0 1 023 03 0 000.0 7 80.0 3 4 31.4 7 11 5.8 00.0 7 80.3 9 20.0 1 0图1 刺槐林地和苹果地再分布过程土壤剖面实测值与模拟值比较图2 刺槐林地和苹果地水分再分布过程土壤水分实测值与模拟值线性分析2.3 模型验证利用标定的水力参数(表2)模拟了2 0 2 1年51 1月在刺槐林地和苹果地土壤水分动态。整体上,模拟值与实测值较接近。其相关性分析如图3所示。刺槐林地和苹果地的R2分别为0.6 3,0.6 1,RM S E分别为0.0 2 6c m3/c m3,0.0 2 9c m3/c m3,误差在允许范围之内3 7。说明本文采用的双孔隙模型可以较好地预测该区域刺槐林地和苹果地的土壤水分动态变化。但是,在2 6 0d后该模型低估了1 6 0c m和3 0 0c m深度的土壤水分(图4),可能是模型忽略了入渗与再分布之间的滞后效应,使得模拟结果出现一定误差3 8。因为在第2 6 0天出现了一次大的降雨过程,剖面水分运动以入渗为主,而模型参数优化以再分布过程为主,滞后效应导致二者水分特征曲线的差异影响了模型精度。2.4 刺槐林地和苹果地水量平衡要素动态变化在干旱半干旱黄土高原区,降雨是唯一的水分来源1 1。水量平衡要素主要为植物蒸腾、土壤蒸发、土壤储水量变化量和深层渗漏。利用上述优化验证后的模型对刺槐林地和苹果地水量平衡各要素进行模拟,结果显示水量平衡平均相对误差()分别为0.1 8%和0.0 4%,误差较小,对模拟结果的影响可忽略不计,表明该模型可以较好地模拟刺槐林地和苹果地水量平衡各要素的动态变化。水量平衡各要素模拟结果(图5)显示,刺槐林地和苹果地的水量平衡各要素的动态变化规律存在明241 水 土 保 持 研 究 第3 0卷显差异。由图5可知,刺槐林地的蒸腾量大于蒸发量,两者总和为(6 2 97 1)mm,占到降雨量的9 5%。其中,蒸腾量为(3 7 27 0)mm,在平水年,占到降雨量的4 1%6 1%,在干旱年,可达到5 8%8 1%;蒸发量为(2 5 64 4)mm,在 平 水 年 占 到 降 雨 量 的2 3%4 3%,干旱年达3 7%6 6%。而苹果地的蒸腾量略大于蒸发量,两者总和是(6 0 85 8)mm,占到降雨量的9 2%。其中蒸腾量为(3 2 03 4)mm,在平水年,占降雨量的3 6%4 8%,在干旱年可达到4 7%6 6%;蒸发量为(2 8 8 3 6)mm,在平水年占到降雨量的3 0%4 4%,干旱年达到4 5%6 2%。相比于刺槐林地,苹果地的蒸腾量较小,主要原因是苹果种植密度低,叶面积指数小,模拟期间苹果地平均最大叶面积指数是2.4 1,刺槐林平均最大叶面积指数是3.6 3。图3 刺槐林地和苹果地土壤水分实测值与模拟值线性分析图4 刺槐林地和苹果地深度土壤水分模拟值与实测值比较 模拟期间,刺槐林地和苹果地的 S的波动趋势与降雨量基本一致。其中,刺槐林地的土壤储水量变化量为(-0.6 19 6)mm,苹果地为(-3.2 68 7)mm。从总趋势来看,两种林地土壤储水量均略有减小,主要由于在干旱年刺槐林和苹果树的蒸散发量大于降雨量,消耗土壤储水量,使其出现亏损,而丰水年降雨量虽然大于蒸散发量,土壤储水量得到了一定的补充,但多年的补充量小于亏损量。相对于苹果地,刺槐林地S值波动更为明显,模拟期间S的方差明显大于苹果地,主要原因是刺槐林种植密度大,叶面积指数大,干旱年土壤水分消耗大,而丰水年土壤水分得到补充3 9。此外,模拟期间,刺槐林地深层渗漏量为(3 44 5)mm,平水年占到降雨量的6%,干旱年达到4%。苹果地的深层渗漏量为(5 94 8)mm,平水年占到降雨量的1 0%,干旱年占到7%。降雨量大的年份深层渗漏均明显增加,尤其是2 0 2 1年,刺槐林地和苹果地的深层渗漏量达到1 6 4mm和1 8 7mm,分别占到降341第5期 郭天崎等:基于双孔隙模型研究刺槐林地和苹果地水量平衡要素的差异性雨量的1 8.1 6%和2 0.6 3%。相比于刺槐林地,模拟期间苹果地的深层渗漏量更大,主要由于苹果地深层(3 03 0 0c m)的大孔隙多,而大孔隙是降雨对深层渗漏主要的补给途径4 0。注:T r是植物蒸腾量;E是土壤蒸发量;D是深层渗漏量;S是土壤储水量变化量;P是降雨量。图5 2 0 1 2-2 0 2 1年刺槐林地和苹果地水量平衡要素动态变化3 结 论(1)利用干旱半干旱黄土高原区的刺槐林地和苹果地再分布过程和实测土壤剖面含水量数据对H y d r u s-1 D中的双孔隙模型参数进行标定和验证,R2在0.6 10.7 4,RM S E在0.0 2 60.0 4 8c m3/c m3;水量平衡在0.0 4%0.1 8%,误差在合理范围之内。因此该模型可以较好地用于模拟黄土高原区刺槐林地和苹果地的水量要素动态变化。(2)刺槐林地和苹果地的水量平衡要素主要受降雨量的影响,但水量平衡要素的动态变化存在明显差异。相比于苹果地,刺槐林地的蒸腾量更大,蒸发量和深层渗漏较小,S的波动更明显。参考文献:1 陈明玉,邵明安,李同川,等.黄土高原典型切沟土壤水分时空分布特征及其影响因素J.土壤学报,2 0 2 1,5 8(2):3 8 1-3 9 0.2 张琨,吕一河,傅伯杰,等.黄土高原植被覆盖变化对生态系统服务影响及其阈值J.地理学报,2 0 2 0,7 5(5):9 4 9-9 6 0.3 郭军庭,张志强,唐丽霞.晋西黄土区小流域典型植被水文功能优化研究J.水土保持通报,2 0 1 0,3 0(3):3 6-4 0.4 王艳萍,王力,韩雪,等.黄土塬区不同土地利用方式土壤水分 消 耗 与 补 给 变 化 特 征 J.生 态 学 报,2 0 1 5,3 5(2 2):7 5 7 1-7 5 7 9.5 白岗栓,邹超煜,杜社妮.渭北旱塬果园自然生草对土壤水分及苹果树生长的影响J.农业工程学报,2 0 1 8,3 4(3):1 5 1-1 5 8.6 安文明,梁海斌,王聪,等.黄土高原阴/阳坡向林草土壤水分随 退 耕 年 限 的 变 化 特 征 J.生 态 学 报,2 0 1 7,3 7(1 8):6 1 2 0-6 1 2 7.7 邵明安,贾小旭,王云强,等.黄土高原土壤干层研究进展与展望J.地球科学进展,2 0 1 6,3 1(1):1 4-2 2.8 于贵瑞,李文华,邵明安,等.生态系统科学研究与生态系统管理J.地理学报,2 0 2 0,7 5(1 2):2 6 2 0-2 6 3 5.9 郭慧敏,张彦军,刘庆芳,等.黄土高原半干旱区土壤呼吸对土地利用变化的响应J.自然资源学报,2 0 1 4,2 9(1 0):1 6 8 6-1 6 9 5.1 0 杨磊,张子豪,李宗善.黄土高原植被建设与土壤干燥化:问题与展望J.生态学报,2 0 1 9,3 9(2 0):7 3 8 2-7 3 8 8.1 1 程立平,王亚萍,齐光.黄土塬区苹果种植对深剖面土壤水分和地下水补给的影响J.地理研究,2 0 2 1,4 0(9):2 6 8 4-2 6 9 4.1 2 杨文治,田均良.黄土高原土壤干燥化问题探源J.土壤学报,2 0 0 4,4 1(1):1-6.1 3 S o i lS c i e n c eG l o s s a r yT e r m sC o mm i t t e e,S o i lS c i e n c eS o c i e t yo fAm e r i c a.G l o s s a r yo f s o i l s c i e n c e t e r m s2 0 0 8R.A S A-C S S A-S S S A,2 0 0 8.1 4 高朝侠,徐学选,赵娇娜,等.土壤大孔隙流研究现状与发展趋势J.生态学报,2 0 1 4,3 4(1 1):2 8 0 1-2 8 1 1.1 5 B a r g u sT o b e l l aA,R e e s eH,A l m a wA,e t a l.T h ee f f e c to f t r e e s o np r e f e r e n t i a l f l o wa n d s o i l i n f i l t r a b i l i t y i na n a g r o-f o r e s t r yp a r k l a n di ns e m i a r i d B u r k i n aF a s oJ.W a t e rR e s o u r c e sR e s e a r c h,2 0 1 4,5 0(4):3 3 4 2-3 3 5 4.1 6 M a r k o v iM,F i l i p o v iV,L e g o v iT,e t a l.E v a l u a t i o no fd i f f e r e n t s o i lw a t e rp o t e n t i a l b y f i e l dc a p a c i t y t h r e s h o l di nc o m b i n a t i o nw i t hat r i g g e r e di r r i g a t i o n m o d u l eJ.S o i l a n dW a t e rR e s e a r c h,2 0 1 5,1 0(3):1 6 4-1 7 1.1 7 A l c n t a r aLRPD,C o u t i n h oAP,N e t oS M DS,e ta l.M o d e l i n go f t h eH y d r o l o g i c a lP r o c e s s e s i nC a a t i n g a441 水 土 保 持 研 究 第3 0卷a n dP a s t u r eA r e a s i nt h eB r a z i l i a nS e m i-A r i dJ.W a-t e r,2 0 2 1(1 3):1 8 7 7.1 8 江胜国.国内土壤容重测定方法综述J.湖北农业科学,2 0 1 9,5 8(S 2):8 2-8 6,9 1.1 9 李映强,曾觉延.关于土壤通气孔隙概念及其测定方法的商榷J.土壤通报,1 9 9 0,2 1(2):8 7-8 8.2 0 邹文安,姜波,顾李华.土壤水分常数的测定J.水文,2 0 1 5,3 5(4):6 2-6 6.2 1 Y a n gF,H u a n g M B,L iC H,e ta l.C h a n g e si ns o i lm o i s t u r ea n do r g a n i cc a r b o nu n d e rd e e p-r o o t e dt r e e so fd i f f e r e n t s t a n da g e so nt h eC h i n e s eL o e s sP l a t e a uJ.A g r i c u l t u r e,E c o s y s t e m&E n v i r o n m e n t,2 0 2 2,3 2 8(1):1-1 4.2 2 张宾,赵明,董志强,等.作物高产群体L A I动态模拟模型的建立与检验J.作物学报,2 0 0 7,3 3(4):6 1 2-6 1 9.2 3 A g a r w a lP,S h a r m aPK.A n a l y s i so fc r i t i c a lp a r a m e-t e r so f f l o wa n ds o l u t et r a n s p o r t i np o r o u sm e d i aJ.W a t e rS c i e n c e&T e c h n o l o g y W a t e rS u p p l y,2 0 2 0,2 0(8):3 4 4 9-3 4 6 3.2 4 R a b b e lI,B o g e n a H,N e u w i r t hB,e ta l.U s i n gS a pF l o w D a t at oP a r a m e t e r i z et h eF e d d e s W a t e rS t r e s sM o d e l f o rN o r w a yS p r u c eJ.W a t e r,2 0 1 8,1 0(3):2 7 9,D o i:o r g/1 0.3 3 9 0/w 1 0 0 3 0 2 7 9.2 5 L iY,Y uY,S u nR,e t a l.S i m u l a t i o no f s o i lw a t e r d y-n a m i c s i nab l a c kl o c u s tp l a n t a t i o no nt h eL o e s sP l a t-e a u,W e s t e r n S h a n x iP r o v i n c e,C h i n aJ.W a t e r,2 0 2 1,1 3(9):1 2 1 3.2 6 L iBB,W a n gY Q,H i l lRL,e ta l.E f f e c t so fa p p l eo r c h a r d sc o n v e r t e df r o mf a r m l a n d so ns o i lw a t e rb a l-a n c e i nt h ed e e pl o e s sd e p o s i t sb a s e do nHY D RU S-1 Dm o d e lJ.A g r i c u l t u r e,E c o s y s t e m sa n dE n v i r o n m e n t,2 0 1 9,2 8 5(C):1 0 6 6 4 5.2 7 李蓝君,宋孝玉,夏露,等.黄土高原沟壑区典型造林树种蒸散发对气候变化的响应J.农业工程学报,2 0 1 8,3 4(2 0):1 4 8-1 5 9.3 8 Z h a oC Y,F e n gZD,C h e nG D.S o i lw a t e rb a l a n c es i m u l a t i o no f a l f a l f a