自组织映射神经网络(SOM)降尺度方法对江淮流域逐日降水量的模拟评估

利用1961~2002年ERA-40逐日再分析资料和江淮流域56个台站逐日观测降水量资料,引入基于自组织映射神经网络(Self-Organizing Maps,简称SOM)的统计降尺度方法,对江淮流域夏季(6~8月)逐日降水量进行统计建模与验证,以考察SOM对中国东部季风降水和极端降水的统计降尺度模拟能力。结果表明,SOM通过建立主要天气型与局地降水的条件转换关系,能够再现与观测一致的日降水量概率分布特征,所有台站基于概率分布函数的Brier评分(Brier Score)均近似为0,显著性评分(Significance Score)全部在0.8以上;模拟的多年平均降水日数、中雨日数、夏季总降水量、日降水强度、极端降水阈值和极端降水贡献率区域平均的偏差都低于11%;并且能够在一定程度上模拟出江淮流域夏季降水的时间变率。进一步将SOM降尺度模型应用到BCCCSM1.1(m)模式当前气候情景下,评估其对耦合模式模拟结果的改善能力。发现降尺度显著改善了模式对极端降水模拟偏弱的缺陷,对不同降水指数的模拟较BCC-CSM1.1(m)模式显著提高,降尺度后所有台站6个降水指数的相对误差百分率基本在20%以内,偏差比降尺度前减小了40%~60%;降尺度后6个降水指数气候场的空间相关系数提高到0.9,相对标准差均接近1.0,并且均方根误差在0.5以下。表明SOM降尺度方法显著提高日降水概率分布,特别是概率分布曲线尾部特征的模拟能力,极大改善了模式对极端降水场的模拟能力,为提高未来预估能力提供了基础。     

全球增暖1.5/2℃下中国区域极端降水的风险变化及其影响因子

利用24个CMIP6全球气候模式的逐日降水模拟资料,基于广义极值分布(GEV)模型,研究了全球增暖1.5/2℃下我国20、50和100 a重现期极端降水的未来风险变化。可以发现,相对于历史时期(1995—2014年),全球升温1.5和2℃下极端降水发生概率风险空间分布相近,总体上呈现增加趋势,但额外增暖0.5℃将导致更高的风险。如50 a重现期极端降水,在增暖1.5/2℃下其重现期将分别变为17/14 a,极端降水将变得更加频繁。不同区域对气候变暖的响应存在区域差异,其中中国西部长江黄河中上游和青藏高原地区、中国东部长江黄河中下游及其以南地区,极端降水发生概率比达到3以上,局部更是达到5以上,为我国极端降水气候变化响应高敏感区域。进一步,基于概率分布函数从理论角度探讨了位置和尺度参数对发生概率风险的影响与贡献度量,并用于探讨极端降水气候平均态和变率变化对极端降水发生风险的影响,结果显示：位置和尺度参数的增量变化、风险变化率存在着显著的东西部差异,从而导致极端降水发生风险的影响因素存在差异。如中国西部尽管极端降水气候平均态和变率变化幅度不大,但因风险变化率较高,从而导致该区域的发生风险大...     

基于转移累计概率分布统计降尺度方法的未来降水预估研究:以湖南省为例

基于最新一代CMIP5（Coupled Model Intercomparison Project Phase 5）模式历史情景和未来RCP4.5情景下的模式逐日降水数据,使用转移累计概率分布（CDF-t）统计降尺度方法,从空间变化和时间变率两个方面评估该降尺度方法对湖南日降水量模拟能力的改善效果,并在此基础上对未来降水量变化进行预估。结果表明, CMIP5气候模式由于分辨率较低,无法细致反映湖南地形变化和大气环流影响导致的区域降水变化特征。经过CDF-t统计降尺度处理之后,模式对湖南降水的时、空分布模拟与实况更为接近,绝大部分模式对降水空间结构的模拟能力都有显著提高。基于CDF-t统计降尺度的多模式集合预估结果表明,21世纪湖南省日降水量呈弱的增多趋势（0.95%/（10 a））。21世纪初、中和末期相对于1986—2005年的气候平均态,湖南省日降水量分别增加了4.6%、5%和5.2%。3个时期湖南省日平均降水变化的空间分布存在较强的一致性,皆表现为湖南西北、东北和东南3个地区降水增幅最为显著,且随着辐射强迫的增大,3个地区降水增幅也呈递增趋势。需要指出的是,预估结果在模式之间存在一定差异,并且这种差异随着辐射强迫的增大而增大。      

