生态学报  2015, Vol. 35 Issue (2): 489-496

文章信息

龚雪文, 李仙岳, 史海滨, 田德龙, 李祯, 彭遵原
GONG Xuewen, LI Xianyue, SHI Haibin, TIAN Delong, LI Zhen, PENG Zunyuan
番茄、玉米套种膜下滴灌条件下农田地温变化特征
The interplanting between tomato and maize of soil temperature dynamics under mulched drip irrigation
生态学报, 2015, 35(2): 489-496
Acta Ecologica Sinica, 2015, 35(2): 489-496
http://dx.doi.org/10.5846/stxb201303240506

文章历史

收稿日期:2013-03-24
网络出版日期:2014-03-25
番茄、玉米套种膜下滴灌条件下农田地温变化特征
龚雪文1, 李仙岳1 , 史海滨1, 田德龙2, 李祯1, 彭遵原1    
1. 内蒙古农业大学水利与土木建筑工程学院, 呼和浩特 010018;
2. 水利部牧区水利科学研究所, 呼和浩特 010020
摘要:土壤水热状况是影响作物生长的重要因素,以膜下滴灌及套种双重效应为背景,针对番茄、玉米套种农田设置了高(番茄:180 mm,玉米:270 mm)、中(番茄:132 mm,玉米:202.5 mm)、低(番茄:84 mm,135 mm)3种不同的灌溉定额,采用定位监测法测定番茄带膜间裸地(A)、覆膜番茄行间(B)、番茄带与玉米带膜间裸地(C)、覆膜玉米行间(D)共4个位置的不同深度主根区含水率和地温动态变化,研究了地膜覆盖、土壤水分以及作物遮阴对套种农田不同行间,不同深度土壤温度的影响特征。结果表明:全生育期内,4个水平位置的地温大小顺序与土壤深度有关,10 cm以上的地温大小顺序为C > D > A > B,10 cm以下的地温大小顺序为D > B > A > C;从土壤温度场的分布情况来看,地温的最大变幅出现在C处,为8.10 ℃;最小变幅出现在B处,为4.71 ℃;含水率与地温的关系表现为:5 cm以上地温与含水率呈反比关系,20 cm以下地温与含水率呈正比关系。可见,套种覆膜滴灌农田合适的土壤水分状况将有利于根区良好水热的形成,从而更有利于作物的生长。研究结果将对西北地区滴灌农田管理有一定的参考价值。
关键词套种    膜下滴灌    地温    含水率    
The interplanting between tomato and maize of soil temperature dynamics under mulched drip irrigation
GONG Xuewen1, LI Xianyue1 , SHI Haibin1, TIAN Delong2, LI Zhen1, PENG Zunyuan1    
1. Water Conservancy and Civil Engineering College of Inner Mongolia Agricultural University, Huhhot 010018, China;
2. Institute of Water Resources for Pastoral Area, Ministry of Water Resources, Huhhot 010020, China
Abstract:It is very significant for maintaining the optimal soil moisture and soil temperature in the crop field, because that is the important foundation of developing an optimal irrigation system. The environment of reasonable soil moisture and temperature in the field will not only enhance the water use efficiency but also benefit crop growth and increase crop yield, especially for interplanting field with complex crop environments. In the experiments, the interplanting of tomato and corn under mulched drip irrigation were tested in Hetao Irrigation District in 2012. We set up three treatments which are the high irrigation quota (tomato: 180 mm, corn: 270 mm), medium irrigation quota (tomato: 132 mm, corn: 202.5 mm) and low irrigation quota (tomato: 84 mm, corn: 135 mm). The moisture content and soil temperature of different depths of root zone were measured in the middle of bare soil between tomatoes belt (A), the middle of mulched tomato lines (B), the middle of bare soil between tomatoes and corn belt (C), the middle of mulched corn lines (D) and crop growth and yield were observed in different growing seasons. The effects of plastic film, soil moisture and crop shade on soil temperature at different rows and depths of inter-planting were studied by the experiment. The results showed that it had apparent effect of mulch on heat preservation under 15 cm soil depths, but the soil temperatures of whole root zone were influenced by crop shade and soil moisture; The soil temperatures of four horizontal positions were significant variation for different the soil depths (5 cm, 10 cm, 15 cm, 20 cm, 25 cm), and the proper sequence of average soil temperature above 10 cm soil depth in different positions was: C > D > A > B, the one under 10 cm soil depth was: D > B > A > C; The study on distribution of the soil temperature in the interplanting field showed that the maximum variational amplitude of soil temperature appeared at the C position, and it would reach 8.10 ℃; the minimum variational amplitude of soil temperature appeared in the B position, and it was only 4.71 ℃. There was definitely some inverse correlation between the soil temperature and soil moisture above 5 cm soil depth, however direct proportion was observed below 20 cm soil depth. The study also found that there was strong temperature effect in the middle of lines between high and low crops (such as C position) during daytime, and the temperature can influence the crop of two sides, so mulched cultivation has obvious effects for restraining the loss of heat. It is thus clear that the optimal soil water and soil temperature are beneficial to form better soil moisture and temperature environment and high crop yield in the intercropping field under mulched drip irrigation. The study will be very meaningful for drip irrigation farm management in the Northwest Region.
Key words: intercropping    mulched drip irrigation    soil temperature    soil water    

