Search NASA⌕ Search

SEARCH · Search NASA

Results for “Return Period”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On the uncertainty of long-period return values of extreme daily precipitation

Methods for calculating return values of extreme precipitation and their uncertainty are compared using daily precipitation rates over the Western U.S. and Southwestern Canada from a large ensemble of climate model simulations. The roles of return-value estimation procedures and sample size in uncertainty are evaluated for various return periods. We compare two different generalized extreme value (GEV) parameter estimation techniques, namely L-moments and maximum likelihood (MLE), as well as empirical techniques. Even for very large datasets, confidence intervals calculated using GEV techniques are narrower than those calculated using empirical methods. Furthermore, the more efficient L-moments parameter estimation techniques result in narrower confidence intervals than MLE parameter estimation techniques at small sample sizes, but similar best estimates. It should be noted that we do not claim that either parameter fitting technique is better calibrated than the other to estimate long period return values. While a non-stationary MLE methodology is readily available to estimate GEV parameters, it is not for the L-moments method. Comparison of uncertainty quantification methods are found to yield significantly different estimates for small sample sizes but converge to similar results as sample size increases. Finally, practical recommendations about the length and size of climate model ensemble simulations and the choice of statistical methods to robustly estimate long period return values of extreme daily precipitation statistics and quantify their uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Return period of extreme rainfall substantially decreases under 1.5 °C and 2.0 °C warming: a case study for Uttarakhand, India

In June 2013, Uttarakhand experienced a hydro-meteorological disaster due to a 4 d extreme precipitation event of return period more than 100 years, claiming thousands of lives and causing enormous damage to infrastructure. Using the weather@home climate modelling system and its Half a degree Additional warming, Prognosis and Projected Impacts simulations, this study investigates the change in the return period of similar events in a 1.5 °C and 2 °C warmer world, compared to current and pre-industrial levels. We find that the likelihood of such extreme precipitation events will significantly increase under both future scenarios. We also estimate the change in extreme river flow at the Ganges; finding a considerable increase in the risk of flood events. Our results also suggest that until now, anthropogenic aerosols may have effectively counterbalanced the otherwise increased meteorological flood risk due to greenhouse gas (GHG) induced warming. Disentangling the response due to GHGs and aerosols is required to analyses the changes in future rainfall in the South Asia monsoon region. More research with other climate models is also necessary to make sure these results are robust.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluation and projection of long period return values of extreme daily precipitation in the CMIP5 and CMIP6 models

Using a non-stationary Generalized Extreme Value statistical method, we calculate selected extreme daily precipitation indices and their 20 year return values from the CMIP5 and CMIP6 climate models over the historical and future periods. We evaluate model performance of these indices and their return values in replicating similar quantities calculated from multiple gridded observational products. Difficulties in interpreting model quality in the context of observational uncertainties are discussed. Projections are framed in terms of specified global warming target temperatures rather than at specific times and under specific emissions scenarios. The change in framing shifts projection uncertainty due to differences in model climate sensitivity from the values of the projections to the timing of the global warming target. At their standard resolutions, we find there are no meaningful differences between the two generations of models in their quality or projections of simulated extreme daily precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Future Intensity‐Duration‐Frequency Curves of Extreme Precipitation in the Midwest United States From Convection‐Permitting Modeling