CDF-T方法在站点尺度日降水预估中的应用

基于1961-2015年上海降水观测数据和8个全球气候模式GCMs模拟的日降水量数据,采用累计概率分布函数构建转换模型CDF-T建立了站点尺度日降水量的统计降尺度模型。结果表明,降尺度模型显著改善了GCMs对降水日数偏多、降水强度偏低和降水量偏少的模拟结果。与利用全年日降水序列建模结果相比,利用汛期日降水序列建模更好地刻画了汛期降水的累计概率分布曲线,同时提高了汛期总降水量、降水强度和年平均暴雨日数、暴雨量、暴雨强度的均值和变化趋势的降尺度效果。模型对较长年份的暴雨重现期订正效果更佳。与当代(2006-2015年)气候相比,2016-2095年上海降水呈现以下特征:全年和汛期总降水量和降水强度增加,降水日数减少,未来可能出现更多的旱涝年;汛期降水极端性增强,暴雨降水均值和极端值均增加;50年以上重现期的年最大日降水量未来呈前40年减少后40年增加的变化。CDF-T模型为站点尺度气候变化影响评估和未来预估提供降尺度技术和基础气候数据。     

非齐次隐马尔可夫降尺度方法对江淮流域夏季逐日降水的模拟及其评估

引入非齐次隐马尔可夫模型(Nonhomogeneous hidden Markov model,NHMM)统计降尺度方法,利用1961—2002年江淮流域夏季逐日降水资料、欧洲中期天气预报中心(ECMWF)的ERA-40再分析资料建立模型,检验其对东部季风区(以江淮流域为代表)夏季日降水的模拟能力,并对比BCC-CSM1.1(m)模式NHMM降尺度前后的模拟效果。结果表明,NHMM降尺度方法通过建立降水概率分布态间转移参数与大尺度环流变量的联系,对江淮流域逐日降水量具有较好的降尺度效果。模拟的各站日降水量概率分布函数(PDF)曲线与观测非常接近,布赖尔评分(Brier Score,S_B)均小于0.11%,显著性评分(Significance Score,Ss)均大于0.84;夏季总降水量、降水日数、中雨日数、降水强度和95%分位降水量指数的多年平均场偏差百分率绝对值低于10%,前3个指数的空间相关系数高于0.9;该方法对各降水指数的年际变率也有一定的模拟能力,模拟得到的各指数的区域平均年际序列与观测序列的相关系数为0.62—0.87。对BCC-CSM1.1(m)模式的模拟结果进行降尺度后,SB较降尺度前平均减小0.57%,Ss平均增大0.23,皆表明降尺度后的概率分布函数曲线更接近于观测;各降水指数在多数台站的偏差百分率绝对值由大于40%降至10%以内,空间相关系数普遍提高至0.8以上。NHMM降尺度方法能够有效提高BCC-CSM1.1(m)模式对江淮流域夏季日降水的模拟能力,相对气候模式具有显著的"增值",未来可进一步利用该方法进行气候变暖背景下的日降水变化预估。     

针对江苏省气温和降水的高分辨率降尺度模拟及预估北大核心

