A feedback mechanism in annual rainfall, Central Sudan

Annual rainfall in many parts of the world is an independent process. Yet annual rainfall series in some regions of Africa show characteristics incompatible with such an hypothesis.

The annual rainfall process in Central Sudan is weakly dependent. The hypothesis that the dependence is due to a “feedback mechanism” is investigated using a mathematical model based on the water balance of the neighbouring region, Bahr Elghazal.     

Pooled flood frequency analysis: a comparison based on peaks-over-threshold and annual maximum series

Despite some theoretical advantages of peaks-over-threshold (POT) series over annual maximum (AMAX) series, some practical aspects of flood frequency analysis using AMAX or POT series are still subject to debate. Only minor attention has been given to the POT method in the context of pooled frequency analysis. The objective of this research is to develop a framework to promote the implementation of pooled frequency modelling based on POT series. The framework benefits from a semi-automated threshold selection method. This study introduces a formalized and effective approach to construct homogeneous pooling groups. The proposed framework also offers means to compare the performance of pooled flood estimation based on AMAX or POT series. An application of the framework is presented for a large collection of Canadian catchments. The proposed POT pooling technique generally improved flood quantile estimation in comparison to the AMAX pooling scheme, and achieved smaller uncertainty associated with the quantile estimates.     

Frequency analyses of annual extreme rainfall series from 5 min to 24 h

The parameter estimation methods of (1) moments, (2) maximum‐likelihood, (3) probability‐weighted moments (PWM) and (4) self‐determined PWM are applied to the probability distributions of Gumbel, general extreme values, three‐parameter log‐normal (LN3), Pearson‐3 and log‐Pearson‐3. The special method of computing parameters so as to make the sample skewness coefficient zero is also applied to LN3, and hence, altogether 21 candidate distributions resulted. The parameters of these distributions are computed first by original sample series of 14 successive‐duration annual extreme rainfalls recorded at a rain‐gauging station. Next, the parameters are scaled by first‐degree semi‐log or log‐log polynomial regressions versus rainfall durations from 5 to 1440 min (24 h). Those distributions satisfying the divergence criterion for frequency curves are selected as potential distributions, whose better‐fit ones are determined by a conjunctive evaluation of three goodness‐of‐fit tests. Frequency tables, frequency curves and intensity–duration–frequency curves are the outcome. Copyright © 2010 John Wiley & Sons, Ltd.     

Regional frequency analysis of rainfall extremes in Southern Malawi using the index rainfall and L-moments approaches

Rainfall extremes often result in the occurrence of flood events with associated loss of life and infrastructure in Malawi. However, an understanding of the frequency of occurrence of such extreme events either for design or disaster planning purposes is often limited by data availability at the desired temporal and spatial scales. Regionalisation, which involves “trading time for space” by pooling together observations for stations with similar behavior, is an alternative approach for more accurate determination of extreme events even at ungauged areas or sites with short records. In this study, regional frequency analysis of rainfall extremes in Southern Malawi, large parts of which are flood prone, was undertaken. Observed 1-, 3-, 5- and 7-day annual maximum rainfall series for the period 1978–2007 at 23 selected rainfall stations in Southern Malawi were analysed. Cluster analysis using scaled at-site characteristics was used to determine homogeneous rainfall regions. L-moments were applied to derive regional index rainfall quantiles. The procedure also validated the three rainfall regions identified through homogeneity and heterogeneity tests based on Monte Carlo simulations with regional average L-moment ratios fitted to the Kappa distribution. Based on assessments of the accuracy of the derived index rainfall quantiles, it was concluded that the performance of this regional approach was satisfactory when validated for sites not included in the sample data. The study provides an estimate of the regional characteristics of rainfall extremes that can be useful in among others flood mitigation and engineering design.     

Non-stationary frequency analysis of heavy rainfall events in southern France

Abstract

Heavy rainfall events often occur in southern French Mediterranean regions during the autumn, leading to catastrophic flood events. A non-stationary peaks-over-threshold (POT) model with climatic covariates for these heavy rainfall events is developed herein. A regional sample of events exceeding the threshold of 100 mm/d is built using daily precipitation data recorded at 44 stations over the period 1958–2008. The POT model combines a Poisson distribution for the occurrence and a generalized Pareto distribution for the magnitude of the heavy rainfall events. The selected covariates are the seasonal occurrence of southern circulation patterns for the Poisson distribution parameter, and monthly air temperature for the generalized Pareto distribution scale parameter. According to the deviance test, the non-stationary model provides a better fit to the data than a classical stationary model. Such a model incorporating climatic covariates instead of time allows one to re-evaluate the risk of extreme precipitation on a monthly and seasonal basis, and can also be used with climate model outputs to produce future scenarios. Existing scenarios of the future changes projected for the covariates included in the model are tested to evaluate the possible future changes on extreme precipitation quantiles in the study area.

