全 文 :书青海省三江源地区退化草地蒸散特征
李婧梅1,2,3,蔡海1,2,程茜4,乔春连1,2,褚晖1,2,陈懂懂1,
徐世晓1,赵新全1,赵亮1
(1.中国科学院西北高原生物研究所,青海 西宁810008;2.中国科学院研究生院,北京100049;3.青海省工程咨询中心,
青海 西宁810001;4.青海师范大学数学系,青海 西宁810001)
摘要:了解三江源地区退化草地的蒸散状况有助于认识该地区退化生态系统水循环,对该地区的生态安全有着重
要的意义。本研究利用涡度相关法观测青海省玛沁县小嵩草退化草地(34.35°N;100.50°E,海拔3963m)的水热
通量,对该地区2006-2008年全年的蒸散情况进行了特征分析。研究表明,2006-2008年三江源退化草地的蒸散
值分别为452.24,474.24,459.57mm,降水量分别为460.7,496.1,480.1mm,其中有75%以上的降水分布在生长
季,全年水平上蒸散量与降水量之比(ET/PPT)为95%以上,蒸散量与降水量基本持平。生长季日蒸散量1.8~
1.9mm/d,而非生长季日蒸散量<0.6mm/d。蒸散日变化和年变化均为单峰型,每年6-7月蒸散量最大,为
2.0~2.2mm/d;日蒸散峰值出现在正午左右,其最大值生长季为0.21mm/h,非生长季为0.10mm/h。退化草地
蒸散的主要环境因子是净辐射,其次为饱和水汽压亏和空气温度。随着生长季的推进,潜热对这些因素的敏感性
逐渐增大。与未退化的矮嵩草草甸相比,退化草地生长季蒸散量较小,而非生长季蒸散量较大,这一结果表明,高
寒草地的植被盖度驱动草地生态系统的蒸散。
关键词:涡度相关法;蒸散;三江源地区;青海省;退化草地
中图分类号:S812.29 文献标识码:A 文章编号:10045759(2012)03022311
蒸散(evapotranspiration,ET)是植被蒸腾与植被间裸露土壤间蒸发之和[1],它作为水循环(降水、径流、蒸
散)中的一个关键环节,在调节区域水分平衡、能量平衡方面起着极为重要的作用。蒸散值可以直接反映区域气
候对陆地生态系统的长期潜在影响。土壤水分含量、生物量、养分状况和水分收支等生态因子和生态过程都受蒸
散的影响[2]。由于近来全球范围内水资源的短缺和分布不平衡,水资源的合理利用和管理定量化日益迫切,蒸散
问题的研究越来越受到人们的重视。
草地是陆地生态系统重要的组成部分之一,天然草地面积占据全球自然植被的32%[3],是目前人类活动影
响较为严重的区域,维持其生态功能的正常对全球及区域的生态平衡有重要意义。因此,在全球变化研究中,草
地生态系统被列为重要的研究对象。它的水分循环过程无疑对全球的水循环过程起着至关重要的作用。三江源
地区是长江、黄河、澜沧江的发源地,也是我国和亚洲最重要的河流上游关键源区[4]。3条江河每年向下游供水
600亿m3,其中长江年平均径流量179.3亿 m3;黄河年平均径流量232亿 m3,占整个黄河流域水资源总量的
49%;澜沧江年平均径流量108.9亿m3[5],因此有“江河源”之称,一直以来被誉为“中华水塔”。同时该区域平均
海拔4000m左右,自然环境脆弱,对气候变化反应强烈,是生态系统最敏感的地区之一,对我国乃至全球的生态
环境有着极为深远的影响[6]。近来,过度放牧,人类活动加剧,加之生态系统脆弱,使三江源区的草地生产力下
降,生态环境恶化,水源涵养能力急剧下降,在这些因素的综合影响下三江源地区草地大面积退化:植物和土壤质
量衰退,生产力、经济潜力和服务功能降低,环境变劣以及生物多样性或复杂程度降低,恢复功能减弱或失去恢复
能力[7]。其中,中度以上退化的草地面积达0.12×108hm2[8],占可利用草地面积的50%~60%,严重地区已沦
为次生裸地或利用价值极低的“黑土滩”(以嵩草属植物为建群种的高寒草地严重退化后,生草土层被破坏,形成
第21卷 第3期
Vol.21,No.3
草 业 学 报
ACTAPRATACULTURAESINICA
223-233
2012年6月
收稿日期:20110401;改回日期:20110721
基金项目:国家973计划(2010CB8335010114),国家自然科学基金项目(41030105,30770419,30970519)和青海省科技项目(2011Z734)
资助。
作者简介:李婧梅(1984),女,青海人,在读硕士。Email:jml107@163.com
通讯作者。Email:lzhao@nwipb.ac.cn
的大面积次生裸地),约占退化草地总面积的40%,并逐年加快增长[9]。据刘纪远等[8]研究发现,三江源地区草
场已呈全面退化的趋势,加之该生态系统脆弱,使草地生产力下降,生态环境恶化,水源涵养能力急剧下降,不仅
使源区居民生活受到极大影响,同时也威胁着长江、黄河流域乃至东南亚诸国的生态安全[10]。在这种背景下,有
必要对这一地区的水分蒸散状况做出科学判断。
近来,涡度相关技术在观测生态系统物质和能量交换中得到了广泛应用[11],在不同的生态系统间对蒸散和
能量平衡的研究已成为国际联网研究的热点之一,而之前的研究多集中于放牧[12,13]、不同水分条件[14]以及正常
状况下不同植被状况[2,15]等对蒸散的研究,利用涡度相关法针对退化生态系统蒸散情况的研究却鲜有报道。
本研究利用涡度相关技术在三江源小嵩草(犓狅犫狉犲狊犻犪狆狔犵犿犪犲犪)退化草地测定的水汽通量数据,对比分析了
2006-2008年该地区蒸散值的变化特征,简单分析了可能的影响因子。为进一步研究退化草地蒸散的变化提供
了科学依据,也为三江源地区退化草地的修复和人工草地的建植工作奠定了理论基础。
1 材料与方法
1.1 研究区概况
本研究位于青海省果洛藏族自治州玛沁县内大武镇东南部15km处格多牧委会草场,地理位置34.35°N,
100.50°E,海拔3963m,多年平均降水量500mm,太阳辐射年总量为5500~6800MJ/m2;无绝对无霜期,牧
草生长期为110~130d[16]。土壤类型以高山草甸土和高山灌丛草甸土为主,土壤表层和亚表层中的有机质含量
丰富,地下0~40cm土层砂石含量约10%。该地区为冬季放牧,建群种为小嵩草,伴生种有矮嵩草(犓狅犫犲狉狊犻犪
犺狌犿犻犾犻狊)、垂穗披碱草(犈犾狔犿狌狊狀狌狋犪狀狊)、早熟禾(犘狅犪spp.)、黄帚橐吾(犔犻犵狌犾犪狉犻犪狏犻狉犵犪狌狉犲犪)、铁棒槌(犃犮狅狀犻狋狌犿
狆犲狀犱狌犾狌犿)、鹅绒委陵菜(犘狅狋犲狀狋犻犾犾犪犪狀狊狉犻狀犪)、二柱头草(犛犮犻狉狆狌狊犱犻狊狋犻犵犿犪狋犻犮狌狊)、异针茅(犛狋犻狆犪犪犾犻犲狀犪)、矮火
绒草(犔犲狅狀狋狅狆狅犱犻狌犿狀犪狀狌犿)、细叶亚菊(犃犼犪狀犻犪狋犲狀狌犻犳狅犾犻犪)、美丽风毛菊(犛犪狌狊狊狌狉犲犪狊狌狆犲狉犫犪)、三裂叶碱毛茛
(犎犪犾犲狉狆犲狊狋犲狊狋狉犻犮狌狊狆犻狊)等。生长期盛期7月植被覆盖度约75%,高度约5cm,根系主要分布于0~10cm的土
壤表层。4月底5月初为植物返青期。该草场中度退化,主要表现为:植被覆盖率低,杂草呈斑块状分布,部分土
壤裸露,地面时有约0.25m2的鼠洞,地带性植被嵩草属植物在群落内的重要值降低。
1.2 观测方法
本研究将每年1月1日-4月20日划分为生长季前期(第1~120天),4月21日-10月27日(第121~300
天)为生长季,其余为非生长季。
涡度相关观测系统安置于地形平坦、开阔,下垫面均一的中心地带,有足够大的“风浪区”,满足安装通量塔的
条件。本研究中对水汽通量观测的主要传感器为Li7500开路红外(CO2/H2O)气体分析仪,距地面220cm。观
测时使用的主要仪器及安装高度见表1。每年4月采用高纯氮标定Li7500红外气体分析仪的CO2 和 H2O的
零点。仪器采样频率为10Hz,每15min输出1组数据。
涡度相关技术通过直接测定大气中物理量的脉动与平均风速的协方差求算湍流通量的方法,求算显热(sen
sibleheatflux,H)和潜热(latentheatflux,LE)通量的公式为:
犎=ρ犆狆ω′犜′ (1)
犔犈=ρλω′狇′ (2)
式中,犎 为显热通量(W/m2),犔犈 为潜热通量(W/m2),ρ为空气密度(kg/m
3),犆狆 为空气的定压比热
(J/kg·k),λ为水的汽化潜热(MJ/kg),ω′为垂直湍流速度脉动(m/s),犜′为气温脉动(K),狇′为比湿脉动,上横线
表示协方差。
1.3 数据处理
原始数据必须进行质量控制后才可用于分析。本研究采用李春等[17]提出的通量数据的处理流程进行数据
处理:主要有野点剔除,WPL校正[13]以及三维坐标旋转。首先对各观测值进行阈值判断,根据当地气候特征,各
观测值的合理范围为LE(-500,700),H(-500,700),土壤热通量G(-100,250),温度Ta(-50,50)。同时由
于降水、电力不稳定、虫禽干扰引起的错误数据也需要剔除;然后使用 WPL方法校正由于空气水汽传输引起的
水汽密度和CO2 波动造成的通量计算误差;最后采用三维坐标旋转修正对于仪器倾斜和地面坡度导致的通量误
差。通量数据经筛选后剩余的数据白天、夜间和全天覆盖率分别为54.9%,68.58%,62.72%。
422 ACTAPRATACULTURAESINICA(2012) Vol.21,No.3
表1 观测要素和设备
犜犪犫犾犲1 犈犾犲犿犲狀狋狊犪狀犱犻狀狊狋狉狌犿犲狀狋狊狅犳犿犲犪狊狌狉犲犿犲狀狋狊
测定要素
Measureelements
传感器和分析仪
Instrumentsandsensors
制造商
Manufacturers
安装高度
Location(cm)
风向 Winddirection 034AL Campbel 110,220
风速 Windspeed 014A Campbel 110,220
温度、湿度Airtemperatureandhumidity HUMP45C Campbel 110,220
降水量Precipitation TE525MM Campbel 50
光合有效辐射Photosyntheticphotonfluxdensity LI190SB LiCor 150
净辐射Netradiation CNR1 KippandNonen 150
地表温度Soilsurfacetemperature 170 Campbel 0
土壤温度Soiltemperature 150T Campbel -5,-20,-30
土壤热通量Soilheatflux HFT3 Campbel -5
土壤含水量Soilwatercontent CS615 Campbel -5,-20,-50
三维风速和温度Verticalwindspeedandairtemperature CSAT3 Campbel 220
CO2和水热通量CO2andH2Oflux Li7500 LiCor 220
数据采集与通信Datalogger CR23X/CR5000 Campbel
另外,使用涡度相关法测量物质通量时,湍流强度对数据质量会有一定程度的影响:夜间大气层结稳定,湍流
较弱,存在通量低估的问题[18]。为解决这一问题,一般剔除摩擦风速U(0.15~0.30m/s)的夜间通量值。本研
究以夜间摩擦风速0.2m/s为临界值,舍去小于0.2的通量观测值[19]。对于各种原因引起的数据缺失和前期剔
除的观测数据,采用Falge等[20]的方法进行数据的插补。其中常规气象数据的缺失,用距离观测点10km左右
的相邻站点的数据回归后,使用预测值代替缺失值。水汽通量的缺失采用 MDV法(日平均法)插补缺失值:运用
滑动平均的方法用临近观测日内同时刻的数据插补缺失的数据,其中生长季时间窗口为7d,非生长季时间窗口
为14d。
能量平衡公式可表示为:
犔犈+犎=犚狀-犌-犛 (3)
式中,犚狀为净辐射(W/m2),犌为土壤热通量(W/m2),由2个土壤热通量板(地下0.05m)测量取平均值,犛为储
存项,由于三江源地区退化草地植被较矮小(小于0.3m),犛可以忽略[19]。
在能量闭合验证测定系统的基础上,利用FAO提供的方法通过潜热通量计算蒸散值[21]:
犈犜=犔犈/λρ (4)
式中,犔犈为潜热通量(W/m2),λ为水的汽化潜热(MJ/kg)。
1.4 数据分析
采用SAS8.0(SASInstituteInc.,Cary,NC,USA)软件进行统计分析。
2 结果与分析
2.1 退化草地环境因子动态
图1与表2为观测点2006-2008年空气温度、土壤温度、土壤含水量、降水量、光合有效辐射、饱和水汽压差
的变化规律。研究期间这些环境因子的变化趋势相似又略有不同。2008年空气温度较2006和2007年低(表
2),生长季与生长季前期尤为明显。冬季空气温度最低值为-26.28℃,夏季的最大值21.77℃。温度日变异较
大,非生长季温差最大可达42.75℃,生长季温差可达31.79℃。2006-2008年生长季平均温度为6.8℃,非生长
季平均温度为-5.27℃。3年内5cm土壤温度约为3℃,略高于空气温度3.2℃(表2)。研究期间各年降水为
420~490mm,年差异较小,其中生长季降水占全年降水的90%以上。大部分总降水量小于1.0mm/d,约有
74.58%降水少于5mm/d。
522第21卷第3期 草业学报2012年
图1 青海三江源退化草地2006-2008年环境因子的变化
犉犻犵.1 犜犺犲犲狀狏犻狉狅狀犿犲狀狋犪犾犳犪犮狋狅狉狊狏犪狉犻犪狋犻狅狀犱狌狉犻狀犵2006狋狅2008
误差线为标准误Thebarsindicatedthestandarderror.下同Thesamebelow.
由于生长季降水丰沛,生长季内土壤含水量在降水事件发生后呈上升趋势,但是在发生中到大雨(18~35
mm)土壤含水量并不随之剧烈上升。生长季内土壤含水量在0.22~0.35cm3/cm3 间波动,非生长季内土壤含
水量显著下降,尤以12月至次年1月明显(<0.1cm3/cm3)(图1)。
研究期间年蒸散量小于年降水量,分别为452.24,474.24,459.57mm,低于当年降水量8.40,15.86,20.53
622 ACTAPRATACULTURAESINICA(2012) Vol.21,No.3
mm(表2)。在整个生长季前期,累计蒸散量均大于累计降水量(图2)。生长季盛期,随着降水量增多,蒸散量与
降水量逐渐持平,然后降水量反超蒸散量。随着生长季的结束,降水量减少,但由于生长季累计降水量保持在一
个较高的水平,始终保持着大于累计蒸散量的状态。由于2007年生长季降水发生时间较早,同时比2006,2008
年多,因此,2007年累计降水量约于6月2日超过累计蒸散量,早于2006和2008年。2006,2008年累计降水量
超过累计蒸散量的时间大致在8月中下旬,比2007年晚2个月。
表2 各年度不同生长阶段各环境因子的变化
犜犪犫犾犲2 犆狅犿狆犪狉犻狊狅狀狅犳犲狀狏犻狉狅狀犿犲狀狋犪犾犮狅狀犱犻狋犻狅狀狊犱狌狉犻狀犵狋犺犲犵狉狅狑犻狀犵犪狀犱狀狅狀犲犵狉狅狑犻狀犵狊犲犪狊狅狀犻狀狊狋狌犱狔狊犻狋犲
生长阶段
Growingphase
年份
Year
空气温度
Ta
(℃)
土壤温度
Ts
(℃)
饱和水汽
压差VPD
(kPa)
光合有效
辐射PPFD
(mol/m2·d)
降水
PPT
(mm)
土壤湿度
SWC
(cm3/cm3)
日蒸散量
Dailyevapotranspiration
(mm/d)
蒸散量
Evapotranspiration
(mm)
生长季前期 2006 -4.58 -2.22 0.28 64.03 18.3 0.11 0.67 72.77
Pregrowing 2007 -6.18 -2.52 0.26 67.67 16.4 0.11 0.73 79.43
season 2008 -6.67 -2.47 0.21 79.89 8.2 0.11 1.01 80.10
生长季 2006 5.65 9.84 0.34 71.37 434.5 0.22 1.86 353.69
Growingseason 2007 6.23 8.28 0.35 72.52 440.7 0.24 1.92 365.03
2008 5.43 7.12 0.34 68.07 454.6 0.23 1.81 344.74
非生长季 2006 -7.72 -2.34 0.20 57.86 7.9 0.11 0.39 25.77
Nongrowing 2007 -5.36 -0.12 0.26 73.28 39.6 0.15 0.60 29.79
season 2008 -7.06 -2.76 0.20 59.22 17.3 0.11 0.53 34.74
全年Annual 2006 0.18 4.04 0.30 67.24 460.7 0.17 1.24 452.24
2007 0.64 3.60 0.31 70.87 496.1 0.19 1.36 474.24
2008 -0.15 2.70 0.28 66.90 480.1 0.17 1.32 459.57
Ta:Airtemperature;Ts:Soiltemperature;VPD:Saturationvaporpressuredifference;PPFD:Photosyntheticphotonfluxdensity;PPT:Precipitati
on;SWC:Soilwatercontent.
2.2 退化草地蒸散的动态特征
研究期间三江源地区退化草地蒸散全年日变化均呈单峰型(图3)。日出前及日落后,蒸散量波动极小,
10:00以后逐渐上升,正午前后到达最大值。生长季0:00-9:00(BST:北京时间)蒸散量小于0.1mm/h,非生长
季小于0.05mm/h,6:00-13:00生长季约增加0.027mm/h,非生长季约增加0.005mm/h,随后又逐渐下降,
日落后趋于平缓。无论生长季或非生长季日蒸散最大值出现时刻一致。不同年度蒸散量日变化的情况是:生长
季内在10:00以后,2006年蒸散量上升速度高于2007与2008年,15:00以后,其下降速度也低于2007与2008
年。生长季前期与非生长季,2006年蒸散量出现最大值后,其变化幅度很小,直至16:00以后才有明显下降,且
全天都小于2007与2008年。而2007和2008年在生长季前期与非生长季出现最大值后,其变化态势与生长季
类似。方差分析表明,无论生长季或非生长季3年蒸散日变化均无显著差异。
三江源地区退化草地蒸散的季节变化变异较大:蒸散量的最大值出现在生长季,期间蒸散平均值约1.8~
1.9mm/d;而非生长季蒸散量趋于0,平均约为0.6mm/d(图4)。早春时分蒸散量相对较小,但从5月开始迅速
增长,6,7月出现年蒸散最大值,为2.0~2.3mm/d,之后随着植被凋落和可利用能量的减少而逐渐减小。至生
长期末期蒸散量减小的速率较快,从10月1.12/1.58/1.32mm/d降至11月0.32/0.50/0.71mm/d(2006/
2007/2008)。月均蒸散日平均值也呈单峰型变化,极值为0.26和2.30mm/d,分别出现于2007年1月与6月。
2.3 潜热与环境因子的相关关系
为了进一步探讨蒸散与各环境因子的关系,采用相关分析与回归分析对不同生长阶段潜热通量与空气温度、
土壤温度、饱和水汽压差以及净辐射关系进行分析(表3)。一方面,通过相关分析发现,无论在哪个生长阶段,净
辐射对该区域潜热值的影响最大,其次是饱和水汽压差,空气温度和土壤温度。另一方面,回归分析的结果表明,
722第21卷第3期 草业学报2012年
随着生长的推进,潜热对各因子的敏感程度(a)有所不同。在生长季,潜热对空气温度、净辐射和饱和水汽压差等
因素的变化较全年其他时间敏感。其中温度的贡献为6.46W/(m2·℃),即温度每升高1℃,潜热变化6.46
W/m2;饱和水汽压差的贡献为126.40W/(m2·kPa),即饱和水汽压差每升高1kPa,潜热变化126.40W/m2;
净辐射的贡献为净辐射每升高1μmol/(m·s),潜热变化0.21W/m
2;而潜热对土壤温度的敏感程度却与其他
因子表现出不同的趋势:生长季前期,土壤温度对潜热变化的贡献为2.86W/(m2·℃)。
图2 2006-2008年累积蒸散量和累积降水量的变化
犉犻犵.2 犆犺犪狀犵犲狅犳犮狌犿狌犾犪狋犻狏犲狆狉犲犮犻狆犻狋犪狋犻狅狀(犘犘犜),
犲狏犪狆狅狋狉犪狀狊狆犻狉犪狋犻狅狀(犈犜)犳狉狅犿2006狋狅2008
图3 三江源退化草地2006-2008年蒸散量日变化
犉犻犵.3 犇犻狌狉狀犪犾狏犪狉犻犪狋犻狅狀狅犳狋犺犲犲狏犪狆狅狋狉犪狀狊狆犻狉犪狋犻狅狀犱狌狉犻狀犵2006狋狅2008
狅犳狋犺犲犱犲犵狉犪犱犲犱犵狉犪狊狊犾犪狀犱
图4 三江源退化草地蒸散月均值(黑色粗线)及逐日蒸散值变化(灰色)
犉犻犵.4 犇犪犻犾狔犲狏犪狆狅狋狉犪狀狊狆犻狉犪狋犻狅狀(犵狉犲狔)犪狀犱狊狌狆犲狉犻犿狆狅狊犲犱犿狅狀狋犺犾狔犿犲犪狀犈犜(犺犲犪狏狔犫犾犪犮犽犾犻狀犲)
犳狅狉狋犺犲狆犲狉犻狅犱犑犪狀狌犪狉狔2006狋狅犇犲犮犲犿犫犲狉2008
822 ACTAPRATACULTURAESINICA(2012) Vol.21,No.3
为了进一步确定所有因子对潜热的综合影响,将各因子逐步回归后得到以下关系式:
犔犈=0.14犚狀+54.58犞犘犇+0.89犜狊-0.18犜犪+69.21犛犠犆-5.52(犚2=0.63,犘<0.0001)
其参数估计见表4。
3 讨论
基于2006-2008年三江源退化草地的涡度相关
系统观测数据,初步研究了其蒸散变化的动态特征,并
对其影响因子作了简单分析。
3.1 蒸散动态特征
蒸散量的日变化呈单峰型变化。无论生长季或非
生长季每日正午12:00-14:00出现最大值,生长季平
均日最大蒸散值为0.21mm/h;非生长季平均日最大
蒸散值为0.10mm/h。月均日变化也为单峰型,表现
为:在生长季前期随着温度回升,土壤消融,植被返青,
蒸散值逐渐增大,每年6,7月蒸散出现全年最大值,月
均蒸散达2.0mm/d以上,此时降水也达到全年的最
大值,说明蒸散和降水是同步的,并不出现时滞。另一
方面,草地生物量和叶面积在每年6,7月达到最大,土
壤含水量同时逐渐增大,也驱动了蒸散的加速。3年
最大日蒸散量出现在2007年6月12日,为3.26
mm/d。非生长季土壤冻结,植被枯萎,土壤蒸发和植
物蒸腾量都非常微弱,导致这个时期的蒸腾量趋于0。
植被特征和降水的季节性分布是决定蒸散值季节变化
的主要因素。三江源地区退化草地季节变化趋势与内
蒙古典型草原的蒸散变化[22]相似(年蒸散最大值3.5,
4.0mm/d),其变化更多的是受土壤水分的影响,还有
部分原因是凋落物增加了水分截留却减少了水分渗
表3 青海三江源退化草地不同生长阶段潜热
与环境因子的关系
犜犪犫犾犲3 犆犺犪狉犪犮狋犲狉犻狊狋犻犮狊狅犳犾犻狀犲犪狉狉犲犵狉犲狊狊犻狅狀犪狀犪犾狔狊犻狊(狔=犪狓+犫)
狅犳犾犪狋犲狀狋犺犲犪狋狏狊.犪犻狉狋犲犿狆犲狉犪狋狌狉犲,狊狅犻犾狋犲犿狆犲狉犪狋狌狉犲犪狋
狋犺犲犱犲狆狋犺狅犳5犮犿,犞犘犇犪狀犱狀犲狋狉犪犱犻犪狋犻狅狀(犚狀)
生长阶段
Period
生长季前期
Pregrowing
season
生长季
Growing
season
非生长季
Nongrowing
season
全年
Annual
LEvs.Ta 狉 0.52 0.55 0.54 0.51
犪 2.70 6.46 2.41 2.87
犫 37.97 14.97 31.03 33.65
LEvs.VPD 狉 0.53 0.52 0.57 0.54
犪 83.00 126.40 85.89 102.34
犫 1.91 8.70 -4.00 3.14
LEvs.Rn 狉 0.67 0.63 0.68 0.61
犪 0.12 0.21 0.10 0.16
犫 13.16 29.09 12.28 20.68
LEvs.Ts 狉 0.29 0.13 0.21 0.32
犪 2.86 1.80 1.81 2.37
犫 28.87 29.63 18.00 24.33
注:狉为两因子间线性相关程度,犢=犪狓+犫为两因子的回归关系式。
Note:狉isthelinearcorrelationcoefficientoftwofactors,犢=犪狓+犫is
theregressionoftwofactors.
透。而在本研究中由于植被盖度低,生态位重叠少,资源利用率低,凋落物量较少,因此凋落物不会影响蒸散量。
由于全球各气候区降水分布迥异,蒸散强弱不同,地球表面草地生态系统的水分平衡状况存在显著差异。表
5为三江源地区退化草地与其他草地生态系统蒸散研究的比较,可以看出本研究区域的蒸散值处于这些研究结
果的中间范围。
表4 多元回归的参数估计
犜犪犫犾犲4 犘犪狉犪犿犲狋犲狉犲狊狋犻犿犪狋犲狅犳狋犺犲犿狌犾狋犻犳犪犮狋狅狉狊狉犲犵狉犲狊狊犻狅狀
变量
Variable
参数估计
Parameterestimate
标准参数估计
Standardizeestimate
标准误差
Standarderr.
狋值
狋value
犘值
犘value
截距Inception -5.52 0 0.52 -10.53 <0.001
净辐射Netration(Rn) 0.14 0.56 0 203.47 <0.001
空气温度 Airtemperature(Ta) -0.18 -0.03 0.03 -5.11 <0.001
土壤温度Solitemperature(Ts) 0.89 0.11 0.04 20.25 <0.001
饱和水汽压差 Airvaporpressdifict(VPD) 54.58 0.27 0.64 84.86 <0.001
土壤水分Soilwatercontent(SWC) 69.21 0.09 2.69 25.69 <0.001
922第21卷第3期 草业学报2012年
表5 不同草地类型蒸散状况的比较
犜犪犫犾犲5 犜犺犲犮狅犿狆犪狉犻狊狅狀狅犳犲狏犪狆狅狋狉犪狀狊狆犻狉犪狋犻狅狀(犈犜)犪犿狅狀犵犱犻犳犳犲狉犲狀狋犵狉犪狊狊犾犪狀犱狋狔狆犲狊
纬度
Latitude
经度
Longitude
海拔
Altitude
(m)
植被类型
Vegetation
type
优势种
Dominant
species
年蒸散量
Annual
ET
(mm)
蒸散/潜热Evapotranspiration/Latent
生长季
Growing
season
(mm/d)
非生长季
Nongrowing
season
(mm/d)
最大值
Maximum
(mm/d)
降水量
Precipitation
(mm)
文献
Reference
37.60°N 101.30°E 3250 高寒草甸
Alpinemeadow
矮嵩草
犓狅犫狉犲狊犻犪犺狌犿犻犾犻狊
341~
426
1.90~
2.23
<0.40 - 554.0~
666.0
[15]
47.21°N 108.74°E 1235 典型草原
Typicalsteppe
针茅、苔草、冷蒿
犛狋犻狆犪犽狉狔犾狅狏犻犻,
犆犪狉犲狓犱狌狉犻狌狊犮狌犾犪,
犃狉狋犲犿犻狊犻犪犳狉犻犵犻犱犪
163 1.70~
2.80
<0.50 2.80 248 [12]
43.53°N 106.67°E 1200 典型草原
Typicalsteppe
羊草
犔犲狔犿狌狊犮犺犻狀犲狀狊犻狊
350 4.22±
0.96a
0.53±
0.40a
3.50~
4.00
350.00 [22]
49.43°N 12.56°E 951 湿生混合草场
Mositmixed
grassland
粗穗冰草、长喙婆罗
门参犃犵狉狅狆狔狉狅狀
犱犪狊狔狊狋犪犮犺狔狌犿,
犜狉犪犵狅狆狅犵狅狀犱狌犫犻狌狊
250 - <0.50 2.20~
4.50
75.10~
295.30b
[2]
46.77°N 100.92°E 518 放牧的混生草场
Grazed mixed
grassprairie
冰草、垂穗草、针茅
犃犵狉狅狆狔狉狅狀spp.,
犅狅狌狋犲犾狅狌犪spp.,
犛狋犻狆犪spp.
489 2.40~
2.70
- 4.40 213.00~
373.00
[23]
不放牧的混生
草场Nongrazed
mixedgrass
prairie
冰草、垂穗草、针茅
犃犵狉狅狆狔狉狅狀spp.,
犅狅狌狋犲犾狅狌犪spp.,
犛狋犻狆犪spp.
455 2.30~
2.90
- 4.70 261.00~
386.00
西部芽草草原
(放 牧)Grazed
westernwheat
grassgrassland
蓝茎冰草
犘犪狊犮狅狆狔狉狌犿狊犿犻狋犺犻犻
497 2.60~
3.00
- 4.40 225.00~
401.00
34.35°N 100.50°E 3963 高寒退化草地
Alpinedegraded
grassland
小嵩草
犓狅犫狉犲狊犻犪狆狔犵犿犪犲犪
454~
474
1.81~
1.92
<0.6 3.20 460.00~
496.00
本研究
Thisstudy
a:MJ/(m2·d)为潜热单位TheunitofLE;b:生长季降水量,其余均为全年降水量Growingseasonprecipitation,therestaretheannualprecipita
tion.
3.2 蒸散影响因子
在整个研究期间三江源退化草地净辐射与潜热的相关性最高,随着生长季的推进,水热条件变好,潜热对各
环境因子的相关性及敏感性逐渐增强,在生长季盛期达到最高。而温度、净辐射、土壤水分以及饱和水汽压差的
综合影响可以解释研究区蒸散变异的63%。
在水分充足的情况下,某一地区的蒸散状况是由可用能量所决定的[24]。而可用能量很大程度上由反照率、
植被高度、植被盖度、植被组成、土壤湿度、太阳辐射所决定[25]。土壤水分作为蒸散的影响因素已在大量文献中
证明[26],在三江源地区,由于地下根系多生长在0~10cm,因此0~5cm的土壤水分含量对蒸散的影响最大。饱
032 ACTAPRATACULTURAESINICA(2012) Vol.21,No.3
和水汽压的增加可以使叶片的气孔导度上升,进而加剧植被蒸腾。此外,对退化草地的反照率研究表明[27],该区
域年均反照率为0.22,生长季(5-9月)的平均反照率为0.18,非生长季的为0.25,生长季反照率略低于未退化
的海北矮嵩草草甸(0.20),可能原因是,海北地区生长季期间处于高原植物的开花期,而其花朵鲜艳的颜色会对
反照率有一定的影响,同时,退化草地土壤颜色较暗,降水多,土壤含水量较高,也会使反照率下降。然而无论是
海北矮嵩草草甸还是退化草地,相比于其他草地生态系统,生长季地表反照率的值都略高[28]。尽管青藏高原地
区太阳辐射强烈但是在地面辐射收入中,射入长波辐射(大气逆辐射)小,仅占平原地区的50%~70%[29],而净长
波辐射却比平原地区高得多[29],因此相对高的反照率和净长波辐射会使净辐射的量相应减少,相应地使净辐射
对蒸散的贡献下降。
3.3 三江源退化草地的水分收支情况
水循环过程中,降水是输入项,而蒸散、渗透、径流和植被截留是输出项,以此来完成水分收支的平衡。在陆
地水循环中,蒸散是生态系统水分散失最大的一个分量,大约有70%归还到地面的降水通过蒸腾和蒸发过程散
失[30]。研究区域3年蒸散量大致与降水量持平,退化草地蒸散量与降水之比为95.6%~98.2%,生长季ET/
PPT为75.8%~82.4%,全年ET/PPT大于青海海北地区矮嵩草草甸[15](57.6%~64.0%),内蒙古克氏针茅
(犛狋犻狆犪犽狉狔犾狅狏犻犻)、寸草(犆犪狉犲狓犱狌狉犻狌狊犮狌犾犪)、冷蒿(犃狉狋犲犿犻狊犻犪犳狉犻犵犻犱犪)占优势的草原[12],可能的原因是由于退化
草地植被退化,植被盖度低,水分涵养能力下降,大部分降水通过植被蒸腾和土壤蒸发散失了。但这个值在已发
表的草地生态系统ET/PPT的范围内(40%~150%)[12]。
退化草地由于地势平坦,每次降水量较小,因此,径流项可以忽略。生长季ET/PPT较低是由于在植被盖
度、根系分布、土壤颗粒结构、土壤透水性等因素的共同作用下,引起水分渗透[31],使土壤含水量上升。此外,在
退化草地由于植物根系多分布在0~10cm土层,根系的存在一定程度上减少土壤表面水分散失。
相比于青藏高原东部矮嵩草草甸[15]、加拿大北部温性草原[2]的蒸散研究,三江源退化草地生长季蒸散量较
低(1.81~1.92mm/d),而非生长季蒸散量却较高(0.39~0.59mm/d)。这是由于生长季降水丰沛,可利用能量
高,三江源退化草地蒸散值在生长季内随着水热条件转好,植被覆盖度逐渐上升,植被蒸腾也随之加剧,而土壤表
面的蒸发降低[15]。但即使在生长盛期,植被盖度也只有75%,因此减少了生长季的蒸腾量,土壤蒸发降低的值也
有限;在非生长季该地区降水量比海北地区小,积雪量较少。在海北矮嵩草草甸、金露梅灌丛草甸、藏嵩草沼泽草
甸积雪对水汽通量影响的研究[32]表明,有积雪时水汽浓度显著高于无积雪。此外,矮嵩草草甸非生长季VPD高
于三江源地区(0.43~0.48kPa)。由这些综合因素可以推出矮嵩草草甸非生长季蒸散应该大于三江源退化草
地,但本研究中却得出相反的结果,这可能是由于草地的退化引起下垫面特征的改变,同时冬季的放牧活动使残
留的枯草因牲畜践踏、风折而发生风蚀[10],减少了植被对水分的截留,使近地面的土壤水分在天气状况较好的情
况下迅速散失[13]。
4 结论
通过3年的蒸散研究,初步表明青海省三江源地区小嵩草草甸退化草地年蒸散量基本与年降水量持平,2006
-2008年年蒸散量差异不大。生长季蒸散占全年蒸散量的75%以上,而非生长季蒸散量较少。相比于未退化的
海北矮嵩草草甸,该地区生长季蒸散略小而非生长季蒸散值略大:生长季退化草地的植被盖度较低,即使是生长
盛期地表也有裸露土壤,而非生长季净辐射对蒸散值影响较大,放牧活动也对这种结果有一定的作用。蒸散日变
化,季节变化都呈单峰型。其中每天正午左右,每年6-7月出现最大值,这种变化表明三江源退化草地蒸散的季
节动态与植被特征和降水的季节性分布相关。
参考文献:
[1] BrutsaertW.EvaporationintotheAtmosphere:Theory,History,andApplications[M].Boston:KluwerAcademicPublic
ers,1982.
[2] WeverLA,FlanaganLB,CarlsonPJ.Seasonalandinterannualvariationinevapotranspiration,energybalanceandsurface
conductanceinanortherntemperategrassland[J].AgriculturalandForestMeteorology,2002,112(1):3149.
132第21卷第3期 草业学报2012年
[3] AdamsJM,FaureH,FaureDenardL,犲狋犪犾.IncreasesinterrestrialcarbonstoragefromtheLastGlacialMaximumtothe
present[J].Nature,1990,348:711714.
[4] 马致远.三江源地区水资源的涵养和保护[J].地球科学进展,2004,19(13):108111.
[5] 陈孝全,苟新京,高延林,等.三江源自然保护区生态环境[M].西宁:青海人民出版社,2002.
[6] 王,洪绂曾,宗锦耀.三江源地区草地资源现状及持续利用途径[J].草地学报,2005,(1):2831.
[7] 李博.中国北方草地退化及其防治对策[J].中国农业科学,1997,30(6):19.
[8] 刘纪远,徐新良,邵全琴.近30年来青海三江源地区草地退化的时空特征[J].地理学报,2008,63(4):364376.
[9] 王启基,来德珍,景增春,等.三江源区资源与生态环境现状及可持续发展[J].兰州大学学报 (自然科学版),2005,41(4):
3137.
[10] 曹广民,龙瑞军.三江源区“黑土滩”型退化草地自然恢复的瓶颈及解决途径[J].草地学报,2009,(1):49.
[11] BaldocchiDD.Assessingtheeddycovariancetechniqueforevaluatingcarbondioxideexchangeratesofecosystems:past,
presentandfuture[J].GlobalChangeBiology,2003,9(4):479492.
[12] LiSG,AsanumaJ,KotaniA,犲狋犪犾.EvapotranspirationfromaMongoliansteppeundergrazinganditsenvironmentalcon
straints[J].JournalofHydrology,2007,333(1):133143.
[13] BremerDJ,AuenLM,HamJM,犲狋犪犾.Evapotranspirationinaprairieecosystem:effectsofgrazingbycattle[J].Agrono
myJournal,2001,93(2):338348.
[14] MeyersTP.AcomparisonofsummertimewaterandCO2fluxesoverrangelandforwelwateredanddroughtconditions[J].
AgriculturalandForestMeteorology,2001,106:205214.
[15] GuS,TangY,CuiX,犲狋犪犾.CharacterizingevapotranspirationoverameadowecosystemontheQinghaiTibetanPlateau[J].
JournalofGeophysicalResearch,2008,113:D08118,doi:10.1029/2007JD009173.
[16] 赵亮,古松,周华坤,等.青海省三江源区人工草地生态系统CO2 通量[J].植物生态学报,2008,32(3):544554.
[17] 李春,何洪林,刘敏,等.ChinaFLUXCO2 通量数据处理系统与应用[J].地球信息科学,2008,10(5):557565.
[18] 温学发,于贵瑞,孙晓敏.基于涡度相关技术估算植被/大气间净CO2 交换量中的不确定性[J].地球科学进展,2004,
19(4):658663.
[19] ZhaoL,LiY,XuS,犲狋犪犾.Diurnal,seasonalandannualvariationinnetecosystemCO2exchangeofanalpineshrublandon
QinghaiTibetanplateau[J].GlobalChangeBiology,2006,12(10):19401953.
[20] FalgeE,BaldocchiD,OlsonR,犲狋犪犾.Gapfilingstrategiesforlongtermenergyfluxdatasets[J].AgriculturalandForest
Meteorology,2001,107(1):7177.
[21] LiY,ZhouL,XuZ,犲狋犪犾.Comparisonofwatervapour,heatandenergyexchangesoveragriculturalandwetlandecosys
tems[J].HydrologicalProcesses,2009,23(14):20692080.
[22] HaoY,WangY,HuangX,犲狋犪犾.Seasonalandinterannualvariationinwatervaporandenergyexchangeoveratypical
steppeinInnerMongolia,China[J].AgriculturalandForestMeteorology,2007,146(12):5769.
[23] FrankAB.Evapotranspirationfromnorthernsemiaridgrasslands[J].AgronomyJournal,2003,95(6):15041509.
[24] HansonRL.Evapotranspirationanddroughts[J].U.S.GeologicalSurveyWaterSupplyPaper,1988,2375:99104.
[25] SongJ.Phenologicalinfluencesonthealbedoofprairiegrasslandandcropfields[J].InternationalJournalofBiometeorology,
1999,42(3):153157.
[26] LafleurPM,RouseWR,CarlsonDW.Energybalancedifferencesandhydrologicimpactsacrossthenortherntreeline[J].
InternationalJournalofClimatology,1992,12(2):193203.
[27] 冯超,古松,赵亮,等.青藏高原三江源区退化草地生态系统的地表反照率特征[J].高原气象,2010,29(1):7077.
[28] GaoF,SchaafCB,StrahlerAH,犲狋犪犾.MODISbidirectionalreflectancedistributionfunctionandalbedoclimatemodeling
gridproductsandthevariabilityofalbedoformajorglobalvegetationtypes[J].JounalofGeophysicalResearch,2005,110:
D01104,doi:10.1029/2004JD005190.
[29] 季国良,江灏,吕兰芝.青藏高原的长波辐射特征[J].高原气象,1995,14(4):451458.
[30] RosenbergNJ,BladBL,VermaSB.Microclimate:TheBiologicalEnvironment[M].2nded.NewYork:JohnWiley&
Sons,1983.
232 ACTAPRATACULTURAESINICA(2012) Vol.21,No.3
[31] 王一博,王根绪,吴青柏,等.植被退化对高寒土壤水文特征的影响[J].冰川冻土,2010,32(5):989998.
[32] 赵亮,徐世晓,伏玉玲,等.积雪对藏北高寒草甸CO2 和水汽通量的影响[J].草地学报,2005,13(3):242247.
犆犺犪狉犪犮狋犲狉犻狕犻狀犵狋犺犲犲狏犪狆狅狋狉犪狀狊狆犻狉犪狋犻狅狀狅犳犪犱犲犵狉犪犱犲犱犵狉犪狊狊犾犪狀犱犻狀狋犺犲
犛犪狀犼犻犪狀犵狔狌犪狀犚犲犵犻狅狀狅犳犙犻狀犵犺犪犻狆狉狅狏犻狀犮犲
LIJingmei1,2,3,CAIHai1,2,CHENGQian4,QIAOChunlian1,2,CHUHui1,2,
CHENDongdong1,XUShixiao1,ZHAOXinquan1,ZHAOLiang1
(1.NorthwestInstituteofPlateauBiology,ChineseAcademyofSciences,Xining810008,
China;2.GraduateUniversityofChineseAcademyofSciences,Beijing100049,China;
3.QinghaiEngineeringConsultingCenter,Xining810001,China;4.Department
ofMathematics,QinghaiNormalUniversity,Xining810001,China)
犃犫狊狋狉犪犮狋:Theevapotranspiration(ET)ofdegradedgrasslandintheSanjiangyuanregionwasinvestigatedto
furtherelucidatethewatercycleinthisregion.Theeddycovariancetechniquewasusedtomeasureevapotrans
pirationandsensibleheatfluxcontinuouslyfor3yearsinthisdegradedgrasslandinQinghaiprovince,China.
Theannualprecipitations(PPT)were460.7,496.1,and480.1mmforthethreeconsecutiveyears.Theannu
alETswere452.24,474.24,and459.57mm,varyinginaccordwiththeannualPPT.About75%ofthean
nualEToccurredduringthegrowingseasonfromDOY121to300.TheratioofETtoPPTwasgreaterthan
95%duringthestudyperiod,indicatingthatthereissomebalancebetweenETandPPTandtheaverageinthe
growingseasonforeachofthethreeyearswere1.86,1.92,1.81mm/d,respectively.However,theaveraged
ETwaslowerthan0.6mm/dduringthenongrowingseason.Theannualanddailydynamicsweresinglepeak
typecurves.ThemaximummonthlyETratewas2.0-2.2mm/d,whichoccurredmainlyinJuneorJuly.The
dailypeakappearedatnoonandwas0.21and0.10mm/hduringthegrowingandnongrowingseasonsrespec
tively.ETwasmostsensitivetonetradiationduringthegrowingseasonwhiletheeffectofairtemperatureand
vapourpressuredeficitweresecondary.Comparedwiththenondegradedareaof犓狅犫狉犲狊犻犪犺狌犿犻犾犻狊meadow,
theETofthedegradedgrasslandwaslowerinthegrowingseasonbuthigherinthenongrowingseason.The
lowervegetationcoverwasattributedtothisresult.
犓犲狔狑狅狉犱狊:eddycovariancemethod;evapotranspiration;Sanjiangyuanregion;Qinghaiprovince;degraded
grassland
332第21卷第3期 草业学报2012年