

Evapotranspiration simulation and its spatio-temporal variation characteristics in Fenqin Region
HANXiang

通讯作者:
收稿日期:2018-03-22
修回日期:2018-05-30
网络出版日期:2018-08-25
版权声明:2018《资源科学》编辑部《资源科学》编辑部
基金资助:
作者简介:
-->
展开
摘要
关键词:
Abstract
Keywords:
-->0
PDF (14228KB)元数据多维度评价相关文章收藏文章
本文引用格式导出EndNoteRisBibtex收藏本文-->
1 引言
蒸散是水循环的重要过程,将全球陆地60%~65%的降水从地表输送到大气[1,2];也是气候系统中连接大气和陆面的重要环节,50%的太阳辐射通过蒸散进行对流交换[3,4]。蒸散将水文循环、能量收支等紧密联系起来,在区域水热平衡中有重要的指标作用[1,2,3,4]。1880—2012年期间全球平均温度升高了0.85℃(0.65℃~1.06℃),增温对蒸散的影响效应被学术界所关注,至今已取得了大量成果[4,5,6,7]。如Jung等利用通量观测、遥感和地面气象观测资料研究发现1998—2010年南半球蒸散以减少趋势为主,1982—2008年北半球以增加或者无明显变化趋势为主[4]。而也有研究表明,气候变暖背景下全球1982—2012年蒸散量呈上升趋势,可能由植被覆盖、辐射、水分等因子导致[6,7]。
区域大尺度蒸散较难直接获取其观测值,目前主要通过模型模拟得到,如Noah陆面模式[8]、Penman-Monteith模型[9]、Priestley-Taylor模型[10]、水量平衡方法(Water Balance Method)[11,12,13]等。结合地表能量平衡的蒸散模型受下垫面不均匀、大气平流等干扰,区域能量可能不完全闭合,而水量平衡方法不受大气复杂变化的影响,对各区域有普遍适用性[11,12,13]。该方法通常在年或更长的时间尺度上假设水储量不变,而在年内尺度或人类活动剧烈的地区,不考虑水储量变化可能对认识和理解区域蒸散时空变化特征产生一定偏差[14,15]。2002年以来,Gravity Recovery and Climate Experiment(GRACE)卫星反演的水储量变化数据的开发与应用,为利用水量平衡方法估算全球陆面蒸散提供了有效途径[14,15,16,17,18,19]。如Zeng等基于该方法研究发现1982—2009年全球年均陆地蒸散为604mm,且以1.10mm/a的趋势增加,其中亚马逊流域和亚洲东南部增加速率最高[14]。Wan等研究2002年4月—2013年4月美国月蒸散量发现,地区蒸散由东南向西北呈大致减少趋势[19]。水量平衡方法数据易获取且计算简便,在大尺度地区特别是水文、气象数据较难获取区域的研究中具有较高应用价值。
位于黄土高原的汾沁地区(汾河流域和沁河流域)受到增温等气候变化、引沁入汾和退耕还林还草等人类活动的影响显著,水热过程复杂[20,21]。贺添等运用MOD16产品分析2001—2010年中国蒸散,其中汾沁地区年均蒸散由北向南增加,大多区域呈减少趋势[22]。而Liu等运用SPS(Simple two-Parameter steady-State)模型估算黄河流域1961—2006年蒸散及Gao等模拟黄土高原1990—2014年蒸散结果均表明,汾沁地区蒸散呈由东南向西北大致呈减少趋势[23,24]。Zhao等运用遥感获取的植被指数和观测的气候数据分别模拟黄土高原2000—2012年植被和土壤的蒸散,研究发现汾沁地区季节蒸散均由北向南大致增加,春季、秋季蒸散多呈增加趋势,夏季蒸散北部多呈增加趋势而南部地区有减少趋势,冬季蒸散多呈减少趋势[25]。Zeng等基于水量平衡模拟全球蒸散结果表明汾沁地区1998—2009年蒸散变化速率为(-1~4)mm/a,大多地区呈增加趋势[14],而Zhang等基于Penman-Monteith模型和互补关系模拟黄河流域蒸散结果表明,汾沁地区1960—2005年蒸散呈减少趋势,其速率为(-2.4~-0.8)mm/a[26]。近年来,汾沁地区不同时间尺度蒸散量数据仍较为缺乏,其空间分布及变化速率亦有待进一步明确。加强人类活动剧烈地区更高精度的蒸散模拟、时空变化特征及区域蒸散变化影响关键因子分析是对已有研究的重要补充。本文通过结合GRACE数据的水量平衡方法,研究和评价2003—2015年汾沁地区蒸散在不同时间尺度的变化特征,对深入理解陆面水循环对气候变化和人类活动的响应具有重要意义,并可为地区水资源规划、生态和环境保护等提供科学支撑。
2 研究区概况、数据来源与研究方法
2.1 研究区概况
汾沁地区主要包括汾河流域和沁河流域(图1)。汾河和沁河是山西省最大的两条河流,也是黄河 的重要支流。汾沁地区介于110.7°E—113.6°E和34.9°N—39.0°N,位于吕梁山与太行山之间,地区总面积为52 538.58km2,占黄河流域面积的6.98%。研究区大部分位于山西省境内(51 218.2km2),少部分位于河南省(1320.38km2)。汾沁地区位于黄土高原东部,地势由西北向东南大致降低,地形较为起伏,海拔91~2799m,平均海拔1095m。属于中纬度半干旱、半湿润型气候过渡区,为暖温带大陆性季风气候,雨热同期,四季分明。2003—2015年均降水为547.00mm,年均温为9.18℃,平均风速为2.28m/s。地区植被状况整体较好,植被覆盖度由西北向东南逐渐增加。大部分地区被落叶阔叶林、农业植被覆盖,也有少量针叶林、针阔混交林、灌木丛等[20,27,28]。
图1汾沁地区海拔及主要水文站点
-->Figure 1The altitude and major hydrological stations in Fenqin Region
-->
该地区水文站点主要有静乐、上静游、汾河水库、寨上、兰村、芦家庄、汾河二坝、文峪河水库、义棠、赵城、柴庄、浍河水库、河津、孔家坡、飞岭、润城、五龙口、山路平、武陟,共19个(图1)。根据DEM(Digital Elevation Model,数字高程模型)提取汾沁地区子流域,其中汾河流域划分为静乐以上区间、岚河流域、静乐-汾河水库区间、汾河水库-寨上区间、寨上-兰村区间、潇河上游芦家庄区间、兰村-汾河二坝区间、文峪河流域、汾河二坝-义棠区间、义棠-赵城区间、赵城-柴庄区间、浍河流域、柴庄-河津区间;沁河流域划分为孔家坡以上区间、孔家坡-飞岭区间、飞岭-润城区间、润城-五龙口区间、丹河流域、五龙口-武陟区间[29]。
2.2 数据来源
2.2.1 遥感数据GRACE数据选用在国际上认可度较高的美国德克萨斯大学奥斯汀分校空间研究中心CSR(Center of Space Research),德国地学研究中心GFZ(Geo Forschungs Zentrum)、美国宇航局喷气推进实验室JPL(Jet Propulsion Laboratory)三种产品,数据来源于http://grace.jpl.nasa.gov,时间分辨率为月,空间分辨率为1°。月水储量用三种产品的均值来表示,少数月份数据缺失,采用滑动平均法插补。GRACE卫星监测水储量变化的原理是由于地球系统的物质质量重新分布导致不同时间尺度的地球重力场变化,从而利用足够精度和时空分辨率的重力观测量分析其物质迁移和交换[30,31]。一般认为,在季节性或其他短时间尺度上(如月、周),陆地水储量变化是引起地球重力场变化的主要原因[16,17]。
GLDAS(Global Land Data Assimilation System,全球陆面数据同化系统) NOAH(National Centers for Environmental Prediction/Oregon State University/Air Force/Hydrologic Research Lab Model)产品基于Noah陆面模式[8,32],来源于https://earthdata.nasa.gov,时间分辨率为月,空间分辨率为0.25°。本文提取汾沁地区0~10cm、10~40cm、40~100cm与100~200cm土壤水、雪、植物冠层水、蒸散数据。
MOD16A2(Moderate Resolution Imaging Spectroradiometer Global Evapotranspiration Data Set,MODIS全球蒸散产品)基于改进的Penman-Monteith模型[9],数据来源于https://lpdaac.usgs.gov,时间分辨率为月,空间分辨率为0.05°。
GLEAM(Global Land Evaporation Amsterdam Model) v3.1b数据基于Priestley-Taylor模型[10,33],来源于https://www.gleam.eu,时间分辨率为日,空间分辨率为0.25°。
高程数据来源于CGIAR-CSI网站http://srtm.csi.cgiar.org/,数据以5°×5°为一图幅组织数据,各图幅之间为无缝拼接,空间分辨率为90m。
2.2.2 实测数据
汾沁地区19个站点2003—2015年径流数据来源于黄河流域水文年鉴[34]。其中2003年、2004年、2006年、2007年部分站点数据缺失,用最小二乘法插补。
2003—2015年月尺度降水、平均温度、风速、相对湿度、水汽压、日照时数数据来源于国家气象科学数据共享服务平台[35]。其中降水数据的空间分辨率为0.5°;平均温度、风速、相对湿度、水汽压、日照时数数据采用五台山、五寨、兴县、原平、离石、太原、阳泉、榆社、隰县、介休、临汾、长治、运城、侯马、阳城、三门峡、孟津、郑州等18个站点数据,地区90mDEM作为辅助资料,运用ANUSPLIN[36]建立汾沁地区平均温度、风速、相对湿度、水汽压、日照时数的年及月值数据集。
2.3 研究方法
2.3.1 水量平衡方法根据水量平衡原理,在一个闭合流域内,如不考虑相邻区域的水量调入与调出,水量收入和支出的差额,等于该地区的水储量变化量[37]。一般的流域水平衡方程可表达为:
式中P为降水量(mm);ET为蒸发量(mm);Q为径流深(mm);ΔS为水储量变化(mm)。蒸散量即等于进入流域的水资源量与损失部分的差,即:
通常在年或更长的时间尺度上水储量变化较小,可将其忽略不计[19,38],蒸散量(ET-NGRACE)即为降水与径流之差。
2.3.2 降尺度方法
本文将GRACE反演的分辨率为1°的水储量 变化(当月水储量与前月水储量之差)降尺度为0.25°[19]。首先,通过GLDAS数据获取研究区分辨率为0.25°的水储量数据(不同深度土壤水、雪、植物冠层水之和),并计算月水储量变化。将GLDAS获取的水储量变化数据合并为1°,其公式如下:
式中ΔSGL为基于GLDAS计算的1°水储量变化量(mm);ΔSgl为基于GLDAS的0.25°水储量变化量(mm);ai为0.25°栅格的面积(m2),A为1°栅格的面积(m2)。
每度GLDAS与GRACE的水储量变化量之间的差值B为:
式中ΔSGR为1°GRACE水储量变化量(mm)。将差值B分配到每0.25°,即差值bi:
式中ΔSglmin表示基于GLDAS的水储量变化量每度的最小值。将bi加到原基于GRACE估算的水储量变化数据,即为降尺度后分辨率为0.25°的水储量变化量ΔS′:
这种降尺度的方法不但保证了计算结果的精度,而且降尺度的GRACE水储量变化与原水储量变化在对应的1°网格的值相同。
2.3.3 数理统计
(1)时间序列数据对比。相关系数是指在消除其他变量影响的条件下单独研究两个变量间的相关程度的方法。本文用该方法表示结合GRACE数据的水量平衡方法模拟的汾沁地区蒸散(ET-GRACE)、与基于GLDAS产品的区域蒸散(ET-GLDAS)、基于MOD16A2产品的区域蒸散(ET-MOD16)、基于GLEAM数据的区域蒸散(ET-GLEAM)的相关性,公式为:

式中yi为第i月蒸散量;
方差S、变异系数CV、极差R来衡量数据波动大小的量:

式中
用偏度与峰度描述蒸散序列的分布状态。偏度大于(小于)0表示数据呈右(左)偏分布,偏度越大(小),右(左)偏程度越高。峰度的绝对值越大,表示其与正态分布的差异程度越大。
(2)趋势分析及显著性检验。线性回归是利用线性回归方程的最小二乘函数对变量之间关系进行建模的回归分析。本文用该方法分析2003—2015年汾沁地区蒸散、气候因子的变化趋势:
式中b为变化速率;x为时间序列;y为时间序列对应蒸散或气候因子的值;j为年份。
用Mann-Kendall检验判断蒸散变化趋势的显著性,序列检验结果的置信度达到0.1(即p<0.1)则变化趋势显著。
(3)偏相关分析。偏相关分析指在控制其他变量的线性影响下分析两个变量间的线性相关性。本文用该方法分析时间序列上气候因子(降水、平均温度、风速、湿度、水汽压、日照时数)对蒸散的影响时,考虑到因子间相互影响,在控制某两个气候因子的基础上,单独分析某一气候因子对蒸散的影响,得到偏相关系数。
3 结果及分析
3.1 蒸散结果对比
将结合GRACE数据的水量平衡方法模拟所得汾沁地区蒸散分别与ET-GLDAS、ET-MOD16、ET-GLEAM对比(图2),相关性分别为0.81(p<0.01)、0.74(p<0.01)、0.78(p<0.01),一致性较好。与ET-NGRACE相比,ET-GRACE的标准差、变异系数、极差均较小,分别为41.42mm、0.94、228.45mm,离散程度较低;偏度和峰度较小,分别为0.90、1.01,更符合正态分布。ET-NGRACE中极端值占7.05%,ET-GRACE中极端值占3.85%,ET-GRACE比ET-NGRACE的极端值减少3.20%(表1)。
图22003—2015年五种常用月蒸散产品对比
-->Figure 2Comparisons of monthly ET products using commonly from 2003 to 2015
-->
Table 1
表1
表12003—2015年汾沁地区ET-GRACE与ET-NGRACE计算结果统计对比
Table 1Statistics comparisons between results of ET-GRACE and ET-NGRACE in Fenqin Region from 2003 to 2015
标准差/mm | 变异系数 | 极差/mm | 偏度 | 峰度 | 极端值/% | |
---|---|---|---|---|---|---|
ET-GRACE | 41.42 | 0.94 | 228.45 | 0.90 | 1.01 | 3.85 |
ET-NGRACE | 46.92 | 1.06 | 265.08 | 1.48 | 2.69 | 7.05 |
新窗口打开
在多年平均尺度上,结合GRACE数据的水量平衡方法所模拟的2003—2015年均蒸散量为530.19mm,不考虑水储量变化的水量平衡方法模拟所得年均蒸散量为529.82mm,两者相差不大。年尺度上,ET-GRACE呈增加趋势(1.13mm/a,p>0.1)而ET-NGRACE呈减少趋势(-1.18mm/a,p>0.1),且前者波动更平稳(图3)。季节尺度上,汾沁地区2003—2015年春季ET-GARCE多高于ET-NGRACE,两者均呈增加趋势,变化速率分别为1.12mm/a(p>0.1)、1.52mm/a(p>0.1)(图4a);夏季,ET-NGRACE常年高于ET-GRACE,时间序列上均呈减少趋势,变化速率差异较大,分别为-1.16mm/a(p>0.1)、-2.85mm/a(p>0.1)(图4b);秋季,ET-GRACE与ET-NGRACE均呈增加趋势,变化速率分别为2.13mm/a(p>0.1)、0.75mm/a(p>0.1),该差异与2003年秋季ET-NGRACE相对较高相关(图4c);冬季,ET-GRACE多高于ET-NGRACE,两者均呈减少趋势,变化速率分别为-0.95mm/a(p>0.1)、-0.60mm/a(p>0.1)(图4d)。可见,年及季节等较短时间尺度上,运用水量平衡方法模拟汾沁地区蒸散时,需考虑水储量变化。