Abstract During the last four decades, global warming has statistically significant intensified extreme precipitation events in the Midwestern United States (defined here as the region covering Illinois, Indiana, Ohio, and Kentucky), leading to increased risks to human life, property, and infrastructure. To enable climate change adaptation and resilience across various economic and social sectors in this region, updated information about future climate changes, specifically at finer spatial scales, is essential. Leveraging a new 150‐year dynamical downscaling data set at convection‐permitting resolution, this study introduces a framework to construct the projected future intensity‐duration‐frequency (IDF) curves of heavy precipitation, which are prominent tools for infrastructure design and water resources management. This framework generates IDF curves at both sub‐daily and multi‐day duration utilizing hourly in situ observations as well as quantile‐based statistical techniques in bias‐correction and return levels selection. The assumption of non‐stationarity in the distribution parameter fitting process is also implemented in this workflow. Compared to historical IDF curves for 1980–2022, future projected IDF curves for 2058–2100 under Representative Concentration Pathway (RCP) 4.5 and RCP 8.5 scenarios indicate an average intensity increase of approximately 15% and 25%, respectively, across 74 stations, considering both annual and seasonal timescales. Future projections suggest that extreme precipitation events may become more severe across six investigated return periods, with longer return periods showing a greater increase. The frequency of future extreme precipitation events in the Midwest region is also projected to double. Furthermore, current results reveal spatial heterogeneity of future trends across stations owing to the high‐resolution input data set. Plain Language Summary This study investigates the evolving nature of extreme precipitation events in the Midwestern United States under a changing climate. By leveraging a high‐resolution dynamical downscaling data set, we construct projected intensity‐duration‐frequency (IDF) curves for future extreme rainfall events. These curves serve as vital tools for infrastructure planning and water resource management. Our analysis reveals a significant increase in both the intensity and frequency of extreme precipitation events in the region. Future projected IDF curves for the late century indicate an average intensity increase of approximately 15%–25% compared to historical values. Moreover, the frequency of such events is expected to double. Spatial heterogeneity in future trends is observed across different stations within the Midwest, highlighting the importance of high‐resolution modeling in capturing localized climate variability. These findings underscore the urgent need for climate adaptation strategies to mitigate the increasing risks associated with extreme precipitation events in the region. Key Points This study introduces a workflow to construct future intensity‐duration‐frequency (IDF) curves over the Midwest United States using a new convection‐permitting modeling data set The current IDF construction workflow reproduces well the historical observed IDF 30 curves in summer months with median relative errors of 2.4% among 74 stations and 6 investigated durations The projected IDF curves show diverse future trends of extreme precipitation across stations, with intensity increases of approximately 15% and 25% under RCP4.5 and RCP8.5 climate scenarios, respectively, and a doubling of frequency on average

Nguyen, Trung↗
DOE OSTI2025

Probabilistic Hazard Assessment for Tornadoes, Straight-Line Wind, and Extreme Precipitation at the Savannah River Site