本文利用WRF模式,以25 km分辨率数值模拟结果作为驱动场,对江苏省现代和未来15 a的气候进行5 km高分辨率降尺度模拟及预估。结果显示,高分辨率降尺度模拟对其驱动场具有显著的提高,降水的负偏差和气温冷偏差均有所降低,其模拟的降水与气温概率分布与观测更为接近;对于极端指数,WRF模式能够模拟出其基本分布,除连续湿润日数CWD和极端高温TXx之外,高分辨率模拟对其他指数的模拟均有显著的提升。在RCP8.5排放情景下对未来气候变化的预估表明,江苏降水在夏季以减少为主,在春季则以增加为主,全年平均降水存在减少趋势;未来0~1 mm·d^(-1)的微弱降水发生概率将增加,小雨、中到大雨以及暴雨发生的概率均降低,而暴雨强度的增强导致极端强降水R95显著增加;气温25℃以上高温发生的概率在未来有所增加,而0℃以下的低温发生概率减小,从而导致暖持续日数显著增加,而冷持续日数减小,另外,极端高温和极端低温都有显著的升高。     

全球变暖对区域极端气温出现概率的影响

在适当的概率分布模式基础上，考察冬夏极端气温出现概率对全球气候变暖的敏感率，以上海、南京、东台、徐州站为试验实例。分析结果表明，无论是冬季逐日最低(极端)或平均最低气温或夏季逐日最高(极端)或平均最高气温的概率分布的敏感率几乎都呈现出气温方差变化的影响大于气温均值变化的影响。     

锦州地区逐日降水概率分布特征

利用锦州地区的逐日降水量观测资料对逐日降水量的概率分布进行了统计分析,采用最大似然估计法得到Gamma函数分布的形状参数α和尺度参数β,通过Gamma概率分布模拟观测站点逐日降水的概率分布。结果表明：锦州地区逐日降水频率整体趋势先上升后下降,基本呈对称式分布,降水概率有一定的振荡,个别日会出现远超相邻日期的降水频率,7月21日降水频率最高,在不计微量降水的情况下,最低逐日降水概率有多个日期为0。各季降水频率偏低是造成义县地区干旱的原因之一;北镇夏季平均降水频率最低,但其夏季平均降水量却为锦州地区最高,说明北镇可能易出现较大量级降水或易出现极端降水天气。清明期间降水频率在50%以上、高考期间降水频率在80%以上,符合大众日常对特殊日期降水情况的认知;逐日降水频率可以为公众气象服务提供新的思路。凌海、北镇更容易出现极端降水天气;锦州地区日降水出现小雨天气概率最高,暴雨以上降水概率较低,锦州地区各站极少出现大暴雨以上量级降水,对锦州降水量级预报,尤其是暴雨或大暴雨以上降水量级的预报起到一定的指示作用。     

21世纪前期长江中下游流域极端降水预估及不确定性分析

在全球变暖背景下,极端降水的频率、强度以及持续时间均在显著增加,尤其是对于气候变化敏感的长江中下游流域。由于模式本身、温室气体排放情景以及自然变率存在较大的不确定性,因此未来预估变化的不确定性一直备受关注。为了能够得到对于未来极端降水更为准确的预估结果,使用NEX-GDDP（NASA Earth Exchange Global Daily Downscaled Projections）提供的19个CMIP5降尺度高分辨率数据（0.25°×0.25°）,给出21世纪前期（2016—2035年）长江中下游流域极端降水的可能变化。根据长江中下游流域178个气象站1981—2005年的逐日降水量数据,计算了能够代表极端降水不同特征的指数,在评估模拟能力的基础上给出了21世纪前期RCP4.5情景下极端降水的变化。结果表明,降尺度结果对长江中下游流域极端降水有很好的模拟能力,除R90N外,所有模式模拟其余指数的空间结构与观测的相关系数均超过了0.6。其中所有模式模拟PRCPTOT和R10的相关系数均超过0.95。21世纪前期,长江中下游地区降水趋于极端化,尤其是在流域的西部地区。极端降水日数的变化在减少,表明对于极端降水的贡献主要来自于极端降水日的较大日降水量,而非极端降水日数。未来预估不确定性的大值区主要位于流域的南部地区,流域的西部地区不确定性较低,西部地区极端降水的增加应该受到更多的重视。      