图32003—2015年汾沁地区年尺度ET-GRACE和ET-NGRACE变化趋势对比
-->Figure 3Comparisons of temporal variations between ET-GRACE and ET-NGRACE at annual scale in Fenqin Region from 2003 to 2015
-->

图42003—2015年汾沁地区季节尺度ET-GRACE和ET-NGRACE变化趋势对比
-->Figure 4Comparisons of temporal variations between ET-GRACE and ET-NGRACE at seasonal scale in Fenqin Region from 2003 to 2015
-->
3.2 多时间尺度上蒸散的空间分布特征
3.2.1 年蒸散空间分布特征结合GRACE数据反演的水储量变化,基于水量平衡方程模拟汾沁地区2003—2015年平均蒸散量为530.19mm(表2)。研究区蒸散量空间差异明显,由北向南大致呈增加趋势。汾河流域静乐以上区间多年平均蒸散量最低为378.82mm,丹河流域多年平均蒸散量最高为641.02mm。沁河流域多年平均蒸散量为577.18mm,高于汾河流域(513.67mm)。
Table 2
表2
表22003—2015年汾沁地区年及季节尺度平均蒸散量
Table 2Annual and seasonal mean evapotranspiration in Fenqin Region from 2003 to 2015 (mm)
全年 | 春季 | 夏季 | 秋季 | 冬季 | |
---|---|---|---|---|---|
汾沁地区均值 | 530.19 | 123.93 | 263.36 | 123.39 | 19.50 |
汾河流域均值 | 513.67 | 113.74 | 258.57 | 124.34 | 17.02 |
静乐以上区间 | 378.82 | 65.21 | 245.31 | 87.04 | 0.00 |
岚河流域 | 503.84 | 88.74 | 276.35 | 125.39 | 13.36 |
静乐-汾河水库区间 | 584.44 | 54.40 | 304.32 | 180.97 | 44.74 |
汾河水库-寨上区间 | 548.68 | 111.05 | 291.67 | 132.68 | 13.27 |
寨上-兰村区间 | 590.31 | 119.09 | 305.88 | 144.27 | 21.07 |
潇河上游芦家庄区间 | 523.93 | 113.00 | 295.97 | 101.74 | 13.23 |
兰村-汾河二坝区间 | 487.79 | 118.17 | 259.39 | 105.50 | 4.74 |
文峪河流域 | 462.62 | 81.18 | 259.50 | 115.35 | 6.58 |
汾河二坝-义棠区间 | 481.55 | 110.89 | 239.20 | 114.31 | 17.15 |
义棠-赵城区间 | 554.13 | 129.47 | 268.38 | 134.90 | 21.38 |
赵城-柴庄区间 | 558.47 | 136.15 | 257.64 | 138.57 | 26.10 |
浍河流域 | 567.54 | 154.55 | 248.67 | 139.16 | 25.16 |
柴庄-河津区间 | 528.89 | 132.35 | 227.70 | 138.43 | 30.42 |
沁河流域均值 | 577.18 | 152.93 | 276.98 | 120.69 | 26.56 |
孔家坡以上区间 | 513.85 | 126.43 | 274.22 | 98.99 | 14.22 |
孔家坡-飞岭区间 | 525.11 | 129.16 | 262.84 | 107.89 | 25.21 |
飞岭-润城区间 | 596.65 | 154.30 | 279.51 | 132.41 | 30.43 |
润城-五龙口区间 | 515.41 | 151.04 | 254.22 | 94.95 | 15.19 |
五龙口-武陟区间 | 580.65 | 141.71 | 277.42 | 116.96 | 44.56 |
丹河流域 | 641.02 | 180.14 | 295.36 | 138.21 | 27.29 |
新窗口打开
3.2.2 季节蒸散空间分布特征
汾沁地区2003—2015年不同季节蒸散量大致均呈现南部高于北部的分布状态,但各季节平均 蒸散差异明显(表2)。春季(3月、4月、5月)地区平均蒸散量为123.93mm,由西北向东南增加的趋势较明显,静乐-汾河水库区间平均蒸散量最低,为54.40mm,丹河流域平均蒸散量最高,为180.14mm。夏季(6月、7月、8月)蒸散最高,为263.36mm。子流域蒸散,量普遍较高,寨上-兰村区间蒸散量最高,为305.88mm,柴庄-河津区间最低,为227.70mm。秋季(9月、10月、11月)为123.39mm,秋季子流域蒸散量为87.04~180.97mm。冬季(12月、1月、2月) 最低,为19.50mm,子流域蒸散普遍偏低,五龙口- 武陟区间最高为44.56mm,静乐以上区间最低为0.00mm。春季、夏季、冬季,汾河流域平均蒸散量(113.74mm、258.57mm、17.02mm)低于沁河流域(152.93mm、276.98mm、26.56mm),而秋季沁河流域平均蒸散(120.69mm)低于汾河流域(124.34mm)。
3.3 蒸散时间变化趋势
3.3.1 年蒸散变化趋势汾沁地区2003—2015年蒸散量呈现不显著增加趋势(1.13mm/a,p>0.1),年际波动较平稳(变异系数为0.08)(图3、图5)。其中,2010年蒸散量最低,为478.22mm,2011年蒸散量最高,为614.57mm。2006年、2007年、2009年、2011年、2013年蒸散量增加,其中2011年增加最多,较2010年增加136.35mm,2009年蒸散仅较2008年增加9.99mm。蒸散量减少的年份中,2012年降低最多,较上一年降低103.59mm,2015降低最少,较2014年降低5.88mm。

