Non-stationary future return levels for extreme rainfall over Extremadura (southwestern Iberian Peninsula)

Based on a previous study for temperature, a new method for the calculation of non-stationary return levels for extreme rainfall is described and applied to Extremadura, a region of southwestern Spain, using the peaks-over-threshold approach. Both all-days and rainy-days-only datasets were considered and the 20-year return levels expected in 2020 were estimated taking different trends into account: first, for all days, considering a time-dependent threshold and the trend in the scale parameter of the generalized Pareto distribution; and second, for rainy days only, considering how the mean, variance, and number of rainy days evolve. Generally, the changes in mean, variance and number of rainy days can explain the observed trends in extremes, and their extrapolation gives more robust estimations. The results point to a decrease of future return levels in 2020 for spring and winter, but an increase for autumn.     

A statistical model for karst spring discharge estimation under extensive groundwater development and extreme climate change

ABSTRACT

To acquire better understanding of spring discharge under extreme climate change and extensive groundwater pumping, this study proposed an extreme value statistical decomposition model, in which the spring discharge was decomposed into three items: a long-term trend; periodic variation; and random fluctuation. The long-term trend was fitted by an exponential function, and the periodic variation was fitted by an exponential function whose index was the sum of two sine functions. A general extreme value (GEV) model was used to obtain the return level of extreme random fluctuation. Parameters of the non-linear long-term trend and periodic variation were estimated by the Levenberg-Marquardt algorithm, and the GEV model was estimated by the maximum likelihood method. The extreme value statistical decomposition model was applied to Niangziguan Springs, China to forecast spring discharge. We showed that the modelled spring discharge fitted the observed data very well. Niangziguan Springs discharge is likely to continue declining with fluctuation, and the risk of cessation by August 2046 is 1%. The extreme value decomposition model is a robust method for analysing the nonstationary karst spring discharge under conditions of extensive groundwater development/pumping, and extreme climate changes.

Editor D. Koutsoyiannis; Associate editor J. Ward     

A multiscaling‐based intensity–duration–frequency model for extreme precipitation

Rainfall intensity–duration–frequency (IDF) curves are a standard tool in urban water resources engineering and management. They express how return levels of extreme rainfall intensity vary with duration. The simple scaling property of extreme rainfall intensity, with respect to duration, determines the form of IDF relationships. It is supposed that the annual maximum intensity follows the generalized extreme value (GEV) distribution. As well known, for simple scaling processes, the location parameter and scale parameter of the GEV distribution obey a power law with the same exponent. Although, the simple scaling hypothesis is commonly used as a suitable working assumption, the multiscaling approach provides a more general framework. We present a new IDF relationship that has been formulated on the basis of the multiscaling property. It turns out that the GEV parameters (location and scale) have a different scaling exponent. Next, we apply a Bayesian framework to estimate the multiscaling GEV model and to choose the most appropriate model. It is shown that the model performance increases when using the multiscaling approach. The new model for IDF curves reproduces the data very well and has a reasonable degree of complexity without overfitting on the data.     

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

A hierarchical Bayesian approach for the analysis of climate change impact on runoff extremes

In this study the impact of climate change on runoff extremes is investigated over the Pacific Northwest (PNW). This paper aims to address the question of how the runoff extremes change in the future compared to the historical time period, investigate the different behaviors of the regional climate models (RCMs) regarding the runoff extremes and assess the seasonal variations of runoff extremes. Hydrologic modeling is performed by the variable infiltration capacity (VIC) model at a 1/8° resolution and the model is driven by climate scenarios provided by the North American Regional Climate Change Assessment Program (NARCCAP) including nine regional climate model (RCM) simulations. Analysis is performed for both the historical (1971–2000) and future (2041–2070) time periods. Downscaling of the climate variables including precipitation, maximum and minimum temperature and wind speed is done using the quantile‐mapping (QM) approach. A spatial hierarchical Bayesian model is then developed to analyse the annual maximum runoff in different seasons for both historical and future time periods. The estimated spatial changes in extreme runoffs over the future period vary depending on the RCM driving the hydrologic model. The hierarchical Bayesian model characterizes the spatial variations in the marginal distributions of the General Extreme Value (GEV) parameters and the corresponding 100‐year return level runoffs. Results show an increase in the 100‐year return level runoffs for most regions in particular over the high elevation areas during winter. The Canadian portions of the study region reflect higher increases during spring. However, reduction of extreme events in several regions is projected during summer. Copyright © 2013 John Wiley & Sons, Ltd.     