套种不仅能提高单位面积的产量,而且能极大提高光、水、热等资源的利用效率[1, 2],据报道,在内蒙古河套灌区套种已超过总种植面积的70%,然而由于套种农田对太阳辐射的有效拦截和利用不同[3],导致套种农田不同行间位置的地温分布特征不同,而且覆膜后,塑料薄膜能够有效将太阳能转化为热能,有效调节作物根区环境[4, 5, 6]。因此将膜下滴灌技术与套种农田结合起来探讨作物耕层温度变化规律,对提高套种农田水热利用效率具有重要的现实意义。近些年针对套种农田采用膜下滴灌技术的研究主要集中在作物系数、水分利用效率、水分迁移过程以及产量品质的影响[7, 8, 9, 10]等方面,而针对水、热特性的研究主要以套种效应、覆膜效应、覆盖效应、土壤水分效应等单项研究为主,如套种可有效利用太阳辐射提高产量[3];覆膜可有效提高土壤耕层温度,促进作物生长[11];覆盖可抑制地温的大幅度波动[12];土壤含水量与土壤温度呈显著负相关[13]。另外,国内外学者还对覆膜条件或套种条件下的土壤热效应及地温的极值变化理论和模拟等方面做了研究[14, 15, 16, 17, 18, 19]。显然,这些研究主要是针对或覆膜或套种或滴灌等单项技术开展的,而对套种农田的不同行间位置地温分布特征尚未做深入研究。故本文主要针对套种及膜下滴灌耦合条件下多作物农田不同行间位置的地温变化规律展开研究,并分析了套种模式下覆膜效应及土壤水分效应对地温影响的综合效应,本研究对于套种农田和膜下滴灌耦合作用下水热理论发展具有一定的意义,并为该技术的进一步推广奠定了基础。

1 材料与方法 1.1 区域概况

2012年在内蒙古河套灌区磴口试验站(40.20 E,107.1 N)开展了番茄套种玉米的试验,该试验田地处干旱半荒漠草原地带,夏季多风干旱,冬季少雨寒冷,属中温带大陆性季风气候,海拔2004 m,年均气温7.2 ℃,年均日照时数2190 h,年均降水量198.3 mm,年均蒸发量2460.3 mm,年均风速2.74 m/s,多年平均无霜132 d。试验田土壤质地为夹砂粘土,1 m内土壤平均容重为1.50 g/cm3,凋萎系数为8%,田持为24.1%。

1.2 试验设计

供试作物为玉米(品种为中地77)和番茄(品种为屯河48),采用“一膜一管两行”的滴灌布设方式,采用4行番茄套2行玉米的种植模式。在当地灌水经验的基础上,试验设置了高、中、低3种不同的灌水处理,本试验考虑了番茄与玉米需水量不同的特点,在每个小区中分别对番茄和玉米设定了不同的灌水定额(表 1),采用精度为0.001 m3的水表(NSD6)控制灌水量,每个处理均设置了3个重复。

表1 全生育期番茄套种玉米的灌溉制度 Table 1 Irrigation scheduling for intercropping between tomato and maize during the whole growing season
处理号 Treatment小区面积/m2 Plot area 灌水次数 Irrigation frequency 灌水定额/mm Irrigation quota 灌溉定额/mm Irrigation quota
番茄Tomato玉米Corn 番茄Tomato玉米Corn 番茄Tomato玉米Corn 番茄Tomato玉米Corn
T1142.547.58922.530180270
T2142.547.58916.522.5132202.5
T3142.547.58910.51584135