图52003—2015年汾沁地区年蒸散变化空间分布示意
注:数字1~19表示子流域代码,详见
-->Figure 5Distribution of trends in annual evapotranspiration in Fenqin Region from 2003 to 2015
-->
汾沁地区子流域2003—2015年蒸散量变化趋势及其显著性分析表明:研究区北部、中部大部分地区及沁河流域入黄口年蒸散量呈现趋势,其中静乐-汾河水库-寨上-兰村区间及文峪河流域增加显著,分别为16.40mm/a(p<0.1),12.15mm/a(p<0.1),13.00mm/a(p<0.1)。汾沁地区南部的赵城-柴庄区间、浍河流域、柴庄-河津区间、飞岭-润城区间、润城-五龙口区间、丹河流域均呈减少趋势,速率分别为-0.49mm/a(p>0.1)、-3.13mm/a(p>0.1)、-0.85mm/a(p>0.1)、-7.63mm/a(p>0.1)、-12.30mm/a(p>0.1)、 -13.51mm/a(p>0.1)。汾河流域年蒸散呈增加趋 势(3.82mm/a(p>0.1)),而沁河流域呈减少趋势 (-6.51mm/a(p>0.1))。从年际波动看,各子流域变异系数为0.12~0.20,年际波动较平稳,沁河流域年际波动(变异系数为0.13)高于汾河流域(变异系数为0.09)。
3.3.2 季节蒸散变化趋势
2003—2015年季节蒸散变化趋势如图4、图6。春季(图4a、图6a),汾沁地区、汾河流域和沁河流域蒸散均呈不显著增加趋势,分别为1.12mm/a(p>0.1)、1.22mm/a(p>0.1)、0.83mm/a(p>0.1)。静乐以上区间、潇河上游芦家庄区间、润城-五龙口区间呈减少趋势,速率分别为-1.93mm/a(p>0.1)、-0.33mm/a(p>0.1)、-0.02mm/a(p>0.1),其他子流域呈增加趋势,静乐-汾河水库区间及五龙口-武陟区间最为显著,分别为4.00mm/a(p<0.1)、6.31mm/a(p<0.1)。

