Hydrological Drought Assessment Based on the Standardized Streamflow Index: A Case Study of the Three Cape Provinces of South Africa

Global impacts of drought conditions pose a major challenge towards the achievement of the 2030 Sustainable Development Goals. As a result, a clarion call for nations to take actions aimed at mitigating the adverse negative effects, managing key natural resources and strengthening socioeconomic development can never be overemphasized. The present study evaluated hydrological drought conditions in three Cape provinces (Eastern, Western and Northern Cape) of South Africa, based on the Standardized Streamflow Index (SSI) calculated at 3- and 6-month accumulation periods from streamflow data spanning over the 3.5 decades. The SSI features were quantified by assessing the corresponding annual trends computed by using the Modified Mann–Kendall test. Drought conditions were also characterized in terms of the duration and severity across the three Cape provinces. The return levels of drought duration (DD) and drought severity (DS) associated with 2-, 5-, 10-, 20- and 50-year periods were estimated based on the generalized extreme value (GEV) distribution. The results indicate that hydrological drought conditions have become more frequent and yet exhibit spatial contrasts throughout the study region during the analyzed period. To this end, there is compelling evidence that DD and DS have increased over time in the three Cape provinces. Return levels analysis across the studied periods also indicate that DD and DS are expected to be predominant across the three Cape provinces, becoming more prolonged and severe during the extended periods (e.g., 20- and 50-year). The results of the present study (a) contribute to the scientific discourse of drought monitoring, forecasting and prediction and (b) provide practical insights on the nature of drought occurrences in the region. Consequently, the study provides the basis for policy- and decision-making in support of preparedness for and adaptation to the drought risks in the water-linked sectors and robust water resource management. Based on the results reported in this study, it is recommended that water agencies and the government should be more proactive in searching for better strategies to improve water resources management and drought mitigation in the region.

Using high-resolution regional climate models to estimate return levels of daily extreme precipitation over Bavaria

Extreme daily rainfall is an important trigger for floods in Bavaria. The dimensioning of water management structures as well as building codes is based on observational rainfall return levels. In this study, three high-resolution regional climate models (RCMs) are employed to produce 10- and 100-year daily rainfall return levels and their performance is evaluated by comparison to observational return levels. The study area is governed by different types of precipitation (stratiform, orographic, convectional) and a complex terrain, with convective precipitation also contributing to daily rainfall levels. The Canadian Regional Climate Model version 5 (CRCM5) at a 12 km spatial resolution and the Weather and Forecasting Research (WRF) model at a 5 km resolution both driven by ERA-Interim reanalysis data use parametrization schemes to simulate convection. WRF at a 1.5 km resolution driven by ERA5 reanalysis data explicitly resolves convectional processes. Applying the generalized extreme value (GEV) distribution, the CRCM5 setup can reproduce the observational 10-year return levels with an areal average bias of +6.6 % and a spatial Spearman rank correlation of ρ=0.72. The higher-resolution 5 km WRF setup is found to improve the performance in terms of bias (+4.7 %) and spatial correlation (ρ=0.82). However, the finer topographic details of the WRF-ERA5 return levels cannot be evaluated with the observation data because their spatial resolution is too low. Hence, this comparison shows no further improvement in the spatial correlation (ρ=0.82) but a small improvement in the bias (2.7 %) compared to the 5 km resolution setup. Uncertainties due to extreme value theory are explored by employing three further approaches. Applied to the WRF-ERA5 data, the GEV distributions with a fixed shape parameter (bias is +2.5 %; ρ=0.79) and the generalized Pareto (GP) distributions (bias is +2.9 %; ρ=0.81) show almost equivalent results for the 10-year return period, whereas the metastatistical extreme value (MEV) distribution leads to a slight underestimation (bias is −7.8 %; ρ=0.84). For the 100-year return level, however, the MEV distribution (bias is +2.7 %; ρ=0.73) outperforms the GEV distribution (bias is +13.3 %; ρ=0.66), the GEV distribution with fixed shape parameter (bias is +12.9 %; ρ=0.70), and the GP distribution (bias is +11.9 %; ρ=0.63). Hence, for applications where the return period is extrapolated, the MEV framework is recommended. From these results, it follows that high-resolution regional climate models are suitable for generating spatially homogeneous rainfall return level products. In regions with a sparse rain gauge density or low spatial representativeness of the stations due to complex topography, RCMs can support the observational data. Further, RCMs driven by global climate models with emission scenarios can project climate-change-induced alterations in rainfall return levels at regional to local scales. This can allow adjustment of structural design and, therefore, adaption to future precipitation conditions.

