免费文献传递   相关文献

Characteristics of surface energy fluxes over a sparse shrubland ecosystem in the farming-pastoral zone of the Loess Plateau, Northwest China.

黄土高原农牧交错带稀疏自然植被生态系统的地表能量特征


基于2011—2012年黄土高原农牧交错带稀疏自然植被生态系统的地表能量通量以及气象数据,对该地区能量平衡各分量(净辐射、感热、潜热和土壤热通量)以及波文比进行日、季节动态的特征分析,研究了潜热通量和感热通量对不同强度降雨事件响应程度的差异,并分析了潜热通量和感热通量的主控因子.结果表明: 该地区净辐射、感热、潜热和土壤热通量的日、季节动态曲线均为单峰型曲线,净辐射、感热通量、潜热通量和土壤热通量的年平均值分别为78.19、33.32、24.91和2.65 W·m-2.在全年能量收支平衡中,感热通量占净辐射的43%,潜热通量占32%,土壤热通量占3%,表明对于黄土高原农牧交错带自然稀疏灌木生态系统,全年能量主要以感热的形式交换.生长季感热和潜热占净辐射的比例相同(36%);而在非生长季,感热占主导,占净辐射的比例高达54%.潜热通量在强、弱降雨事件发生后明显升高,感热通量则明显下降.潜热通量与净辐射、水汽压差及植被参数均显著相关,感热通量与净辐射及空气温度梯度显著相关.

Based on the energy flux and meteorological data during 2011-2012 over a sparse shrubland ecosystem in the farming-pastoral zone of the Loess Plateau, this study investigated the diurnal and seasonal variations of the energy balance components, and discussed the responses of the latent and sensible heat fluxes to different intensities of rainfall events. In addition, we identified the major environmental controlling factors on latent and sensible heat fluxes via correlation analysis. The results showed that the diurnal and seasonal variations of net radiation (Rn), sensible heat flux (H), latent heat flux (LE) and soil heat flux (G) all showed single-peak curves. The annual mean values of Rn, H, LE and G were 78.19, 33.32, 24.91 and 2.65 W·m-2, respectively. The ratios of energy budget components to net radiation were 43% (H/Rn), 32% (LE/Rn), and 3% (G/Rn), which indicated that sensible heat flux was the major form of energy loss at this site. In the growing season, the ratios of sensible heat flux and latent heat flux to net radiation were nearly the same (36%); while in the non-growing season, sensible heat flux accounted for 54% of net radiation. Latent heat flux increased sharply after heavy and weak rainfall events, while sensible heat flux decreased sharply at the same time. Continuous rainfall events would lead to fluctuations in latent and sensible heat fluxes. There were significant correlations between latent heat flux and net radiation, vapor pressure deficit and vegetation parameter, while remarkable correlations were found between sensible heat flux and net radiation, and air temperature gradient.


全 文 :黄土高原农牧交错带稀疏自然植被生态系统
的地表能量特征∗
龚婷婷  雷慧闽∗∗  焦  阳  杨汉波  杨大文
(清华大学水利水电工程系水沙科学与水利水电工程国家重点实验室, 北京 100084)
摘  要  基于 2011—2012 年黄土高原农牧交错带稀疏自然植被生态系统的地表能量通量以
及气象数据,对该地区能量平衡各分量(净辐射、感热、潜热和土壤热通量)以及波文比进行
日、季节动态的特征分析,研究了潜热通量和感热通量对不同强度降雨事件响应程度的差异,
并分析了潜热通量和感热通量的主控因子.结果表明: 该地区净辐射、感热、潜热和土壤热通
量的日、季节动态曲线均为单峰型曲线,净辐射、感热通量、潜热通量和土壤热通量的年平均
值分别为 78.19、33.32、24.91和 2.65 W·m-2 .在全年能量收支平衡中,感热通量占净辐射的
43%,潜热通量占 32%,土壤热通量占 3%,表明对于黄土高原农牧交错带自然稀疏灌木生态
系统,全年能量主要以感热的形式交换.生长季感热和潜热占净辐射的比例相同(36%);而在
非生长季,感热占主导,占净辐射的比例高达 54%.潜热通量在强、弱降雨事件发生后明显升
高,感热通量则明显下降.潜热通量与净辐射、水汽压差及植被参数均显著相关,感热通量与
净辐射及空气温度梯度显著相关.
关键词  能量平衡; 潜热通量; 感热通量; 稀疏灌木; 黄土高原
文章编号  1001-9332(2015)06-1625-09  中图分类号  P404, P464  文献标识码  A
Characteristics of surface energy fluxes over a sparse shrubland ecosystem in the farming⁃
pastoral zone of the Loess Plateau, Northwest China. GONG Ting⁃ting, LEI Hui⁃min, JIAO
Yang, YANG Han⁃bo, YANG Da⁃wen (State Key Laboratory of Hydroscience and Engineering, De⁃
partment of Hydraulic Engineering, Tsinghua University, Beijing 100084, China) . ⁃Chin. J. Appl.
Ecol., 2015, 26(6): 1625-1633.
Abstract: Based on the energy flux and meteorological data during 2011-2012 over a sparse shrub⁃
land ecosystem in the farming⁃pastoral zone of the Loess Plateau, this study investigated the diurnal
and seasonal variations of the energy balance components, and discussed the responses of the latent
and sensible heat fluxes to different intensities of rainfall events. In addition, we identified the major
environmental controlling factors on latent and sensible heat fluxes via correlation analysis. The re⁃
sults showed that the diurnal and seasonal variations of net radiation (Rn), sensible heat flux (H),
latent heat flux (LE) and soil heat flux (G) all showed single⁃peak curves. The annual mean values
of Rn, H, LE and G were 78.19, 33.32, 24.91 and 2.65 W·m
-2, respectively. The ratios of ener⁃
gy budget components to net radiation were 43% (H / Rn), 32% ( LE / Rn), and 3% (G / Rn),
which indicated that sensible heat flux was the major form of energy loss at this site. In the growing
season, the ratios of sensible heat flux and latent heat flux to net radiation were nearly the same
(36%); while in the non⁃growing season, sensible heat flux accounted for 54% of net radiation.
Latent heat flux increased sharply after heavy and weak rainfall events, while sensible heat flux de⁃
creased sharply at the same time. Continuous rainfall events would lead to fluctuations in latent and
sensible heat fluxes. There were significant correlations between latent heat flux and net radiation,
vapor pressure deficit and vegetation parameter, while remarkable correlations were found between
sensible heat flux and net radiation, and air temperature gradient.
Key words: energy balance; sensible heat flux; latent heat flux; sparse shrubland; Loess Plateau.
∗国家自然科学基金项目(51209117,51139002)和水沙科学与水利水电工程国家重点实验室自主研究课题(2014⁃KY⁃04)资助.
∗∗通讯作者. E⁃mail: leihm@ mail.tsinghua.edu.cn
2014⁃09⁃05收稿,2015⁃03⁃09接受.
应 用 生 态 学 报  2015年 6月  第 26卷  第 6期                                                         
Chinese Journal of Applied Ecology, Jun. 2015, 26(6): 1625-1633
    陆面与大气之间的能量交换表征着大气物理气