图62003—2015年汾沁地区季节蒸散变化趋势空间分布示意
注:数字1~19表示子流域代码,详见
-->Figure 6Distribution of trends in seasonal evapotranspiration in Fenqin Region from 2003 to 2015
-->
夏季(图4b、图6b),汾沁地区蒸散呈减少趋势,-1.16mm/a(p>0.1),波动较平稳,变异系数为0.18。汾河流域呈增加趋势,0.42mm/a(p>0.1),沁河流域呈减少趋势,-5.68mm/a(p>0.1)。研究区北部及中部蒸散量多呈增加趋势,汾河水库-寨上-兰村区间及文峪河水库增加最为显著,分别为7.63mm/a(p<0.1)、6.16mm/a(p<0.1)、7.64mm/a(p<0.1);南部地区多呈减少趋势,润城-五龙口-武陟区间最显著,分别为-9.41mm/a(p<0.1)、-8.05mm/a(p<0.1)。
秋季(图4c、图6c),汾沁地区、汾河流域、沁河流域蒸散均呈不显著增加趋势,分别为2.13mm/a(p>0.1)、2.69mm/a(p>0.1)、0.52mm/a(p>0.1)。浍河流域、飞岭-润城区间、丹河流域呈减少趋势,其他子流域均呈增加趋势,静乐以上区间、潇河上游芦家庄区间、五龙口-武陟区间增加显著,分别为4.54mm/a(p<0.1)、6.73mm/a(p<0.1)、9.17mm/a(p<0.1)。
冬季(图4d、图6d),汾沁地区、汾河流域、沁河流域蒸散均呈不显著减少趋势,分别为-0.95mm/a(p>0.1)、-0.52mm/a(p>0.1)、-2.17mm/a(p>0.1),蒸散量离散程度均较高。汾河流域上游的岚河流域、静乐-汾河水库-寨上-兰村区间呈增加趋势,其他子流域均呈减少趋势,其中润城-五龙口区间减少显著,为-3.92mm/a(p<0.1)。
3.4 气候因子对蒸散的影响分析
蒸散作为复杂的水热传输过程,受到降水、温度、风速、湿度、水汽压、日照时数等气候因子的综合影响[39,40]。年尺度上(图7,见1667页),汾沁地区2003—2015年降水、风速、相对湿度、水汽压均呈不显著 减少趋势(-1.17mm/a,p>0.1;-0.004m/(s·a),p>0.1;-0.16%/a,p>0.1;-0.001hPa/a,p>0.1),年均温呈不显著增加趋势(0.03℃/a,p<0.1),日照时数呈显著增加趋势(-16.00h/a,p<0.1),其中降水、日照时数波动程度较大(变异系数分别为0.15、0.10)。气候因子变化趋势及波动程度受到个别极值年份影响较大,如2003年降水偏高,该年后降水呈增加趋势。蒸散与降水、平均温度、相对湿度呈正相关,与风速、水汽压、日照时数呈负相关(表3,见1667页)。综合来看,研究区蒸散增加受温度增加、水汽压减少影响较大;2011年蒸散量最高与该年降水较高相关,2010年蒸散量最低与该年降水较低相关。
图72003—2015年汾沁地区气候因子变化趋势
-->Figure 7Trends of climate factors in Fenqin Region from 2003 to 2015
-->
Table 3
表3
表32003—2015年汾沁地区主要气候因子与蒸散的偏相关分析
Table 3Partial correlation between climate factors and evapotranspiration in Fenqin Region from 2003 to 2015
降水 | 平均温度 | 风速 | 相对湿度 | 水汽压 | 日照时数 | |
---|---|---|---|---|---|---|
全年 | 0.28 | 0.55 | -0.55 | 0.41 | -0.24 | -0.23 |
春季 | 0.50 | 0.23 | -0.58 | -0.09 | -0.13 | -0.14 |
夏季 | 0.97* | 0.44 | 0.82* | -0.22 | 0.44 | 0.20 |
秋季 | 0.83* | 0.45 | 0.30 | 0.27 | -0.64* | -0.63* |
冬季 | -0.10 | 0.18 | 0.42 | 0.27 | -0.27 | -0.16 |
新窗口打开
季节尺度上(图7,见1667页),汾沁地区2003—2015年春、秋降水呈增加趋势(1.68mm/a,0.24mm/a),夏、冬降水呈减少趋势(-2.51mm/a,-0.51mm/a),夏季平均降水较高;温度均呈不显著增加趋势,夏季 平均温度最高而冬季最低;春、秋风速呈减少趋势(-0.01m/(s·a),-0.01m/(s·a)),而夏、冬呈增加趋势(0.01m/(s·a),0.01m/(s·a)),春季平均风速较高;春、夏、秋季相对湿度呈减少趋势,冬季呈增加趋势;春、夏、冬水汽压呈减少趋势,而秋季呈增加趋势;春、秋、冬日照时数呈减少趋势(-3.15h/a,-9.74h/a, -3.57h/a),其中秋季减少最显著(p<0.1),而夏季呈增加趋势(0.46h/a)。汾沁地区蒸散受年内水热差异影响较大,由偏相关分析可知(表3):夏季蒸散减少与降水减少、风速增加显著正相关(0.97,p<0.1;0.82,p<0.1);而秋季蒸散增加与降水增加显著正相关(0.83,p<0.1),与水汽压增加、日照时数减少显著负相关(-0.64,p<0.1;-0.63,p<0.1)。
4 讨论
本文模拟结果中存在一些负值,可能与水储量数据的获取及处理过程中有一定误差相关[41]。GRACE数据基于卫星观测地球重力场变化,具有一定的信号误差,其反演得到的水储量数据,除灌溉、跨流域调水、城市化等人类活动引起外,还包括了地壳运动等大尺度非水因素引起的质量变化[16]。目前尚没有统一有效的水文模型用于GRACE数据反演结果的验证[17,42]。一些研究中用水文模型模拟结果与之相比,但相关模型尚不能完全准确地模拟水文变化,且模型应用存在区域差异,验证某一特定区域的水文要素时可能存在争议[43,44,45]。研究GRACE数据的反演结果与实测数据、水文模型模拟数据和遥感数据之间的对比及造成差异的原因,并选择有效的验证方法来降低GRACE数据反演结果的不确定性会是未来研究的重要方向之一。文中结合GLDAS NOAH产品将GRACE数据反演的1°水储量变化降尺度为0.25°。降尺度前后汾沁地区月水储量变化在时间序列上相关性均较高(0.91~1.00),但空间格局上仍有较大差异。该空间分布的不确定性可能与子流域空间尺度较小、降尺度方法有待完善相关[19,46]。尚需进一步考虑基于多源数据及相关模型的降尺度方法研究,以提高降尺度后水储量变化的精度[19]。
本文将运用水量平衡方法的模拟结果与国际广泛应用的GLDAS、MOD16、GLEAM产品对比表明模拟结果有一定的可信度,但模拟结果之间尚有差异。与已有研究结果相比,如Zhang等研究表明1960—2005年汾沁地区蒸散呈减少趋势,可能与该研究时段内降水、风速、日照时数减少等气候因子综合作用相关[26];贺添等研究表明2001—2010年汾沁地区北部蒸散多在300~450mm,低于本文模拟结果,可能由于不同遥感蒸散发模型相关假设及控制因子等不同,导致模拟结果具有差异[22]。总之,目前大尺度蒸散模拟结果仍有一定不确定性,运用通量观测数据验证模拟结果非常重要[47,48]。由于目前汾沁地区相关实测数据相对缺乏,其地面验证部分的研究有待加强。
蒸散时空变化具有一定的不确定性,与其受气候变化和人类活动的共同影响相关[39]。汾沁地区2003年蒸散量较高,特别是忽略水储量变化的水量平衡方法模拟结果蒸散高达698.79mm,可能与该年降水量高达720.42mm,冷空气活动频繁,气温偏低等气候要素的共同影响相关[49]。另外,气候因子与蒸散间并非线性关系,用偏相关研究所得结果不一定是气候因子对蒸散的影响,不能消除蒸散对气候因子的反作用。人类活动引起的土地利用变化、土地水源涵养能力变化、跨流域调水等下垫面条件改变也可作用于蒸散[23,50]。近年来汾沁地区蒸散变化可能与其积极采取退耕还林还草、修建林草等措施改善水源涵养功能[20],及“引沁入汾”实施期间修建张峰水库、河口村水库、五马水库等水利工程有关[51]。蒸散作为复杂、连续的地理过程,受到多重因子的共同作用,今后需进一步明确汾沁地区蒸散的影响机制,为理解及应对区域气候变化提供科学依据。
5 结论
蒸散是全球能量循环、水循环、碳循环的关键环节,研究其时空分布对地区水资源的合理利用及气候变化的深入理解具有重要意义。本文结合GARCE数据反演的水储量变化,基于水量平衡方法模拟了汾沁地区2003—2015年蒸散量并分析其时空变化,主要得到以下结论:(1)年及季节尺度上,水储量变化对水量平衡方法模拟汾沁地区蒸散影响较大。该模拟结果较忽略水储量变化的水量平衡方法模拟结果在时间序列上波动更平稳,可靠性更高(变异系数、标准差、极差分别减少0.12、5.50mm、36.63mm,极端值减少3.20%)。
(2)汾沁地区2003—2015年均蒸散量为530.19mm,大致呈现由北向南增加的趋势;年际波动较平稳(变异系数为0.08),其中2010年蒸散量最低(478.22mm),2011年蒸散量最高(614.57mm)。季节尺度上,汾沁地区夏季平均蒸散量最高(263.36mm)且波动较平稳,冬季蒸散量最低(19.50mm)且离散程度较大。
(3)汾沁地区2003—2015年蒸散变化的主要因素为降水和温度,蒸散年际波动的主要因素为降水。季节尺度上,夏季、秋季蒸散变化与降水相关性最显著(p<0.1)。
The authors have declared that no competing interests exist.
参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
[47] | [J]. [J]. |
[48] | [J]. |
[49] | |
[50] | [J]. |
[51] | [J]. [J]. |
[1] | |
[2] | [J]. |
[3] | [J]. |
[4] | [J]. |
[5] | . |
[6] | [J]. |
[7] | [J]. |
[8] | [J]. |
[9] | [J]. |
[10] | [J]. |
[11] | [J]. |
[12] | [J]. |
[13] | [J]. [J]. |
[14] | [J]. |
[15] | [J]. |
[16] | [J]. |
[17] | [J]. |
[18] | [J]. |
[19] | [J]. |
[20] | [J]. |
[21] | [J]. [J]. |
[22] | [J]. [J]. |
[23] | [J]. |
[24] | [J]. |
[25] | [J]. |
[26] | [J]. |
[27] | [J]. |
[28] | [J]. |
[29] | [J]. [J]. |
[30] | [J]. |
[31] | [J]. |
[32] | [EB/OL]. ( |
[33] | [J]. |
[34] | [M]. [M]. |
[35] | [EB/OL]. ( [EB/OL]. ( |
[36] | [J]. [J]. |
[37] | |
[38] | [J]. |
[39] | [J]. |
[40] | [J]. |
[41] | [J]. |
[42] | [J]. |
[43] | [J]. |
[44] | [J]. |
[45] | [J]. |
[46] | [J]. |