Application of physical scaling towards downscaling climate model precipitation data

Physical scaling (SP) method downscales climate model data to local or regional scales taking into consideration physical characteristics of the area under analysis. In this study, multiple SP method based models are tested for their effectiveness towards downscaling North American regional reanalysis (NARR) daily precipitation data. Model performance is compared with two state-of-the-art downscaling methods: statistical downscaling model (SDSM) and generalized linear modeling (GLM). The downscaled precipitation is evaluated with reference to recorded precipitation at 57 gauging stations located within the study region. The spatial and temporal robustness of the downscaling methods is evaluated using seven precipitation based indices. Results indicate that SP method-based models perform best in downscaling precipitation followed by GLM, followed by the SDSM model. Best performing models are thereafter used to downscale future precipitations made by three global circulation models (GCMs) following two emission scenarios: representative concentration pathway (RCP) 2.6 and RCP 8.5 over the twenty-first century. The downscaled future precipitation projections indicate an increase in mean and maximum precipitation intensity as well as a decrease in the total number of dry days. Further an increase in the frequency of short (1-day), moderately long (2–4 day), and long (more than 5-day) precipitation events is projected.     

Analysis and mapping of present and future drought conditions over Greek areas with different climate conditions

Estimation of drought in a certain temporal and spatial scale is crucial in climate change studies. The current study targets on three agricultural areas widespread in Greece, Ardas River Basin in Northeastern Greece, Sperchios River Basin in Central Greece, and Geropotamos River Basin in Crete Island in South Greece that are characterized by diverse climates as they are located in various regions. The objective is to assess the spatiotemporal variation of drought conditions prevailing in these areas. The Standardized Precipitation Index (SPI) was used to identify and assess the present and future drought conditions. Future simulated data were derived from a number of Regional Climatic Models (RCMs) from the ENSEMBLES European Project. The analysis was performed for the future periods of 2021–2050 and 2071–2100, implementing A1B and B1 scenarios. The spatial analysis of the drought conditions was performed using a combined downscaling technique and the Ordinary Kriging. The Mann-Kendall test was implemented for trend investigation. During both periods and scenarios, drought conditions will tend to be more severe in the upcoming years. The decrease of the SPI values in the Sperchios River Basin is expected to be the strongest, as it is the only study area that will show a negative balance (in SPI values), regarding the drought conditions. For the Ardas and the Geropotamos River Basins, a great increase of the drought conditions will occur during the 2021–2050 period, while for 2071–2100 period, the decrease will continue but it will be tempered. Nevertheless, the situation in all study areas according to the SPI classification is characterized as “Near-normal”, in terms of drought conditions.     

Tibetan Plateau summer precipitation: covariability with circulation indices

Relations between Tibetan Plateau precipitation and large-scale climate indices are studied based on the Standardized Precipitation Index (SPI) and the boreal summer season. The focus is on the decadal variability of links between the large-scale circulation and the plateau drought and wetness. Analysis of teleconnectivity of the continental northern hemisphere standardized summer precipitation reveals the Tibetan Plateau as a major SPI teleconnectivity center in south-eastern Asia connecting remote correlation patterns over Eurasia. Employing a moving window approach, changes in covariability and synchronizations between Tibetan Plateau summer SPI and climate indices are analyzed on decadal time scales. Decadal variability in the relationships between Tibetan Plateau summer SPI and the large-scale climate system is characterized by three shifts related to changes in the North Atlantic, the Indian Ocean, and the tropical Pacific. Changes in the North Atlantic variability (North Atlantic Oscillation) result in a stable level of Tibetan Plateau summer SPI variability; the response to changes in tropical Pacific variability is prominent in various indices such as Asian monsoon, Pacific/North America, and East Atlantic/Western Russia pattern.

     

An Experiment of a Statistical Downscaling Forecast Model for Summer Precipitation over China

A combination of the optimal subset regression (OSR) approach, the coupled general circulation model of the National Climate Center (NCC-CGCM) and precipitation observations from 160 stations over China is used to construct a statistical downscaling forecast model for precipitation in summer. Retroactive forecasts are performed to assess the skill of statistical downscaling during the period from 2003 to 2009. The results show a poor simulation for summer precipitation by the NCC- CGCM for China, and the average spatial anomaly correlation coefficient (ACC) is 0.01 in the forecast period. The forecast skill can be improved by OSR statistical downscaling, and the OSR forecast performs better than the NCC-CGCM in most years except 2003. The spatial ACC is more than 0.2 in the years 2008 and 2009, which proves to be relatively skillful. Moreover, the statistical downscaling forecast performs relatively well for the main rain belt of the summer precipitation in some years, including 2005, 2006, 2008, and 2009. However, the forecast skill of statistical downscaling is restricted to some extent by the relatively low skill of the NCC- CGCM.     

