气象学报  2015, Vol. 73 Issue (5): 895-909   PDF    
http://dx.doi.org/10.11676/qxxb2015.063
中国气象学会主办。
0

文章信息

朱连华, 江志红, 刘征宇. 2015.
ZHU Lianhua, JIANG Zhihong, LIU Zhengyu. 2015.
春季陆面植被对长江流域夏季降水可预报性的影响分析
The impact analysis of spring vegetation on the summer precipitation predictability over the Yangtze River Basin
气象学报, 73(5): 895-909
Acta Meteorologica Sinica, 73(5): 895-909.
http://dx.doi.org/10.11676/qxxb2015.063

文章历史

收稿日期: 2014-07-31
改回日期: 2015-04-25
春季陆面植被对长江流域夏季降水可预报性的影响分析
朱连华1,2, 江志红1 , 刘征宇3    
1. 气象灾害教育部重点实验室, 气象灾害预报预警与评估协同创新中心, 南京信息工程大学, 南京, 210044;
2. 南京信息工程大学数学与统计学院统计系, 南京, 210044;
3. 大气海洋科学系气候研究中心, 美国威斯康星大学麦迪逊分校, 麦迪逊, 53706
摘要: 利用美国全球监测与模型研究中心(GIMMS)1982—2006年逐月归一化植被指数(NDVI)、美国国家海洋和大气局(NOAA)1854—2008年海温资料以及中国国家气候中心(NCC)1951—2006年160站月降水资料,通过旋转经验正交函数分解(REOF)和相关分析获得了长江流域夏季降水预报序列和植被、海温预报因子集。基于最优子集回归方法(OSR),并借助交叉验证(CV)以及空间重建等手段,构建了单独以前期春季海温为预报因子和同时引入前期春季海温与归一化植被指数为因子的两类预报模型,对比分析引入陆面植被因子前后长江流域夏季降水预报效果改善状况,评估春季陆面植被对长江流域夏季降水可预报性的影响及预报效果的稳健性。结果表明:(1)相对于海温因子,春季陆面植被因子对长江流域夏季降水预报具有同样重要性,引入春季归一化植被指数后,长江流域夏季降水预报得到明显改善,相关系数平均由0.49提升到0.66,提高0.17左右,模型解释方差提升平均60%左右,其中单纯海温因子预报效果较差的汉江—淮河地区和淮河流域地区,相关系数更是提高了0.20—0.30,模型解释方差提升1倍左右;(2)交叉验证预报表明,相对于仅考虑海温因子模拟情形,交叉预报相关系数下降较多,模型稳健性较低,引入归一化植被指数后,长江流域夏季降水预报稳健性得到明显提升,长江中下游及其以南的长江三角洲地区、洞庭湖—鄱阳湖地区改善尤为明显;(3)长江流域降水可预报性存在明显的区域差异,嘉陵江流域地区、汉江—洞庭湖地区预报效果最好,汉江—淮河地区、淮河流域地区、长江三角洲地区预报效果最差,但引入归一化植被指数后预报效果提高最明显,而洞庭湖—鄱阳湖地区虽然模拟效果较好,但预报稳健性较低,交叉验证相关系数降幅达到0.27,这也从侧面说明了长江流域夏季降水分区预报的重要性。
关键词: 归一化植被指数     长江流域     夏季降水     最优子集回归     交叉验证    
The impact analysis of spring vegetation on the summer precipitation predictability over the Yangtze River Basin
ZHU Lianhua1,2, JIANG Zhihong1 , LIU Zhengyu3    
1. Key Laboratory of Meteorological Disaster Ministry of Education, Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science and Technology, Nanjing 210044, China;
2. School of Mathematics and Statistics, Nanjing University of Information Science and Technology, Nanjing 210044, China;
3. Center for Climatic Research, Department Atmosphere and Oceanic Science, University of Wisconsin-Madison, Madison 53706, USA
Abstract: The Rotated Empirical Orthogonal Functions (REOF) and correlation analysis are used to get the summer precipitation forecasting objects, the vegetation and ocean factors over the Yangtze River Basin, based on the 1982-2006 monthly Normalized Difference Vegetation Index (NDVI) from the Global Inventory Monitoring and Modeling Studies (GIMMS) in USA, the 1854-2008 sea surface temperature (SST) data supported by the American National Oceanic and Atmospheric Administration (NOAA) and the 1951-2006 monthly precipitation dataset of 160 stations from the National Climate Center (NCC) of China. The optimal subset regression models (OSR), the cross validation (CV) tests and space reconstruction methods, are respectively introduced to analyze the improvement owing to the incorporation of the vegetation factors, the spring vegetation impacts on the summer precipitation predictability and robustness over the Yangtze River Basin, under the two different situations of ocean factors alone and both ocean and land vegetation included. The results show that: (1) The spring land vegetation is at least as important as the ocean temperature indices. Compared with the pure pre-spring SST forecast models, the predictability of the Yangtze River Basin is obviously improved after the introduction of the pre-spring vegetation factors. The average forecast correlations are increased by about 0.17 from 0.49 to 0.66 (the explained variance is increased by about 60%), especially for the poor prediction areas using SST factors alone such as the Han River-Huai River subarea and the Huai River Basin subarea in which the correlations are raised by about 0.20-0.30 (the explained variance is increased by about 100%). (2) The cross validation results indicate that the pure SST forecast models of the Yangtze River Basin have the low prediction robustness in which the cross validation forecast correlations have the great drop. The introduction of the vegetation NDVI factors can receive the better performance in which the forecast correlations are significantly improved over the middle-lower reaches of the Yangtze River Basin and its south such as the Yangtze River Delta subarea and the Dongting Lake-Poyang Lake subarea. And, (3) there are the obvious predictability difference among the subareas over the Yangtze River Basin, which show the necessity of the subareas forecasts one by one. The best performance areas are the Jialing River Basin subarea and the Han River-Dongting Lake subarea, and the worst ones are the Han River-Huai River subarea, the Huai River Basin subarea and the Yangtze River Delta subarea. In spite of the well simulation results, the Dongting Lake-Poyang Lake subarea shows the weaker robustness for the cross validation tests in which the forecast correlation coefficients are reduced by 0.27.
Key words: NDVI;Yangtze River Basin;Summer precipitation;Optimal subset regression;Cross validation    
1 引言

长江流域夏季降水是东亚夏季风降水的主要部分,也是中国夏季降水型态预报的关键区域,该区域的夏季降水预报一直是中国汛期降水季节预测的重点与难点(钱维宏等,2010)。目前,中国夏季汛期降水统计预报主要考虑海温(SST)、大气环流场(如副热带高压)、积雪等因子(Ding et al,2009Wu et al,2009魏凤英等,2010冯蕾等,2011郭玲等,2012霍飞等,2014祁莉等,2014周放等,2014),但预测结果表明仅考虑海温和副热带高压(副高)等影响,该区夏季降水的月预报技巧一直相对较低。