Editor Z.W. Kundzewicz; Associate editor K. Hamed

Citation Tramblay, Y., Neppel, L., Carreau, J., and Najib, K., 2013. Non-stationary frequency analysis of heavy rainfall events in southern France. Hydrological Sciences Journal, 58 (2), 280–294.     

The use of the exponent K(q) function to delimit homogeneous regions in regional frequency analysis of extreme annual daily rainfall

The regional frequency analysis of extreme annual rainfall data is a useful methodology in hydrology to obtain certain quantile values when no long data series are available. The most crucial step in the analysis is the grouping of sites into homogeneous regions. This work presents a new grouping criterion based on some multifractal properties of rainfall data. For this purpose, a regional frequency analysis of extreme annual rainfall data from the Maule Region (Chile) has been performed. Daily rainfall data series of 53 available stations have been studied, and their empirical moments scaling exponent functions K(q) have been obtained. Two characteristics parameters of the K(q) functions (γ_max and K(0)) have been used to group the stations into three homogeneous regions. Only five sites have not been possible to include into any homogenous regions, being the local frequency analysis of extreme daily rainfall the most appropriate method to be used at these locations. Copyright © 2014 John Wiley & Sons, Ltd.     

Non-stationary frequency analysis of extreme precipitation in South Korea using peaks-over-threshold and annual maxima

The conventional approach to the frequency analysis of extreme precipitation is complicated by non-stationarity resulting from climate variability and change. This study utilized a non-stationary frequency analysis to better understand the time-varying behavior of short-duration (1-, 6-, 12-, and 24-h) precipitation extremes at 65 weather stations scattered across South Korea. Trends in precipitation extremes were diagnosed with respect to both annual maximum precipitation (AMP) and peaks-over-threshold (POT) extremes. Non-stationary generalized extreme value (GEV) and generalized Pareto distribution (GPD) models with model parameters made a linear function of time were applied to AMP and POT respectively. Trends detected using the Mann–Kendall test revealed that the stations showing an increasing trend in AMP extremes were concentrated in the mountainous areas (the northeast and southwest regions) of South Korea. Trend tests on POT extremes provided fairly different results, with a significantly reduced number of stations showing an increasing trend and with some stations showing a decreasing trend. For most of stations showing a statistically significant trend, non-stationary GEV and GPD models significantly outperformed their stationary counterparts, particularly for precipitation extremes with shorter durations. Due to a significant-increasing trend in the POT frequency found at a considerable number of stations (about 10 stations for each rainfall duration), the performance of modeling POT extremes was further improved with a non-homogeneous Poisson model. The large differences in design storm estimates between stationary and non-stationary models (design storm estimates from stationary models were significantly lower than the estimates of non-stationary models) demonstrated the challenges in relying on the stationary assumption when planning the design and management of water facilities. This study also highlighted the need of caution when quantifying design storms from POT and AMP extremes by showing a large discrepancy between the estimates from those two approaches.     

Complexity analysis of rainfall and runoff time series based on sample entropy in different temporal scales

This study applied sample entropy (SampEn) to rainfall and runoff time series to investigate the complexity of different temporal scales. Rainfall and runoff time series with intervals of 1, 10, 30, 90, and 365 days for the Wu-Tu upstream watershed were used. Thereafter, SampEn was computed for the five rainfall and runoff time series. The results show that for the various temporal scales, comparisons of the complexity between the rainfall and runoff time series based on the SampEn are inconsistent. Calculating the dynamic SampEn further elucidated variations of the complexity in the rainfall and runoff time series. In addition, the results show that SampEn measures of the rainfall and runoff time series are typically higher than the approximate entropy measures of the rainfall and runoff time series for a specific temporal scale. The complexity increases when the sample size increases for a specific temporal scale. Furthermore, temporal scales with low complexity and high predictability are obtained from the variations of SampEn for the rainfall and runoff time series with different temporal scales, thereby providing a reference for determining the appropriate temporal scale for rainfall and runoff time series forecasting.     