候系统的下边界条件.地表的能量通量是陆面与大
气能量交换的关键内容[1],并以边界层的湍流运动
形式进行传输.它决定了边界层湍流及扩散的强度
和稳定性[2],另一方面也影响着地表温度、水分的
传输,以及植被的生长状况[3] .因此,有关陆面⁃大气
间的能量通量研究一直是生态学、大气科学和水文
学等多学科领域的重点研究内容[4] .研究陆面⁃大气
之间能量通量动态的方法很多,如采用涡度相关技
术进行观测分析,采用模型(如 SVAT 模型、SHAW
模型)进行模拟分析[2,5] .涡度相关技术因其具有采
样频率高、干扰小、测量精度高以及全天候连续、自
动采集数据等优势,已经成为国际上公认的研究陆
地生态系统与大气之间能量、水分交换的方法[5-8] .
目前已经发展出成熟的全球通量观测网络(Ameri⁃
FLUX,Euro⁃FLUX,Asia⁃FLUX),对研究不同下垫面
条件下陆面⁃大气之间的能量、水分的交换过程起到
了强有力的促进作用.
黄土高原是华北季风区向西北干旱区的过渡
带,是典型的生态脆弱区[5],受周边地区风沙威胁
严重,植被退化、土壤干化严重[9] .近几十年来,国内
对北方地区能量通量的特征研究进行了大量的野外
观测试验,但是多集中于干旱、极端干旱地区的荒
漠[10-11] 和青藏高原等地区[12-13]、 农田生态系
统[8,14-15]和草地生态系统[1,7,16],而对黄土高原典型
的稀疏灌木自然植被生态系统的能量通量特征研究
较少[5] .由于各类陆地生态系统的差异显著,不同下
垫面的陆面⁃大气能量通量的交换强度以及不同时
间尺度上的通量变化动态都存在着很大的差别,为
此,本文选择黄土高原农牧交错带自然稀疏灌木生
态系统为研究对象,基于 2011年 7月 1日—2012年
6月 30日的涡度相关系统观测资料和气象观测资
料,分析该生态系统能量平衡各分量的日、季节动态
特征,以及全年能量收支平衡特征,以增进对黄土高
原典型自然植被生态系统陆面⁃大气之间水热交换
机制的理解.
1  研究地区与研究方法
1􀆰 1  研究区概况
陕西榆林生态水文观测站 ( 38° 26′ 43. 6″ N,
109°28′2.7″ E)位于黄土高原的农牧交错带,海拔
1233 m,于 2011年 6 月中旬建成并开始观测.该地
区为毛乌素沙地向陕北黄土高原的过渡地带,气候
属于半湿润气候向半干旱气候的过渡[17] . 1951—
2011年该地区多年平均降雨量为 402 mm,多年平
均气温 8􀆰 4 ℃ .主要植被类型为沙柳( Salix psam⁃
mophila)和沙蒿(Artemisia ordosica).另外,当地还有
地锦(Parthenocissus tricuspidata)、茅草( Imperata cy⁃
lindrica)、旋覆花( Inula japonica)等多种物种,植被
稀疏低矮,呈斑块状分布.土壤类型为砂土(含有
98%的砂粒).
1􀆰 2  仪器布设
观测仪器安装在一个高 10 m 的观测塔上和 2
个土壤剖面(深度为 2 m)中.涡度相关系统安设在
距离地表 7.53 m高度的悬臂上,包括三维超声风速
仪(CSAT3,Campbell Scientific, Inc.,USA)和开路式
红外 H2O / CO2气体分析仪(LI7500A,LI⁃COR,Inc.,
USA),并通过数据采集器(CR5000,Campbell Scien⁃
tific, Inc.,USA)采集数据,同时在此高度上安装了
净辐射传感器 ( CNR4, KIPP&ZONEN, Delft, the
Netherlands).由于该研究区域植被分布稀疏,考虑
到植被覆盖区和裸土区水分⁃能量的传输有显著不
同,因此在植被区和裸土区各设置一处土壤剖面,土
壤热通量板 ( HFP01SC, Campbell Scientific, Inc.,
USA)安装在距地表 3 cm 深度处,并分别在距离地
表 8 个不同深度(5、10、20、40、60、80、120 和 160
cm)埋设探头,观测土壤体积含水率(CS616,Camp⁃
bell Scientific, Inc.,USA).在距离地表 0.7 m高度处
安装翻斗式雨量计(TE525MM,Campbell Scientific,
Inc.,USA)获取降雨信息.
1􀆰 3  数据处理
数据的质量控制对于数据分析十分关键,外界
环境的影响以及仪器本身的问题都会导致采集到的
原始数据存在不可避免的误差,因此需要对原始数
据进行处理.气象数据(辐射,空气温、湿度)的日动
态过程近似于正弦曲线[18],因此需要去除数据序列
中严重偏离正弦曲线的异常点.涡度相关系统受到
天气因素(如降雨、强风等)以及仪器本身的影响
(如仪器电压低于额定电压,仪器本身受强风影响
发生晃动等)较大,因此采用 LI⁃COR 公司研制开发
的 Eddypro软件进行质量控制并得到 30 min的平均
数据.
由于对数据进行了质量控制,导致数据序列出
现空值,因此需要对数据进行插补.气象数据(如辐
射、气温等)缺测较少,因此直接采用线性插值得到
连续数据.涡度相关系统的数据采用 Falge 等[19]提
出的方法,对连续缺测 1 h 及以下的数据采用线性
插值,对连续缺测 1 h 以上的数据采用平均日变化
6261 应  用  生  态  学  报                                      26卷
(MDV, mean diurnal variation)方法进行插值,即采
用邻近 7 d相同时段的观测值的平均值进行插补.
采用南(裸土区)、北(植被区)不同深度的土壤
含水率加权平均得到该深度的平均土壤含水率
(θ),采用南、北土壤热通量加权得到平均的土壤热
通量(G):
θ=Mθc+(1-M)θb (1)
G=MGc+ (1-M)Gb (2)
式中:θc、Gc和 θb、Gb分别表示植被区和裸土区的土
壤含水率、土壤热通量;M 为植被覆盖度,由月平均
的归一化植被指数(NDVI)计算得到,采用广泛应用
的 Gutman等[20]提出的算式:
M=(NDVI-NDVImin) / (NDVImax-NDVImin) (3)
式中:NDVImax和 NDVImin为全植被和裸土覆盖像元
的 NDVI值,参考文献[21]分别取 0.8 和 0.05,计算
出来的植被覆盖度能够较好地反映该地区的下垫面
覆盖情况.月平均 NDVI数据是由日尺度的 NDVI 数
据平均得到,日尺度的 NDVI 是选择涡度相关的源
区范围内 30个像元(250 m×250 m)的 NDVI数据平
均得到.每个像元的 NDVI由 MODIS⁃Terra 卫星测得
的日尺度的空间分辨率为 250 m 的地表反射率
(http: / / reverb.echo.nasa.gov / reverb)计算得到,采用
近红外(NIR)和可见光(VIS)波段的数据,算式为:
NDVI=(RNIR-RVIS) / (RNIR+RVIS) (4)
能量平衡方程简化为:
Rn =LE+H+G (5)
式中:Rn为净辐射通量;H为感热通量;LE 为潜热通
量;G为地表热通量.将一年中相同时刻的净辐射、
感热、潜热和土壤热通量数值进行平均,分析其日动
态过程.
2  结果与分析
2􀆰 1  能量通量的日动态特征
由图 1 可以看出,黄土高原农牧交错带稀疏自
然植被生态系统中各能量通量的日动态曲线均为单
峰曲线,净辐射的日动态曲线在 13:00 达到峰值
(379.7 W·m-2),随后开始下降,并在 19:00 之后
基本维持不变.潜热通量和感热通量的日动态曲线
的振幅比净辐射小,感热通量的日动态曲线在13:30
达到峰值(141.3 W·m-2),而潜热通量则相对延
迟,在 14:00 达到峰值(67.4 W·m-2).土壤热通量
的日动态曲线峰值(88.7 W·m-2)出现的时间为
13:00.
在9:00—17:00,能量分配以感热为主,且在
图 1  黄土高原农牧交错带稀疏自然植被能量通量的日动
态变化
Fig.1  Diurnal dynamics of energy flux over the sparse shrub⁃
land ecosystem in the farming⁃pastoral zone of the Loess Plateau.
9:00—14:30表现为 G / Rn > LE / Rn,而在 14: 30—
18:00表现为 LE / Rn>G / Rn .夜间时段,能量分配比例
表现为 G / Rn>H / Rn >LE / Rn .在日出(6:30—8:30)、
日落(17:30—18:30)时段内,能量分配比例存在较
大的波动.出现这种现象主要是由于该生态系统的
热源在大气与土壤之间相互转换导致[4] .各能量分
量的方向或大小在这些时段内会发生变化,导致能
量分配比例曲线出现了跳跃式的变化,这种剧烈变
化主要体现在 LE / Rn和 G / Rn,且两者的峰值出现互
补现象.这主要是由于在日出时段,净辐射和土壤热
通量的数值都经历由负值转为正值的过程,在日落
时段,净辐射和土壤热通量的数值都经历由正值转
为负值的过程,而潜热通量在日出、日落时段均为正
值,因此出现了 LE / Rn和 G / Rn曲线互补的现象.
将一年中相同时刻的感热通量和潜热通量数值
进行平均,得到波文比的年均日动态过程曲线,其日
动态曲线过程与能量平衡各分量的日动态曲线过程
相似,夜间时段的数值均<0,直到 7:30 开始波文比
>0,并在 12:00 达到峰值(2.25),一天中波文比数
值>1的时段为 9:00—17:00,这表明该研究区域白
天的净辐射主要分配给感热通量.
2􀆰 2  能量通量的季节动态特征
由图2可以看出,净辐射、感热、潜热和土壤热
72616期                    龚婷婷等: 黄土高原农牧交错带稀疏自然植被生态系统的地表能量特征           
图 2  黄土高原农牧交错带稀疏自然植被能量通量的季节
动态
Fig.2  Seasonal dynamics of energy flux over the sparse shrub⁃
land ecosystem in the farming⁃pastoral zone of the Loess Plateau.
通量 4个分量在季节变化上都呈单峰曲线形式,遵
循着夏季高、冬季低的变化趋势.净辐射的年平均值
为 78.19 W·m-2,感热通量的年平均值为 33􀆰 32
W·m-2,潜热通量的年平均值为 24.91 W·m-2,土
壤热通量的全年平均值为 2.65 W·m-2 .夏季(6—8
月)各能量通量的平均值均较大,净辐射、感热、潜
热和土壤热通量的平均值分别为 126. 46、45􀆰 15、
45􀆰 72和 9.80 W·m-2,而冬季(12 月—次年 2 月)
的净辐射、感热、潜热和土壤热通量的平均值只有
29.00、12.27、8.31 和-6.46 W·m-2 .潜热通量在 11
月末至次年 3 月初数值偏低且变化小.主要是因为
在这段时间内,该地区的自然植被已经凋萎,植被与
大气之间不存在由水汽传输产生的热量[1],即不存
在植被蒸腾;另一方面,这段时间内土壤进入封冻
期,土壤与大气之间由于水汽传输产生的热量少,即
土壤蒸发很少,因此导致了该时段内潜热通量数值
的偏低.
不同生态系统的潜热通量和感热通量的季节动
态存在差异.本研究区内一年生自然植被在该研究
时段内未受到人为因素的干扰,潜热通量和感热通
量的季节动态呈现单峰曲线形式.但对于受到人为
因素干扰的生态系统或者轮作式的农田生态系统,
潜热通量和感热通量的季节动态会出现双峰曲线形
式.例如,李辉东等[7]在科尔沁温带草甸对能量平衡
分量进行分析,发现感热通量的季节动态曲线为双
峰曲线,峰值分别出现在 4和 9月,由于人为割草的
原因,导致感热通量在 9月出现第 2个峰值;雷慧闽
等[22-23]对我国华北平原位山引黄灌区典型农田生
态系统(夏玉米、冬小麦轮作)进行研究,发现潜热
通量的季节动态呈现双峰曲线形式,峰值分别出现
在玉米和小麦生长的旺季.
本研究区能量平衡各分量的月平均最小值都出
现在冬季 12月(表 1),感热通量的月均最大值出现
在 4月,而净辐射、潜热通量和土壤热通量的月均最
大值都出现在 6月.这主要是由于在 4 月,净辐射迅
速增加,而地表植被还未进入生长期,下垫面的比热
较大[7],空气和土壤的升温较大,导致感热通量增
大较快,从而出现峰值.
    波文比是反映地表吸收净辐射能量在潜热通量