近年来,越来越多理论分析和数值模拟表明,陆-气尤其是大气和植被存在复杂的相互作用(Liu et al,2006; 崔林丽等,2011Notaro et al,20112012; Barriopedro et al,2012; Peng et al,2012; 张少波等,2013Wang et al,2014)。一方面陆面植被(主要以归一化植被指数(NDVI)表征)对气候因子存在响应,植被生长发育状况及分布格局会受到如降水、温度等气候条件影响;另一方面陆面植被对气候因子也存在反馈,植被通过与大气进行水分以及能量、动量交换等方式影响降水、温度等气候要素,而植被自身延续性决定植被会月度、季节乃至更长时间尺度局地、非局地区域范围影响着大气,这为利用前期植被因子进行短期气候预报提供了可行性(Lee et al,2008)。

不少观测研究表明,陆面植被对中国区域的降水存在直接的影响。如在局地反馈分析方面,张井勇等(2003)研究指出,中国多数地区前期春季归一化植被指数增大,当年夏季降水也相应增多,而青藏高原、东部干旱-半干旱区和华中更是植被影响降水的3个高相关敏感区;Li等(2005)通过模式对比研究显示引入植被覆盖度(FVC)和叶面积指数(LAI)后有利于模式更好地刻画降水的年际变率,并指出植被覆盖度、叶面积指数增大,地面反照率减小,地面净辐射增大,从而导致局地蒸发和降水增强。而在非局地反馈分析方面,Xue(1996)符淙斌等(2001)Zhang等(2005)研究显示内蒙古草原沙漠化可能会导致季风环流减弱、局地降水减少,从而导致中国北部、南部降水减少,而中部区域降水增多;丁一汇等(2005)进一步指出中国西北地区绿化则有利于黄河流域降水增多,而长江流域和江南地区降水却不同程度减少;郑益群等(2002a2002b)研究指出长江—淮河流域降水频率增大、华北洪涝强度增大,可能是北方草原和南方常绿阔叶林改变的综合影响;范广洲等(2009)Wang等(2010)Zuo等(2011)指出青藏高原植被增加,中国华南、华北汛期降水增多,长江流域及其以南降水减少;周锁铨等(1995)张少波等(2013)认为青藏高原无植被后大气湿度和地面潜热输送均减小,削弱了扰动发展和高原上东西方向波的传播,很可能影响青藏高原低涡或西南涡的东移,进而影响中国夏季降水。

可见,围绕陆面植被与中国夏季降水关系的分析已开展了不少研究工作,然而需要指出的是,上述研究多集中在诊断分析方面,在短期气候预测方面直接将陆面植被作为预报因子引入预报模型中尚不多见。本研究拟选取长江流域夏季降水为预报对象,通过构建单独考虑前期春季海温为预报因子,以及同时引入前期春季海温与归一化植被指数为因子,建立最优子集回归预报模型,并采用交叉验证以及空间重建等统计手段,对比分析引入陆面植被因子前、后长江流域夏季降水预报效果改善的情况,评估春季陆面植被对长江流域夏季降水可预报性的影响及预报效果的稳健性

2 资料选取与方法介绍2.1 资料选取

所用资料有:(1)1951—2006年中国国家气候中心(NCC)160站降水资料,这里选取(28°—34°N,105°—122°E)区域内降水表示长江流域,共33个站,覆盖长江中上游及下游区域,具体空间分布和站点见图 1表 1;(2)美国国家海洋和大气局(NOAA)提供的1854—2008年海温资料,空间分辨率为2°×2°,文中仅选取太平洋和印度洋区域(40°S—60°N,40°E—80°W);(3)美国马里兰大学全球监测与模型研究中心(GIMMS)提供的归一化植被指数卫星观测资料,时间长度为1982年1月—2006年12月,分辨率为8 km×8 km(Tucker et al,19862005),这里将归一化植被指数资料处理为1°×1°月平均资料,区域限定为亚洲区域(10°S—60° N,60°—140° E)。考虑到与归一化植被指数资料时间一致,降水、海温建模分析时均截取1982—2006年时段,另外计算6—8月(JJA)平均降水量表征长江流域夏季降水,而前期春季海温、归一化植被指数则计算3—5月(MAM)平均值。

图 1 长江流域(28°-34°N,105°-122°E)33个站点空间分布以及 夏季降水旋转经验正交函数分解的分区Fig. 1 Distributions of the 33 stations of the Yangtze River Basin and the corresponding subdivisions of the summer precipitation based on REOF
表 1 长江流域夏季降水分区及其与旋转经验正交函数分解主成分(PC)时间系数的相关分析Table 1 The subdivisions of the summer precipitation and correlations with REOF PCs over the Yangtze River Basin
区域 名称 站点 时间系数 相关系数
YZ-Ⅰ 嘉陵江流域地区 达县 重庆 南充 内江 汉中 6 0.77
YZ-Ⅱ 汉江—淮河地区 南阳 信阳 安康 郧县 4 0.80
YZ-Ⅲ 淮河流域地区 清江 蚌埠 阜阳 东台 南京 合肥 3 0.82
YZ-Ⅳ 汉江—洞庭湖地区 汉口 钟祥 岳阳 宜昌 常德 恩施 酉阳 2 0.77
YZ-Ⅴ 长江三角洲地区 上海 杭州 安庆 屯溪 九江 宁波 5 0.78
YZ-Ⅵ 洞庭湖—鄱阳湖地区 衢县 贵溪 南昌 长沙 温州 1 0.89
2.2 方法介绍

(1)降水预报序列构建:利用旋转经验正交函数对长江流域夏季降水资料(Z)进行时空分解:

式中,A=(a1a2,…,ak),F=(f1f2,…,fk);ai、fi (i=1,2,…,k)分别为提取的第i个旋转经验正交函数分解的空间模态及其对应的时间序列。利用上述结果对长江流域进行降水区划,从图 2表 1不难看出,旋转经验正交函数分解的空间模态能够较好地体现所在区域特征,各模态时间序列fi(i=1,2,…,k)与对应区域降水平均序列相关系数超过0.75,为便于空间预报效果分析,这里选取fi代表所在区域降水预报序列,同时依据地理特征从西向东、由北向南依次将序列重新记为Y=(y1y2,…,yk)。

图 2 长江流域夏季降水旋转经验正交函数分解(REOF)模态的空间分布
((a)—(f)分别对应前6个旋转经验正交函数分解模态,阴影区的旋转载荷向量绝对值大于等于0.5)
Fig. 2 Spatial patterns of the REOF modes of the summer precipitation over the Yangtze River Basin
((a)-(f) are corresponding to the first six REOF modes, and areas with values larger than 0.5 are shaded)

(2)预报因子序列选取:通过相关普查,分析不同降水区域 yi(i=1,2,…,k)与前期春季海温和归一化植被指数的相关关系,选取相关系数较大、置信度超过0.9且数量超过5个以上的同号格点区域,然后以矩形盒子覆盖上述区域(具体选取区域见图 3、4),并计算区域平均值序列。考虑到预报量yi为标准化时间系数,为了保持一致,这里对上述所得预报因子序列进行相应标准化处理,由此分别获得mi个海温、ni个归一化植被指数标准化预报因子集,合记为Xi=(x1SSTx2SST,…,xmiSST,xh1x2NDVI,…,xniNDVI)。