Temporal and spatial variability of annual and seasonal rainfall over Ethiopia

Abstract

Characterization of the seasonal and inter-annual spatial and temporal variability of rainfall in a changing climate is vital to assess climate-induced changes and suggest adequate future water resources management strategies. Trends in annual, seasonal and maximum 30-day extreme rainfall over Ethiopia are investigated using 0.5° latitude?×?0.5° longitude gridded monthly precipitation data. The spatial coherence of annual rainfall among contiguous rainfall grid points is also assessed for possible spatial similarity across the country. The correlation between temporally coinciding North Atlantic Multidecadal Oscillation (AMO) index and annual rainfall variability is examined to understand the underlying coherence. In total 381 precipitation grid points covering the whole of Ethiopia with five decades (1951–2000) of precipitation data are analysed using the Mann-Kendall test and Moran spatial autocorrelation method. Summer (July–September) seasonal and annual rainfall data exhibit significant decreasing trends in northern, northwestern and western parts of the country, whereas a few grid points in eastern areas show increasing annual rainfall trends. Most other parts of the country exhibit statistically insignificant trends. Regions with high annual and seasonal rainfall distribution exhibit high temporal and spatial correlation indices. Finally, the country is sub-divided into four zones based on annual rainfall similarity. The association of the AMO index with annual rainfall is modestly good for northern and northeastern parts of the country; however, it is weak over the southern region.

Editor Z.W. Kundzewicz; Associate editor S. Uhlenbrook

Citation Wagesho, N., Goel, N.K., and Jain, M.K. 2013. Temporal and spatial variability of annual and seasonal rainfall over Ethiopia. Hydrological Sciences Journal, 58 (2), 354–373.     

Time series analysis for evaluating hydrological responses of pore-water pressure to rainfall in a slope

Characterization of pore-water pressure at the soil–atmosphere interface is a major requirement in relation to slope instability. A rain-gauge and five piezometers (BH1–BH5) were installed on a slope located in a Korean military base. The upper slope (BH4, BH5) was covered with plastic sheets to prevent rainwater from percolating into the slope due to safety issues. Rainfall is matched by prompt changes in the pore-water pressure except at BH5. Due to the plastic cover, the pore water does not show any significant change in the early period by evapotranspiration. From correlation analysis, two wells (BH3, BH5) have longer memory effects due to matrix flow of past precipitation. Two principal components show hydrological responses of pore water to rainfall during intense rainfall, but PC2 does not indicate any important changes in low or no rainfall. This study suggests that correlation analysis with PCA can be a valuable tool for interpreting datasets consisting only of pore-water pressure.     

Comparison between the peaks-over-threshold method and the annual maximum method for flood frequency analysis

Abstract

Flood frequency analysis can be made by using two types of flood peak series, i.e. the annual maximum (AM) and peaks-over-threshold (POT) series. This study presents a comparison of the results of both methods for data from the Litija 1 gauging station on the Sava River in Slovenia. Six commonly used distribution functions and three different parameter estimation techniques were considered in the AM analyses. The results showed a better performance for the method of L-moments (ML) when compared with the conventional moments and maximum likelihood estimation. The combination of the ML and the log-Pearson type 3 distribution gave the best results of all the considered AM cases. The POT method gave better results than the AM method. The binomial distribution did not offer any noticeable improvement over the Poisson distribution for modelling the annual number of exceedences above the threshold.

Editor D. Koutsoyiannis

Citation Bezak, N., Brilly, M., and ?raj, M., 2014. Comparison between the peaks-over-threshold method and the annual maximum method for flood frequency analysis. Hydrological Sciences Journal, 59 (5), 959–977.     

Assessment of GCM simulations of annual and seasonal rainfall and daily rainfall distribution across south‐east Australia