Uncertainty analysis of strong earthquake hazard estimation based on the generalized Pareto distribution model: A case study of the northeastern Tibetan Plateau

A statistical method to analyze the uncertainties of strong earthquake hazard estimation is proposed from Generalized Pareto distribution (GPD) model using the northeastern Tibetan Plateau earthquake catalogue (1885–2017) data. For magnitude threshold of 5.5, the magnitude return levels in 20, 50, 100, 200, and 500 years are 7.19, 7.70, 7.99, 8.22, and 8.45, respectively. The corresponding 95% confidence intervals are [6.77, 7.60], [7.27, 8.12], [7.53, 8.44], [7.71, 8.72], and [7.84, 9.06], respectively. The upper magnitude limit obtained from this GPD model is 9.07 and its 80% confidence interval is [8.16, 9.98]. The sensitivity analysis by the Morris method indicates that the input magnitude threshold has a relatively large influence on the estimation results. Thus, threshold selection is important for the GPD model construction. The sensitivity characteristic ranking of input factors become increasingly stable with the increasing of return period, which implies that GPD model is more suitable for estimating strong earthquakes magnitude return levels and upper magnitude limit. The GPD modeling approach and qualitative uncertainty analysis methods for strong earthquake hazard estimations proposed in this paper can be applied to seismic hazard analysis elsewhere.

An Analysis of Precipitation Extreme Events Based on the SPI and EDI Values in the Free State Province, South Africa

Recognizing that, over the last several years, extreme rainfall has led to hazardous stress in humans, animals, plants, and even infrastructure, in the present study, we aimed to investigate the characteristics of droughts over the Free State (FS) Province of South Africa in order to determine the future likelihood of reoccurrences of precipitation extremes using the generalized extreme value distribution (GEV) and extreme frequency analysis (EFA). In this regard, daily rainfall datasets from nine South African weather service homogenous climatic districts, spanning from 1980 to 2019, were used to compute: (a) the total annual rainfall, (b) the Effective Drought Index (EDI), and (c) the Standard Precipitation Index (SPI). The SPI was calculated for 3, 6, and 12 month accumulation periods (hereafter SPI-3, SPI-6, and SPI-12, respectively). The trend analysis results of the EDI and SPI-3, -6, and -12 showed that the Free State Province is generally negative, illustrating persistent drought. An analysis of the GEV parameters across the EDI and SPI-3, -6, and -12 values illustrated that the location, scale, and shape parameters exhibited a noticeable spatial variability across the Free State Province with the location parameter largely negative, the scale parameter largely positive, while the shape parameter pointed to an inherent Type III (Weibull) GEV distribution. In addition, the return levels for the drought/wet duration and severity of the EDI and SPI-3, -6, and -12 values generally showed increasing patterns across the corresponding return periods; the spatial contrasts were only noticeable in the return levels derived from the wet/drought duration and severity derived from SPI-3, -6, and -12 values (and not in the EDI). Further, the EFA results pointed to a noticeable spatial contrast in the return periods derived from the EDI and SPI-3, -6, and -12 values for each of the extreme precipitation categories: moderately wet, severely wet, extremely wet to moderately dry, and severely dry. Over four decades, the FS Province has generally experienced a suite of extreme precipitation categories ranging from moderately wet, severely wet, extremely wet to moderately dry, severely dry, and extremely dry conditions. Overall, the present study contributes towards implementation of effective drought early warning systems and can be used to enhance drought related policy and decision making in support of water resource management and planning in the FS Province.

Expected Recurrence of Extreme Winds in Northwestern Sahara and Associated Uncertainties