图 3 长江流域6个区域夏季降水与前期太平洋、印度洋区域春季海温的相关分析
(a. 嘉陵江流域地区(YZ-Ⅰ),b. 汉江—淮河地区(YZ-Ⅱ), c. 淮河流域地区(YZ-Ⅲ), d. 汉江—洞庭湖地区(YZ-Ⅳ), e. 长江三角洲地区(YZ-Ⅴ), f. 洞庭湖—鄱阳湖地区(YZ-Ⅵ); 浅色、深色阴影区分别表示通过0.10、0.05显著性检验,红色方框区域为最优子集回归模型入选预报因子)
Fig. 3 Correlation patterns of the summer precipitation over the six subareas of the Yangtze River Basin and the previous spring SST of the Pacific Ocean and Indian Ocean
(a. YZ-Ⅰ, b. YZ-Ⅱ, c. YZ-Ⅲ, d. YZ-Ⅳ, e. YZ-Ⅴ, f. YZ-Ⅵ; The deep and light shaded areas are respectively significant at the 0.10 and 0.05 significant level, and the regions for selected predictors are shown in red bold boxes)

(3)预报模型构建与评价:对于入选的预报因子集Xi(i=1,2,…,k),分别建立基于海温和同时引入海温与归一化植被指数两种情形下第i个降水预报区yi最优子集回归模型(OSR):

式中,xi1xi2,…,xip为预报因子集Xi中入选的p个因子(这里因子个数p最大限定为3个,事实上后续分析也可以看出预报因子个数增加到3个时拟合效果已经增幅趋缓)。

考虑到资料长度有限,这里引入Bootstrap抽样(王松桂,1985)思想用来评估模型(2)参数估计的稳健性,通过1000次Bootstrap抽样重建回归获得模型参数95%置信区间估计值和显著性水平;其次,通过对比两种情形下模型拟合优度R2和观测值yi与模型(2)估计值i相关系数R的改善情况,评估陆面植被对长江流域各区域夏季降水预报的贡献;而为了进一步分析长江流域各站点预报情况,这里借助空间重建用预估值yi替代对应区域的旋转经验正交函数时间序列fj,并代入式(1)可得预估值,从而获得长江流域各站点的预报序列。通过计算并对比两种情形下站点观测值与预报值的相关系数R,评估春季陆面植被对长江流域各站点夏季降水可预报性的影响。

(4)交叉验证(CV)预报:由于植被资料时段为1982—2006年共25年,长度有限,为进一步验证模型预报效果和稳健性,这里引入交叉验证预报方法,具体做法是:首先将观测值yi依据时间等分为5段,即1982—1986、1987—1991、1992—1996、1997—2001和2002—2006年;其次逐个保留其中一段5年资料作为独立检验样本,即用剩余资料建立模型(2),然后对保留的5年进行预报,由此得到25年的预报值序列i,类似于(3)中的方法,借助于空间重建和相关系数R来评估长江流域夏季降水区域、站点交叉验证预报效果。

3 长江流域分区预报模型的构建3.1 长江流域夏季降水区划

考虑到长江流域夏季降水存在明显的区域性差异(张永领等,2006张文等,2007),利用旋转经验正交函数分解对长江流域(28°—34°N,105°—122°E)33个站点1951—2006年夏季平均降水量进行时空分解,从而获得长江流域的基本空间型和降水区划。图 2给出了长江流域夏季降水前7个经验正交函数分解主分量(累计方差为74%)进行极大方差旋转分解后前6个旋转经验正交函数分解空间模态的分布(其中由于第7个旋转经验正交函数分解主分量只有温州一个站载荷较大,未给出其模态分布)。根据旋转经验正交函数分解空间模态的高载荷分布(等值线绝对值不小于0.5)将长江流域分为6个区,按照从西向东、由北向南的原则依次命名为:嘉陵江流域地区(YZ-Ⅰ)、汉江—淮河地区(YZ-Ⅱ)、淮河流域地区(YZ-Ⅲ)、汉江—洞庭湖地区(YZ-Ⅳ)、长江三角洲地区(YZ-Ⅴ)和洞庭湖—鄱阳湖地区(YZ-Ⅵ),其中将温州并到第6区中(各分区中包含站点和对应模态详见图 1表 1)。表 1中最后一列给出了6个区域站点平均值与旋转经验正交函数分解对应时间系数(PC值)的相关系数,相关系数均超过0.75,说明6个旋转经验正交函数分解的PC系数代表了6个区域的特征。而以6个旋转经验正交函数分解模态及其对应PC时间系数进行重建观测序列,对比重建序列与原始序列的曲线(图略),可见总体拟合效果较好(解释程度达到0.85)。因此,这里选取PC时间系数代表上述6个区域的降水特征作为分析、预报对象,所得模拟、预测序列代入到旋转经验正交函数分解关系模型中进行空间重建,用于分析空间各站模拟、预测效果,这与冯蕾等(2011)中所用方法相似。

3.2 长江流域夏季降水与前期春季海温、 归一化植被指数的相关分析

由于植被指数资料仅为1982—2006年,这里截取上述代表 6个区域的PC时间系数1982—2006年的值,并与前期春季海温以及归一化植被指数进行相关分析(图 34),为后续构建最优子集回归模型提供预报因子集合。

图 4 长江流域6个区域夏季降水与前期亚洲区域春季归一化植被指数的相关分析
(a. 嘉陵江流域地区(YZ-Ⅰ),b. 汉江—淮河地区(YZ-Ⅱ), c. 淮河流域地区(YZ-Ⅲ), d. 汉江—洞庭湖地区(YZ-Ⅳ), e. 长江三角洲地区(YZ-Ⅴ), f. 洞庭湖—鄱阳湖地区(YZ-Ⅵ); 浅色、深色阴影区分别表示通过0.10、0.05显著性检验,红色方框区域为最优子集回归模型入选预报因子)
Fig. 4 Correlation patterns of the summer precipitation over the six subareas of the Yangtze River Basin and the previous spring NDVI over Asian areas
(a. YZ-Ⅰ,b. YZ-Ⅱ, c. YZ-Ⅲ, d. YZ-Ⅳ, e. YZ-Ⅴ, f. YZ-Ⅵ; The deep and light shaded areas are respectively significant at the 0.10 and 0.05 significant level, and the regions for selected predictors are shown in red bold boxes)
3.2.1 长江流域夏季降水与前期春季海温的相关