每个处理均布设了4组精度为0.1 ℃的水银直角地温计(RM-004),作物种植方式及地温计的布设方式如图 1所示。

图 1 番茄、玉米套种方式及地温计与滴灌带布置图 Fig. 1 Latout of intercropping between tomato and maize and geothermometer and drip irrigation system
1.3 试验方法 1.3.1 气象资料测量

离试验区300 m的空旷地中设有自动气象站(HOBO),用于观测降水、气温、湿度、太阳辐射、大气压等。

1.3.2 土壤含水率测量

采用TRIME-PICO TDR和土钻相结合的方法,以TRIME为主,每隔3 d测量1次,定期用土钻取土烘干的方法加以校正,TDR与地温计位置相同(见图 1),采用定位监测法测量含水率,测量深度为0—100 cm。

1.3.3 地温观测时间与方法

在6—9月,利用直角地温计(RM-004)对5、10、15、20、25 cm处的土壤温度进行定点观测,测量时间与土壤含水率一致,每隔3 d测量1次,每天的观测时间为8:00—20:00,2 h测1次。

1.3.4 作物生长指标测量

用直尺测量玉米和番茄的株高;叶面积的测量:玉米叶面积=叶长×叶宽×0.75,番茄采用折算系数法计算[20],以上指标均在每个生育期内测3次,取平均值。

1.3.5 试验观测日期

玉米4月20日播种,9月10日收获,番茄5月5日移苗,8月10日采摘第1茬,8月25日采摘结束。试验观测时间从6月开始到9月结束,试验观测期内气温变化见图 2

图 2 试验观测期内的气温变化 Fig. 2 The temperature change during the experimental period
1.4 数据分析

本试验数据采用Spss 18.0和Excel 2003进行统计分析,采用Excel 2003和Surfer 8.0 软件绘制图形。

2 结果与分析 2.1 套种农田不同位置及不同深度地温日变化规律

在套种农田不同行间的不同位置土壤温度日变化亦有所不同。通过对比分析6—9月份地温数据可以看出,高水分处理T1各行间位置的地温变化过程(图 3)均是随着土层深度的增加地温日波动呈减小趋势,其中距地表5 cm处地温日波动最剧烈,距地表25 cm处地温日波动最平缓,距地表5、10、15、20、25 cm的日最高与最低地温差值分别为:12.72、8.78、6.09、4.44、2.96 ℃,存在极显著性差异(F=12.74005,P=0.007301)。

图 3 整个观测期内T1的日平均地温变化图 Fig. 3 The daily variation law of soil temperature during the whole observation period of T1

套种农田不同行间位置的地温受地膜覆盖和作物遮阴的影响较大。透明塑料薄膜对土壤的保温效应主要体现在土层15 cm以下,从表 2中可以看出,B、D位置的地温明显高于A、C位置,随着土层的加深,膜内与膜外平均地温的差值变化呈抛物线形,且均差在20 cm处达到最大值,可见在套种农田中,地膜可隔绝土壤与外界的水分交换和显热交换,使土壤的热通量在20 cm处达到最大。地面增温的唯一来源是太阳辐射,套种农田不同行间位置地温受作物遮阴的影响程度也不相同,这在5 cm以上土层表现的较为明显,从图 3中可以看出,距地表5 cm处的地温大小关系表现为C>D>A>B。这是由于C处位于高矮作物的行间,遮阴面积最小,日照时间最长,因此表层(5 cm处)的地温也最高;而B位于低矮作物行间,遮阴面积较大,地表所受太阳直射时间最短,因而地温也最低。可见,在立体种植农田中,高杆作物遮阴区的地温要远远高于矮杆作物遮阴区的地温;在15 cm以下,膜内的日平均地温要高于膜外。

表2 15 cm以下膜内与膜外的日平均地温 Table 2 Average day soil temperature of covered and uncovered treatments below 15 cm
层深/cm Layer depth膜内平均地温/℃ Average temperature in film mulching 膜外平均地温/℃ Average temperature in no mulching 膜内与膜外地温的均差/℃ Mean temperature difference of film mulching and no mulching
BD AC
1524.7625.0124.2323.630.96
2023.9324.2622.7222.261.61
2522.1022.0421.4921.580.53
2.2 套种农田不同处理及不同位置地温显著性分析