Robust intensification of hydroclimatic intensity over East Asia from multi-model ensemble regional projections

This study assesses the hydroclimatic response to global warming over East Asia from multi-model ensemble regional projections. Four different regional climate models (RCMs), namely, WRF, HadGEM3-RA, RegCM4, and GRIMs, are used for dynamical downscaling of the Hadley Centre Global Environmental Model version 2–Atmosphere and Ocean (HadGEM2-AO) global projections forced by the representative concentration pathway (RCP4.5 and RCP8.5) scenarios. Annual mean precipitation, hydroclimatic intensity index (HY-INT), and wet and dry extreme indices are analyzed to identify the robust behavior of hydroclimatic change in response to enhanced emission scenarios using high-resolution (12.5 km) and long-term (1981–2100) daily precipitation. Ensemble projections exhibit increased hydroclimatic intensity across the entire domain and under both the RCP scenarios. However, a geographical pattern with predominantly intensified HY-INT does not fully emerge in the mean precipitation change because HY-INT is tied to the changes in the precipitation characteristics rather than to those in the precipitation amount. All projections show an enhancement of high intensity precipitation and a reduction of weak intensity precipitation, which lead to a possible shift in hydroclimatic regime prone to an increase of both wet and dry extremes. In general, projections forced by the RCP8.5 scenario tend to produce a much stronger response than do those by the RCP4.5 scenario. However, the temperature increase under the RCP4.5 scenario is sufficiently large to induce significant changes in hydroclimatic intensity, despite the relatively uncertain change in mean precipitation. Likewise, the forced responses of HY-INT and the two extreme indices are more robust than that of mean precipitation, in terms of the statistical significance and model agreement.

     

Comparison of regional climate model and statistical downscaling simulations of different winter precipitation change scenarios over Romania

Summary Regional climate model and statistical downscaling procedures are used to generate winter precipitation changes over Romania for the period 2071–2100 (compared to 1961–1990), under the IPCC A2 and B2 emission scenarios. For this purpose, the ICTP regional climate model RegCM is nested within the Hadley Centre global atmospheric model HadAM3H. The statistical downscaling method is based on the use of canonical correlation analysis (CCA) to construct climate change scenarios for winter precipitation over Romania from two predictors, sea level pressure and specific humidity (either used individually or together). A technique to select the most skillful model separately for each station is proposed to optimise the statistical downscaling signal. Climate fields from the A2 and B2 scenario simulations with the HadAM3H and RegCM models are used as input to the statistical downscaling model. First, the capability of the climate models to reproduce the observed link between winter precipitation over Romania and atmospheric circulation at the European scale is analysed, showing that the RegCM is more accurate than HadAM3H in the simulation of Romanian precipitation variability and its connection with large-scale circulations. Both models overestimate winter precipitation in the eastern regions of Romania due to an overestimation of the intensity and frequency of cyclonic systems over Europe. Climate changes derived directly from the RegCM and HadAM3H show an increase of precipitation during the 2071–2100 period compared to 1961–1990, especially over northwest and northeast Romania. Similar climate change patterns are obtained through the statistical downscaling method when the technique of optimum model selected separately for each station is used. This adds confidence to the simulated climate change signal over this region. The uncertainty of results is higher for the eastern and southeastern regions of Romania due to the lower HadAM3H and RegCM performance in simulating winter precipitation variability there as well as the reduced skill of the statistical downscaling model.     

Simulation of extreme precipitation indices in the Yangtze River basin by using statistical downscaling method (SDSM)