从长江流域6个区域夏季降水与前期春季海温的相关分析场(图 3)不难看出,太平洋区域是中国长江流域夏季降水预报的主要信号区,在影响中国东部夏季降水中起主导作用(郭玲等,2012祁莉等,2014)。由图 3可见,长江中上游嘉陵江流域地区(图 3a)在太平洋区域相关信号区最大,其与赤道东太平洋(TEP)海温成显著正相关,相关系数均在0.5左右,与赤道西太平洋(TWP)以及澳大利亚东部南太平洋(SP)则成负相关,相关系数在-0.5左右;汉江—洞庭湖地区与嘉陵江流域地区相关分布相似,但相关系数和显著区域大小均减小较多,相关系数仅有0.2—0.3(图 3d);赤道东太平洋海温与洞庭湖—鄱阳湖地区则为弱正相关(图 3f),与长江下游长江三角洲地区则为弱负相关(图 3e),而在长江中下游以北汉江—淮河流域地区则基本无信号(图 3bc)。

印度洋也是影响中国长江流域降水的关键海区,前期春季印度洋海温与长江中上游嘉陵江流域地区、汉江—洞庭湖地区(图 3ad)及长江下游以北淮河流域地区相关程度较高(图 3c),均成正相关,淮河流域地区显著区域最大,信号最强,相关达到0.4左右;而在汉江—淮河地区、长江三角洲地区前期春季印度洋相关信号减弱,成弱负相关(图 3be),仅-0.2左右,且相关显著区域偏西或偏东;在洞庭湖—鄱阳湖地区前期春季印度洋则基本无显著相关信号(图 3f)。Yang等(20072010)系列研究表明印度洋在ENSO对中国东亚夏季降水影响中扮演着充电器作用:赤道东太平洋厄尔尼诺(拉尼娜)现象出现后,通过大气桥和海洋动力过程对印度洋充电,形成热带印度洋海盆一致模(IOM);该模态一直持续到夏季,通过南亚高压、西太平洋副高对中国东亚夏季降水异常产生影响,导致印度洋放电。

3.2.2 长江流域夏季降水与前期春季归一化植被指数的相关

与海温相关分析类似,在长江流域6个降水区与前期春季归一化植被指数相关分析场(图 4)中,长江中上游嘉陵江流域地区显著区域大小依然最大,而长江中游以北汉江—淮河地区、长江下游三角洲地区植被信号依然较弱。影响长江中上游嘉陵江流域地区降水的主要陆面植被区为亚洲中西部南北两个区域,均为一致负相关,南部显著区域为印度半岛(IND)、中南半岛(ICP),北部主要为蒙古草原(MO)及其以北区域,随着该区域前期春季植被的增加,嘉陵江流域地区降水会呈现减少趋势。长江中游汉江—洞庭湖地区降水(图 4d)则主要与亚洲北部(NA)俄罗斯中部区域以及西部青藏高原区域前期春季植被成负相关。长江中游以北区域汉江—淮河地区(图 4b)最大负相关显著区位于相邻长江中上游片区植被区域,而长江下游淮河流域地区降水则与前期春季中国中、南部区域植被成正相关(图 4c),长江三角洲地区则与前期春季青藏高原区域植被成正相关(图 4e),长江中下游以南洞庭湖—鄱阳湖地区则与亚洲西部区域(WA)前期春季植被成正相关,与东南亚区域(SEA)植被成负相关(图 4f)。

不难看出,无论是与前期春季海温还是归一化植被指数,长江流域6个降雨区相关显著区域大小上存在着明显区域差异,嘉陵江流域地区、汉江—洞庭湖地区通过检验的显著区域最大,洞庭湖—鄱阳湖地区居中,汉江—淮河地区、淮河流域地区以及长江三角洲地区则显著区域最小,这从客观上说明了降水分区预报的必要性。相对于传统太平洋海温预报因子,前期春季印度洋海温也是影响中国长江流域各区域夏季降水预报的关键预报因子;而影响中国长江流域各区域夏季降水的前期春季陆面归一化植被指数相关显著区主要是青藏高原植被,其次为北方蒙古草原及其以北区域,南方印度半岛、中南半岛以及东南亚等区域。范广洲等(2009)Lee等(20082009)在研究春季植被对夏季降水影响时,通过分析潜热、感热通量、地表温度和850 hPa风场等变化,也指出青藏高原、亚洲北方草原、印度半岛以及东南亚等区域归一化植被指数变化,导致地面热源发生变化,从而引起中国及相关地区大气环流变化,进而造成中国夏季降水发生相应改变。

3.3 长江流域夏季降水分区预报模型构建

为了构建长江流域夏季降水各分区预报模型,首先从上述与前期春季海温、归一化植被指数相关分析场中选取高值相关区域来构建海陆下垫面预报因子集,具体选取区域见图 3、4中矩形区(实线(虚线)为正(负)相关区域),通过计算区域平均值来获得预报因子序列。表 2给出了各分区海陆下垫面预报因子及其相关系数情况,代入式(2),可得长江流域6个降水区单纯海温和同时引入海温与归一化植被指数因子两种情形下1、2、3个因子的最优子集回归预报模型(限于篇幅,详细建模过程不再给出)。

表 2 江流域各分区海陆下垫面预报因子及其相关系数Table 2 The areas and the correlation values of the potential predictors of the 6 different subareas of the Yangtze River Basin
响应变量 预测因子 区域 相关系数 p
嘉陵江流域地区(YZ-Ⅰ) 海温 TEP 10°S—5°N, 135°—80°W 0.562 0.003
SP 40°—25°S,165°—80°W -0.652 0.000
TWP 28°S—20°N,147°E—148°W -0.691 0.000
TIO 21°—12°S, 84°—105°E 0.495 0.012
NP 39°—48°N, 165°E—165°W -0.512 0.009
植被 NA 50°—60°N, 60°—90°E -0.551 0.004
MO 47°—50°N, 85°—96°E -0.541 0.005
IND 15°—21°N, 74°—84°E -0.623 0.001
ICP 10°—20°N, 98°—107°E -0.607 0.001
汉江—淮河地区(YZ-Ⅱ) 海温 SIO 32°—25°S, 40°—52°E -0.390 0.054
SNEP 25°—35°N, 128°—110°W 0.351 0.085
植被 TP 31°—35°N, 95°—104°E -0.577 0.003
淮河流域地区(YZ-Ⅲ) 海温 TIO 11°S—16°N,50°—75°E 0.495 0.012
植被 SA 20°-30°N,93°—103°E 0.577 0.003
MEC 32°—36°N, 112°—120°E 0.458 0.021
汉江—洞庭湖地区(YZ-Ⅳ) 海温 SNEP 29°—39°N,170°—150°W -0.424 0.035
TIO 21°—10°S, 65°—95°E 0.465 0.019
植被 NA 53°—60°N, 85°—116°E -0.524 0.007
TP 32°—40°N, 70°—79°E -0.487 0.014
长江三角洲地区(YZ-Ⅴ) 海温 TIO 15°—9°S, 120°—140°E -0.478 0.016
植被 TP 32°—35°N, 80°—92°E 0.497 0.012
洞庭湖—鄱阳湖地区(YZ-Ⅵ) 海温 TEP 18°S—0°,105°—88°W 0.424 0.035
SNWP 20°—40°N,110°—125°E 0.572 0.003
植被 WA 25°—30°N, 68°—80°E 0.535 0.006
SEA 3°S—3°N,110°—116°E -0.538 0.006
注:预报因子名称以因子所在区域来命名,海温方面:TEP、TWP、SP、NP、SNWP、SNEP依次表示热带东太平洋、热 带西太平洋、南太平洋、北太平洋、西北太平洋以及副热带东北太平洋,TIO、SIO分别表示热带、副热带印度洋。 植被方面:MO、IND、ICP、MEC、TP依次表示蒙古草原、印度半岛、中南半岛、中国中东部、青藏高原区域; NA、SA、WA、SEA分别表示亚洲北部、南部、西部和东南部区域。