6—9月份是番茄、玉米进行营养生长的重要阶段,地温的变化直接影响其根系竞争水分和养分的能力。本文通过该时期地温数据进行平均,分析不同位置、不同深度的地温变化规律。结果显示:套种农田不同位置地温的差异性受到地膜覆盖、作物遮阴及土壤水分的综合影响,从不同位置地温平均值变化来看,T1、T2的番茄一侧(A、B)平均地温无差异,而玉米一侧(C、D)的平均地温存在极显著差异(表 3),T3的4个位置平均地温均在0.05水平上显著。说明作物遮阴和土壤水分是导致套种农田不同位置地温出现差异的主要原因。从番茄一侧的不同深度地温变化来看,T1在25 cm范围内均无差异,而T2和T3分别在15 cm处和15 cm以下出现了显著差异性;从玉米一侧的不同深度地温变化来看,T1、T2、T3的地温均在15 cm处出现了差异性。这是因为在覆膜滴灌条件下,不同土壤水分所形成的独特湿润体不同,导致土壤水分在10—15 cm之间出现了差异性,从而改变了土壤的热容量,加之套种农田不同行间位置的遮阴率不同,热量在土层中分配时,其量值有所改变,因而出现了以上差异[21]。综上可知,套种农田在营养生长的共生期,作物遮阴和水分差异是影响地温的主要因素,而地膜覆盖的影响变小。

表3 不同处理及不同位置地温显著性分析/℃ Table 3 Analyze the significance of soil temperature at different processing and positions
同列数值后不同大小写字母表示差异达1%和5%显著水平
处理号 Treatment位置 Position5 cm10 cm15 cm20 cm25 cm平均值 Average value
T1A127.10bB24.98cB24.51bB23.78bB22.48aA24.57bB
B127.57bB25.05cB24.11bB23.89bB22.52aA24.63bB
C129.69aA27.41bA24.63bB22.78cB22.61aA25.42bAB
D129.84aA28.41aA26.34aA25.59aA22.93aA26.62aA
T2A227.55bB26.42cC24.42bAB24.22aA23.85aA25.29bB
B227.62bB26.15cC24.19bB23.84aA22.47aA24.85bB
C230.17aA27.59bB25.14aA23.72aA23.44aA26.01aAB
D230.23aA28.53aA26.37bAB25.14aA23.49aA26.75aA
T3A328.02bB26.97cB24.97bBC24.62aA24.19aA25.75bB
B327.84bB26.87cB24.32cC23.80bB22.31cB25.03cB
C330.29aA27.88bA25.30bB24.10abAB22.38cB25.99bAB
D330.40aA28.61aA26.45aA24.49aAB23.58bA26.71aA
2.3 套种农田不同土壤水分对地温变化的影响

由于7月份套种农田两种作物生长发育趋于稳定,本文为了分析不同高低作物套种不同位置地温的差异,以及不同土壤水分对地温的影响,在7月份对土壤水分、不同位置地温进行了连续监测。结果显示:

距地表5 cm处,土壤受外界环境的影响较大,高水分T1的平均体积含水率明显高于T2和T3,其中T1、T2相差3.8%,T1、T3相差7.2%(图 4),T1、T2、T3之间的体积含水率呈显著性差异(F=6.223587,P=0.020088)。而地温则相反,高水分T1各位置的平均地温明显低于T2和T3,其中T1、T2的平均地温相差0.81 ℃,T1、T3相差1.65 ℃(图 5),T1、T2、T3之间的地温无显著性差异(F=1.44877,P=0.282652)。3个处理的不同位置含水率均表现为B>A>D>C,地温则为D>C>B>A的规律。不同处理的膜内位置B(或D),含水率从T1减小到T2和T3依次降低2.73%、6.04%(或2.71%、6.63%),地温则依次上升0.95 ℃、2.02 ℃(或0.47 ℃、0.97 ℃);不同处理的膜外位置A(或C),含水率从T1减小到T2和T3依次降低3.16%、6.44%(或6.59%、9.71%),地温则依次上升1.16 ℃、1.97 ℃(或0.64 ℃、1.65 ℃)。说明在距地表5 cm处,地温随含水率的降低而上升,其中膜外的上升比率是膜内的1.2倍,表明在土壤表层膜内与膜外的温度差较小。

