skip to main content

Title: Probabilistic Evaluation of Drought in CMIP6 Simulations

As droughts have widespread social and ecological impacts, it is critical to develop long‐term adaptation and mitigation strategies to reduce drought vulnerability. Climate models are important in quantifying drought changes. Here, we assess the ability of 285 CMIP6 historical simulations, from 17 models, to reproduce drought duration and severity in three observational data sets using the Standardized Precipitation Index (SPI). We used summary statistics beyond the mean and standard deviation, and devised a novel probabilistic framework, based on the Hellinger distance, to quantify the difference between observed and simulated drought characteristics. Results show that many simulations have less thanerror in reproducing the observed drought summary statistics. The hypothesis that simulations and observations are described by the same distribution cannot be rejected for more thanof the grids based on ourdistance framework. No single model stood out as demonstrating consistently better performance over large regions of the globe. The variance in drought statistics among the simulations is higher in the tropics compared to other latitudinal zones. Though the models capture the characteristics of dry spells well, there is considerable bias in low precipitation values. Good model performance in terms of SPI does not imply good more » performance in simulating low precipitation. Our study emphasizes the need to probabilistically evaluate climate model simulations in order to both pinpoint model weaknesses and identify a subset of best‐performing models that are useful for impact assessments.

« less
 ;  ;  ;  ;  ;  