为了进一步说明因子选取和回归系数估计情况,同时便于后期预报效果对比分析,这里选取回归拟合、交叉验证预报相对较好的情形作为代表进行分析(表 3,出于方便,用R(Reg)、R(CV)分别表示两者情形下模型模拟、交叉预报相关系数,R2(Reg)为模拟模型拟合优度)。从表 3可见,仅前期春季海温因子时,入选海洋海温因子区域最多的为热带印度洋,其次为副热带东北太平洋,前者依次在嘉陵江流域地区、淮河流域地区、汉江—洞庭湖地区以及长江三角洲地区入选为预报因子,后者依次在汉江—淮河地区、汉江—洞庭湖地区入选;其他因子主要是赤道东太平洋、赤道西太平洋、澳大利亚东部南太平洋、副热带西北太平洋和副热带印度洋。引入前期春季归一化植被指数因子后,印度洋海温因子TIO、SIO除嘉陵江流域地区外其他区域均得到保留。可见,相对于太平洋,前期春季印度洋海温也是影响中国长江流域各区域夏季降水预报的关键因子。相对于海温因子,入选前期春季归一化植被指数因子最多的为青藏高原地区,其他入选植被因子基本与相关分析场类似,主要以亚洲南北两个区域的因子为主。

表 3 长江流域各分区两种预报模型最终入选海陆下垫面预报因子及其结果Table 3 The introduced predictors and forecast results of the 6 subareas of the Yangtze River Basin for the two different forecasting models
区域 模型1: 海温
海温 R (Reg) R 2(Reg) R (CV)
嘉陵江流域地区 SP-0.25(-0.59,0.12) TWP **-0.44(-0.78,-0.15) TIO **0.38(0.09,0.65) 0.82 0.67 0.68
汉江—淮河地区 SIO-0.34(-0.73,-0.07) SNEP *0.29(-0.00,0.67) 0.48 0.23 0.32
淮河流域地区 TIO *0.55(0.06,1.17) 0.50 0.25 0.34
汉江—洞庭湖地区 SNEP-0.24(-0.58,0.21) TIO *0.30(-0.02,0.65) 0.52 0.27 0.36
长江三角洲地区 TIO ***-0.53(-0.80,-0.29) 0.48 0.23 0.35
洞庭湖—鄱阳湖地区 TEP0.25(-0.30,0.90) SNWP **0.55(0.00,0.88) 0.61 0.37 0.38
区域 模型2:海温+归一化植被指数
海温 归一化植被指数 R (Reg) R 2(Reg) R(CV)
嘉陵江流域地区 TEP ***0.46(0.16,0.78) IND ***-0.45(-0.65,-0.07) ICP *-0.26(-0.63,-0.05) 0.85 0.72 0.73
汉江—淮河地区 SIO-0.20(-0.56,0.18) SNEP ** 0.32(0.06,0.66) TP **-0.53(-1.00,-0.20) 0.70 0.49 0.61
淮河流域地区 TIO0.26(-0.30,1.12) SA *0.45(0.05,0.88) MEC0.13(-0.34,0.53) 0.64 0.41 0.50
汉江-洞庭湖地区 TIO **0.32(0.06,0.59) NA ***-0.36(-0.62,-0.13) TP-0.22(-0.68,0.01) 0.71 0.50 0.63
长江三角洲地区 TIO *-0.41(-0.79,0.11) TP0.44(0.03,1.10) 0.61 0.37 0.53
洞庭湖—鄱阳湖地区 SNWP **0.45(0.03,0.73) WA0.17(-0.31,0.62) SEA *-0.47(-1.34,-0.20) 0.72 0.52 0.42
注:因子下面为相应参数估计值和bootstrap回归95%置信区间估计值;*,**,***分别表示bootstrap回归该因子通过0.10、0.05、0.01显著性水平检验。

对于上述入选预报因子,表 3给出了两种情形下最优子集回归模型参数估计值。对比表 2,可见各因子估计值符号总体上与前面相关分析一致,能够较好地反映预报因子对相应区域夏季降水预报的贡献。值得关注的是,由于植被资料所限,样本长度仅为25年,这影响了参数估计的可靠性。因此,这里引入bootstrap抽样检验思想,通过1000次bootstrap抽样重建回归获得预报因子参数95%置信区域估计值和显著性水平。从表 3可见,两种情形预报因子大部分通过了bootstrap回归显著性水平检验:单纯海温因子情形,6个分区预报模型11个因子中有7个因子通过了0.05或0.10显著性水平检验,其中长江三角洲地区TIO更是通过0.01显著性水平检验;引入陆面归一化植被指数因子后,6个分区预报模型17个因子中有11个因子通过了0.05或0.10显著性水平检验,其中嘉陵江流域地区TEP、IND以及汉江—洞庭湖地区NA因子通过了0.01显著性水平检验。而各因子Bootstrap回归95%置信区间估计值所在区间也基本与上述参数估计符号一致,验证了模型参数估计的有效性。但也不难看出,未通过bootstrap显著性水平检验的预报因子,其置信区间估计值虽然总体覆盖相应参数估计值,但置信上下限估计值符号存在异号情况,这与最优子集回归模型选取因子有关,也与模型样本有限,估计值稳健性不够有着较大关系。

4 陆面植被引入前后夏季降水预测效果对比4.1 预测随因子数的变化

从各个分区预报模型在不同因子个数以及不同因子组合情形下长江流域各分区降水拟合和交叉验证预报结果(图 5)可见,随着入选预报因子数的增多,长江流域各分区拟合效果相应增加,而增加幅度方面则存在差异:只有1个因子的时候,相关系数普遍维持在0.4—0.6;当预报因子由1个增加到2个时,相关系数提升到了0.5—0.7,增幅明显;而2个因子增加到3个因子时相关系数增加幅度则不太明显(图 5a)。相对于上述拟合结果,交叉预报虽然相关系数平均降低0.10—0.15,但交叉预报特征总体上与上述结果类似,随着预报因子个数增多,交叉预报相关系数也随着增大,但在增大幅度上不如前者。不难看出,无论是拟合还是交叉验证,引入第3个因子后对降水预报提高贡献减小,后续增加因子意义已经不太明显,这也是本文分析中入选因子最终限定不超过3个的主要原因。