Recent data sets for three meteorological phenomena with the potential to inflict damage on SRS facilities – tornadoes, straight-line winds, and heavy precipitation – are analyzed using appropriate statistical techniques to estimate the occurrence probabilities for these events in the future. Summaries of the results for DOE-mandated return periods and comparisons to similar calculations performed in 2013 by Werth et al. (W2013) are given. Using tornado statistics for i) the combined states of Georgia and South Carolina, and ii) a 2⁰ square area surrounding SRS, we calculated the probability per year of any location at SRS being struck by a tornado (the ‘strike’ probability) and the probability that any point will experience winds above set thresholds. The strike probability was calculated to be 7.04E-4 (1 chance in 1420) per year and tornadic wind speeds for DOE mandated return periods of 50,000 years (corresponding to wind design category 3 (WDC-3), and 125,000 years (meeting WDC-4) (USDOE, 2016) were estimated to be 132 mph and 147 mph, respectively. By contrast, default tornado wind speeds taken from ANSI/ANS-2.3-2011 are somewhat higher: 161 mph for return periods of 50,000 years and 173 mph every 125,000 years (ANS, 2011). Although the ANS and the SRS evaluation used the same basic model (Ramsdell and Rishel, 2007), the region defined in ANS 2.3 that encompasses the SRS also includes areas of the Great Plains and lower Midwest, regions with much higher occurrence frequencies of strong tornadoes. The SRS straight-line wind values associated with various return periods were calculated by fitting existing wind data to a GEV1 distribution and extrapolating the values for any return period from the tail of that function. For the DOE mandated return periods, we expect straight-line winds of 117 mph every 2500 years (the required WDC-3 standard) and 125 mph every 6250 years (WDC-4) at any point within the SRS. These values are similar to those from the ANS-2.3-2011 report, which has wind speeds of 125mph and 133 mph for return periods of 2500 years and 6250 years, respectively. For extreme precipitation, we compared the fits of two different theoretical extreme-value distributions and applied the one that fit the data best for each of several accumulation periods. The DOE mandated 6-hr accumulated rainfall for return periods of 10,000 years (corresponding to precipitation design category 3 (PDC-3) and 25,000 years (PDC-4) were estimated as 9.1 inches and 10.1 inches, respectively. For the 24-hr rainfall return periods of 10,000 years and 25,000 years, total rainfall estimates were 12.02 inches and 13.17 inches, respectively, higher than comparable values provided in the W2013 report.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Climatology and Extreme Value Analysis of Large Hail in China

In this paper, hail size records are analyzed at 2254 stations in China and a hail size climatology is developed based on gridded hail observations for the period 1960–2015. It is found that the annual percentiles of hail size records changed sharply and national-wide after 1980, therefore two periods, 1960–79 and 1980–2015, are studied. There are some similarities between the two periods in terms of the characteristics of hail size such as the spatial distribution patterns of mean annual maximum hail size and occurrence week of annual maximum hail size. The 1980–2015 period had higher observation density than the 1960–79 period, but showed smaller mean annual maximum hail size, especially in northern China. In the majority of grid boxes, the annual maximum hail size experienced a decreasing trend during the 1980–2015 period. A Gumbel extreme value model is fitted to each grid box to estimate the return periods of maximum hail size. The scale and location parameter of the fitted Gumbel distributions are higher in eastern China than in western China, thereby reflecting a greater likelihood of large hail in eastern China. In southern China, the maximum hail size exceeds 127 mm for a 10-yr return period, whereas in northern China maximum hail size exceeds this threshold for a 50-yr return period. Furthermore, the Gumbel model is found to potentially underestimate the maximum hail size for certain return periods, but provides a more informed picture of the spatial distribution of extreme hail size and the regional differences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An Open‐Source, Physics‐Based, Tropical Cyclone Downscaling Model With Intensity‐Dependent Steering

Abstract An open‐source, physics‐based tropical cyclone (TC) downscaling model is developed, in order to generate a large climatology of TCs. The model is composed of three primary components: (a) a random seeding process that determines genesis, (b) an intensity‐dependent beta‐advection model that determines the track, and (c) a non‐linear differential equation set that determines the intensification rate. The model is entirely forced by the large‐scale environment. Downscaling ERA5 reanalysis data shows that the model is generally able to reproduce observed TC climatology, such as the global seasonal cycle, genesis locations, track density, and lifetime maximum intensity distributions. Inter‐annual variability in TC count and power‐dissipation is also well captured, on both basin‐wide and global scales. Regional TC hazard estimated by this model is also analyzed using return period maps and curves. In particular, the model is able to reasonably capture the observed return period curves of landfall intensity in various sub‐basins around the globe. The incorporation of an intensity‐dependent steering flow is shown to lead to regionally dependent changes in power dissipation and return periods. Advantages and disadvantages of this model, compared to other downscaling models, are also discussed.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

A Non‐Intrusive Machine Learning Framework for Debiasing Long‐Time Coarse Resolution Climate Simulations and Quantifying Rare Events Statistics

Abstract Due to the rapidly changing climate, the frequency and severity of extreme weather is expected to increase over the coming decades. As fully‐resolved climate simulations remain computationally intractable, policy makers must rely on coarse‐models to quantify risk for extremes. However, coarse models suffer from inherent bias due to the ignored “sub‐grid” scales. We propose a framework to non‐intrusively debias coarse‐resolution climate predictions using neural‐network (NN) correction operators. Previous efforts have attempted to train such operators using loss functions that match statistics. However, this approach falls short with events that have longer return period than that of the training data, since the reference statistics have not converged. Here, the scope is to formulate a learning method that allows for correction of dynamics and quantification of extreme events with longer return period than the training data. The key obstacle is the chaotic nature of the underlying dynamics. To overcome this challenge, we introduce a dynamical systems approach where the correction operator is trained using reference data and a coarse model simulation nudged toward that reference. The method is demonstrated on debiasing an under‐resolved quasi‐geostrophic model and the Energy Exascale Earth System Model (E3SM). For the former, our method enables the quantification of events that have return period two orders longer than the training data. For the latter, when trained on 8 years of ERA5 data, our approach is able to correct the coarse E3SM output to closely reflect the 36‐year ERA5 statistics for all prognostic variables and significantly reduce their spatial biases.

Barthel Sorensen, B.↗
DOE OSTI2023

Machine Learning of Key Variables Impacting Extreme Precipitation in Various Regions of the Contiguous United States

Abstract Amplification in extreme precipitation intensity and frequency can cause severe flooding and impose significant social and economic consequences. Variations in extreme precipitation intensity, frequencies, and return periods can be attributed to many physical variables across spatial and temporal scales. Here we employ ensemble machine learning (ML) methods, namely random forest (RF), eXtreme Gradient Boosting (XGB), and artificial neural networks (ANN), to explore key contributing variables to monthly extreme precipitation intensity and frequency in six regions over the United States. We further establish emulators for return periods. Results show that the ML models for intensity perform better in regions with obvious seasonality (i.e., Northern Great Plains, Southern Great Plains, and West Coast) than the other three regions (Northeast, Southwest, and Rocky Mountains), while for frequency the models perform well for most regions. The Shapley additive explanation is used to help explain the relationships between extreme precipitation characteristics and identify top variables for RF and XGB. We find that latent heat flux, relative humidity, soil moisture, and large‐scale subsidence are key common variables across the regions for both monthly intensity and frequency, and their compound effects are non‐negligible. The developed ML models capture the probability and return period of extreme precipitation well for all regions and may be used for decision making (e.g., infrastructure planning and design).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Technoeconomic Assessment of Phosphoric Acid and Rare Earth Element Recovery from Phosphoric Acid Sludge

Sustainability faces many challenges, including the availability of materials necessary for technological advancement. Rare earth elements (REEs), for example, are key materials for several manufacturing industries that can unlock renewable energy and sustainable development. In this study, a decanter centrifuge has been employed to successfully separated phosphoric acid and REE-containing particles from phosphoric acid sludge with concentrations ranging from 1000 to 2200 ppm REEs. Operating efficiently with up to 35 wt.% solids, the centrifuge was demonstrated to achieve approximately 95% phosphoric acid recovery and 90% REE recovery in a single pass, eliminating the need for additional processing steps. This breakthrough supports a proposed rare earth oxide (REO) recovery process integrating phosphoric acid (PA), elemental phosphorus (P4), and REO into two potential pathways: PA-REO and PA-P4-REO. These processes aim to reintroduce recovered phosphoric acid into the main product to significantly increase output and revenue. Post-separation, phosphorus-rich particles can be converted to P4, while REE-containing solids undergo further treatment including acid leaching, extraction/stripping, precipitation, and calcination to produce a marketable REO material. Technoeconomic analysis indicates promising profitability, with the PA-REO process showing a delta net present value (ΔNPV) of USD 441.8 million over a 12-year period and expected return within a year of construction, while the PA-P4-REO process yields a ΔNPV of USD 178.7 million over a 12-year return period. Both pathways offer robust financial prospects and demonstrate the feasibility of commercial-scale REO recovery from phosphoric acid sludge, offering an economically feasible approach to produce REEs for future sustainable development challenges related to sustainability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantifying the compound hazard of freezing rain and wind gusts across CONUS

Abstract The co-occurrence of freezing rain, ice accumulation and wind gusts (FZG) poses a significant hazard to infrastructure and transportation. However, quantification of the frequency and intensity of FZG is challenged by the lack of direct icing measurements. In this work, we evaluate and then apply an energy balance model to high-frequency data collected during 2005–2022 to derive hourly ice accumulation at 883 stations across the contiguous USA. These estimates are combined with wind gust observations to compute time series of hourly FZG hazard magnitude using the Sperry–Piltz Ice Accumulation (SPIA) index. Results are evaluated using Storm Reports of damage and economic disruption. The hourly SPIA estimates are also used to (i) derive a geospatial atlas of the hazard including the 50 yr return period event intensities for each US state derived using superstations, and (ii) describe storylines of significant events in terms of meteorological drivers and socioeconomic impacts. The highest values of SPIA during the 18 yr study period occur in a region extending from the Southern Great Plains over the Midwest into the densely populated Northeast. States in these regions also have high 50 yr return period maximum radial ice accumulation of 3–5 cm and co-occurring wind gusts >30 ms −1 . These values are comparable to past estimates for the 500 yr event which may imply this hazard has been previously underestimated. This atlas can be used to inform optimal FZG hazard mitigation strategies for each state/region.

58 GEOSCIENCES↗
DOE OSTI2021

Application of Point Precipitation Frequency Estimates to Watersheds

This report documents work sponsored by the U.S. Nuclear Regulatory Commission (NRC) at the Oak Ridge National Laboratory (ORNL) as part of the RES project, “Application of Point Precipitation Frequency Estimates to Watersheds.” This project was implemented as part of the Probabilistic Flood Hazard Assessment (PFHA) Research Program. The objective of the PFHA Research Program is to develop tools and guidance on the use of PFHA methods to risk-inform NRC’s licensing of new facilities as well as licensing and oversight of currently operating facilities as they relate to flooding hazards. Many nuclear power plants (NPPs) are located on or near rivers so riverine flooding hazards need to be considered in their design and operation. Probabilistic riverine flood models are important tools for realistic assessment of flooding risks. However, these models require areal estimates of the depth, duration, and frequency of rainfall distributed over the watershed, which are not often available. Point precipitation frequency estimates are more widely available. For example, the National Oceanic and Atmospheric Administration (NOAA) has published NOAA Atlas 14, which provides point precipitation frequency estimates for 5-minute through 60-day durations at average recurrence intervals of 1-year through 1,000-year. The research documented in this report addresses areal reduction factors (ARFs), which can be used to convert the widely available point precipitation frequency estimates, to estimates of areal precipitation frequency over a watershed. The most widely used ARF source is Technical Paper 29 (TP-29) published by the then U.S. Weather Bureau in 1958. However, both the methods and the underlying precipitation data used to produce TP-29 are seriously out of date. For example, due to the small gauge network available at the time of TP-29’s compilation, ARF estimates developed are only for watersheds smaller than about 400 square miles. Due to the relatively short record lengths of precipitation data available, frequency considerations could not be accurately determined. Other factors such as regional climate and seasonality were not addressed. Several newer methods have been published since TP-29 was developed and both the type and quantity of precipitation data have increased significantly, along with computational resources and analytical tools such as geographic information systems. This report reviewed and assessed the available precipitation products and methods for conducting ARF analysis. The work applied up-to-date precipitation data products and analysis methods with a novel watershed-based approach to investigate how ARF estimates vary across different methods, data sources, geographical locations, return periods, and seasons. The overall findings reported here regarding basic ARF trends are in line with other recent studies showing that ARFs decrease with increasing area, increase with increasing duration, and decrease with increasing return period. This study found significant differences among the available ARF methods. This work also found a strong geographical variability across different US hydrologic regions, suggesting that the ARF are specific to regional climate patterns and geographical characteristics and should not be applied arbitrarily to other locations. The results also reveal the importance of data record length, especially for high return level ARFs. The work reported in NUREG/CR-7271 will assist NRC staff in assessing different classes of ARF methods in conjunction with available rainfall data sets. It will also support the development of guidance for application of point precipitation data in PFHAs. It should be noted that the ARF values presented in this report for any location or region were developed for the purposes of comparing methods and investigating the factors that influence ARFs. They should not be considered official and should not be used in leu of a site-specific analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Gulf of Mexico Offshore Wind Energy Hurricane Risk Assessment

NREL's feasibility assessment of offshore wind in the Gulf of Mexico concluded that hurricane risk was one of the major challenges that would need to be overcome for a mature offshore wind industry to develop in the Gulf of Mexico. To ensure the robust design of wind turbines in the Gulf of Mexico, it is critical to understand the added risk posed by the threat of major hurricanes, as those affecting the Gulf of Mexico region have a significant potential to exceed design limits prescribed by the International Electrotechnical Commission (IEC) wind design standards. To satisfy this charge, this project defines the wind hazard for the Gulf of Mexico Offshore Wind Energy area using the hurricane hazard model develop by Applied Research Associates and published extensively in the open literature. In doing so, the return periods associated with the IEC Class 1A and Typhoon Class limit-state hurricanes are estimated on a grid with nominal resolution of 10 km to determine where hurricane risk results in the exceedance of the IEC design criteria. On the same grid, wind speeds hazard contours associated with return periods varying from 50 to 1,000 years are also estimated. An additional challenge in assessing hurricane wind speed risk in the Gulf of Mexico arises from inconsistent terminology across the Saffir-Simpson hurricane scale and the IEC design criteria. Saffir-Simpson definitions are based on 1-minute sustained wind speeds estimated at 10-m height over marine terrain, while the IEC uses a different averaging period (3-second versus 1-minute) and reference height (assumed herein a hub height of 150 m versus 10 m). Employing the latest research on turbulence characteristics of the hurricane boundary layer, conversions between various durations (e.g., 3-seconds, 1-minute, 10-minutes, 1-hour) and between elevations near the surface (10 m) to near hub height (assumed herein 150 m) are developed. IEC Class 1A and Typhoon Class limit states are also provided in terms of an equivalent Saffir-Simpson hurricane wind speed category.

17 WIND ENERGY↗
DOE OSTI2026

CONUS-wide Projected Flood Frequency and Uncertainty Estimates, Version 1.0

This dataset presents a large-ensemble of CONUS-wide projected flood frequency and uncertainty estimates across ~2.7 million NHDPlusV2 river reaches over the CONUS. The framework producing this dataset leverages a multi-model, uncertainty-aware modeling framework that allows evaluating shifts in flood frequences at the stream reach level across the CONUS. CONUS-wide ensemble streamflow projections generated from hydrologic simulations driven by downscaled and bias-corrected Coupled Model Intercomparison Project Phase 6 (CMIP6) outputs are used to derive these flood frequency and uncertainty estimates over the period 1980 - 2099. A spatially consistent regional L-moment algorithm is applied across clusters defined by the US Hydrologic Unit Code Subregions (HUC4s and HUC8s) and NHDPlusV2 stream orders to estimate flood frequencies. The dataset also includes at-site based flood estimates that allow for the comparison between local and regional approach-based estimates, assess projected changes, and characterize their uncertainties. For more reliable estimation of rare flood frequencies such as 500 and 1000-year return periods, super-ensemble based estimates are also included in the dataset. This dataset is derived to support the "Impact-Informed Dam Safety Risk Assessment for Securing Hydropower Assests" project for the US Department of Energy (DOE) Hydropower and Hydrokinetic Office (H2O). For further details, refer to Kao et al. (2022), Ghimire et al. (2023), Ghimire et al. (2025), and Hosking and Wallis (1997).

Ghimire, Ganesh [ORNL] (ORCID:0000000242843941)↗
DOE OSTI2021

Extreme Wind and Waves in U.S. East Coast Offshore Wind Energy Lease Areas

The Outer Continental Shelf along the U.S. east coast exhibits abundant wind resources and is now a geographic focus for offshore wind deployments. This analysis derives and presents expected extreme wind and wave conditions for the sixteen lease areas that are currently being developed. Using the homogeneous ERA5 reanalysis dataset it is shown that the fifty-year return period wind speed (U50) at 100 m a.s.l. in the lease areas ranges from 29.2 to 39.7 ms -1 . After applying corrections to account for spectral smoothing and averaging period, the associated pseudo-point U50 estimates are 34 to 46 ms -1 . The derived uncertainty in U50 estimates due to different distributional fitting is smaller than the uncertainty associated with under-sampling of the interannual variability in annual maximum wind speeds. It is shown that, in the northern lease areas, annual maximum wind speeds are generally associated with intense extratropical cyclones rather than cyclones of tropical origin. Extreme wave statistics are also presented and indicate that the 50-year return period maximum wave height may substantially exceed 15 m. From this analysis, there is evidence that annual maximum wind speeds and waves frequently derive from the same cyclone source and often occur within a 6 h time interval.

16 TIDAL AND WAVE POWER↗
DOE OSTI2019

The El Niño event of 2015–2016: climate anomalies and their impact on groundwater resources in East and Southern Africa

The impact of climate variability on groundwater storage has received limited attention despite widespread dependence on groundwater as a resource for drinking water, agriculture and industry. Here, we assess the climate anomalies that occurred over Southern Africa (SA) and East Africa, south of the Equator (EASE), during the major El Niño event of 2015–2016, and their associated impacts on groundwater storage, across scales, through analysis of in situ groundwater piezometry and Gravity Recovery and Climate Experiment (GRACE) satellite data. At the continental scale, the El Niño of 2015–2016 was associated with a pronounced dipole of opposing rainfall anomalies over EASE and Southern Africa, north–south of ~12°S, a characteristic pattern of the El Niño–Southern Oscillation (ENSO). Over Southern Africa the most intense drought event in the historical record occurred, based on an analysis of the cross-scale areal intensity of surface water balance anomalies (as represented by the standardised precipitation evapotranspiration index – SPEI), with an estimated return period of at least 200 years and a best estimate of 260 years. Climate risks are changing, and we estimate that anthropogenic warming only (ignoring changes to other climate variables, e.g. precipitation) has approximately doubled the risk of such an extreme SPEI drought event. These surface water balance deficits suppressed groundwater recharge, leading to a substantial groundwater storage decline indicated by both GRACE satellite and piezometric data in the Limpopo basin. Conversely, over EASE during the 2015–2016 El Niño event, anomalously wet conditions were observed with an estimated return period of ~10 years, likely moderated by the absence of a strongly positive Indian Ocean zonal mode phase. The strong but not extreme rainy season increased groundwater storage, as shown by satellite GRACE data and rising groundwater levels observed at a site in central Tanzania. We note substantial uncertainties in separating groundwater from total water storage in GRACE data and show that consistency between GRACE and piezometric estimates of groundwater storage is apparent when spatial averaging scales are comparable. These results have implications for sustainable and climate-resilient groundwater resource management, including the potential for adaptive strategies, such as managed aquifer recharge during episodic recharge events.

54 ENVIRONMENTAL SCIENCES↗