Future changes in precipitation extremes over China using the NEX-GDDP high-resolution daily downscaled data-set

最近,NASA发布了一套基于CMIP5 21个耦合模式输出的高分辨率降尺度逐日数据集,简称NEX-GDDP。本文评估了NEX-GDDP对中国极端降水的模拟性能。研究发现:(1)相比CMIP5直接输出结果,NEX-GDDP能够更好刻画中国极端降水的空间分布;(2)未来中国极端降水事件明显增多、强度增强,NEX-GDDP在区域尺度上给出了更多的气候变化信息;(3)NEXGDDP预估的中国未来极端降水变化的不确定性范围相比CMIP5直接输出结果明显减少,使得预估结果更加可靠.     

Homogeneous clusters over India using probability density function of daily rainfall

The Indian landmass has been divided into homogeneous clusters by applying the cluster analysis to the probability density function of a century-long time series of daily summer monsoon (June through September) rainfall at 357 grids over India, each of approximately 100 km × 100 km. The analysis gives five clusters over Indian landmass; only cluster 5 happened to be the contiguous region and all other clusters are dispersed away which confirms the erratic behavior of daily rainfall over India. The area averaged seasonal rainfall over cluster 5 has a very strong relationship with Indian summer monsoon rainfall; also, the rainfall variability over this region is modulated by the most important mode of climate system, i.e., El Nino Southern Oscillation (ENSO). This cluster could be considered as the representative of the entire Indian landmass to examine monsoon variability. The two-sample Kolmogorov-Smirnov test supports that the cumulative distribution functions of daily rainfall over cluster 5 and India as a whole do not differ significantly. The clustering algorithm is also applied to two time epochs 1901–1975 and 1976–2010 to examine the possible changes in clusters in a recent warming period. The clusters are drastically different in two time periods. They are more dispersed in recent period implying the more erroneous distribution of daily rainfall in recent period.

     

Regional analysis of annual maximum rainfall using TL-moments method

Information related to distributions of rainfall amounts are of great importance for designs of water-related structures. One of the concerns of hydrologists and engineers is the probability distribution for modeling of regional data. In this study, a novel approach to regional frequency analysis using L-moments is revisited. Subsequently, an alternative regional frequency analysis using the TL-moments method is employed. The results from both methods were then compared. The analysis was based on daily annual maximum rainfall data from 40 stations in Selangor Malaysia. TL-moments for the generalized extreme value (GEV) and generalized logistic (GLO) distributions were derived and used to develop the regional frequency analysis procedure. TL-moment ratio diagram and Z-test were employed in determining the best-fit distribution. Comparison between the two approaches showed that the L-moments and TL-moments produced equivalent results. GLO and GEV distributions were identified as the most suitable distributions for representing the statistical properties of extreme rainfall in Selangor. Monte Carlo simulation was used for performance evaluation, and it showed that the method of TL-moments was more efficient for lower quantile estimation compared with the L-moments.     

Downscaling daily precipitation time series using a combined circulation- and regression-based approach

The aim of this paper is to introduce a new conditional statistical model for generating daily precipitation time series. The generated daily precipitation can thus be used for climate change impact studies, e.g., crop production, rainfall–runoff, and other water-related processes. It is a stochastic model that links local rainfall events to a continuous atmospheric predictor, moisture flux, in addition to classified atmospheric circulation patterns. The coupled moisture flux is proved to be capable of capturing continuous property of climate system and providing extra information to determine rainfall probability and rainfall amount. The application was made to simultaneously downscale daily precipitation at multiple sites within the Rhine River basin. The results show that the model can well reproduce statistical properties of daily precipitation time series. Especially for extreme rainfall events, the model is thought to better reflect rainfall variability compared to the pure CP-based downscaling approach.     

Fitting the distribution of dry and wet spells with alternative probability models