图 5 长江流域各分区夏季降水预报拟合结果
(a. 对应最优子集回归模型拟合结果,b. 交叉验证预报模型拟合结果;P1、P2、P3表示预报因子数为1、2、3个)
Fig. 5 Forecasting model results for the 6 different subareas of the Yangtze River Basin
(a. the modeling result of the optimal subset regression models, b. the modeling result of the cross validating models; P1, P2, and P3 indicate the different predictors number being 1, 2, and 3, respectively)
4.2 引入植被因子前后模拟效果的对比

为了进一步评估陆面植被对长江流域夏季降水可预报性的影响,这里选取上述各个因子数下拟合相关系数最高情形为代表(表 3),从各个分区模拟、空间模拟两个方面对比分析仅考虑海温和引入归一化植被指数因子后两种情形下的模拟效果。

4.2.1 长江流域夏季降水各分区模拟效果对比

表 3(图 5a)中可见,单纯海温情形长江以南及长江中上游降水拟合效果要好于长江中下游及其以北区域。嘉陵江流域地区拟合效果最好,相关系数达到0.82,拟合优度为0.67,其次为洞庭湖—鄱阳湖地区,相关系数为0.61,拟合优度为0.37,而淮河流域地区、汉江—洞庭湖地区相关系数只有0.51左右,拟合优度约为0.26;汉江—淮河地区以及长江三角洲地区则拟合效果最差,相关系数为0.48,拟合优度为0.23。

引入陆面归一化植被指数因子后,长江流域各分区拟合相关系数平均提高约0.14,提升幅度达26%,拟合优度平均提高约0.17,模型解释方差提升平均约60%(表 3图 5b):汉江—淮河地区、汉江—洞庭湖地区拟合效果更是得到大幅提升,相关系数增加了约0.20,拟合优度提升约0.25,模型解释方差达到或超过原来的2倍;其次淮河流域地区、长江三角洲地区也提升了约0.13,拟合优度提升约0.15;洞庭湖—鄱阳湖地区增加幅度略小,相关系数提升约0.11,拟合优度提升约0.15;而嘉陵江流域地区则增幅不明显。

4.2.2 长江流域夏季降水空间模拟效果对比

利用表 3最终入选模型模拟结果,结合长江流域旋转经验正交函数分解的空间模态,通过空间重建可得1982—2006年长江流域33个站点降水模拟序列。计算区域平均值,可得长江流域整个区域年际变化拟合曲线及其相关系数分布(图 6),结果显示:单纯考虑海温因子时长江流域整个区域模拟相关系数只有0.46,这与前面分区预报中长江中下游及其以北片区模拟效果较差有关(图 6a);引入归一化植被指数因子后,长江流域整个区域模型拟合效果得到显著提升,模拟相关系数增加0.16达到0.62(图 6b)。

图 6 长江流域夏季降水年际变化预报结果及其相关系数分布
(a.单纯海温因子预报模型, b. 海温和归一化植被指数组合因子预报模型; 实折线、破折线及点虚线分别对应观测、回归拟合和交叉预报序列)
Fig. 6 Inter-annual forecasting and correlations results of the Yangtze River Basin
(a. the only SST factors, b. both SST and NDVI factors; solid lines, dash lines and dotted lines are the observed values, forecasting time series of optimal subset regression models and cross validating models, respectively)

图 7则给出了长江流域夏季降水 33个站点模型拟合相关系数的空间分布。类似于区域预报,考虑单纯海温因子时(图 7a),拟合效果最好的为嘉陵江流域地区,站点相关系数平均约为0.68,最高达到0.75;而长江中下游则基本以长江为界,长江以南汉江—洞庭湖地区、洞庭湖—鄱阳湖地区预报效果要优于长江以北汉江—淮河地区、淮河流域地区,前者站点拟合相关系数在0.5—0.6,后者则在0.36—0.47;长江中游以北区域更是约在0.37,拟合效果最低,而长江下游长江三角洲地区虽然预报因子较少,但拟合相关系数也维持在0.43左右。

图 7 长江流域夏季降水33个站点最优子集回归预报相关系数的空间分布
(a. 单纯海温因子预报模型, b. 海温和归一化植被指数组合因子预报模型)
Fig. 7 Forecasting correlation distributions of the 33 stations over the Yangtze River Basin for the optimal subset regression models
(a. the only SST factors, b. both SST and NDVI factors)

引入植被因子后(图 7b),长江流域各个站点与对应各个分区预报一致,模拟相关系数得到大幅度提高,相关系数平均由0.49提升到了0.66,提高0.17。长江中下游以北区域汉江—淮河地区和淮河流域地区模拟效果得到显著改善,相关系数由仅考虑海温因子的0.36—0.47提高到0.54—0.72,平均提高0.23,局部站点更是提高0.30左右。长江三角洲地区引入归一化植被指数因子后相关系数也普遍提升0.15左右,由原来的约0.43提升到约0.58;汉江—洞庭湖地区也普遍提升约0.19,洞庭湖—鄱阳湖地区则提升约0.13;嘉陵江流域地区提升幅度不明显,只有0.09左右。

4.3 引入植被覆盖因子预测效果的对比

由于植被资料时段为1982—2006年,长度有限,这里采用每5年资料作为独立样本,利用表 3中构建的最优子集回归模型进行交叉验证预报,并类似于前述分析利用空间重建得到空间33个站点交叉验证预报值。

仅考虑海温因子,相对于前面各个分区模拟,各个分区交叉预报相关系数高低分布相似,但预测相关系数则由约0.57下滑到约0.41,说明模型稳健性较差(表 3),然而经过空间重建后长江流域整体预报则基本维持不变(图 6a);引入归一化植被指数因子后,各个分区交叉验证预测总体上得到提高,比仅考虑海温因子的情况提高了约0.17,达到了0.57,且模型稳健性得到小幅度提高,相对于其拟合情形降幅略微减小,只为0.14左右(表 3),长江流域整个区域交叉验证预测相关系数也由0.44提升到0.57,提升0.13(图 6b)。

图 8给出了长江流域33个观测站点交叉验证预报相关系数的空间分布。对比前面站点模拟情形(图 7),仅考虑海温因子时(图 8a),嘉陵江流域地区、汉江—洞庭湖地区站点交叉预报相关系数改变较小,仅降低约0.08;而长江中下游区域站点除拟合效果本来就不好的长江以北区域外,原来拟合效果较好的洞庭湖—鄱阳湖地区也下降明显,相关系数降幅最高达到0.27。引入归一化植被指数因子后(图 8b),长江中下游及其以北地区站点交叉预报效果得到明显提高,相对于仅考虑海温情形(图 8a),汉江—淮河地区交叉验证相关系数平均提升0.27,淮河流域地区提升0.17,汉江—洞庭湖地区也提升了约0.19,而交叉预报效果并不太好的长江三角洲地区、洞庭湖—鄱阳湖地区则在引进陆面植被后,站点交叉验证预报效果得到大幅度改进,相关系数普遍达到0.42、0.45,分别提高了约0.19、0.12。可见,相对于仅考虑海温因子情形,引入归一化植被指数因子后,陆面植被改进了长江流域中下游及其以北区域站点预报效果,提升了长江中下游及其以南区域站点预报稳健性。