A flexible nonlinear modelling framework for nonstationary generalized extreme value analysis in hydroclimatology

Parameters in a generalized extreme value (GEV) distribution are specified as a function of covariates using a conditional density network (CDN), which is a probabilistic extension of the multilayer perceptron neural network. If the covariate is time or is dependent on time, then the GEV‐CDN model can be used to perform nonlinear, nonstationary GEV analysis of hydrological or climatological time series. Owing to the flexibility of the neural network architecture, the model is capable of representing a wide range of nonstationary relationships. Model parameters are estimated by generalized maximum likelihood, an approach that is tailored to the estimation of GEV parameters from geophysical time series. Model complexity is identified using the Bayesian information criterion and the Akaike information criterion with small sample size correction. Monte Carlo simulations are used to validate GEV‐CDN performance on four simple synthetic problems. The model is then demonstrated on precipitation data from southern California, a series that exhibits nonstationarity due to interannual/interdecadal climatic variability. Copyright © 2009 Her Majesty the Queen in right of Canada. Published by John Wiley & Sons, Ltd.     

Climate change impacts on extreme precipitation of water supply area in Istanbul: use of ensemble climate modelling and geo-statistical downscaling

ABSTRACT

Numerous statistical downscaling models have been applied to impact studies, but none clearly recommended the most appropriate one for a particular application. This study uses the geographically weighted regression (GWR) method, based on local implications from physical geographical variables, to downscale climate change impacts to a small-scale catchment. The ensembles of daily precipitation time series from 15 different regional climate models (RCMs) driven by five different general circulation models (GCMs), obtained through the European Union (EU)-ENSEMBLES project for reference (1960–1990) and future (2071–2100) scenarios are generated for the Omerli catchment, in the east of Istanbul city, Turkey, under scenario A1B climate change projections. Special focus is given to changes in extreme precipitation, since such information is needed to assess the changes in the frequency and intensity of flooding for future climate. The mean daily precipitation from all RCMs is under-represented in the summer, autumn and early winter, but it is overestimated in late winter and spring. The results point to an increase in extreme precipitation in winter, spring and summer, and a decrease in autumn in the future, compared to the current period. The GWR method provides significant modifications (up to 35%) to these changes and agrees on the direction of change from RCMs. The GWR method improves the representation of mean and extreme precipitation compared to RCM outputs and this is more significant, particularly for extreme cases of each season. The return period of extreme events decreases in the future, resulting in higher precipitation depths for a given return period from most of the RCMs. This feature is more significant with downscaling. According to the analysis presented, a new adaption for regulating excessive water under climate change in the Omerli basin may be recommended.     

Nonstationary extreme value analysis of temperature extremes in China

In a changing climate, the common assumption of stationarity of climate extremes has been increasingly challenged, raising the need to incorporate non-stationarity in extreme value modeling. In this study, quantile regression is used to identify the trends of annual temperature extremes and their correlations with two large climate patterns, the western Pacific subtropical high (WPSH) and the Arctic Oscillation (AO) at 357 stations in China. Statistical significant positive trends and correlations between warm (or cold) temperature extremes and WPSH (or AO) have been detected at most stations. The influence of WPSH on warm extremes is significant in southern China, while the AO mainly affects the cold extremes in northern and eastern China. Then, annual temperature extremes are fitted to generalized extreme value (GEV) distributions with time-varying parameters. The summer (or winter) mean daily maximum (or minimum) temperatures and two climate indices, the WPSH index and the AO index, are chosen as covariates. In total, 16 candidate GEV distribution models are constructed, and the best fitting model with the lowest Bayesian information criterion (BIC) is selected. The 20-year return levels of annual warm (or cold) extremes in the period 1961–1980 and 1991–2010 are computed and compared. The changes of 20-year return levels of annual warm and cold extremes are jointly determined by trend and distributional changes of annual temperature extremes. Analysis of large scale atmospheric circulation changes indicate that a strengthening anticyclonic circulation and increasing geopotential height in recent decades may have contributed to the changes in temperature extremes in China.     

