湖泊科学   2019, Vol. 31 Issue (6): 1685-1697.  DOI: 10.18307/2019.0608.
0

研究论文

引用本文 [复制中英文]

刘子豪, 陆建忠, 黄建武, 陈晓玲, 张玲, 盛颖东, 基于CMIP5模式鄱阳湖流域未来参考作物蒸散量预估. 湖泊科学, 2019, 31(6): 1685-1697. DOI: 10.18307/2019.0608.
[复制中文]
LIU Zihao, LU Jianzhong, HUANG Jianwu, CHEN Xiaoling, ZHANG Ling, SHENG Yingdong. Prediction and trend of future reference crop evapotranspiration in the Poyang Lake Basin based on CMIP5 Models. Journal of Lake Sciences, 2019, 31(6): 1685-1697. DOI: 10.18307/2019.0608.
[复制英文]

基金项目

武汉市应用基础前沿专项(2019020701011502)、国家自然科学基金项目(41971402)、湖北省自然科学基金项目(2019CFB736)、中央高校基本科研业务费专项资金(2042018kf0220)、江西省水工程安全与资源高效利用工程研究中心开放基金重点项目(OF201601)和测绘遥感信息工程国家重点实验室专项科研经费联合资助

作者简介

刘子豪(1995~), 男, 硕士研究生; E-mail:18827029266@163.com

通信作者

陆建忠, E-mail:lujzhong@whu.edu.cn

文章历史

2019-04-11 收稿
2019-04-17 收修改稿

码上扫一扫

基于CMIP5模式鄱阳湖流域未来参考作物蒸散量预估
刘子豪1 , 陆建忠2 , 黄建武1 , 陈晓玲2 , 张玲2 , 盛颖东1     
(1: 华中师范大学地理过程分析与模拟湖北省重点实验室, 武汉 430079)
(2: 武汉大学测绘遥感信息工程国家重点实验室, 武汉 430079)
摘要:预测未来气候情境下鄱阳湖流域参考作物蒸散量(Reference crop Evapotranspiration,ET0)的时空分布可为流域水资源的优化管理,为科学应对气候变化对农业生产的影响提供基础数据支撑.利用鄱阳湖流域14个气象站点1961-2014年逐日气象数据,采用Penman-Monteith公式计算出历史ET0;基于同期美国环境中心(NCEP)再分析数据及2006-2100年CMIP5中CNRM-CM5模式在RCP4.5和RCP8.5情景下的预测数据,经统计降尺度模型(statistical downscaling model,SDSM)模拟和偏差校正,预测流域未来ET0;通过Mann-Kendall检验、普通克里金插值和空间自相关法分析了流域1961-2100年ET0的时空演变特征.结果表明:NCEP再分析资料与流域ET0建立的逐步回归降尺度模型模拟效果较好,CNRM-CM5模式降尺度模拟结果经偏差校正后,精度明显提高,适宜流域未来ET0的预估.鄱阳湖流域在基准期1961-2010年ET0整体上呈减小趋势,空间分布上呈南北高、中间低的特点,表现出明显的空间差异性.RCP4.5、RCP8.5情景下未来3个时期鄱阳湖流域ET0较基准期均呈不同程度的增加趋势,其空间分布整体表现为东高西低、局地略有突出;无论是在基准期或是未来情景下的3个时期,ET0均具有较强的空间自相关性.在RCP8.5情景下,鄱阳湖1961-2100年干旱指数呈现出较为明显的上升趋势,流域的干旱状况随时间加剧,2011-2100年间流域绝大部分地区由湿润区转为半湿润区,干旱指数自南向北递减,赣江流域将是鄱阳湖流域未来干旱风险的重点防范区.
关键词蒸散    CMIP5    统计降尺度    未来气候    干旱指数    参考作物蒸散量    鄱阳湖流域    
Prediction and trend of future reference crop evapotranspiration in the Poyang Lake Basin based on CMIP5 Models
LIU Zihao1 , LU Jianzhong2 , HUANG Jianwu1 , CHEN Xiaoling2 , ZHANG Ling2 , SHENG Yingdong1     
(1: Key Laboratory of Geographic Process Analysis and Simulation of Central China Normal University, Wuhan 430079, P. R. China)
(2: State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan 430079, P. R. China)
Abstract: Predicting the temporal and spatial distribution of Reference crop Evapotranspiration(ET0) in the Poyang Lake Basin in the future climate scenario can provide optimal management of water resource in the watershed and provide basic data support for scientific response to the impact of climate change on agricultural production. In this paper, the daily meteorological data of 14 meteorological stations in Poyang Lake Basin from 1961 to 2014 were used to calculate the historical ET0 by the Penman-Monteith formula; Based on the reanalysis data of the US Environmental Center(NCEP) and the CNRM-CM5 model in CMIP5 from 2006 to 2100. The prediction data in the RCP4.5 and RCP8.5 scenarios are predicted by the statistical downscaling model(SDSM) and the bias correction to predict the future ET0 of the basin; Through M-K test, ordinary kriging interpolation and spatial autocorrelation, the temporal and spatial evolution characteristics of ET0 in the basin from 1961 to 2100 were analyzed. The results show that:the simulation effect of NCEP reanalysis data and the stepwise regression scale reduction model established by ET0 of the basin is good, and the accuracy of the scale reduction simulation results of CNRM-CM5 model is significantly improved after the deviation correction, which is suitable for the estimation of future ET0 of the basin. The overall ET0 of Poyang Lake Basin during the base period from 1961 to 2010 shows a decreasing trend, and the spatial distribution is high in the north and south, and low in the middle, showing obvious spatial differences. Under the scenario of RCP4.5 and RCP8.5, the ET0 of Poyang Lake Basin in the next three periods shows an increasing trend in different degrees compared with the base period. ET0 has strong spatial autocorrelation in both benchmark period and future scenario. In RCP8.5 scenario, the aridity index from 1961 to 2100 in Poyang Lake Basin, showing a more obvious rise and increased with time, the valley of drought conditions in most of river basin from 2011 to 2100 by the humid to semi-humid areas, aridity index decreases from south to north, the Gan river basin will be the focus of the Poyang Lake Basin in the future drought risk prevention area.
Keywords: Evapotranspiration    CMIP5    statistical downscaling    future climate    aridity index    reference crop evapotranspiration(ET0)    Poyang Lake Basin    