In this study, the applicability of the statistical downscaling model (SDSM) in modeling five extreme precipitation indices including R10 (no. of days with precipitation ≥10?mm?day^?1), SDI (simple daily intensity), CDD (maximum number of consecutive dry days), R1d (maximum 1-day precipitation total) and R5d (maximum 5-day precipitation total) in the Yangtze River basin, China was investigated. The investigation mainly includes the calibration and validation of SDSM model on downscaling daily precipitation, the validation of modeling extreme precipitation indices using independent period of the NCEP reanalysis data, and the projection of future regional scenarios of extreme precipitation indices. The results showed that: (1) there existed good relationship between the observed and simulated extreme precipitation indices during validation period of 1991–2000, the amount and the change pattern of extreme precipitation indices could be reasonably simulated by SDSM. (2) Under both scenarios A2 and B2, during the projection period of 2010–2099, the changes of annual mean extreme precipitation indices in the Yangtze River basin would be not obvious in 2020s; while slightly increase in the 2050s; and significant increase in the 2080s as compared to the mean values of the base period. The summer might be the more distinct season with more projected increase of each extreme precipitation indices than in other seasons. And (3) there would be distinctive spatial distribution differences for the change of annual mean extreme precipitation indices in the river basin, but the most of Yangtze River basin would be dominated by the increasing trend.     

Validating a regional climate model’s downscaling ability for East Asian summer monsoonal interannual variability

Performance of a regional climate model (RCM), WRF, for downscaling East Asian summer season climate is investigated based on 11-summer integrations associated with different climate conditions with reanalysis data as the lateral boundary conditions. It is found that while the RCM is essentially unable to improve large-scale circulation patterns in the upper troposphere for most years, it is able to simulate better lower-level meridional moisture transport in the East Asian summer monsoon. For precipitation downscaling, the RCM produces more realistic magnitude of the interannual variation in most areas of East Asia than that in the reanalysis. Furthermore, the RCM significantly improves the spatial pattern of summer rainfall over dry inland areas and mountainous areas, such as Mongolia and the Tibetan Plateau. Meanwhile, it reduces the wet bias over southeast China. Over Mongolia, however, the performance of precipitation downscaling strongly depends on the year: the WRF is skillful for normal and wet years, but not for dry years, which suggests that land surface processes play an important role in downscaling ability. Over the dry area of North China, the WRF shows the worst performance. Additional sensitivity experiments testing land effects in downscaling suggest the initial soil moisture condition and representation of land surface processes with different schemes are sources of uncertainty for precipitation downscaling. Correction of initial soil moisture using the climatology dataset from GSWP-2 is a useful approach to robustly reducing wet bias in inland areas as well as to improve spatial distribution of precipitation. Despite the improvement on RCM downscaling, regional analyses reveal that accurate simulation of precipitation over East China, where the precipitation pattern is strongly influenced by the activity of the Meiyu/Baiu rainfall band, is difficult. Since the location of the rainfall band is closely associated with both lower-level meridional moisture transport and upper-level circulation structures, it is necessary to have realistic upper-air circulation patterns in the RCM as well as lower-level moisture transport in order to improve the circulation-associated convective rainfall band in East Asia.     

基于SPI的1961－2020年昌吉地区作物生长季气象干旱时空特征研究

基于昌吉地区7个气象站1961－2020年降水量资料，计算昌吉地区作物生长季标准化降水指数（SPI-7）。运用趋势分析法、M-K突变检验法和小波分析法探究了昌吉地区作物生长季SPI-7指数的年际和年代变化特征；在此基础上分析了作物生长季干旱的站次比和干旱强度的年际变化,并结合该区实际发生的旱灾对SPI进行了验证。结果表明：1961—2020年昌吉地区作物生长季标准化降水指数以0.08/10 a的速率呈微弱的正趋势（变湿），在年代变化趋势中呈现出变干-变湿-变干的变化波动, 1981年标准化降水指数由低到高突变；干旱强度呈增加趋势，干旱发生的区域面积有轻微减少的趋势；干旱强度在全区范围内主要为轻旱和中旱等级，并表现为全域性干旱和区域性干旱；空间分布上看干旱率最高区域在东部地区，轻旱主要集中在东部，中旱、重旱和特旱集中在西部地区，干旱强度大的区域大致分布在西部地区；在周期性变化方面，SPI指数存在着6年、9年、16年周期震荡；历史旱灾与SPI指数干旱评价结果吻合率较高，SPI指数在昌吉地区作物生长季的干旱监测与分析中具有较好的实用性。     

基于标准降水指标的新疆干旱特征演变