Award ID(s):
1839336 1839441
Publication Date:
Journal Name:
Earth's Future
DOI PREFIX: 10.1029
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Weakly nonlinear, bi‐periodic patterns of waves that propagate in the‐direction with amplitude variation in the‐direction are generated in a laboratory. The amplitude variation in the‐direction is studied within the framework of the vector (vNLSE) and scalar (sNLSE) nonlinear Schrödinger equations using the uniform‐amplitude, Stokes‐like solution of the vNLSE and the Jacobi elliptic sine function solution of the sNLSE. The wavetrains are generated using the Stokes‐like solution of vNLSE; however, a comparison of both predictions shows that while they both do a reasonably good job of predicting the observed amplitude variation in, the comparison with the elliptic function solution of the sNLSE has significantly less error when the ratio of‐wavenumber to the two‐dimensional wavenumber is less than about 0.25. For ratios between about 0.25 and 0.30 (the limit of the experiments), the two models have comparable errors. When the ratio is less than about 0.17, agreement with the vNLSE solution requires a third‐harmonic term in the‐direction, obtained from a Stokes‐type expansion of interacting, symmetric wavetrains. There is no evidence of instability growth in the‐direction, consistent with the work of Segur and colleagues, who showed that dissipationmore »stabilizes the modulational instability. Finally, there is some extra amplitude variation in, which is examined via a qualitative stability calculation that allows symmetry breaking in that direction.

    « less
  2. Abstract

    Silicon stable isotope ratios (30Si) of over 150 stream water samples were measured during seven storm events in six small critical zone observatory (CZO) catchments spanning a wide range in climate (sub‐humid to wet, tropical) and lithology (granite, volcanic, and mixed sedimentary). Here we report a cross‐site analysis of this dataset to gain insight into stream30Si variability across low‐order catchments and to identify potential climate (i.e., runoff), hydrologic, lithologic, and biogeochemical controls on observed stream Si chemical and isotopic signatures. Event‐based30Si exhibit variability both within and across sites (−0.22‰ to +2.27‰) on the scale of what is observed globally in both small catchments and large rivers. Notably, each site shows distinct30Si signatures that are preserved even after normalization for bedrock composition. Successful characterization of observed cross‐site behavior requires the merging of two distinct frameworks in a novel combined model describing both non‐uniform fluid transit time distributions and multiple fractionating pathways in application to low‐order catchments. The combined model reveals that site‐specific architecture (i.e., biogeochemical reaction pathways and hydrologic routing) regulates stream silicon export signatures even when subject to extreme precipitation events.

  3. Abstract

    Meteorological drought indices like the Standardized Precipitation Evaporation Index (SPEI) are frequently used to diagnose “ecological drought,” despite the fact that they were not explicitly designed for this purpose. More recently developed indices like the Evaporative Stress Index (ESI), which is based on the degree of coupling between actual to potential evapotranspiration, may better capture dynamic plant response to moisture limitations. However, the skill of these indices at describing plant water stress is rarely evaluated at sub‐seasonal timescales over which drought evolves. Moreover, it remains unclear how variability in phenological timing impacts and complicates early drought detection. Here, we compared the ability of ESI and SPEI to reflect the dynamics of ecological drought in forests and grasslands, based on anomalies of Gross Primary Productivity (GPP), surface conductance (Gs, a proxy for stomatal conductance), soil moisture, and vapor pressure deficit. ESI performed better than SPEI in capturing the dynamics of GPP andGs, but still missed early ecological drought signals due to biases linked to earlier onset of spring leaf development. Thus, we developed a modified variant of the ESI () that accounts for the complicating effects of phenological shifts in leaf area index (LAI). Thedetected droughtmore »onset up to 7–10 weeks earlier than SPEI and ESI. Additionally, drought onset dates determined fromare close to (±2 weeks) the dates determined from LAI‐corrected anomalies ofGs, and GPP, as well as the onset dates of soil water deficit and atmospheric aridity.

    « less
  4. Abstract

    Terrestrial lightning frequently serves as a loss mechanism for energetic electrons in the Van Allen radiation belts, leading to lightning‐induced electron precipitation (LEP). Regardless of the specific causes, energetic electron precipitation from the radiation belts in general has a significant influence on the ozone concentration in the stratosphere and mesosphere. The atmospheric chemical effects induced by LEP have been previously investigated using subionospheric VLF measurements at Faraday station, Antarctica (65.25°S, 64.27°W,L= 2.45). However, there exist large variations in the precipitation flux, ionization production, and occurrence rate of LEP events depending on the peak current of the parent lightning discharge, as well as the season, location, and intensity of the thunderstorm activity. These uncertainties motivate us to revisit the calculation of atmospheric chemical changes produced by LEP. In this study, we combine a well‐validated LEP model and first‐principles atmospheric chemical simulation, and investigate three intense storms in the year of 2013, 2015, and 2017 at the magnetic latitude of 50., 32., and 35., respectively. Modeling results show that the LEP events in these storms can cumulatively drive significant changes in the,, andconcentration in the mesosphere. These changes aremore »as high as,, andat 75–85 km altitude, respectively, and comparable to the effects typically induced by other types of radiation belt electron precipitation events. Considering the high occurrence rate of thunderstorms around the globe, the long‐term global chemical effects produced by LEP events need to be properly quantified.

    « less
  5. Abstract

    We present a statistical investigation of the effects of interplanetary magnetic field (IMF) on hemispheric asymmetry in auroral currents. Nearly 6 years of magnetic field measurements from Swarm A and C satellites are analyzed. Bootstrap resampling is used to remove the difference in the number of samples and IMF conditions between the local seasons and the hemispheres. Currents are stronger in Northern Hemisphere (NH) than Southern Hemisphere (SH) for IMF Bin NH (Bin SH) in most local seasons under both signs of IMF B. For Bin NH (Bin SH), the hemispheric difference in currents is small except in local winter when currents in NH are stronger than in SH. During Band Bin NH (Band Bin SH), the largest hemispheric asymmetry occurs in local winter and autumn, when the NH/SH ratio of field aligned current (FAC) is 1.180.09 in winter and 1.170.09 in autumn. During Band Bin NH (Band Bin SH), the largest asymmetry is observed in local autumn with NH/SH ratio of 1.160.07 for FAC. We alsomore »find an explicit Beffect on auroral currents in a given hemisphere: on average Bin NH and Bin SH causes larger currents than vice versa. The explicit Beffect on divergence‐free current during IMF Bis in very good agreement with the Beffect on the cross polar cap potential from the Super Dual Auroral Radar Network dynamic model except at SH equinox and NH summer.

    « less