图 4 不同灌溉定额水平下土壤的平均含水率 Fig. 4 Comparison of average soil moisture content for different irrigation quota
图 5 不同处理在5cm和20cm处的平均温度日变化 Fig. 5 The average temperature of the diurnal variation (soil depth 5cm and 20cm)

2 0 cm处,土壤受外界环境的影响变小,而受水分的影响变大,高水分T1的平均含水率仍高于T2和T3,其中T1、T2相差3.61%,T1、T3相差6.9%(图 4),无显著性差异(F=3.338497,P=0.082302)。地温亦为T1高于T2和T3,其中T1、T2的平均地温相差0.87 ℃,T1、T3相差1.48 ℃(图 5),T1、T2、T3之间的地温无显著性差异(F=1.981115,P=0.193657)。3个处理的不同位置体积含水率均表现为B>A>D>C,地温则为D>B>C>A的规律。不同处理的膜内位置B(或D),含水率从T1减小到T2和T3依次降低2.91%、6.58%(或2.47%、5.77%),地温则依次降低1.04 ℃、1.63 ℃(或1.03 ℃、1.52 ℃);不同处理的膜外位置A(或C),含水率从T1减小到T2和T3依次降低3.6%、6.82%(或5.47%、8.4%),地温则依次降低1 ℃、0.9 ℃(或0.43 ℃、1.52 ℃)。说明在20 cm处,地温随含水率的降低而降低,其中膜外的降低比率是膜内的2.1倍,表明较深层(20 cm以下)膜内地膜波动幅度小于膜外地温波动幅度。

2.4 套种农田不同处理及不同位置土壤温度场的分布特征

为进一步研究套种农田不同行间位置之间地温的传递效应,本文对作物生长旺期,两种作物都发育充分条件下选择7、8月份观测数据进行平均后绘成套种农田地温二维分布图(图 6)。由于监测的是农田日地温的平均值,从图中可以看出,A、C膜外位置在土体内温度变化幅度大,B、D膜内位置在土体内温度变化幅度小,3个处理在C位置处在土壤上层地温明显高于其它位置,这主要是C位置两侧遮阴小,吸热多的原故,而B、D位置尽管覆膜,但白天这两个位置受两侧作物的遮阴影响,地温反而没有C位置大。在土壤下层,C位置的地温明显低于其他位置地温,这是由于C位置无覆膜,保温效果差,到了夜间,上层热量难以向下传递。可见,在15 cm处C位置地温低于B、D位置地温,无覆膜情况下白天太阳辐射增温效应仅对土壤上层有效,而覆膜情况下则0—25 cm都能升温。另外对于不同水分处理温度场的分布主要是土体内含水量不同,对于高水分处理的T1土体内平均含水量较大,因而比热大,所以最高地温与最低地温的差值变小,3个处理的最高地温与最低地温的差值分别为6.02、7.26、8.04 ℃。

图 6 不同处理土壤剖面的温度分布特征 Fig. 6 Comparison of soil profile characteristics of temperature distribution
3 结论

3.1     番茄套种玉米能有效利用太阳辐射提高作物根区温度,其不同行间位置的地温变化主要受地膜覆盖、作物遮阴和土壤水分含量的影响。在套种农田中,地膜的保温作用主要体现在15 cm以下,而作物遮阴和土壤水分含量则对整个耕层都有影响。可见,合理的套种模式和适宜的土壤水分含量有利于作物根区形成良好的生长环境,从而促进作物的生长。

3.2      套种农田不同深度土层土壤含水率与地温的关系不同。本试验通过分析番茄套种玉米表层(5 cm)及深层(20 cm)的地温与含水率的关系得出:5 cm以上的地温与含水率呈反比关系,20 cm以下地温与含水率呈正比关系,可见,在套种农田中,地温与含水率在不同土层呈现的关系也不相同,而并不像単种农田一样,两者呈负相关关系。

3.3      从套种农田不同行间位置地温的传递规律可以看出,高矮作物行间位置(C处)在白天有明显的增温效应,且温度可有效传递到作物两侧,而覆膜可有效防止土壤温度的散失,因此高矮作物套种应选择恰当的行距,从而有效提高作物根区的温度。