运用标准降水指标 (Standardized Precipitation Index, SPI) 对新疆地区53个雨量站1957—2009年日降水量资料进行全面分析,研究了不同干旱等级发生概率的空间分布变化规律。同时,采用Mann-Kendall趋势检验法分析了各月份标准降水指标值、干旱强度和干旱历时的变化趋势,探讨了新疆地区干旱时空演变特征。结果表明:北疆易发生中等及以上干旱,而南疆易发生轻度干旱。总的来说,北疆干旱强度有下降的趋势,干旱历时趋于缩短,南疆南部干旱强度和干旱历时有轻微上升,东疆中部干旱情况显著恶化。具体来说,北疆冬季干旱程度有减弱趋势,而对于农业生产较重要的春、夏、秋季,西部干旱加剧;南疆夏季干旱有减弱趋势;东疆中部四季干旱程度有轻微加剧趋势。     

Drought in the Southern United States over the 20th century: variability and its impacts on terrestrial ecosystem productivity and carbon storage

Drought is one of the most devastating natural hazards faced by the Southern United States (SUS). Drought events and their adverse impacts on the economy, society and environment have been extensively reported during 1895?C2007. Our aim is thus to characterize drought conditions in the SUS and explore the impacts on terrestrial ecosystem function (i.e., net primary productivity (NPP) and net carbon exchange (NCE)). Standard precipitation index (SPI) was used to characterize drought intensity and duration, and a process-based ecosystem model was used to explore the relationship between drought and ecosystem function. Combining overall information on growing-season SPI, drought area and duration, we concluded there was no significant change in drought conditions for the SUS during 1895?C2007. However, increased drought intensity was found for many areas in the east, resulting in significant decreases in NPP for these areas, with the largest decrease up to 40% during extreme droughts. Changes in precipitation patterns increased C emissions of 0.16 Pg (1 Pg?=?10¹⁵?g) in the SUS during 1895?C2007. The west (dry region) acted as a C sink due to increased precipitation, while the east (water-rich region) acted as a C source due to increased drought intensity. Both NPP and NCE significantly increased along a gradient of declining drought intensity. Changes in precipitation resulted in C sources in forest, wetland, and cropland ecosystems, while C sinks in shrubland and grassland ecosystems. Changes in air temperature could either enhance or reduce drought impacts on NPP and NCE across different vegetation types.     

Interannual summer air temperature variability over Greece and its connection to the large-scale atmospheric circulation and Mediterranean SSTs 1950–1999

The interannual and decadal variability of summer (June to September) air temperature in the northeastern Mediterranean is analysed for the period 1950 to 1999. Extremely hot and cool summers are illustrated by means of composite analysis. The combined influence of the large-scale atmospheric circulation and thermic predictors on local temperature is assessed by means of an objective approach based on empirical orthogonal functions and canonical correlation analysis. Monthly values of sea level pressure, geopotential heights, atmospheric thickness and Mediterranean sea surface temperatures are used as predictor fields and air temperature from 24 observational sites spread over Greece and western Turkey constitute the predictand variable. Results indicate that more than 50% of the total summer temperature variability can be explained linearly by the combination of eight large-scale predictor fields on two canonical correlation modes. The first canonical mode is related to a more meridional circulation at the upper tropospheric levels, which favours local land-sea contrasts in the associated local temperature pattern. Variations of this mode are found to be responsible for the occurrence of extreme events and decadal trends in regional temperature, the latter being characterized by a cooling in the early 1960s and a warming in the early 1990s. The second canonical mode pictures variations in the intensity of the zonal circulation over the Atlantic area that drive temperature anomalies affecting mainly the Aegean Sea and the west of Greece. Our results suggest the potential of statistical downscaling for Greek summer temperature with reliable climate forecasts for planetary-scale anomalies.     

An Improved Statistical-Dynamical Downscaling Scheme and its Application to the Alpine Precipitation Climatology

Summary An improved statistical-dynamical downscaling method for the regionalization of large-scale climate analyses or simulations is introduced. The method is based on the disaggregation of a multi-year time-series of large-scale meteorological data into multi-day episodes of quasi-stationary circulation. The episodes are subsequently grouped into a defined number of classes. A regional model is used to simulate the evolution of weather during the most typical episode of each class. These simulations consider the effects of the regional topography. Finally, the regional model results are statistically weighted with the climatological frequencies of the respective circulation classes in order to provide regional climate patterns. The statistical-dynamical downscaling procedure is applied to large-scale analyses for a 12-year climate period 1981–1992. The performance of the new method is demonstrated for winter precipitation in the Alpine region. With the help of daily precipitation analyses it was possible to validate the results and to assess the different sources of errors. It appeared that the main error originates from the regional model, whereas the error of the procedure itself was relatively unimportant. This new statistical-dynamical downscaling method turned out to be an efficient alternative to the commonly used method of nesting a regional model continuously within a general circulation model (dynamical downscaling). Received April 8, 1999 Revised July 30, 1999     