在全球气候变化的背景下, 人类赖以生存的粮食安全、水土资源、生态环境等诸多方面受到巨大挑战[1-2].政府间气候变化专门委员会(Intergovernmental Panel on Climate Change, IPCC)第5次评估报告指出[3], 自1880 -2012年, 全球表面平均气温上升约0.85℃, 对全球的气候和环境尤其是水循环产生重要影响.参考作物蒸散量(Reference crop Evapotranspiration, ET0)作为连接气象要素和水循环的重要物理量, 是水文模型研究区域水量平衡的一项重要参数, 是评价区域内水资源配置和计算作物生产潜力的重要指标[4], ET0的变化直接关系到流域的干湿状况和植被用水, 因此研究及预估参考作物蒸散量的时空变化对于深入了解流域水分循环规律和农田水分管理具有重要意义.

目前, 大多数学者就鄱阳湖流域的蒸散研究主要基于历史气候数据和遥感数据, 如涂安国等[5]利用M-K突变检验和归因分析研究了鄱阳湖流域近50年潜在蒸散量的变化特征及其影响因子.刘健等[6]运用互补相关蒸发模型和小波分析探究了流域实际蒸发量及变化规律.朱婧瑄等[7]利用MODIS数据分析了2000 -2013年鄱阳湖流域土地利用变化对蒸散量的影响.赵晓松等[8]基于MODIS数据, 应用地面温度-植被指数三角关系法反演2000 -2009年鄱阳湖湖区蒸散的时空分布及主要气候因子的影响.相关研究主要集中在对鄱阳湖流域蒸散量的时空变化规律和影响因素的分析, 对流域未来蒸散量的分析还尚未有研究提及.预估未来气候特征的研究主要是基于全球气候模式(Global Climate Models, GCMs)和降尺度模型[9-10], 应用降尺度技术可以将大尺度、低空间分辨率的全球气候模式数据转换为小尺度高空间分辨率的区域地面信息[11], 这种方法在国外的比利时[12]、英国[13]及我国的黑河流域[14]和泾河流域[15]等地区已取得相关研究成果.本文基于第五阶段全球气候耦合模式比较计划(CMIP5)模式并结合统计降尺度模型预估我国鄱阳湖流域未来参考作物蒸散量的时空变化特征, 以期深入了解流域水分循环与水量平衡规律, 为旱涝灾害的成因分析提供科学的依据, 填补鄱阳湖流域关于未来参考作物蒸散量研究的空白.

1 材料与方法 1.1 研究区概况

鄱阳湖流域位于江西省北部, 总面积16.22×104 km2, 由修水、饶河、信江、抚河、赣江五大子流域组成, 地形三面环山, 中部平原、丘陵与盆地交错分布, 属典型的亚热带湿润性季风气候, 降水丰富, 日照充足, 蒸发旺盛, 年平均气温为17.1℃, 年平均降水量为1632 mm, 日照百分率达43 % ~47 %, 年均蒸发量800~1100 mm[16-17].受五大支流和长江水系影响, 近年来, 鄱阳湖流域旱涝灾害频发, 原有的水量平衡被打破.

1.2 数据来源 1.2.1 气象站点实测数据

本文选取鄱阳湖流域共14个气象站点(图 1)1961 -2014年逐日监测数据, 气象指标分别为最高气温、最低气温、平均气温、平均风速、日照时数、平均相对湿度以及小型蒸发皿蒸发量等, 数据来源于中国气象数据共享服务中心.

图 1 鄱阳湖流域范围及站点分布 Fig.1 Poyang Lake Basin and sites distribution
1.2.2 NCEP再分析数据

NCEP/NCAR再分析数据是由美国国家环境预测中心(NCEP)和大气研究中心(NCAR)联合发布的全球26个预报因子1961 -2014年的日序列再分析资料, 网格分辨率是2.5°×2.5°.借鉴现有研究基础[18-19], 本文采用的NCEP再分析资料备选预报因子分别为地面温度(air)、地面气压(pres)、地面相对湿度(rhum)、地面风速(wspd)、地面经向风速(vwnd)、地面纬向风速(uwnd)、海平面压力(slp), 以及在500/850 hPa不同高度处的温度、气压、相对湿度、风速、经向风速、纬向风速以及垂直风速因子(omega)和位势高度(hgt).这些备选因子通过对地面和不同高度处的水热传输及扩散过程的影响, 进而影响整个流域的蒸散量.

1.2.3 CMIP5全球气候模式

综合考虑全球气候模式(GCMs)在我国各地区的适用性、模式输出数据的连续性、完整性以及模型的分辨率高低, 本文GCMs数据选取法国气象研究中心推出的CNRM-CM5模式[14, 18]所输出的1961 -2100年逐日序列, 空间网格分辨率为1.406°×1.401°, 预报因子与NCEP再分析资料一致. CMIP5气候模式包含4种“典型浓度路径”(Representative Concentration Pathways, RCPs), 按低至高排列分别为RCP2.6、RCP4.5、RCP6.0和RCP8.5, 每一种情景都指示温室气体、活性化学气体和气溶胶的不同排放浓度, 本文在此选取应用最为广泛的RCP4.5和RCP8.5浓度路径.

1.3 研究方法 1.3.1 Penman-Monteith公式