Rainfall frequency analysis using a mixed GEV distribution: a case study for annual maximum rainfalls in South Korea

Extreme rainfalls in South Korea result mainly from convective storms and typhoon storms during the summer. A proper way for dealing with the extreme rainfalls in hydrologic design is to consider the statistical characteristics of the annual maximum rainfall from two different storms when determining design rainfalls. Therefore, this study introduced a mixed generalized extreme value (GEV) distribution to estimate the rainfall quantile for 57 gauge stations across South Korea and compared the rainfall quantiles with those from conventional rainfall frequency analysis using a single GEV distribution. Overall, these results show that the mixed GEV distribution allows probability behavior to be taken into account during rainfall frequency analysis through the process of parameter estimation. The resulting rainfall quantile estimates were found to be significantly smaller than those determined using a single GEV distribution. The difference of rainfall quantiles was found to be closely correlated with the occurrence probability of typhoon and the distribution parameters.     

Trivariate generalized extreme value distribution in flood frequency analysis

Abstract

The multivariate extension of the logistic model with generalized extreme value (GEV) marginals is applied to provide a regional at-site flood estimate. The maximum likelihood estimators of the parameters were obtained numerically by using a multivariable constrained optimization algorithm. The asymptotic results were checked by distribution sampling techniques in order to establish whether or not those results can be utilized for small samples. A region in northern Mexico with 21 gauging stations was selected to apply the model. Results were compared with those obtained by the most popular univariate distributions, the bivariate approach of the logistic model and three regional methods: station-year, index flood and L-moments. These show that there is a reduction in the standard error of fit when estimating the parameters of the marginal distribution with the trivariate distribution instead of its univariate and bivariate counterpart, and differences between at-site and regional at-site design events can be significant as return period increases.     

Wavelet combined innovative trend analysis for precipitation data in the Euphrates-Tigris basin,Turkey

ABSTRACT

In this study, annual and seasonal precipitation trend analysis was performed in the Euphrates-Tigris basin, Turkey, using innovative trend analysis (ITA) and discrete wavelet transformation. In this context, it was seen that there is a downward trend in winter, spring and annual precipitation, whereas precipitation has an increasing tendency in summer and autumn seasons, in the greater part of the basin. When annual and seasonal data were decomposed into wavelet components, the most significant trends were observed for high-periodic wavelet components, such as D3 (8-year), D4 (16-year) and D5 (32-year), where these components represent the periods of the precipitation data. Then, the relationship between North Atlantic Oscillation (NAO) and trend in precipitation was investigated. In this regard, it was found that there could be a significant relationship between the NAO and precipitation trends of the Euphrates-Tigris basin, especially in winter, based on the wavelet ITA.     

The dependence of UK extreme sea levels and storm surges on the North Atlantic Oscillation

The role of the North Atlantic Oscillation (NAO) in effecting changes in winter extreme high and low waters and storm surges in UK waters has been investigated with the use of a depth-averaged tide+surge numerical model. Spatial patterns of correlation of extreme high and low waters (extreme still water sea levels) with the NAO index are similar to those of median or mean sea level studied previously. Explanations for the similarities, and for differences where they occur, are proposed. Spatial patterns of correlations of extreme high and low and median surge with the NAO index are similar to the corresponding extreme sea-level patterns. Suggestions are made as to which properties of surges (frequency, duration, magnitude) are linked most closely to NAO variability. Several climate models suggest higher (more positive) average values of NAO index during the next 100 years. However, the impact on the UK coastline in terms of increased flood risk should be low (aside from other consequences of climate change such as a global sea-level rise) if the existing relationships between extreme high waters and NAO index are maintained.     

Climate change and modelling of extreme temperatures in Switzerland