和感热通量上的分配比例(β =H / LE) [4] .在 2011—
2012年,该地区波文比出现 2 个峰值(图 3),分别
为 2011年 11 月 24 日(β = 15.69)和 2012 年 4 月 6
日(β= 9.64).第 1个峰值出现的时间正值土壤表层
水分开始进入封冻阶段,土壤与大气之间由于水汽
传输而产生的热量逐渐减少,即潜热通量减少,因此
导致波文比出现峰值.第 2 个峰值主要由于感热通
量的增加导致,感热通量在 4月到达峰值,而此时潜
热通量并没有发生较大变化,因此波文比在 4 月出
现第 2个峰值.对于土壤水在冬季不经历封冻期的
生态系统,波文比的季节变化呈现“U”型分布,冬季
的波文比高,夏季的波文比低[7,24] .研究区波文比的
全年平均值为 1.90,生长季的波文比平均值为 1.11,
非生长季的波文比平均值为 2.44,表明非生长季平
均感热通量明显大于潜热通量,且波动范围较大
表 1  能量通量的月均最大、最小值及出现的月份
Table 1  Maximum and minimum values of energy fluxes
and their occurred months
参数
Parameter
月均最大值
Maximum
monthly
average value
(W·m-2)
最大值月份
Month of
maximum
月均最小值
Minimum
monthly
average value
(W·m-2)
最小值月份
Month of
minimum
净辐射
Net radiation
139.4 6 13.4 12
感热
Sensible heat flux
62.8 4 4.2 12
潜热
Latent heart flux
50.3 6 5.6 12
土壤热通量
Soil heat flux
28.1 6 -9.5 12
8261 应  用  生  态  学  报                                      26卷
图 3  黄土高原农牧交错带稀疏自然植被波文比的季节动态
Fig.3  Seasonal dynamics of Bowen ratio over the sparse shrub⁃
land ecosystem in the farming⁃pastoral zone of the Loess Plateau.
表 2  全年、生长季和非生长季能量平衡分量的收支比例
Table 2  Ratios of energy budget components of the entire
year, growing season and non⁃growing season
时期
Period
潜热通量 /
净辐射
Latent heat flux /
net radiation
感热通量 /
净辐射
Sensible heat flux /
net radiation
土壤热通量 /
净辐射
Soil heat flux /
net radiation
全年
Entire year
31.8 42.6 3.3
生长季
Growing season
36.2 36.4 6.7
非生长季
Non⁃growing season
24.9 53.7 -2.3
(-0.6~15.69).
在该地区全年能量收支平衡(表 2)中,感热通
量占净辐射的 43%,潜热通量占 32%,土壤热通量
占 3%,表明该生态系统能量主要以感热的形式交
换.生长季内,感热通量和潜热通量占净辐射的比例
十分接近(36%),而在非生长季内,感热通量占净
辐射的比例高达 54%,潜热通量只占 25%.
2􀆰 3  影响潜热通量和感热通量的主控因子
2􀆰 3􀆰 1潜热通量和感热通量与降雨的关系  2011—
2012年有 55%的降雨出现在 7—9 月,当降雨发生
时,表层(5和 10 cm)土壤含水率迅速发生变化(图
4),且两深度处的土壤含水率十分接近.20 cm 深度
处土壤含水率对降雨的响应较表层土壤有所延迟,
而这种延迟在 60 cm 深度处表现得更加明显.黄土
高原地区是典型的季节性冻土带,每年冬季都有 2
个月以上的封冻期(本研究区域封冻期为 11 月末
至次年 3月初).该地区表层土壤水分在 11 月末开
始进入冻结阶段,中间层的土壤水随着时间推移也
陆续冻结,120 cm厚度及其以下深度的土壤水没有
经历冻结过程.
选取 2011年 8 月 31 日—2011 年 9 月 10 日的
降雨事件(这段时间的前后3 d未发生降雨事件)和
图 4  不同深度土壤平均含水量值与降雨日累计值的关系
Fig.4  Relationships between daily precipitation and daily mean
soil water content in different depths.
潜热、感热通量进行分析,发现感热通量的低谷对应
于降雨时段(图 5),因为降雨较多时,土壤和空气温
度差降低,导致感热通量较少.2011年 9月 1日发生
强降雨事件(35.1 mm),当天的潜热通量也减少,但
降雨事件发生后的几天内会出现明显升高,主要因
为降雨使得土壤含水量升高,土壤蒸发量和植被蒸
腾增加,因此导致潜热通量升高.另外,连续的降雨
事件也会导致感热通量和潜热通量呈现波动式的变
化[1,8] .
选取 3个典型降雨事件(强降雨事件、弱降雨
事件和微弱降雨事件),分析降雨事件前后 1 d 的潜
热通量和感热通量变化特征(表3) .日累计降雨量
图 5  日尺度潜热、感热通量与降雨量的关系(2011)
Fig.5  Relationships between daily total precipitation (P) and
latent heat flux (LE), sensible heat flux (H) in 2011.
92616期                    龚婷婷等: 黄土高原农牧交错带稀疏自然植被生态系统的地表能量特征           
表 3  潜热通量和感热通量对不同强度降雨过程的响应
Table 3  Latent heat flux and sensible heat flux in response to different intensity rainfall events
参数
Parameter
强降雨事件
Heavy rainfall event
降雨前 1 d
1 d before
rain event
降雨后 1 d
1 d after
rain event
弱降雨事件
Weak rainfall event
降雨前 1 d
1 d before
rain event
降雨后 1 d
1 d after
rain event
微弱降雨事件
Small rainfall event
降雨前 1 d
1 d before
rain event
降雨后 1 d
1 d after
rain event
日照时数
Sunshine duration (h)
12.4 13.3 9.5 8.8 8.9 8.6
净辐射平均值
Average net radiation (W·m-2)
159.2 167.6 48.9 40.7 10.4 8.3
净辐射最大值
Maximum net radiation (W·m-2)
591.8 619.2 362.4 358.8 244.9 236.3
水汽压差平均值
Average vapor pressure deficit (kPa)
1.72 1.62 0.27 0.33 0.10 0.12
水汽压差最大值
Maximum vapor pressure deficit (kPa)
4.33 3.76 0.73 0.67 0.20 0.23
潜热通量平均值
Average latent heat flux (W·m-2)
63.0 77.9 12.7 22.0 7.6 7.9
潜热通量最大值
Maximum latent heat flux (W·m-2)
194.6 240.4 49.0 80.5 40.5 43.4
感热通量平均值
Average sensible heat flux (W·m-2)
50.7 32.6 25.3 18.8 0.85 0.41
感热通量最大值
Maximum sensible heat flux (W·m-2)
181.9 152.6 143.8 141 68.7 62.7
≤1 mm 的降雨事件定义为微弱降雨事件,≥5 mm
的降雨事件为强降雨事件,>1 mm且<5 mm的降雨
事件为弱降雨事件[25] .这 3个典型降雨事件前、后 1
d的日照时数均大于榆林地区日均的日照时数
(6.8~7.9 h) [26],降雨前、后 1 d的净辐射和 VPD条
件相似.对于强降雨事件(2011 年 7 月 28—29 日,
16.2 mm)和弱降雨事件(2011 年 11 月 16—17 日,
3.7 mm),降雨事件次日的潜热通量的平均值和日
动态峰值均高于降雨前 1 d,而降雨事件次日的感热
通量的平均值和日动态峰值则低于降雨前 1 d,并且
这种变化在强降雨事件中表现得更加明显.强降雨
事件的发生导致潜热通量的日平均值增加 14􀆰 9
W·m-2,峰值增加 45.8 W·m-2,感热通量的日平均
值减少 18.1 W·m-2,峰值减少 29.3 W·m-2;弱降
雨事件的发生导致潜热通量的日平均值增加 9􀆰 3
W·m-2,峰值增加 31.5 W·m-2,感热通量的日平均
值减少 6.5 W·m-2,峰值减少 2.8 W·m-2 .微弱降
雨事件(2011年 12 月 22 日, 0.1 mm)的发生对潜
热通量和感热通量的影响没有规律.
2􀆰 3􀆰 2潜热通量和感热通量与环境因子的相关性 
以净辐射 ( Rn )、空气温度梯度 ( ΔT)、土壤温度
(Ts),水汽压差(VPD)和植被参数(NDVI)作为环
境因子代表,分析潜热通量、感热通量与环境因子在
全年、生长季和非生长季的相关关系(表 4),以揭示
不同时期影响潜热通量和感热通量的主控因子.其
中 ,空气温度梯度用2和4 m高度处的温度差值表
表 4  感热通量和潜热通量与环境因子的相关系数
Table 4  Correlation coefficients of latent heat flux, sensible heat flux with environmental factors
参数
Parameter
全年
Entire year
潜热通量
Latent
heat flux
感热通量
Sensible
heat flux
生长季
Growing season
潜热通量
Latent
heat flux
感热通量
Sensible
heat flux
非生长季
Non⁃growing season
潜热通量
Latent
heat flux
感热通量
Sensible
heat flux
净辐射
Net radiation
0.768∗∗ 0.799∗∗ 0.447 0.648∗ 0.640∗ 0.894∗∗
空气温度梯度
Air temperature gradient
0.500 0.721∗∗ - 0.648∗ 0.458 0.781∗∗
土壤温度
Soil temperature
0.748∗∗ 0.548 0.332 0.374 0.458 0.557
水汽压差
Vapor pressure deficit
0.748∗∗ 0.592 0.616∗ - 0.819∗∗ 0.490
植被参数 NDVI 0.656∗ 0.361 0.748∗∗ 0.173 0.346 0.173
∗P<0.05; ∗∗P<0.01.
0361 应  用  生  态  学  报                                      26卷
示,地表温度用 5 cm 深度处的土壤温度表示,选择
白天(9:00—16:00)的潜热通量、感热通量和各环
境因子的平均值进行相关性分析.
    从全年平均尺度分析,潜热通量与净辐射、土壤