Penman-Monteith(P-M)公式1948年由英国科学家彭曼提出, 该公式以能量动态平衡、水汽扩散原理以及空气的热导定律为理论依据, 经联合国粮农组织(FAO)修正, 推荐P-M公式作为计算ET0的标准方法, 参考作物蒸散发(ET0)是指假设作物高度为0.12 m, 有固定为70 s/m的表面阻力, 反照率为0.23的参考冠层的蒸散发[20].其公式为:

$ \mathrm{ET}_{0}=\frac{0.408 \Delta\left(R_{\mathrm{n}}-G\right)+\gamma \frac{900}{T+273} u_{2}\left(e, -e_{\mathrm{a}}\right)}{\Delta+\gamma\left(1+0.34 u_{2}\right)} $ (1)

式中, ET0表示参考作物蒸散量(mm);Rn为太阳净辐射(MJ/(m2·d));T为平均气温(℃);G为土壤热通量(MJ/(m2·d));u2为2 m高度处平均风速;γ为干湿表常数((kPa/℃);Δ为饱和水汽压曲线斜率(kPa/℃);es为饱和水汽压(kPa);ea为实际水汽压(kPa);es-ea表示为饱和气压差. P-M公式作为计算ET0的有效方法, 但在不同地区计算效果差异较大, 本文基于鄱阳湖流域14个站点数据合并计算出的ET0年值序列, 将其与1961 -2014年实测蒸发皿(ETpan)数据进行对比(图 2), 二者间相关系数高达0.89, 所以P-M公式计算出的ET0在鄱阳湖流域具有较好的适用性.

图 2 P-M公式计算ET0值与ETpan观测值比较 Fig.2 Comparison between ET0 value calculated by P-M formula and ETpan observation value
1.3.2 逐步回归分析

本文引入逐步回归分析(stepwise regression)来筛选NCEP备选因子[21], 逐步回归的基本思想是:将变量逐个引入模型并对其逐个进行检验, 如果先引入的变量因后面变量的引入而不显著(P>0.05)时, 则将其剔除.这个过程反复进行, 直到既无显著的变量选入方程, 也无不显著自变量从回归方程中剔除为止.

1.3.3 Mann-Kendall检验法

M-K突变检验法是世界气象组织(WMO)推荐并已广泛使用的一种非参数统计方法[22], 能有效区分长时间数据序列是存在确定的变化趋势或处于自然波动状态, 对于非正态分布的长时序气象水文数据, M-K秩次相关检验具有更加突出的适用性.

1.3.4 空间自相关

根据“地理学第一定律”指出任何事物与其周围事物都存在相关性, 且距离越近的事物相关性越显著.空间自相关(spatial autocorrelation)就是指在一个给定的区域内, 各变量的观测值之间存在不同程度的潜在依赖性.本文选取经典的全局空间自相关指数(Moran's I)来定量地描述鄱阳湖流域ET0的空间依赖性, 计算公式如下[23]

$ {\rm{Moran}}'{\rm{s}}\;I = \frac{{n\sum\limits_{i = 1}^n {\sum\limits_{j = 1}^n {{W_{ij}}} } \left( {{x_i} - \bar x} \right)\left( {{x_j} - \bar x} \right)}}{{\left( {\sum\limits_{i = 1}^n {\sum\limits_{j = 1}^n {{W_{ij}}} } } \right)\sum\limits_{i = 1}^n {{{\left( {{x_i} - \bar x} \right)}^2}} }}(i \ne j) $ (2)

式中, n为空间位置总和;xixj分别表示在空间ij的观测数值;Wij为空间权重矩阵;为所有空间位置的观测均值. Moran's I的取值范围在-1~1之间, 负值表示呈空间负相关关系, 0表示不相关, 正值表示空间正相关关系.通常以标准统计量Z对研究区域的Moran's I进行显著性检验, Z值的计算公式如下:

$ Z=\frac{I-E(I)}{\sqrt{V A R(I)}} $ (3)

其中:

$ E(I)=-\frac{I}{n-1} $ (4)
$ V A R(I)=E(I 2)-E(I) 2 $ (5)

式中, E(I)为I的均值;VAR(I)为I的方差;当Z值为负时, 表示观测数值在空间上趋于分散;当Z值为0时, 观测值在空间上处于独立随机分布状态;当Z值为正时, 表示观测值在空间趋于集聚分布.当Z值满足Z < 1.96时, 即拒绝零假设, 即在95 %的概率下, 各空间单元之间存在空间自相关关系.

2 降尺度模型构建

统计降尺度模型(statistical downscaling model, SDSM)的本质是借助长时间观测数据, 建立大尺度气象环流因素与区域气象要素间的显著性关系.该模式是基于随机天气发生器和多元回归相耦合的原理构建而成, 其主要包括两个内容:(1)基于区域预报量(ET0)与全球大气环流NCEP预报因子之间的统计关系, 构建SDSM模型;(2)将构建好的SDSM模型结合全球气候模型(CNRM-CM5)输出的未来数据, 生成对区域未来ET0的预估.具体流程如图 3.

图 3 基于统计降尺度的CMIP5模式未来气候变化预估 Fig.3 Prediction of future climate change in CMIP5 based on statistical downscaling model
2.1 选择预报因子

参照先前学者选取预报因子的原则[24], 本文选择ET0作为预报量, 用逐步回归分析法初步筛选所有备选预报因子, 依据所选预报因子与ET0回归方程的标准化系数与共线性检验, 进一步剔除影响不显著和存在显著共线性的预报因子, 对各站点重复上述操作, 最终选取的统计降尺度预报因子如表 1.

表 1 统计降尺度预报因子* Tab. 1 Predictors of statistical downscaling models
2.2 模型的交叉验证

基于鄱阳湖流域14个站点1961 -2014年逐日实测数据, 采用P-M公式计算出流域ET0月均值, 建立预报量ET0与所选预报因子NCEP的统计降尺度模型.为了提高未来ET0的预估精度, 以更长的时间段建立降尺度模型, 因此本文的模型率定和验证不划分时间尺度, 即模型率定期同样也是模型的验证期, 以1961 -2014年作为模型的交叉验证期, 以均值相对误差、标准差相对误差和月值拟合度作为降尺度模拟效果的评价指标, 结果见表 2.

表 2 统计降尺度模型交叉验证 Tab. 2 Calibration and verification of statistical downscaling models

表 2所示, 各气象站点的ET0公式计算值与模型模拟值十分接近, 均值相对误差在-0.112 % ~0.000 %之间, 存在略微低估现象.从标准差相对误差和月值拟合度看, 其分别在-5.609 % ~-1.432 %和0.812~0.956之间, 各站点的评价指标值相差不大, 以上表明降尺度模型对鄱阳湖流域的ET0值模拟效果较好.

2.3 CMIP5模式适用性评估

全球气候模式(GCMs)因初始边界条件不同, 情景、模型观测参数、内部结构等设置的差异, 在不同地区对气候因子的模拟效果差异显著, 因此, 本文在将CMIP5模拟数据输入降尺度模型预测流域未来ET0之前, 首先对其进行区域适用性评估. CMIP5中的CNRM-CM5模式的预测起始时间是2006年, 分辨率是1.406°×1.401°, 本文借助反距离权重插值法(IDW)将其调整为与NCEP数据相同的网格分辨率, 采用CNRM-CM5模式2006 -2014年的预测因子, 利用上文建立的统计降尺度模型, 模拟CNRM-CM5模式在RCP4.5和RCP8.5情景下各站点的ET0值(表 3), 在不同情景下ET0均值、相对误差和标准差相对误差较大, 可见直接使用CNRM-CM5数据模拟鄱阳湖流域ET0值的效果较差, 因此本文采用偏差校正法, 通过二次多项式对流域各站点的模式输出值进行校正, 结果如表 3所示.由表 3可以看出, 基于CNRM-CM5模式直接经降尺度模拟所得ET0值与P-M公式计算值差值较大, 经过偏差校正后的差值明显减小.在RCP4.5情景下, 均值相对误差范围为-17.45 % ~1.33 %、标准差相对误差范围为-39.43 % ~28.33 %、月值拟合度为0.623~0.965, 经偏差校正后的结果分别为-2.104 % ~1.022 %、-7.672 % ~-2.064 %和0.848~0.972.在RCP8.5情景下, 均值相对误差范围为-17.79 % ~1.06 %、标准差相对误差范围为-39.85 % ~31.87 %、月值拟合度为0.638~0.933, 经偏差校正后的结果分别为-1.038 % ~1.031 %、-8.422 % ~-2.093 %和0.705~0.934.综上, 通过对模式输出结果进行偏差校正后, 评价指标均值相对误差、标准差相对误差显著减小, 月值拟合度明显升高, 模型的模拟效果得到很大改善, 即在偏差校正基础上的CNRM-CM5模式的降尺度模拟值适宜于鄱阳湖流域未来参考作物蒸散量的预估.

表 3 CMIP5模式适用性评估 Tab. 3 Evaluation of CMIP5 model applicability
3 CMIP5未来参考作物蒸散量(ET0)预估 3.1 参考作物蒸散量的时间变化特征 3.1.1 参考作物蒸散量年内变化特征

依据偏差校正后的统计降尺度模型, 对RCP4.5和RCP8.5两种情景下鄱阳湖流域未来ET0的变化进行预估, 将1961 -2010年作为基准期, 未来按30年时段划分为未来初期(2011 -2040年)、未来中期(2041 -2070年)、未来远期(2071 -2100年). 图 4a图 4b分别为在RCP4.5和RCP8.5情景下鄱阳湖流域ET0的年内变化趋势, 图 4a表明在RCP4.5情景下ET0的年内变化趋势较为一致, 均呈倒“U”字型, 夏季(6、7、8月)为高值区, 7月份时, 基准期和未来3个时期ET0均达到最大值, 冬季(12、1、2月)为低值区, 1月份为年内最低值. RCP4.5情景下未来3个不同时期ET0月均值表现为未来远期>未来中期>未来初期>基准期, 即在RCP4.5情景下未来ET0月均值呈不断增大的趋势. 图 4b为在RCP8.5情景下不同时期ET0的年内变化, 总体趋势与图 4a相似, 呈倒“U”字型分布, 与之不同的是, ET0年内高值时间有后移的趋势, 基准期和未来初期在7月达到最高值, 分别为144.6和151.3 mm, 未来中期和未来远期则在8月时达到最高值, 分别为152.5和156.1 mm.在RCP8.5情景下未来ET0月均值同样呈不断增大的趋势.

图 4 鄱阳湖流域RCP4.5(a)、RCP8.5(b)情景下1961 -2100年参考作物蒸散量年内变化趋势 Fig.4 Change trend of intra-annual reference crop evapotranspiration in Poyang Lake Basin during 1961 to 2100 under the scenarios of RCP4.5(a) and RCP8.5(b)
3.1.2 参考作物蒸散量年际变化特征

以1961-2100年鄱阳湖流域14个站点ET0年均值构成的时间序列进行M-K突变检验及年际变化分析, 结果图 5所示.

图 5 鄱阳湖流域RCP4.5、RCP8.5情景下1961 -2100年ET0 M-K检验和年际变化特征 (UF、UB分别表示顺序变化曲线和逆序变化曲线;CL表示置信水平) Fig.5 Change trend of inter-annual ET0 and M-K test in Poyang Lake Basin during 1961 to 2100 under the scenarios of RCP4.5 and RCP8.5

图 5a~5b可以看出, 鄱阳湖流域在基准期(1961 -2010年), UF曲线呈明显的下降态势, 即ET0整体减小, 而在21世纪以后ET0则是以波动上升趋势为主.通过改变时间序列长度的基础上反复检验结果表明, RCP4.5情景下, 2062年为ET0的突变年, 2083年以后ET0增大的趋势达到显著性水平(P < 0.05);RCP8.5情景下, 2058年为ET0的突变年, 2079年以后ET0增大的趋势达到显著性水平(P < 0.05), 相较于RCP4.5情景, 突变年年份有所提前、上升趋势更为显著. 图 5c~5d为鄱阳湖流域RCP4.5、RCP8.5情景下1961 -2100年ET0年际变化特征, 可以看出, RCP4.5情景下未来3个时期ET0倾向率分别为9.2、14.5和5.6 mm/10 a, 即在RCP4.5情景下未来中期ET0增幅最大、未来初期次之、未来远期相对最小;RCP8.5情景下3个时期ET0的倾向率分别为14.5、16.4和12.9 mm/10 a, 与RCP4.5情景相似, 在RCP8.5情景下ET0增幅同样表现为未来中期最大、未来初期次之、未来远期相对最小.

以泰森多边形法(TPM)计算各气象台站面积权重, 得出全流域未来参考作物蒸散量, 如表 4所示为鄱阳湖流域RCP4.5、RCP8.5情景下不同时间段内ET0的平均值, 可以看出, 在RCP4.5情景下未来3个时期ET0平均值较基准期1961 -2010年的平均值分别增加了3.8 %、5.3 %和6.3 %, 在鄱阳湖流域一共14个站点中, 未来3个时间段ET0均保持持续上升的站点有9个, 占比64.3 %, 即在RCP4.5情景下, 鄱阳湖流域未来ET0整体上保持上升趋势, 但增幅将逐渐减小. RCP8.5情景下3个时间段ET0平均值相较基准期分别增加了4.2 %、6.9 %和8.6 %, 3个时间段ET0均保持持续上升的站点有11个, 占比78.6 %, 这表明鄱阳湖流域在RCP8.5情景下ET0同样保持明显的上升趋势, 且增幅要大于情景RCP4.5.

表 4 鄱阳湖流域为未来不同时段平均ET0(mm) Tab. 4 Mean value of ET0 under different scenarios in Poyang Lake Basin
3.2 参考作物蒸散量空间分布特征 3.2.1 基准期参考作物蒸散量的空间分布

通过对由统计降尺度模型模拟出CNRM-CM5模式的未来参考作物蒸散数据, 以及P-M公式估算出鄱阳湖流域14个站点历史参考作物蒸散数据, 采用普通克里金插值对其进行空间化处理, 得到鄱阳湖流域在基准期和未来不同情景下ET0的空间分布情况.

鄱阳湖流域基准期(1961 -2010年)ET0的空间分布可以看出, 鄱阳湖流域ET0的空间分布总体上表现出南北高、中间低的特点, 南部是以赣州、遂川为中心的高值中心, ET0的变化范围在1111~1160 mm之间, 东北部则是以波阳和南昌为中心的次高值区, ET0的变化范围在1093~1130 mm之间, 而中西部是以修水和宜春为典型代表的低值区, ET0的变化范围在985~1080 mm之间(图 6).总的来说, 鄱阳湖流域在1961 -2010年间, ET0值空间分布上表现出明显的空间差异性.

图 6 鄱阳湖流域1961 -2010年参考作物蒸散量的空间分布 Fig.6 Spatial distribution of reference crop evapotranspiration during 1961 to 2010 in Poyang Lake Basin
3.2.2 未来情景下参考作物蒸散量的空间分布

鄱阳湖流域在RCP4.5情景下ET0的空间分布可以看出(图 7a~7c), 鄱阳湖流域ET0在未来的3个不同时期变化特征非常相似, 从未来初期到未来远期, 南部以赣州、遂川为高值中心的赣江流域逐渐与东北部以波阳、南昌、贵溪和南城为中心的次高值区由原先相互独立的块状到最后逐渐连成一个整体, 高值区范围明显扩大, 而以修水、宜春和樟树三地为中心的低值区范围则逐渐缩小, 其中广昌的ET0增幅最为突出达63.2 mm, 修水的ET0增幅最小仅为12.5 mm.总的来说, 随着时间的推移, 鄱阳湖流域ET0在RCP4.5情景下呈不断增加的态势, 且以未来中期增幅最为显著, 表现出东高西低的空间分布特征.

图 7 RCP4.5、RCP8.5情景下鄱阳湖流域未来ET0的空间分布 Fig.7 Spatial distribution of future ET0 in Poyang Lake Basin under the scenarios of RCP4.5 and RCP8.5

鄱阳湖流域在RCP8.5情景下ET0的空间分布可以看出(图 7d~7f), 鄱阳湖流域在未来初期和未来中期ET0的空间分布特征较为相似, 以修水、宜春、吉安和广昌为中心的低值区与中值区相连, 横切鄱阳湖流域南北部的两个高值中心, 直至未来远期, 以赣州为中心的高值区向北扩张, 广昌和吉安两地ET0值明显增加, 使得南北的高值中心在鄱阳湖流域东部地区连成一个整体, 流域ET0的最大值赣州由未来初期的1198.8 mm到未来远期的1221.9 mm, 相较增加了23.1 mm, 流域ET0的最小值修水由未来初期的1044.7 mm到未来远期的1068.7 mm, 相较增加了24.0 mm, 可见ET0的最高值、最低值的增幅并不显著, 广昌依然是增幅最大的地区(高达93.1 mm).总的来说, 鄱阳湖流域在RCP8.5情景下未来ET0的空间分布整体上呈东高西低, 局地略有突出分布的特征, 未来中期是ET0增幅最大的一个时期, 且流域ET0在RCP8.5情景下增幅要明显高于RCP4.5情景.

为了科学直观地探究鄱阳湖流域未来时期ET0在不同情景下的空间集聚性, 利用GeoDA软件构建基于距离关系的空间权重矩阵, 并分别计算出基准期和RCP4.5、RCP8.5情景下未来不同时期ET0的Moran's I值, 结果如表 5.可以看出, 基准期和未来不同情景下的标准统计量Z(I)值均大于1.96, 且Moran's I值均大于0, 也就意味着鄱阳湖流域无论是在基准期或是未来情景下的不同时期, 其ET0值均具有较强的空间自相关性, 存在一定的空间集聚现象, 根据Moran's I的变化可以看出, 在RCP4.5情景下ET0的空间自相关的强度表现为未来初期 < 未来中期 < 未来远期, 在RCP8.5情景下ET0的空间自相关的强度表现为未来中期 < 未来初期 < 未来远期, 在RCP4.5、RCP8.5情景下的未来远期时期, Moran's I值最高分别达0.309和0.336, 这表明RCP4.5、RCP8.5情景下的未来远期不仅仅是鄱阳湖流域ET0值最高的时期, 同时也是ET0空间相关性、空间集聚最为显著的时期.总的来说, 鄱阳湖流域的ET0的空间分布存在较强的自相关性, 空间集聚强度RCP8.5情景高于RCP4.5情景.

表 5 鄱阳湖流域RCP4.5、RCP8.5情景下ET0的Moran's I Tab. 5 Moran's I value of ET0 in Poyang Lake Basin under the scenarios of RCP4.5 and RCP8.5
3.3 未来情景下干旱指数的时空分布

本文以参考作物蒸散量为基础, 结合统计降尺度模型所输出的未来降水数据, 基于budyko的水热平衡耦合理论[25], 计算鄱阳湖流域在RCP8.5情景下未来干旱指数(Aridity index, AI=ET0/P)时空分布特征, AI对流域干湿状况的划分界限为(干旱:5≤AI < 12;半干旱:2≤AI < 5;半湿润0.75≤AI < 2;湿润:0.35≤AI < 0.75)[26].在RCP8.5情景下, 鄱阳湖流域在1961 -2100年干旱指数呈现出较为明显的上升趋势, 在2045年以前, 流域平均干旱指数低于0.75界限值, 以湿润气候为主要特征, 2045年以后, 干旱指数逐渐上升高于0.75, 表明在21世纪中期, 鄱阳湖流域气候特征将从湿润向半湿润过渡, 干旱风险有所提升(图 8).

图 8 鄱阳湖流域1961 -2100年干旱指数年际变化 Fig.8 Change trend of inter-annualaridity index in Poyang Lake Basin during 1961 to 2100

鄱阳湖流域1961 -2010年以及在RCP8.5情景下2011 -2100年干旱指数的空间分布可以看出(图 9a), 鄱阳湖流域在基准期(1961 -2010年), 除南部的遂川、赣州两地, 流域其余地区AI值均低于0.75, 属典型的湿润区;然而在RCP8.5情景下, 鄱阳湖流域2011 -2100年间其AI值>0.75的范围明显扩大(图 9b), 流域绝大部分地区由湿润区转为半湿润区, 空间分布表现为以南部赣州和吉安为AI高值中心, 向北逐渐减小, 可见, 在未来时期, 鄱阳湖流域的干旱状况随时间加剧, 赣江流域是未来干旱风险的重点防范区.

图 9 鄱阳湖流域1961 -2010年(a)和RCP8.5情景下2011 -2100年(b)干旱指数空间分布 Fig.9 Spatial distribution of aridity index in Poyang Lake Basin from 1961 to 2010(a) and RCP8.5 from 2011 to 2100(b)
4 讨论

全球气候变暖的大背景下, 鄱阳湖流域在1961 -2010年(基准期)这50 a来, 参考作物蒸散量整体上呈现出减小态势, 气温的升高并未导致流域ET0的升高, 与不少学者基于历史时期长时间序列的研究结论相一致, 参考作物蒸散量的减少表明在鄱阳湖流域同样存在“蒸发悖论”现象[14, 27].目前一般认为导致这种现象的原因主要包括太阳辐射、气温日较差、相对湿度和风速等, 此外, 随着社会经济的发展, 人类活动和下垫面状况的改变对ET0的间接影响越来越明显.本研究发现在基准期的后期2000 -2010年, 鄱阳湖流域ET0表现出上升的趋势, 这说明“蒸发悖论”现象存在时间上的不一致性, 由于ET0的影响因子众多, 导致区域ET0存在不同程度上的时空差异, 这对未来时期ET0的预测提出巨大挑战.

本研究预测在未来RCP4.5、RCP8.5情景下预测鄱阳湖流域ET0整体上呈增加趋势, 且后者大于前者, 没有出现类似于历史时期的“蒸发悖论”现象, 与学者对我国其他地区未来ET0预测结果相符[14-15], 研究结果表明, 气温、辐射水平与ET0呈正比. RCP4.5、RCP8.5情景下未来3个不同时期鄱阳湖流域ET0空间分布整体上呈东高西低、局地略有突出分布的特征, 较基准期ET0南北高、中间低的分布特征有明显的改变, 原因可能与全球气候变化和人类活动和下垫面状况的改变密切相关, 在未来的研究中将结合多重数据做进一步分析. RCP4.5作为中间稳定排放情景, 在RCP4.5气候情景下, 气温相较工业革命升高约4.8℃, 预计在2100年以后辐射强迫水平达4.5 W/m2, 更为接近实际发展水平, 具有很大的参考价值. RCP8.5作为按照当前水平不采取减排措施的高浓度排放路径, RCP8.5气候情景下, 气温相较工业革命升高约6.9℃, 预计在2100年以后辐射强迫水平达8.5 W/m2, 则具有更强的警示意义.

此外, 刘昌明等[28]指出, 现今的统计降尺度模型在实际研究应用中存在较大的不确定性, 其中GCMs是最重要的不确定因子, 预报因子的选取和统计降尺度方法的不确定性也会对最终的模拟结果产生一定误差. Baguis P等[12]研究结果表明, 由于气候情景和模式的不确定性, 对蒸散量的影响要比计算公式中的参数校正带来的影响大得多, 从而使得ET0的估算产生不同程度上的差异.因此, 在以后的研究中要加强对特定区域气候特征与大尺度环流之间的数理分析, 从原理上论证二者间的统计关系在时空上的稳态性, 通过对不同尺度的环流进行多尺度多时空的转换, 结合多情景、多种统计降尺度方法和多模式集合的分析方法, 定量地描述不确定性的概率分布, 以更科学地预测和评估未来时期的某一地区的蒸散量.

5 结论

本文基于鄱阳湖流域14个气象站点1961 -2014年的历史数据, 利用P-M公式、大气环流模式CNRM-CM5和统计降尺度模型SDSM, 分析了在RCP4.5、RCP8.5情景下未来3个时期鄱阳湖流域ET0的时空分布特征, 并估算鄱阳湖流域在RCP8.5情景下未来干旱指数的分布特征.主要结论如下:

1) 基于1961 -2014年NECP再分析资料和鄱阳湖流域ET0的降尺度模型交叉验证结果表明模型模拟效果较好;基于2006 -2014年CMIP5的CNRM-CM5模式的降尺度模拟结果经偏差校正后拟合精度明显提高, 校正后模型输出值适宜流域未来ET0的预估.

2) 鄱阳湖流域在1961 -2010年(基准期)ET0整体上呈减小趋势, 倾向率达46.1 mm/10 a, 年均值为1086.38 mm;ET0在空间分布上呈南北高、中间低的特点, 表现出明显的空间差异性.