Global warming can potentially lead to changes in future rainfall and runoff and can significantly impact the regional hydrology and future availability of water resources. All the large‐scale climate impact studies use the future climate projections from global climate models (GCMs) to estimate the impact on future water availability. This paper presents results from a detailed assessment to investigate the capability of 15 GCMs to reproduce the observed historical annual and seasonal mean rainfalls, the observed annual rainfall series and the observed daily rainfall distribution across south‐east Australia. The assessment shows that the GCMs can generally reproduce the spatial patterns of mean seasonal and annual rainfalls. However, there can be considerable differences between the mean rainfalls simulated by the GCMs and the observed rainfall. The results clearly show that none of the GCMs can simulate the actual annual rainfall time series or the trend in the annual rainfall. The GCMs can also generally reproduce the observed daily (ranked) rainfall distribution at the GCM scale. The GCMs are ranked against their abilities to reproduce the observed historical mean annual rainfall and daily rainfall distribution, and, based on the combined score, the better GCMs include MPI‐ECHAM5, MIUB, CCCMA_T47, INMCM, CSIRO‐MK3·0, CNRM, CCCMA_T63 and GFDL 2·0 and those with poorer performances are MRI, IPSL, GISS‐AOM, MIROC‐M, NCAR‐PCM1, IAP and NCAR‐CCSM. However, the reduction in the combined score as we move from the best‐ to the worst‐performing GCMs is gradual, and there is no evident cut‐off point or threshold to remove GCMs from climate impact studies. There is some agreement between the results here and many similar studies comparing the performance of GCMs in Australia, but the results are not always consistent and do significantly disagree with several of the studies. Copyright © 2010 John Wiley & Sons, Ltd.     

Crustal structure beneath the central and western North China from receiver function analysis

The North China Craton (NCC) is one of the oldest cratons on earth. Several important tectonic transformations of Mesozoic-Cenozoic tectonic regime led to the destruction of the North China craton. The knowledge of crustal structure can provide important constraints for the formation and evolution of cratons. New maps of sediment thickness, crustal thickness (H) and v_P/v_S (κ) in the central and western NCC were obtained using sequential H-κ stacking. P-wave receiver functions are calculated using teleseismic waveform data recorded by 405 stations from ChinArray project. Benefiting from the densely distribution of temporary seismic stations, our results reveal details of the crustal structure in the study area. The thickness of sedimentary layer in North China ranges from 0–6.4 km, and the thickest sedimentary layer is in Ordos block and its surroundings (about 2.8–6 km); The thickness of sedimentary layer in the Mongolia fold belt and Yinshan orogenic belt is relatively thin (less than 1 km). The crustal thickness of the study area varies between 27–48 km, of which the crust of the North China Plain is about 30–33 km, the central NCC is about 33–40 km, and the Ordos block is 40–48 km thick. The average v_P/v_S ratios in the study area is mostly between 1.66 and 1.90, and that in the Yanshan-Taihang mountain fold belt is between 1.70 and 1.85, and that in the Ordos block is between 1.65 and 1.90, with an average value of 1.77, indicating the absence of a thick basaltic lower crust. The obvious negative correlation between crustal thickness and average v_P/v_S ratio within Ordos and Central Asia orogenic belt may be related to magmatic underplating during the crustal formation. There is no significant correlation between the crustal thickness and the v_P/v_S ratio in the Lüliang-Taihang mountain fold belt, which may be related to the multiple geological processes such as underplating and crustal extension and thinning in this area. The lack of correlation between crust thickness and topography in the central orogenic belt and the North China Basin indicates the topography of these areas are controlled not only by crustal isostatic adjustment but also by the lithospheric mantle processes.     

Update of likelihood-based ground-motion model selection for seismic hazard analysis in western central Europe