Estimating the probability of the occurrence of hazardous winds is crucial for their impact in human activities; however, this is inherently affected by the shortage of observations. This becomes critical in poorly sampled regions, such as the northwestern Sahara, where this work is focused. The selection of any single methodological variant contributes with additional uncertainty. To gain robustness in the estimates, we expand the uncertainty space by applying a large body of methodologies. The methodological uncertainty is constrained afterward by keeping only the reliable experiments. In doing so, we considerably narrow the uncertainty associated with the wind return levels. The analysis suggest that not necessarily all methodologies are equally robust. The highest 10-min speed (wind gust) for a return period of 50 years is about 45 ms−1 (56 ms−1). The intensity of the expected extreme winds is closely related to orography. The study is based on wind and wind gust observations that were collected and quality controlled for the specific purposes herein. We also make use of a 12-year high-resolution regional simulation to provide simulation-based wind return level maps that endorse the observation-based results. Such an exhaustive methodological sensitivity analysis with a long high-resolution simulation over this region was lacking in the literature.

Local-scale variability of seasonal mean and extreme values of in situ snow depth and snowfall measurements

Daily measurements of snow depth and snowfall can vary strongly over short distances. However, it is not clear if there is a seasonal dependence in these variations and how they impact common snow climate indicators based on mean values, as well as estimated return levels of extreme events based on maximum values. To analyse the impacts of local-scale variations we compiled a unique set of parallel snow measurements from the Swiss Alps consisting of 30 station pairs with up to 77 years of parallel data. Station pairs are usually located in the same villages (or within 3 km horizontal and 150 m vertical distances). Investigated snow climate indicators include average snow depth, maximum snow depth, sum of new snow, days with snow on the ground, days with snowfall, and snow onset and disappearance dates, which are calculated for various seasons (December to February (DJF), November to April (NDJFMA), and March to April (MA)). We computed relative and absolute error metrics for all these indicators at each station pair to demonstrate the potential variability. We found the largest relative inter-pair differences for all indicators in spring (MA) and the smallest in DJF. Furthermore, there is hardly any difference between DJF and NDJFMA, which show median variations of less than 5 % for all indicators. Local-scale variability ranges between less than 24 % (DJF) and less than 43 % (MA) for all indicators and 75 % of all station pairs. The highest percentage (90 %) of station pairs with variability of less than 15 % is observed for days with snow on the ground. The lowest percentage (30 %) of station pairs with variability of less than 15 % is observed for average snow depth. Median differences of snow disappearance dates are rather small (3 d) and similar to the ones found for snow onset dates (2 d). An analysis of potential sunshine duration could not explain the higher variabilities in spring. To analyse the impact of local-scale variations on the estimation of extreme events, 50-year return levels were quantified for maximum snow depth and maximum 3 d new snow sum, which are often used for avalanche prevention measures. The found return levels are within each other's 95 % confidence intervals for all (but three) station pairs, revealing no striking differences. The findings serve as an important basis for our understanding of variabilities of commonly used snow indicators and extremal indices. Knowledge about such variabilities in combination with break-detection methods is the groundwork in view of any homogenization efforts regarding snow time series.

Elevation-dependent trends in extreme snowfall in the French Alps from 1959 to 2019

Climate change projections indicate that extreme snowfall is expected to increase in cold areas, i.e., at high latitudes and/or high elevation, and to decrease in warmer areas, i.e., at mid-latitudes and low elevation. However, the magnitude of these contrasting patterns of change and their precise relations to elevation at the scale of a given mountain range remain poorly known. This study analyzes annual maxima of daily snowfall based on the SAFRAN reanalysis spanning the time period 1959–2019 and provided within 23 massifs in the French Alps every 300 m of elevation. We estimate temporal trends in 100-year return levels with non-stationary extreme value models that depend on both elevation and time. Specifically, for each massif and four elevation ranges (below 1000, 1000–2000, 2000–3000, and above 3000 m), temporal trends are estimated with the best extreme value models selected on the basis of the Akaike information criterion. Our results show that a majority of trends are decreasing below 2000 m and increasing above 2000 m. Quantitatively, we find an increase in 100-year return levels between 1959 and 2019 equal to +23 % (+32kgm-2) on average at 3500 m and a decrease of −10 % (-7kgm-2) on average at 500 m. However, for the four elevation ranges, we find both decreasing and increasing trends depending on location. In particular, we observe a spatially contrasting pattern, exemplified at 2500 m: 100-year return levels have decreased in the north of the French Alps while they have increased in the south, which may result from interactions between the overall warming trend and circulation patterns. This study has implications for natural hazard management in mountain regions.