3) RCP4.5、RCP8.5情景下鄱阳湖流域ET0均呈不同程度的增长趋势, 突变时间分别是2062年和2058年;在RCP4.5和RCP8.5情景下, 鄱阳湖流域未来3个时期年均ET0倾向率分别相较于基准期ET0值均有所提高;ET0空间分布整体上呈东高西低、局地略有突出的分布特征, 未来中期是ET0增幅最大的一个时期;鄱阳湖流域无论是在基准期或是未来情景下的3个时期, ET0均具有较强的空间自相关性, 未来远期是鄱阳湖流域ET0值最高的时期, 也是ET0空间相关性、空间集聚最为显著的时期.

4) 在RCP8.5情景下, 鄱阳湖流域在1961-2100年干旱指数呈现出较为明显的上升趋势, 且在21世纪中期, 鄱阳湖流域气候特征将从湿润向半湿润过渡;在1961-2010年间, 鄱阳湖流域绝大部分属典型的湿润区, 在RCP8.5情景下的2011-2100年, 流域的干旱状况随时间加剧, 流域绝大部分地区由湿润区转为半湿润区, AI值空间分布自南向北递减, 赣江流域将是鄱阳湖流域未来干旱风险的重点防范区.

6 参考文献

[1]
Cai YL. Sensitivity and adaptation of Chinese agriculture under global climate change. Acta Geographica Sinica, 1996(3): 202-212. [蔡运龙. 全球气候变化下中国农业的脆弱性与适应对策. 地理学报, 1996(3): 202-212. DOI:10.3321/j.issn:0375-5444.1996.03.002]
[2]
Coquard J, Duffy PB, Taylor KE et al. Present and future surface climate in the western USA as simulated by 15 global climate models. Climate Dynamics, 2004, 23(5): 455-472. DOI:10.1007/s00382-004-0437-6
[3]
IPCC. Climate Change 2013: The physical science basis. In: Stocker TF, Qin D, Plattner GK et al eds. Final draft underlying scientific-technical assessment. Cambridge, UK and New York, NY, USA: Cambridge University Press, 2012: 5.
[4]
Xin WQ, Wang WG, Shao QX et al. Projection of future reference evapotranspiration change across the haihe river basin. Journal of Basic Science and Engineering, 2014(2): 239-251. [邢万秋, 王卫光, 邵全喜等. 未来气候情景下海河流域参考蒸发蒸腾量预估. 应用基础与工程科学学报, 2014(2): 239-251.]
[5]
Tu AG, Li Y, Nie XF et al. Spatial and temporal variations of reference crop evapotranspiration and its influencing factors in Poyang lake watershed. Ecology and Environmental Sciences, 2017, 26(2): 211-218. [涂安国, 李英, 聂小飞等. 鄱阳湖流域参考作物蒸散量变化特征及其归因分析. 生态环境学报, 2017, 26(2): 211-218.]
[6]
Liu J, Zhang Q, Xu CY et al. Change of actual evapotranspiration of Poyang lake watershed and associated influencing factors in the past 50 years. Resources and Environment in the Yangtze Basin, 2010, 19(2): 139-145. [刘健, 张奇, 许崇育等. 近50年鄱阳湖流域实际蒸发量的变化及影响因素. 长江流域资源与环境, 2010, 19(2): 139-145.]
[7]
Zhu JX, Qi SH, Liu GH et al. Spatiotemporal variation of terrestrial evapotranspiration in Poyang lake watershed from 2001 to 2013. Resources and Environment in the Yangtze Basin, 2016, 25(S1): 9-16. [朱婧瑄, 齐述华, 刘贵花等. 2000-2013年鄱阳湖流域蒸散量时空变化. 长江流域资源与环境, 2016, 25(S1): 9-16.]
[8]
Zhao XS, Liu YB, Wu GP et al. A remote-sensing-based study on evapotranspiration and the environmental factors over the Lake Poyang region. J Lake Sci, 2013, 25(3): 428-436. [赵晓松, 刘元波, 吴桂平. 基于遥感的鄱阳湖湖区蒸散特征及环境要素影响. 湖泊科学, 2013, 25(3): 428-436. DOI:10.18307/2013.0318]
[9]
Hellstrm C, Chen D. Statistical downscaling based on dynamically downscaled predictors:Application to monthly precipitation in Sweden. Advances in Atmospheric Sciences, 2003, 20(6): 951-958. DOI:10.1007/BF02915518
[10]
Murphy J. Predictions of climate change over Europe using statistical and dynamical downscaling techniques. International Journal of Climatology, 2000, 20(5): 489-501. DOI:10.1002/(SICI)1097-0088(200004)20:5<489::AID-JOC484>3.0.CO;2-6
[11]
Liu YH, Guo WD, Feng JM et al. A summary of methods for statistical downscaling of meteorological data. Advances in Earth Science, 2011, 26(8): 837-847. [刘永和, 郭维栋, 冯锦明等. 气象资料的统计降尺度方法综述. 地球科学进展, 2011, 26(8): 837-847.]
[12]
Baguis P, Roulin E, Willems P et al. Climate change scenarios for precipitation and potential evapotranspiration over central Belgium. Theoretical and Applied Climatology, 2010, 99: 273-286. DOI:10.1007/s00704-009-0146-5
[13]
Chun KP, Wheater HS, Onof C. Projecting and hindcasting potential evaporation for the UK between 1950 and 2099. Climate Change, 2012, 113: 639-661. DOI:10.1007/s10584-011-0375-3
[14]
Qi XF, Li WP, Li HT et al. Prediction and trend of future potential evapotranspiration in the Heihe river basin based on CMIP5 models. Journal of Desert Research, 2018, 38(4): 849-857. [祁晓凡, 李文鹏, 李海涛. 基于CMIP5模式的黑河流域潜在蒸散量预估. 中国沙漠, 2018, 38(4): 849-857.]
[15]
Zhao C, Li Z, Liu WZ et al. Spatiotemporal characteristics of potential evapotranspiration during current and future period in the Jinghe watershed. Acta Ecologica Sinica, 2014, 34(19): 5600-5608. [赵姹, 李志, 刘文兆等. 泾河流域潜在蒸散量的时空变异. 生态学报, 2014, 34(19): 5600-5608.]
[16]
Liu C, Qi SH, Tang LL et al. Spatiotemporal pattern of evapotranspiration in Poyang Lake Basin under the joint influences of vegetation restoration and climate change. Geographical Research, 2016, 35(12): 2373-2383. [刘冲, 齐述华, 汤林玲等. 植被恢复与气候变化影响下的鄱阳湖流域蒸散时空特征. 地理研究, 2016, 35(12): 2373-2383.]
[17]
Liu ZH, Lu JZ, Chen XL et al. Attribution of runoff variation in the Fu river basin of Poyang lake based on Budyko hypothesisbased. Henan Science, 2019(8): 1303-1310. [刘子豪, 陆建忠, 陈晓玲等. 基于Budyko假设的鄱阳湖抚河流域径流变化归因分析. 河南科学, 2019(8): 1303-1310. DOI:10.3969/j.issn.1004-3918.2019.08.017]
[18]
Wu HM, Huang AN, He Qing et al. Projection of the spatial and temporal variation characteristics of precipitation over Central Asia of 10 CMIP5 models in the next 50 years. Arid Land Geography, 2013, 36(4): 669-679.
[19]
Moss R, Babiker M, Brinkman S et al. Towards new scenarios for analysis of emissions, climate change, impacts, and response strategies. Geneva, Switzerland: IPCC, 2009: 25.
[20]
Allen RG, Pereira LS, Raes D et al. Crop evapotranspiration-guidelines for computing crop water requirement. Rome, Italy: United Nations Food and Agriculture Organization, 1998.
[21]
Liu SY, Huang Q, Wang YM et al. Simulation and prediction of climatic elements in Jinghe River basin using downscaling method and CMIP5. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(23): 138-144. [刘赛艳, 黄强, 王义民等. 基于统计降尺度和CMIP5模式的泾河流域气候要素模拟与预估. 农业工程学报, 2015, 31(23): 138-144. DOI:10.11975/j.issn.1002-6819.2015.23.018]
[22]
World Meteorological Organization. Analyzing long time series of hydrological data with respect to climate variability. Geneva, Switzerland: Tech. Rep.Wcap-3, 1998.
[23]
Anselin L, Syabri I, Kho Y. Geo Da:An introduction to spatial data analysis. Geographical Analysis, 2005, 38(1): 5-22.
[24]
Wilby RL, Hay LE, Leavesley GH. A comparison of downscaled and raw GCM output:implications for climate change scenarios in the San Juan River basin, Colorado. Journal of Hydrology, 1999, 225(1/2): 67-91.
[25]
Budyko MI. Climate and life. 1974.
[26]
Arora VK. The use of the aridity index to assess climate change effect on annual runoff. Journal of Hydrology, 2002, 265(1/2/3/4): 164-177.
[27]
Ma XN, Zhang MJ, Wang SJ et al. Evaporation paradox in the Yellow river basin. Acta Geographica Sinica, 2012, 67(5): 645-656. [马雪宁, 张明军, 王圣杰等. "蒸发悖论"在黄河流域的探讨. 地理学报, 2012, 67(5): 645-656.]
[28]
Liu CM, Liu WB, Fu GB et al. A discussion of some aspects of statistical downscaling in climate impacts assessment. Advances in Water Science, 2012, 23(3): 427-437. [刘昌明, 刘文彬, 傅国斌等. 气候影响评价中统计降尺度若干问题的探讨. 水科学进展, 2012, 23(3): 427-437.]