The development of the rainfall occurrence model is greatly important not only for data-generation purposes, but also in providing informative resources for future advancements in water-related sectors, such as water resource management and the hydrological and agricultural sectors. Various kinds of probability models had been introduced to a sequence of dry (wet) days by previous researchers in the field. Based on the probability models developed previously, the present study is aimed to propose three types of mixture distributions, namely, the mixture of two log series distributions (LSD), the mixture of the log series Poisson distribution (MLPD), and the mixture of the log series and geometric distributions (MLGD), as the alternative probability models to describe the distribution of dry (wet) spells in daily rainfall events. In order to test the performance of the proposed new models with the other nine existing probability models, 54 data sets which had been published by several authors were reanalyzed in this study. Also, the new data sets of daily observations from the six selected rainfall stations in Peninsular Malaysia for the period 1975–2004 were used. In determining the best fitting distribution to describe the observed distribution of dry (wet) spells, a Chi-square goodness-of-fit test was considered. The results revealed that the new method proposed that MLGD and MLPD showed a better fit as more than half of the data sets successfully fitted the distribution of dry and wet spells. However, the existing models, such as the truncated negative binomial and the modified LSD, were also among the successful probability models to represent the sequence of dry (wet) days in daily rainfall occurrence.     

AN OPERATIONAL STATISTICAL SCHEME FOR TROPICAL CYCLONE INDUCED RAINFALL FORECAST

A non-parametric method is used in this study to analyze and predict short-term rainfall due to tropical cyclones(TCs) in a coastal meteorological station. All 427 TCs during 1953-2011 which made landfall along the Southeast China coast with a distance less than 700 km to a certain meteorological station- Shenzhen are analyzed and grouped according to their landfalling direction, distance and intensity. The corresponding daily rainfall records at Shenzhen Meteorological Station(SMS) during TCs landfalling period(a couple of days before and after TC landfall) are collected. The maximum daily rainfall(R-24) and maximum 3-day accumulative rainfall(R-72) records at SMS for each TC category are analyzed by a non-parametric statistical method, percentile estimation. The results are plotted by statistical boxplots, expressing in probability of precipitation. The performance of the statistical boxplots is evaluated to forecast the short-term rainfall at SMS during the TC seasons in 2012 and 2013. Results show that the boxplot scheme can be used as a valuable reference to predict the short-term rainfall at SMS due to TCs landfalling along the Southeast China coast.     

基于多状态Markov链模式的极端降水模拟试验

文中建立了基于多状态一阶Markov链的逐日降水量随机模式式结合广义帕雷托分布(GPD)产生夏季逐日极端降水量的模拟资料,结果所显示的各种气候特征表明,绝大多数站点(尤其是中国东部多雨地区)都达到较高的精度.分析表明,该模式对中国东部极端降水特征的模拟能力在某些方面优于两状态一阶Markov链模式.对东部6个代表站模拟试验结果表明,月降水均方差、日降水极大值、月半均降水日数、日降水均方差、日平均降水量等指标与实况比较,均证明该模式对逐日降水量的模拟效果较好,基本模拟出降水量的各种特征.对中同东部78个代表站采用的两种模式模拟结果对比发现,除日平均降水量以外,月半均降水日数、日降水平均极大值都与实际观测结果较为一致,总体上优于两状态模式,说明用该模式在全国范围内模拟逐日降水特征尤其是极端降水特征有较高的町行性.例如,由其中6个代表站模拟资料所拟合的极端降水GPD模式具有较高的拟合优度.无论从门限值或重现期值来看都可发现模拟与实测结果有较好的相似性,且两者门限值的误差越小,重现期极值的差距也越小.证明Markov链模式对极端降水的模拟有广泛的适用性.     

Do weak AR4 models bias projections of future climate changes over Australia?