图 8 长江流域夏季降水33个站点交叉验证预报相关系数的空间分布
(a. 单纯海温因子预报模型, b. 海温和归一化植被指数组合因子预报模型)
Fig. 8 Forecasting correlation distributions of the 33 stations over the Yangtze River Basin for the cross validating models
(a. the only SST factors, b. both SST and NDVI factors)
5 结论

利用最优子集回归、交叉验证以及空间重建等方法分别构建了单独考虑前期春季海温因子以及同时引入前期春季海温与归一化植被指数因子两种情形下的预报模型,通过对比分析引入春季归一化植被指数因子前、后长江流域夏季降水区域预报以及空间预报改善状况,系统评估了春季陆面归一化植被指数对长江流域夏季降水可预报性的影响以及预报效果的稳健性,结果表明:

(1)相对于传统的海温因子,春季陆面植被因子对长江流域夏季降水预报同样重要。引入春季归一化植被指数因子后,长江流域夏季降水预报得到明显改善,相关系数平均提升0.17,由原来的0.49提高到0.66,模型解释方差的提升幅度平均60%左右,尤其在单纯海温因子预报效果较差的长江中下游以北汉江—淮河地区和淮河流域地区,站点相关系数由仅考虑海温因子的0.36—0.47提高到0.54—0.72,平均提高0.23,局部站点更是提高0.30左右,模型解释方差提升约1倍。

(2)交叉验证预报表明,相对于仅考虑海温因子模拟情形,交叉预报相关系数下降较多,模型稳健性较低,引入归一化植被指数因子后,长江流域夏季降水预报稳健性得到明显提升,长江中下游及其以南的长江三角洲地区、洞庭湖—鄱阳湖地区改善尤为明显,交叉预报相关系数普遍达到0.42、0.45,分别提高了约0.19、0.12。

(3)长江流域降水可预报性存在明显的区域差异,嘉陵江流域地区、汉江—洞庭湖地区预报效果最好,汉江—淮河地区、淮河流域地区、长江三角洲地区预报效果最差,引入归一化植被指数因子后预报效果提升最明显,而洞庭湖—鄱阳湖地区虽然模拟效果较好,但预报稳健性较低,交叉验证相关系数降幅达到0.27,这也从侧面说明了长江流域夏季降水分区预报的重要性。

值得注意的是,本研究主要是从统计角度探讨了春季陆面植被对长江流域夏季降水可预报性的影响,而有关前期春季陆面植被是如何影响到中国长江流域夏季降水并没有进行过多探讨,这方面工作目前依然是个难点。关于陆面植被对气候的反馈认识主要是基于区域模式模拟获得,但由于陆面过程的复杂性,以及相关试验多采用简单的两极化敏感试验,存在着较大的不确定性和模式偏差。另外,相关分析多为局地或者同期反馈分析,关于非局地反馈以及季节或者更长时间尺度陆面植被对降水的反馈影响机制分析依然较少,这方面的工作有待后续进一步的研究分析。

致 谢:感谢美国威斯康星大学麦迪逊分校Michael Notaro教授提供了GIMMS逐月归一化植被指数资料,感谢中国国家气候中心和美国国家海洋大气局(NOAA)提供的降水和海温资料,感谢王志福和陈威霖老师在研究工作和论文撰写中提供的帮助和指导。