This study models maximum temperatures in Switzerland monitored in twelve locations using the generalised extreme value (GEV) distribution. The parameters of the GEV distribution are determined within a Bayesian framework. We find that the parameters of the underlying distribution underwent a substantial change in the beginning of the 1980s. This change is characterised by an increase both in the level and the variability. We assess the likelihood of the heat wave of the summer 2003 using the fitted GEV distribution by accounting for the presence of a structural break. The estimation results do suggest that the heat wave of 2003 is not that statistically improbable if an appropriate methodology is used for dealing with nonstationarity.     

淮河流域焦岗湖水质参数时空变化及影响因素

焦岗湖是淮河左岸一个天然湖泊,集防洪、灌溉、养殖、旅游等多种功能于一体.利用焦岗湖4个季节水质监测数据,运用Kriging方法,分析焦岗湖水质参数的时空变化及影响因素.结果表明:由于受水文季节变化过程及人类活动等综合影响,焦岗湖水质参数在时间及空间上均存在一定差异.从时间变化来看,夏季透明度较低、秋季较高;溶解氧浓度在春、冬季显著高于夏、秋季;总氮、总磷浓度与高锰酸盐指数均表现为夏季最高、秋季最低.从空间变化来看,4个季节的透明度空间差异较为显著;溶解氧浓度在春、冬季空间分布较为均匀,夏季呈现中心高周围低的变化趋势,秋季则表现为西高东低;总磷浓度春季分布较为均匀,夏、秋及冬季则呈西高东低之势;高锰酸盐指数在春、秋季节呈现东高西低之势,夏季高浓度主要集中在湖区北部,冬季浓度变化不大.     

Teleconnections between oceanic–atmospheric indices and drought over Iran using quantile regressions

ABSTRACT

In this research, the Bayesian quantile regression model is applied to investigate the teleconnections between large oceanic–atmospheric indices and drought standardized precipitation index (SPI) in Iran. The 12-month SPI time series from 138 synoptic stations for 1952–2014 were selected as the drought index. Three oceanic–atmospheric indices, the North Atlantic Oscillation (NAO), the Southern Oscillation Index (SOI) and the Multivariate El Niño/Southern Oscillation Index (MEI), were selected as covariates. The results show that NAO has the weakest impact on drought in different quantiles and different regions in Iran. La Niña conditions amplified droughts through all SPI quantiles in western, Caspian Sea coastal regions and southern regions. The positive phase of MEI significantly modulates low SPI quantiles (i.e. drought conditions) throughout the Zagros region, Caspian Sea coastal regions and southern regions. The study shows that the effect of large oceanic–atmospheric indices have heterogeneous impacts on extreme dry and wet conditions.     

Effect of the North Atlantic Oscillation on the Thermal Characteristics of Lakes in Poland

This paper presents the effect of the North Atlantic Oscillation (NAO) on the thermal characteristics of lakes in Poland. In the analysis, the use was made of monthly air temperatures recorded at fifteen meteorological stations, water temperatures of twelve lakes, and Hurrell’s winter NAO indices. Over the study period (1971–2010), there was a marked increase in the temperatures of both, air and lake waters. Depending on the NAO phase, water temperatures were observed to depart from mean values, being markedly higher than average (even by 1°C) in the positive winter NAO phase. The differences in water temperatures were statistically significant in the winter-spring season. In turn, in the negative NAODJFM phase lake water temperatures in winter and spring were markedly lower than average (in March even by 1.0°C). The unique response of some lakes depends on their morphometric parameters, including their mean depth.     

Time variations of the relationships between the North Atlantic Oscillation and European winter temperature and precipitation

The degree of stationarity of relationships between the NAO index and long European temperature and precipitation series in winter is quantified by running correlations with a time window of 31 years at 29 and 27 stations in Europe, respectively. They indicate major nonstationarities in the NAO-to-surface climate relationships at most stations. The temporal course most common for correlations with temperature is a slight change prior to about 1950, followed by an increase; for precipitation, a typical course is a decrease in the first half of the 20^th century, followed by an increase. The temporal variations in correlations do not result from the presence of trends in the time series. The periods of high correlations with temperature are accompanied with an eastward shift of both NAO action centres; the eastward shift is thus at least partly responsible for the time variations in correlations. huth@ufa.cas.cz     

Analysis of dry/wet conditions using the standardized precipitation index and its potential usefulness for drought/flood monitoring in Hunan Province, China