Regional climate projections using climate models commonly use an “all-model” ensemble based on data sets such as the Intergovernmental Panel on Climate Change’s (IPCC) 4th Assessment (AR4). Some regional assessments have omitted models based on specific criteria. We use a criteria based on the capacity of climate models to simulate the observed probability density function calculated using daily data, model-by-model and region-by-region for each of the AR4 models over Australia. We demonstrate that by omitting those climate models with relatively weak skill in simulating the observed probability density functions of maximum and minimum temperature and precipitation, different regional projections are obtained. Differences include: larger increases in the mean maximum and mean minimum temperatures, but smaller increases in the annual maximum and minimum temperatures. There is little impact on mean precipitation but the better models simulate a larger increase in the annual rainfall event combined with a larger decrease in the number of rain days. The weaker models bias the amount of mean warming towards lower increases, bias annual maximum temperatures to excessive warming and bias precipitation such that the amount of the annual rainfall event is under-estimated. We suggest that omitting weak models from regional scale estimates of future climate change helps clarify the nature and scale of the projected impacts of global warming.     

Increasing Flash Floods in a Drying Climate over Southwest China

In a globally warming world,subtropical regions are generally expected to become drier while the tropics and mid–high latitudes become wetter.In line with this,Southwest China,close to 25?N,is expected to become increasingly prone to drought if annual mean precipitation decreases.However,despite this trend,changes in the temporal distribution of moisture supply might actually result in increased extreme rainfall in the region,whose climate is characterized by distinct dry and wet seasons.Using hourly and daily gauge observations,rainfall intensity changes since 1971 are examined for a network of 142 locations in the region.From the analysis,dry season changes are negligible but wet season changes exhibit a significantly strong downward trend [-2.4%(10 yr)~(-1)],particularly during the past 15 years [-17.7%(10 yr)~(-1)].However,the intensity of events during the wettest of 5% hours appears to steadily increase during the whole period [1.4%(10 yr)~(-1)],tying in with government statistical reports of recent droughts and flooding.If the opposing trends are a consequence of a warming climate,it is reasonable to expect the contradictory trend to continue with an enhanced risk of flash flooding in coming decades in the region concerned.     

Projection and possible causes of summer precipitation in eastern China using self-organizing map

Self-organizing map (SOM) is used to simulate summer daily precipitation over the Yangtze–Huaihe river basin in Eastern China, including future projections. SOM shows good behaviors in terms of probability distribution of daily rainfall and spatial distribution of rainfall indices, as well as consistency of multi-model simulations. Under RCP4.5 Scenario, daily rainfall at most sites (63%) is projected to shift towards larger values. For the early 21st century (2016–2035), precipitation in the central basin increases, yet decreases occur over the middle reaches of the Yangtze River as well as a part of its southeast area. For the late 21st century (2081–2100), the mean precipitation and extreme indices experience an overall increase except for a few southeast stations. The total precipitation in the lower reaches of the Yangtze River and in its south area is projected to increase from 7% at 1.5 °C global warming to 11% at 2 °C, while the intensity enhancement is more significant in southern and western sites of the domain. A clustering allows to regroup all SOM nodes into four distinct regimes. Such regional synoptic regimes show remarkable stability for future climate. The overall intensification of precipitation in future climate is linked to the occurrence-frequency rise of a wet regime which brings longitudinally closer the South Asia High (eastward extended) and the Western Pacific Subtropical High (westward extended), as well as the reduction of a dry pattern which makes the two atmospheric centers of action move away from each other.     

南四湖流域暴雨分布特征及可能日最大降水量计算

利用南四湖流域11县市1971～2007年的暴雨资料,分析南四湖流域首次和末次暴雨的开始和结束时间以及暴雨的时空分布特征,发现南四湖流域暴雨的时空分布差异较大,但日降水极值的概率分布却有一定规律,呈Λ(x)型渐进分布。利用耿贝尔分布计算南四湖流域多年一遇的日最大降水量极值,计算的未来10年、20年、40年的日最大降水量与历史上10年、20年、40年的日最大降水量重现期基本一致,对未来50～200年的估算值也具有一定的参考价值。