Variability and trends in daily precipitation extremes on the northern and southern slopes of the central Himalaya

This study focuses on the precipitation extremes recorded on the northern and southern slopes of the central Himalaya, especially those documented at higher altitudes. Daily precipitation data recorded over a 35-year period at nine meteorological stations in the region were studied. We used the precipitation extreme indices delineated by the Expert Team on Climate Change Detection and Indices (ETCCDI). The spatial and temporal variations in these precipitation extremes were calculated. When regional patterns were investigated to detect any anomalies, only 1 of the 10 precipitation extreme indices from the southern slopes of the central Himalaya showed a statistically significant trend; none from the northern slopes of the central Himalaya highlighted a statistically significant trend. On the southern slopes, all indices increased, apart from the maximum 1-day precipitation (RX1) and simple daily precipitation intensity (SDII) indices. Indices such as the consecutive dry days (CDDs) and RX1 indices exhibited similar increases on both the northern and southern slopes of the central Himalaya. These results suggest that increases in precipitation have been accompanied by an increasing frequency of extremes over the southern central Himalaya. Nonetheless, no relation could be established between the precipitation extreme indices and circulation indices for higher altitudes.

     

Future Changes of Drought and Flood Events in China under a Global Warming Scenario

This study investigates the impact of global warming on drought/flood patterns in China at the end of the 21st century based on the simulations of 22 global climate models and a regional climate model(RegCM3) under the SRES(Special Report on Emissions Scenarios) A1B scenario.The standardized precipitation index(SPI),which has well performance in monitoring the drought/flood characteristics(in terms of their intensity,duration,and spatial extent) in China,is used in this study.The projected results of 22 coupled models and the RegCM3 simulation are consistent.These models project a decrease in the frequency of droughts in most parts of northern China and a slight increase in the frequency in some parts of southern China.Considering China as a whole,the spatial extents of droughts are projected to be significantly reduced.In contrast,future flood events over most parts of China are projected to occur more frequently with stronger intensity and longer duration than those prevalent currently.Additionally,the spatial extents of flood events are projected to significantly increase.     

Impact of drought on wildland fires in Greece: implications of climatic change?

An increasing trend and a statistically significant positive correlation between wildfire occurrence, area burned and drought (as expressed by the Standardized Precipitation Index, SPI) have been observed all over Greece, during the period 1961?C1997. In the more humid and colder regions (Northern and Western Greece) the number of fires and area burned were positively correlated to both summer (SPI6_October) and annual drought (SPI12_September), whereas in the relatively more dry and hot regions (Southern and Central Greece) the number of fires and area burned were correlated only to summer drought. In 1978, Greece entered a period of prolonged drought, possibly as a result of the global climatic change. Data analysis of the period 1978?C1997 revealed a statistically significant increase in the mean annual number of fires, the area burned and the summer and annual drought episodes in the relatively more humid and colder regions (Northern and Western) of Greece (which in the past were characterized by less fires and area burned) compared to the more dry and hot regions (Southern and Eastern Greece), which always presented high fire activity. Additionally, analyzing the two sub-periods (1961?C1977, 1978?C1997) separately, drought was significantly correlated only to fire occurrence during the years 1961?C1977, whereas during 1978?C1997 drought was significantly correlated mainly to area burned. It became obvious that drought episodes, although they are not solely responsible for fire occurrence and area burned, they exert an increasingly significant impact on wildfire activity in Greece.     

Regional climate model projections for the State of Washington