Local dry/wet conditions and extreme rainfall events are of great concern in regional water resource and disaster risk management. Extensive studies have been carried out to investigate the change of dry/wet conditions and the adaptive responses to extreme rainfall events within the context of climate change. However, applicable tools and their usefulness are still not sufficiently studied, and in Hunan Province, a major grain-producing area in China that has been frequently hit by flood and drought, relevant research is even more limited. This paper investigates the spatiotemporal variation of dry/wet conditions and their annual/seasonal trends in Hunan with the standardized precipitation index (SPI) at various time scales. Furthermore, to verify the potential usefulness of SPI for drought/flood monitoring, the correlation between river discharge and SPI at multiple time scales was examined, and the relation between extreme SPI and the occurrence of historical drought/flood events is explored. The results indicate that the upper reaches of the major rivers in Hunan Province have experienced more dry years than the middle and lower reaches over the past 57 years, and the region shows a trend of becoming drier in the spring and autumn seasons and wetter in the summer and winter seasons. We also found a strong correlation between river discharge and SPI series, with the maximum correlation coefficient occurred at the time scale of 2 months. SPI at different time scales may vary in its usefulness in drought/flood monitoring, and this highlights the need for a comprehensive consideration of various time scales when SPI is employed to monitor droughts and floods.     

Recent trends of extreme precipitation indices in the Iberian Peninsula using observations and WRF model results

Spatial and temporal distributions of the trends of extreme precipitation indices were analysed between 1986 and 2005, over the Iberian Peninsula (IP). The knowledge of the patterns of extreme precipitation is important for impacts assessment, development of adaptation and mitigation strategies. As such, there is a growing need for a more detailed knowledge of precipitation climate change.This analysis was performed for Portuguese and Spanish observational datasets and results performed by the Weather Research and Forecast (WRF) model forced by the ERA-Interim reanalysis. Extreme precipitation indices recommended by the Expert Team for Climate Change Detection Monitoring and Indices were computed, by year and season. Then, annual and seasonal trends of the indices were estimated by Theil-Sen method and their significance was tested by the Mann-Kendal test. Additionally, a second simulation forced by the Max Planck Institute Earth System Model (MPI-ESM), was considered. This second modelling configuration was created in order to assess its performance when simulating extremes of precipitation.The annual trends estimated for the 1986–2005, from the observational datasets and from the ERA-driven simulation reveal: 1) negative statistically significant trends of the CWD index in the Galicia and in the centre of the IP; 2) positive statistically significant trends of the CDD index over the south of the IP and negative statistically significant trends in Galicia, north and centre of Portugal; 3) positive statistically significant trends of the R75p index in some regions of the north of the IP; 4) positive statistically significant trends in the R95pTOT index in the Central Mountains Chain, Leon Mountains and in the north of Portugal.Seasonally, negative statistically significant trends of the CWD index were found in Galicia, in winter and in the south of the IP, in summer. Positive statistically significant trends of the CWD index were identified in the Leon Mountains, in spring, and in Galicia, in autumn. For the CDD index, negative statistically significant trends were seen in Valencia, in the spring, and, in Galicia and Portugal (north and centre), in summer. Positive statistically significant trends of the CDD index were found: in the east of the IP, in the winter; in the Cantabrian Mountain, in the spring; and, in the south of the IP, in summer. Regarding to the R75p index, negative statistically significant trends were found in Galicia, in winter and positive statistically significant trends in the north of Portugal, in spring and in the Central Mountains Chain and north of Portugal, in autumn. For the R95pTOT index, negative statistically significant trends were found over the Sierra Cuenca and Sierra Cazorla, in winter and positive statistically significant trends were found over the Sierra Cebollera, in winter and in Castile-la Mancha region, in spring.The results of the annual and seasonal trends of the extreme precipitation indices performed for observational datasets and the simulation forced by ERA-Interim, are similar. The results obtained for the simulation forced by MPI-ESM are not satisfactory, and can be a source of criticism for the use of simulation forced by MPI-ESM in this type of climate change studies. Even for the relatively short period used, the WRF model, when properly forced is a useful tool due to the similar results of Portuguese and Spanish observational datasets and the simulation forced by ERA-Interim.