温度、水汽压差和植被参数均有显著的相关关系,而
在生长季内,潜热通量与净辐射的相关关系不显著,
与植被参数的决定系数较全年有所增加,表明在生
长季,潜热通量受植被生长状况的影响最大;在非生
长季内,潜热通量与水汽压差的相关关系最为显著,
但是与植被参数的相关性不显著,主要因为在非生
长季,植被已经凋萎,潜热通量主要是以土壤蒸发的
形式进行传输.感热通量在全年平均尺度上与净辐
射的相关关系最为显著,其次是空气温度梯度,而在
生长季内,这种显著的相关关系仍然存在,但是确定
性系数较全年平均尺度有所下降,主要是由于在生
长期内植被生长使潜热通量增加,其占净辐射的比
重增加,感热通量占净辐射的比重减少,导致相关关
系减弱.这也合理地解释了在非生长期内感热通量
与净辐射的相关关系最为显著的现象,同时在非生
长季内感热通量也受到空气温度梯度的影响.
3  讨    论
3􀆰 1  不同生态系统的能量平衡收支比较
不同生态系统的能量收支平衡存在显著差异,
对于植被覆盖度较大的生态系统(如农田、森林、浓
密草甸等),植被蒸腾相对强烈,因此潜热输送占据
主导.有研究表明,在我国华北平原位山引黄灌区典
型农田生态系统中,潜热通量占净辐射的比例在小
麦和玉米全盛期分别为 83%和 57%,且在全年平均
尺度上,这一比例为 59%[18];石羊河玉米覆盖区农
田的潜热通量是玉米农田的主要消耗项[14];长白山
阔叶红松林全年能量收支平衡中,潜热通量占净辐
射的比例为 52%[27];内蒙古科尔沁温带草甸全年潜
热通量是最主要的能量支出项,占净辐射的
50%[7];我国北方半干旱地区草地生态系统潜热通
量占净辐射的 56%[28] .对于植被覆盖度较小的生态
系统(如荒漠草地、稀疏灌丛等),地表植被稀疏,其
生态系统持水和滞水的能力较差[4],因此感热通量
占主导.例如,内蒙古荒漠草原全年能量收支平衡中
感热通量占净辐射 58%,潜热通量仅占 26%[29];敦
煌地区荒漠戈壁的感热通量占净辐射的比例高达
84%[10];张掖绿洲荒漠过渡带感热占据主导优势,
占净辐射的 67%,这表明下垫面不同的植被覆盖状
况对地表能量的分配过程起到了重要的调节作
用[30] .
3􀆰 2  环境因素对潜热、感热通量的影响
降水过程对水、热通量的交换有很大的影
响[29] .降雨事件的发生会影响大气的透明度、地表
反射率等,导致辐射平衡发生变化[31] .同时,降雨事
件还会导致温度降低,土壤含水率增加,土壤蒸发和
植被蒸腾增加,使得潜热通量增加,感热通量减
少[32-33] .阳伏林等[29]研究了降雨事件对内蒙古温带
荒漠草原能量平衡的影响,发现潜热通量和感热通
量的日动态峰值在降雨事件发生后较降雨事件发生
前分别有明显的增加和减少.青藏高原高寒草甸[33]
及敦煌戈壁[34]的研究中均发现,降雨事件发生后潜
热通量较降雨事件发生前有明显的增加,感热通量
则有明显的减少.
在黄土高原农牧交错带稀疏自然植被生态系统
中,水热通量受多个环境因素的影响,其中净辐射是
影响潜热通量、感热通量的最重要因素.阳伏林
等[29]和严坤等[35]在不同生态系统中也得到了相同
的结论.下垫面温度是影响土壤和植被中水分的活
跃程度的重要因素,对土壤蒸发和植被蒸腾影响较
大,同时下垫面温度的日变化还直接影响空气温度
及空气温度梯度[36],从而影响湍流输送热量.感热
通量是由于温差导致的地表与大气之间的热交
换[1],因此其与空气温度梯度的相关性很强.水汽压
差是反映空气湿度的参数,也是土壤蒸发和植被蒸
腾的重要驱动因素之一,直接影响土壤蒸发和植被
蒸腾的水汽输送速率[28],因此其与潜热通量的相关
程度高于与感热通量的相关程度.下垫面植被参数
(NDVI)对植被蒸腾产生显著影响,NDVI越大,植被
蒸腾所占蒸散发的比例越大, NDVI 与潜热通量的
相关关系越显著.
4  结    论
黄土高原农牧交错带稀疏自然植被生态系统中
各能量分量的日动态曲线均为单峰曲线.能量分配
的年平均日动态曲线过程,白天以感热为主,夜间能
量分配比例的顺序为 G / Rn>H / Rn>LE / Rn .波文比的
日动态过程也呈现单峰曲线形式,在 12:00 达到峰
值(2.25).能量通量在季节变化上也呈现单峰曲线
形式,夏季高、冬季低.在全年能量收支平衡中,感热
通量占净辐射的 43%,潜热通量占 32%,土壤热通
量占 3%,表明对于半干旱稀疏植被覆盖地区,能量
主要以感热的形式交换.连续的降雨事件会导致感
热通量和潜热通量呈现波动式变化.强、弱降雨事件
13616期                    龚婷婷等: 黄土高原农牧交错带稀疏自然植被生态系统的地表能量特征           
的发生会导致潜热通量在降雨事件发生后有明显的
升高,感热通量则有明显的下降.潜热通量与净辐
射、水汽压差及植被参数均显著相关,而感热通量则
与净辐射、空气温度梯度显著相关.
参考文献
[1]  Zhang G (张   果), Zhou G⁃S (周广胜), Yang F⁃L
(阳伏林). Dynamics of sensible and latent heat fluxes
over a temperate desert steppe ecosystem in Inner Mongo⁃
lia. Chinese Journal of Applied Ecology (应用生态学
报), 2010, 21(3): 597-603 (in Chinese)
[2]  Sun R (孙  睿), Liu C⁃M (刘昌明). A review on re⁃
search of land surface water and heat fluxes. Chinese
Journal of Applied Ecology (应用生态学报), 2003, 14
(3): 434-438 (in Chinese)
[3]   Yang J (杨   娟), Zhou G⁃S (周广胜), Wang Y⁃L
(王云龙), et al. Characteristics of net ecosystem flux
exchanges over Stipa krylovvii steppe in Inner Mongolia.
Chinese Journal of Applied Ecology (应用生态学报),
2008, 19(3): 533-538 (in Chinese)
[4]  Yan R⁃H (闫人华), Xiong H⁃G (熊黑钢), Chen X⁃F
(陈肖飞). Characteristics of land surface energy over:
Achnatherum splendens grassland in the oasis⁃desert eco⁃
tone of northern piedmont of Tianshan Mountains. Acta
Ecologica Sinica (生态学报), 2015, 35(5): 1 - 13
(in Chinese)
[5]  Cheng X⁃R (成向荣), Yu M⁃K (虞木奎), Huang M⁃
B (黄明斌), et al. Simulation of energy balance of
SVAT system in the farming⁃pastoral zone of the Loess
Plateau. Acta Prataculturae Sinica (草业学报), 2011,
20(5): 160-168 (in Chinese)
[6]  Shen Y (沈  艳), Liu Y⁃F (刘允芬), Wang Y (王
堰). Advances in applying the eddy⁃covariance tech⁃
nique to calculate heat, moisture and CO2 flux. Journal
of Nanjing Institute of Meteorology (南京气象学院学
报), 2005, 28(4): 559-566 (in Chinese)
[7]  Li H⁃D (李辉东), Guan D⁃X (关德新), Yuan F⁃H
(袁凤辉), et al. Diurnal and seasonal variations of
energy balance over Horqin meadow. Chinese Journal of
Applied Ecology (应用生态学报), 2014, 25(1): 69-
76 (in Chinese)
[8]  Li Y⁃J (李翊君), Xu Z⁃Z (许振柱),Wang Y⁃L (王
云龙), et al. Latent and sensible heat fluxes and energy
balance in a maize agroecosystem. Chinese Journal of
Plant Ecology (植物生态学报), 2007, 31(6): 1132-
1144 (in Chinese)
[9]  Wang X⁃W (王锡稳), Wang Y⁃R (王毅荣), Zhang
C⁃J (张存杰). Climate evolvement and response of soil
moisture in typical semi⁃arid region on Loess Plateau.
Journal of Desert Research (中国沙漠), 2007, 27(1):
123-129 (in Chinese)
[10]  Zhang Q (张  强), Cao X⁃Y (曹晓彦). The influence
of synoptic conditions on the averaged surface heat and
radiation budget energy over desert or Gobi. Chinese
Journal of Atmospheric Sciences (大气科学), 2003, 27
(2): 245-254 (in Chinese)
[11]  Ji G⁃L (季国良), Ma X⁃Y (马晓燕), Zou J⁃L (邹基
玲), et al. Characteristics of the radiation budget over
oases in arid region. Arid Meteorology (干旱气象),
2003, 21(3): 29-33 (in Chinese)
[12]  Ji G⁃L (季国良), Gu B⁃W (顾本文), Lv L⁃Z (吕兰
芝). Characteristics of atmospheric heating field over
northern Qinghai⁃Xizang Plateau. Plateau Meteorology
(高原气象), 2002, 21(3): 238-242 (in Chinese)
[13]  Bian L⁃G (卞林根), Lu L⁃H (陆龙骅), Lu C⁃G (逯
昌贵), et al. The characteristics of radiation balance
components of the Tibetan Plateau in the summer of
1998. Chinese Journal of Atmosphere Sciences (大气科
学), 2001, 25(5): 577-588 (in Chinese)
[14]  Ding R⁃S (丁日升), Kang S⁃Z (康绍忠), Zhang Y⁃Q
(张彦群), et al. Characteristics of water vapor and heat
fluxes and the controlling factors over a maize field in the
arid inland region. Journal of Hydraulic Engineering (水
利学报), 2014, 45(3): 312-319 (in Chinese)
[15]  Lei H⁃M (雷慧闽), Yang D⁃W (杨大文), Shen Y⁃J
(沈艳君), et al. Energy and water fluxes in the yellow
river irrigation region. Journal of Tsinghua University
(Science & Technology) (清华大学学报: 自然科学
版), 2007, 47(6): 801-804 (in Chinese)
[16]  Wang W⁃Z (王维真), Xu Z⁃W (徐自为), Liu S⁃M
(刘绍民), et al. The characteristics of heat and water
vapor fluxes over different surfaces in the Heihe River
Basin. Advances in Earth Science (地球科学进展),
2009, 24(7): 714-723 (in Chinese)
[17]  Yang S⁃H (杨述河), Yan H⁃L (闫海利), Guo L⁃Y
(郭丽英). The land use change and its eco⁃environmen⁃
tal effects in transitional agro⁃pastoral region: A case
study of Yulin City in northern Shaanxi Province. Pro⁃
gress in Geography (地理科学进展), 2004, 23(6):
49-55 (in Chinese)
[18]  Lei H⁃M (雷慧闽). Ecohydrological Processes in a
Large Irrigated Area of the North China Plain: Field Ob⁃
servation and Modeling. PhD Thesis. Beijing: Tsinghua
University, 2011 (in Chinese)
[19]  Falge E, Baldocchi D, Olson R, et al. Gap filling stra⁃
tegies for defensible annual sums of net ecosystem ex⁃
change. Agricultural and Forest Meteorology, 2001,
107: 43-69
[20]  Gutman G, Ignatov A. The derivation of the green vege⁃
tation fraction from NOAA / AVHRR data for use in nu⁃
merical weather prediction models. International Journal
of Remote Sensing, 1998, 19: 1533-1543
[21]  Montandon LM, Small EE. The impact of soil reflectance
on the quantification of the green vegetation fraction from
NDVI. Remote Sensing of Environment, 2008, 112:
1835-1845
[22]  Lei HM, Yang DW. Seasonal and inter⁃annual variations
in carbon dioxide exchange over a cropland in the North
China Plain. Global Change Biology, 2010, 16: 2944-
2957
[23]  Lei H⁃M (雷慧闽), Yang D⁃W (杨大文), Liu Y (刘
钰). Comparison of the sensible heat flux measurements
2361 应  用  生  态  学  报                                      26卷
on different spatial scales in an irrigated cropland. Jour⁃
nal of Hydraulic Engineering (水利学报), 2011, 42
(2): 136-142 (in Chinese)
[24]  Sun C (孙  成), Jiang H (江  洪), Chen J (陈 
健), et al. Energy flux and balance analysis of Phyl⁃
lostachys edulis forest ecosystem subtropical China. Acta
Ecologica Sinica (生态学报), 2015, 35(12): 1- 13
(in Chinese)
[25]  Long T (龙  桃), Xiong H⁃G (熊黑钢), Zhang J⁃B
(张建兵), et al. Experimental study on grassland soil
evaporation with different rainfall intensity. Journal of
Soil and Water Conservation (水土保持学报), 2010,
24(6): 240-245 (in Chinese)
[26]  Wu L⁃R (吴林荣), Jiang Z⁃H (江志红), Lu Y⁃P (鲁
渊平), et al. Characteristics of the distribution and evo⁃
lution of sunshine duration and percentage in Shaanxi
Province. Shaanxi Meteorological (陕西气象), 2008
(1): 18-21 (in Chinese)
[27]  Wu J, Duan D, Han S, et al. Energy budget above a
temperate mixed forest in northeastern China. Hydrologi⁃
cal Processes, 2007, 21: 2425-2434
[28]  Wang X⁃X (王修信), Liu X (刘   馨), Chen S⁃H
(陈声海), et al. Correlation analysis of water and heat
fluxes with environmental variables over lawn in semi⁃
arid area. Journal of Basic Science and Engineering (应
用基础与工程科学学报), 2008, 16(5): 770- 777
(in Chinese)
[29]  Yang F⁃L (阳伏林), Zhou G⁃S (周广胜). Characte⁃
ristics and driving factors of energy budget over a tempe⁃
rate desert steppe in Inner Mongolia. Acta Ecologica
Sinica (生态学报), 2010, 30(21): 5769-5780 ( in
Chinese)
[30]  Zhang K (张   凯), Wang R⁃Y (王润元), Zhang Q
(张  强), et al. Characteristics of surface radiation and
energy balance as well microclimate within oasis⁃desert
ecotone of Zhangye on typical clear days in summer.
Journal of Desert Research (中国沙漠), 2007, 27(6):
1055-1061 (in Chinese)
[31]  Wang Y⁃F (王研峰). The Response of Land⁃atmos⁃
phere Atmosphere Interaction to Precipitation in the
Semi⁃arid Area. Master thesis. Lanzhou: Lanzhou Uni⁃
versity, 2013 (in Chinese)
[32]  Xu S⁃X (徐世晓), Zhao L (赵  亮), Li Y⁃N (李英
年), et al. Effect of precipitation on cold season CO2
flux of alpine shrub on Qinghai⁃Tibetan Plateau. Journal
of Soil and Water Conservation (水土保持学报),
2007, 21(3): 193-195 (in Chinese)
[33]  Zhang F⁃W (张法伟), Li Y⁃N (李英年), Zhao X⁃Q
(赵新全), et al. Effects of one precipitation process on
CO2 flux and thermal transpiration in alpine meadow of
Qinghai⁃Tibetan Plateau. Chinese Journal of Ecology (生
态学杂志), 2008, 27(10): 1685-1691 (in Chinese)
[34]  Liu H⁃Y (刘宏谊), Yang X⁃G (杨兴国), Zhang Q
(张  强), et al. Contrast of ground surface radiation⁃
energy balance in summer and winter on Dunhuang Gobi.
Journal of Desert Research (中国沙漠), 2009, 29(3):
558-565 (in Chinese)
[35]  Yan K (严  坤), Li S⁃Y (李生宇), Lei J⁃Q (雷加
强), et al. Characteristics of surface energy exchange in
the artificial shelter forest land of the hinterland of Takli⁃
makan Desert. Arid Land Geography (干旱区地理),
2013, 36(3): 433-440 (in Chinese)
[36]  Yokoharis M, Brown RD, Kato Y, et al. The cooling
effect of paddy fields on summertime air temperature in
residential Tokyo. Landscape and Urban Planning,
2001, 53: 17-27
作者简介  龚婷婷,女,1990 年生,博士研究生.主要从事农
业生态学研究. E⁃mail: gtt12@ mails.tsinghua.edu.cn
责任编辑  孙  菊
33616期                    龚婷婷等: 黄土高原农牧交错带稀疏自然植被生态系统的地表能量特征