Scherbaum et al. [(2004) Bull Seismolo Soc Am 94(6): 2164–2185] proposed a likelihood-based approach to select and rank ground-motion models for seismic hazard analysis in regions of low-seismicity. The results of their analysis were first used within the PEGASOS project [Abrahamson et al. (2002), In Proceedings of the 12 ECEE, London, 2002, Paper no. 633] so far the only application of a probabilistic seismic hazard analysis (PSHA) in Europe which was based on a SSHAC Level 4 procedure [(Budnitz et al. 1997, Recommendations for PSHA: guidance on uncertainty and use of experts. No. NUREG/CR-6372-V1). The outcome of this project have generated considerable discussion (Klügel 2005, Eng Geol 78:285–307, 2005b) Eng Geol 78: 285–307, (2005c) Eng Geol 82: 79–85 Musson et al. (2005) Eng Geol 82(1): 43–55]; Budnitz et al. (2005), Eng Geol 78(3–4): 285–307], a central part of which is related to the issue of ground-motion model selection and ranking. Since at the time of the study by Scherbaum et al. [(2004.) Bull Seismolo Soc Am 94(6): 2164–2185], only records from one earthquake were available for the study area, here we test the stability of their results using more recent data. Increasing the data set from 12 records of one earthquake in Scherbaum et al. [(2004) Bull Seismolo Soc Am 94(6): 2164–2185] to 61 records of 5 earthquakes, which have mainly occurred since the publication of the original study, does not change the set of the three top-ranked ground-motion models [Abrahamson and Silva (1997) Seismolo Res Latt 68(1): 94–127; Lussou et al. (2001) J Earthquake Eng 5(1):13–33; Berge-Thierry et al. (2003) Bull Seismolog Soc Am 95(2): 377–389. Only for the lower-ranked models do we obtain modifications in the ranking order. Furthermore, the records from the Waldkirch earthquake (Dec, 5th, 2004, M _w = 4.9) enabled us to develop a new stochastic model parameter set for the application of Campbell’s [(2003) Bull Seismolo Soc Am 93(3): 1012–1033] hybrid empirical model to SW Germany and neighbouring regions.     

Drought hotspot analysis and risk assessment using probabilistic drought monitoring and severity–duration–frequency analysis

Drought hotspot identification requires continuous drought monitoring and spatial risk assessment. The present study analysed drought events in the agriculture‐dominated mid‐Mahanadi River Basin in Odisha, India, using crop water stress as a drought indicator. This drought index incorporated different factors that affect crop water deficit such as the cropping pattern, soil characteristics, and surface soil moisture. The drought monitoring framework utilized a relevance vector machine model‐based classification that provided the uncertainty associated with drought categorization. Using the proposed framework, drought hotspots are identified in the study region and compared with indices based on precipitation and soil moisture. Further, a bivariate copula is employed to model the agricultural drought characteristics and develop the drought severity–duration–frequency (S–D–F) relationships. The drought hotspot maps and S–D–F curves are developed for different locations in the region. These provided useful information on the site‐specific drought patterns and the characteristics of the devastating droughts of 2002 and 2012, characterized by an average drought duration of 7 months at several locations. The site‐specific risk of short‐ and long‐term agricultural droughts are then investigated using the conditional copula. The results suggest that the conditional return periods and the S–D–F curves are valuable tools to assess the spatial variability of drought risk in the region.     

Application of Fourier series in managing the seasonality of convective and monsoon rainfall

ABSTRACT

This study assesses the performance of Fourier series in representing seasonal variations of the tropical rainfall process in Malaysia. Fourier series are incorporated into a spatial-temporal stochastic model in an attempt to make the model parsimonious and, at the same time, capture the annual variation of rainfall distribution. In view of Malaysia’s main rainfall regime, the model is individually fitted for two regions with distinctive rainfall profiles: one being an urban area receiving rainfall from convective activities whilst the other receives rainfall from monsoonal activities. Since both regions are susceptible to floods, the study focuses on the rainfall process at fine resolution. Fourier series equations are developed to represent the model’s parameters to describe their annual periodicity. The number of significant harmonics for each parameter is determined by inspecting the cumulative fraction of total variance explained by the significant harmonics. Results reveal that the number of significant harmonics assigned for the parameters is slightly higher in the region with monsoonal rains. The overall simulation results show that the proposed model is capable of generating tropical rainfall series from convective and monsoonal activities.

Editor D. Koutsoyiannis Associate editor K. Hamed     

Segmentation-based approach for trend analysis and structural breaks in rainfall time series (1851–2006) over India

ABSTRACT

This study analysed long-term rainfall data (1851–2006) over seven climatic zones of India at seasonal and annual scales based on three techniques: (i) linear regression, (ii) multifractal detrended fluctuation analysis (MFDFA) and (iii) Bayesian algorithm. The linear regression technique was used for trend analysis of short-term (30 years) and long-term (156 years) rainfall data. The MFDFA revealed small- and large-scale fluctuations, whereas the Bayesian algorithm helped in quantifying the uncertainty in break-point detection from the rainfall time series. Major break points years identified through Bayesian algorithm were 1888, 1904 and 1976. The MFDFA technique identified that high fluctuation years were between 1871–1890, 1891–1910 and 1951–1970. Linear regression-based analysis revealed 1881–1910 and 1971–2006 as break-point periods in the North Mountainous Indian region. A similar analysis was carried out for India as a whole, as well as its seven climatic zones.