参考文献
[1] Willey R W. Resource use in intercropping systems. Agricultural Water Management, 1990, 17(1): 215-231.
[2] 王仰仁, 李明思, 康绍忠. 立体种植条件下作物需水规律研究. 水利学报, 2003, 7: 90-95.
[3] Awal M A, Koshi H, Ikeda T. Radiation interception and use by maize/peanut intercrop canopy. Agricultural and Forest Meteorology, 2006, 139(1): 74-83.
[4] Li F M, Guo A H, Wei H. Effects of clear plastic film mulch on yield of spring wheat. Field Crops Research, 1999, 63(1): 79-86.
[5] Choudhary V K, Bhambri M C, Pandey N, Sharma H G. Effect of drip irrigation and mulches on physiological parameters, soil temperature, picking patterns and yield in capsicum (Capsicum annuum L.). Archives of Agronomy and Soil Science, 2012, 58(3): 277-292.
[6] 潘渝, 郭谨, 李毅, 张江辉, 丁新利. 地膜覆盖条件下的土壤增温特性. 水土保持研究, 2002, 9(2): 130-134.
[7] Gao Y, Duan A W, Sun J S, Li F S, Liu Z G, Liu H, Liu Z D. Crop coefficient and water-use efficiency of winter wheat/spring maize strip intercropping. Field Crops Research, 2009, 111(1): 65-73.
[8] Sampathkumar T, Pandian B J, Mahimairaja S. Soil moisture distribution and root characters as influenced by deficit irrigation through drip system in cotton-maize cropping sequence. Agricultural Water Management, 2012, 103: 43-53.
[9] Zhang L, van der Wert W, Zhang S, Li B, Spiertz J H J Spiertz. Growth, yield and quality of wheat and cotton in relay strip intercropping systems. Field Crops Research, 2007, 103(3): 178-188.
[10] Makumba W, Janssen B, Oenema O, Akinnifesi F K, Mweta D, Kwesiga F. The long-term effects of a gliricidia-maize intercropping system in Southern Malawi, on gliricidia and maize yields, and soil properties. Agriculture, Ecosystems and Environment, 2007, 118(1): 237-243.
[11] 张治, 田富强, 钟瑞森, 胡和平. 新疆膜下滴灌棉田生育期地温变化规律. 农业工程学报, 2011, 27(1): 44-51.
[12] 孙建, 刘苗, 李立军, 刘景辉, Acharya S N. 不同耕作方式对内蒙古旱作农田土壤水热状况的影响. 生态学报, 2010, 30(6): 1539-1547.
[13] 李全起, 陈雨海, 于舜章, 吴巍, 周勋波, 董庆裕, 余松烈. 灌溉与秸秆覆盖条件下冬小麦农田小气候特征. 作物学报, 2006, 32(2): 306-309
[14] Mahrer Y, Naot O, Rawitz E, Katan J. Temperature and moisture regimes in soils mulched with transparent polyethylene. Soil Science Society of America Journal, 1984, 48(2): 362-367.
[15] 王建东, 龚时宏, 许迪, 马晓鹏, 于颖多. 地表滴灌条件下水热耦合迁移数值模拟与验证. 农业工程学报, 2010, 26(12): 66-71.
[16] 吴丛林, 黄介生, 沈荣开. 地膜覆盖条件下SPAC系统水热耦合运移模型的研究. 水利学报, 2000, (11): 89-96.
[17] Hunt H W, Fountain A G, Doran P T, Basagic H. A dynamic physical model for soil temperature and water in Taylor Valley, Antarctica. Antarctic Science, 2010, 22(4): 419-434.
[18] 李毅, 王文焰, 门旗, 钟新才, 谢香文. 宽地膜覆盖条件下土壤温度场特征. 农业工程学报, 2001, 17(3): 32-36.
[19] 李毅, 绍明安. 新疆农田作物覆膜地温极值的时空变化. 应用生态学报, 2004, 15(11): 2039-2044.
[20] 刘浩, 孙景生, 段爱旺, 孙磊, 梁媛媛. 基于AutoCAD软件确定番茄与青椒叶面积的简易方法. 中国农学通报, 2009, 25(5): 287-293.
[21] 李慧星, 夏自强, 马广慧. 含水量变化对土壤温度和水分交换的影响研究. 河海大学学报(自然科学版), 2007, 35(2): 172-175.