Global climate models do not have sufficient spatial resolution to represent the atmospheric and land surface processes that determine the unique regional climate of the State of Washington. Regional climate models explicitly simulate the interactions between the large-scale weather patterns simulated by a global model and the local terrain. We have performed two 100-year regional climate simulations using the Weather Research and Forecasting (WRF) model developed at the National Center for Atmospheric Research (NCAR). One simulation is forced by the NCAR Community Climate System Model version 3 (CCSM3) and the second is forced by a simulation of the Max Plank Institute, Hamburg, global model (ECHAM5). The mesoscale simulations produce regional changes in snow cover, cloudiness, and circulation patterns associated with interactions between the large-scale climate change and the regional topography and land-water contrasts. These changes substantially alter the temperature and precipitation trends over the region relative to the global model result or statistical downscaling. To illustrate this effect, we analyze the changes from the current climate (1970–1999) to the mid twenty-first century (2030–2059). Changes in seasonal-mean temperature, precipitation, and snowpack are presented. Several climatological indices of extreme daily weather are also presented: precipitation intensity, fraction of precipitation occurring in extreme daily events, heat wave frequency, growing season length, and frequency of warm nights. Despite somewhat different changes in seasonal precipitation and temperature from the two regional simulations, consistent results for changes in snowpack and extreme precipitation are found in both simulations.     

利用条件植被指数评价西藏植被对气象干旱的响应

基于2000—2014年4—10月西藏气象站遥感干旱指数 (条件植被指数,VCI) 和气象干旱指数 (标准降水指数,SPI) 之间的相关性,评估植被对气象干旱的响应特征,通过分析气候环境要素对响应特征的影响并归纳相应规则,获取西藏地区植被对气象干旱有明显响应的区域分布。结果显示:VCI与12周时间尺度的SPI具有较强相关性,说明西藏地区植被生长对降水的响应大约滞后12周;植被对气象干旱响应不敏感的原因主要包括气候极度干燥或极度湿润、土地覆盖类型为森林、年平均归一化植被指数 (NDVI) 值过小、多年NDVI变化标准差过小、有降水之外的其他水源补给等;基于对区域气候环境要素特征的分析,可以得出西藏中部偏南地区植被对气象干旱有明显响应,主要包括拉萨地区、山南地区北部、日喀则地区东部、那曲地区中部和西南部、阿里地区的东南部。     

Uncertainty analysis of statistical downscaling models using Hadley Centre Coupled Model

Regression-based statistical downscaling is a method broadly used to resolve the coarse spatial resolution of general circulation models. Nevertheless, the assessment of uncertainties linked with climatic variables is essential to climate impact studies. This study presents a procedure to characterize the uncertainty in regression-based statistical downscaling of daily precipitation and temperature over a highly vulnerable area (semiarid catchment) in the west of Iran, based on two downscaling models: a statistical downscaling model (SDSM) and an artificial neural network (ANN) model. Biases in mean, variance, and wet/dry spells are estimated for downscaled data using vigorous statistical tests for 30 years of observed and downscaled daily precipitation and temperature data taken from the National Center for Environmental Prediction reanalysis predictors for the years of 1961 to 1990. In the case of daily temperature, uncertainty is estimated by comparing monthly mean and variance of downscaled and observed daily data at a 95 % confidence level. In daily precipitation, downscaling uncertainties were evaluated from comparing monthly mean dry and wet spell lengths and their confidence intervals, cumulative frequency distributions of monthly mean of daily precipitation, and the distributions of monthly wet and dry days for observed and modeled daily precipitation. Results showed that uncertainty in downscaled precipitation is high, but simulation of daily temperature can reproduce extreme events accurately. Finally, this study shows that the SDSM is the most proficient model at reproducing various statistical characteristics of observed data at a 95 % confidence level, while the ANN model is the least capable in this respect. This study attempts to test uncertainties of regression-based statistical downscaling techniques in a semiarid area and therefore contributes to an improvement of the quality of predictions of climate change impact assessment in regions of this type.     

Quantile regression and clustering analysis of standardized precipitation index in the Tarim River Basin,Xinjiang, China

Because drought is a very common and widespread natural disaster, it has attracted a great deal of academic interest. Based on 12-month time scale standardized precipitation indices (SPI12) calculated from precipitation data recorded between 1960 and 2015 at 22 weather stations in the Tarim River Basin (TRB), this study aims to identify the trends of SPI and drought duration, severity, and frequency at various quantiles and to perform cluster analysis of drought events in the TRB. The results indicated that (1) both precipitation and temperature at most stations in the TRB exhibited significant positive trends during 1960–2015; (2) multiple scales of SPIs changed significantly around 1986; (3) based on quantile regression analysis of temporal drought changes, the positive SPI slopes indicated less severe and less frequent droughts at lower quantiles, but clear variation was detected in the drought frequency; and (4) significantly different trends were found in drought frequency probably between severe droughts and drought frequency.