参考文献
崔林丽, 史军. 2011. 中国华东及其周边地区NDVI对气温和降水的月际响应特征. 自然资源学报, 26(12): 2121-2130. Cui L L, Shi J. 2011. Inter-monthly response characteristics of NDVI to the variation of temperature and precipitation in East China and its surrounding areas. J Nat Res, 26(12): 2121-2130 (in Chinese)
丁一汇, 李巧萍, 董文杰. 2005. 植被变化对中国区域气候影响的数值模拟研究. 气象学报, 63(5): 613-621. Ding Y H, Li Q P, Dong W J. 2005. A numerical simulation study of the impacts of vegetation changes on regional climate in China. Acta Meteor Sinica, 63(5): 613-621 (in Chinese)
范广洲, 周定文, 黄荣辉等. 2009. 我国夏季降水与青藏高原春季NDVI的关系. 应用气象学报, 20(4): 492-496. Fan G Z, Zhou D W, Huang R H, et al. 2009. The relationship between summer rainfall over China and spring NDVI variability in Tibetan Plateau. J Appl Meteor Sci, 20(4): 492-496 (in Chinese)
冯蕾, 魏凤英, 朱艳峰. 2011. 基于前春对流层温度和北大西洋涛动的中国夏季降水统计预测模型. 大气科学, 35(5): 963-976. Feng L, Wei F Y, Zhu Y F. 2011. A predictive model for summer precipitation over China based on upper tropospheric temperature and North Atlantic Oscillation in the preceding spring. Chin J Atmos Sci, 35(5): 963-976 (in Chinese)
符淙斌, 袁慧玲. 2001. 恢复自然植被对东亚夏季气候和环境影响的一个虚拟试验. 科学通报, 46(8): 691-695. Fu C B, Yuan H L. 2001. An virtual numerical experiment to understand the impacts of recovering natural vegetation on the summer climate and environmental conditions in East Asia. Chin Sci Bull, 46(14): 1199-1203
郭玲, 何金海, 祝从文. 2012. 影响长江中下游夏季降水的前期潜在预报因子评估. 大气科学, 36(2): 337-349. Guo L, He J H, Zhu C W. 2012. Evaluations of prior potential predictors for the summer rainfall in the middle and lower reaches of the Yangtze River. Chinese J Atmos Sci, 36(2): 337-349 (in Chinese)
霍飞, 江志红, 刘征宇. 2014. 春夏季青藏高原积雪对中国夏末秋初降水的影响及其可能机制. 大气科学, 38(2): 352-362. Huo F, Jiang Z H, Liu Z Y. 2014. Impacts of late spring Tibetan Plateau snow cover on early autumn precipitation. Chinese J Atmos Sci, 38(2): 352-362 (in Chinese)
祁莉, 王晓芳, 何金海等. 2014. 前期西太平洋暖池热含量异常影响长江中下游夏季降水的可能途经. 地球物理学报, 57(6): 1769-1781. Qi L, Wang X F, He J H, et al. 2014. The approach of the previous anomalous heat content in the Western Pacific warm pool affecting the summer rainfall over the middle and lower reaches of the Yangtze River. Chinese J Geophys, 57(6): 1769-1781 (in Chinese)
钱维宏, 陆波. 2010. 我国汛期季度降水预报得分和预报技巧. 气象, 36(10): 1-7. Qian W H, Lu B. 2010. Score and skill of seasonal forecasts of summer precipitation in China. Meteor Mon, 36(10): 1-7 (in Chinese)
王松桂. 1985. 线性模型参数估计的新进展. 数学进展, 14(3): 193-204. Wang S G. 1985. New progress on the parameter estimation in linear model. Adv Math, 14(3): 193-204 (in Chinese)
魏凤英, 黄嘉佑. 2010. 大气环流降尺度因子在中国东部夏季降水预测中的作用. 大气科学, 34(1): 202-212. Wei F Y, Huang J Y. 2010. A study of downscaling factors of atmospheric circulations in the prediction model of summer precipitation in Eastern China. Chinese J Atmos Sci, 34(1): 202-212 (in Chinese)
张井勇, 董文杰, 叶笃正等. 2003. 中国植被覆盖对夏季气候影响的新证据. 科学通报, 48(1): 91-95. Zhang J Y, Dong W J, Ye D Z, et al. 2003. New evidence for effects of land cover in China on summer climate. Chinese Sci Bull, 48(4): 401-405
张少波, 陈玉春, 吕世华等. 2013. 青藏高原植被变化对中国东部夏季降水影响的模拟研究. 高原气象, 32(5): 1236-1245. Zhang S B, Chen Y C, Lv S H, et al. 2013. Numerical simulation of impact of vegetation in Qinghai-Xizang Plateau on summertime precipitation in Eastern China. Plateau Meteor, 32(5): 1236-1245 (in Chinese)
张文, 寿绍文, 杨金虎. 2007. 长江中下游地区汛期极端降水量的异常特征分析. 气象, 33(3): 61-67. Zhang W, Shou S W, Yang J H. 2007. Analysis of anomaly characteristics of extreme precipitation in flood season in the mid-lower reaches of the Yangtze River. Meteor Mon, 33(3): 61-67 (in Chinese)
张永领, 高全洲, 丁裕国等. 2006. 长江流域夏季降水的时空特征及演变趋势分析. 热带气象学报, 22(2): 161-168. Zhang Y L, Gao Q Z, Ding Y G, et al. 2006. Analysis of time-spatial characteristics and evolutional trends of summer precipitation in the Yangtze River catchment. J Trop Meteor, 22(2): 161-168 (in Chinese)
郑益群, 钱永甫, 苗曼倩等. 2002a. 植被变化对中国区域气候的影响Ⅰ: 初步模拟结果. 气象学报, 60(1): 1-16. Zheng Y Q, Qian Y F, Miao M Q, et al. 2002a. The effects of vegetation change on regional climate I: Simulation results. Acta Meteor Sinica, 60(1): 1-16 (in Chinese)
郑益群, 钱永甫, 苗曼倩等. 2002b. 植被变化对中国区域气候的影响Ⅱ: 机理分析. 气象学报, 60(1): 17-30. Zheng Y Q, Qian Y F, Miao M Q, et al. 2002b. The effects of vegetation change on regional climate Ⅱ: Mechanisms. Acta Meteor Sinica, 60(1): 17-30 (in Chinese)
周放, 孙照渤, 许小峰等. 2014. 中国东部夏季暴雨日数的分布特征及其与大气环流和海温的关系. 气象学报, 72(3): 447-464. Zhou F, Sun Z B, Xu X F, et al. 2014. Spatiotemporal characteristics of summer rainstorm days in eastern China and their relationships with atmospheric circulation and SST. Acta Meteor Sinica, 72(3): 447-464 (in Chinese)
周锁铨, 陈万隆. 1995. 青藏高原植被下垫面对东亚大气环流影响的数值试验. 南京气象学院学报, 18(4): 536-542. Zhou S Q, Chen W L. 1995. Numerical experiments with effect of Tibetan vegetation on East Asian atmospheric circulations. J Nanjing Inst Meteor, 18(4): 536-542 (in Chinese)
Barriopedro D, Gouveia C M, Trigo R M, et al. 2012. The 2009/10 drought in China: Possible causes and impacts on vegetation. J Hydrometeor, 13(4): 1251-1267
Ding Y H, Sun Y, Wang Z Y, et al. 2009. Inter-decadal variation of the summer precipitation in China and its association with decreasing Asian summer monsoon Part II: Possible causes. Int J Climatol, 29(13): 1926-1944
Lee E, Chase T N, Rajagopalan B. 2008. Highly improved predictive skill in the forecasting of the East Asian summer monsoon. Water Resources Res, 44(10): 285-297
Lee E, Chase T N, Rajagopalan B, et al. 2009. Effects of irrigation and vegetation activity on early Indian summer monsoon variability. Int J Climatol, 29(4): 573-581
Li W P, Xue Y K. 2005. Numerical simulation of the impact of vegetation index on the interannual variation of summer precipitation in the Yellow River Basin. Adv Atmos Sci, 22(6): 865-876
Liu Z Y, Notaro M, Kutzbach J, et al. 2006. Assessing global vegetation-climate feedbacks from observations. J Climate, 19(5): 787-814
Notaro M, Chen G S, Liu Z Y. 2011. Vegetation feedbacks to climate in the global monsoon regions. J Climate, 24(22): 5740-5756
Notaro M, Gutzler D. 2012. Simulated impact of vegetation on climate across the North American monsoon region in CCSM3.5. Climate Dyn, 38(3-4): 795-814
Peng J, Dong W J, Yuan W P, et al. 2012. Responses of grassland and forest to temperature and precipitation changes in Northeast China. Adv Atmos Sci, 29(5): 1063-1077
Tucker C J, Sellers P J. 1986. Satellite remote sensing of primary production. Int J Remote Sens, 7(11): 1395-1416
Tucker C J, Pinzon J E, Brown M E, et al. 2005. An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data. Int J Remote Sens, 26(20): 4485-4498
Wang F Y, Notaro M, Liu Z Y, et al. 2014. Observed local and remote influences of vegetation on the atmosphere across North America using a model-validated statistical technique that first excludes oceanic forcings. J Climate, 27(1): 362-382
Wang Y X, Zhao P, Yu R C, et al. 2010. Inter-decadal variability of Tibetan spring vegetation and its associations with eastern China spring rainfall. Int J Climatol, 30(6): 856-865
Wu Z W, Wang B, Li J P, et al. 2009. An empirical seasonal prediction model of the East Asian summer monsoon using ENSO and NAO. J Geophys Res, 114(D18): D18120, doi: 10.1029/2009JD011733
Xue Y K. 1996. The impact of desertification in the Mongolian and inner Mongolian grassland on the regional climate. J Climate, 9(9): 2173-2189
Yang J L, Liu Q Y, Xie S P, et al. 2007. Impact of the Indian Ocean SST basin mode on the Asian summer monsoon. Geophys Res Lett, 34(2): L02708, doi: 10.1029/2006GL028571
Yang J L, Liu Q Y, Liu Z Y. 2010. Linking observations of the Asian monsoon to the Indian ocean SST: Possible roles of Indian ocean basin mode and dipole mode. J Climate, 23(21): 5889-5902
Zhang J Y, Dong W J, Fu C B. 2005. Impact of land surface degradation in northern China and southern Mongolia on regional climate. Chinese Sci Bull, 50(1): 75-81
Zuo Z Y, Zhang R H, Zhao P. 2011. The relation of vegetation over the Tibetan Plateau to rainfall in China during the boreal summer. Climate Dyn, 36(5-6): 1207-1219