Search

African easterly waves (AEWs) exert significant influence on local and downstream high-impact weather including tropical cyclone (TC) genesis over the Atlantic. Accurate representation of AEWs in climate and weather prediction models therefore is necessary for skillful predictions. In this study, we examine simulated AEWs, including their evolution, vertical structure, and linkage to tropical cyclone genesis, in the NASA Goddard Earth Observing System Model, version 5 (GEOS-5) atmospheric global climate model. Identified by the leading empirical orthogonal function mode of time-filtered precipitation, the observed westward propagating AEWs along the southern track over the Atlantic are largely captured in GEOS-5, but with a slower phase speed and significantly weaker amplitude downstream off the West Africa coast. The weak downstream development of AEWs in GEOS-5 is accompanied by much reduced TC genesis over the main development region. Further analyses suggest that the slow westward propagation and weaker AEW amplitude downstream can be ascribed to a weak African easterly jet, while overestimated negative (positive) meridional potential vorticity (PV) gradients appear to the north (south) of 10°N in GEOS-5. The greatly overestimated positive meridional PV gradient to the south of 10°N is expected to generate strong horizontal stretching in the AEW wave pattern in the model, which hinders organization of convection and its feedback to sustain the AEW development. Persistent and vigorous AEW precipitation in the Guinea Highlands of the West Africa coast could also be responsible for reduced westward propagation of AEWs in the model.

Since the infamous extreme drought of the 1960s, the climate of the northeastern United States (NEUS) has generally trended toward warmer and wetter conditions. Nonetheless, there is mounting evidence that short-term droughts will continue to pose a significant risk for this region. To better explore the processes governing events such as these, climate models have adopted more complex representations of the fully coupled atmosphere–land–ocean–sea ice system; however, large uncertainties in future projections still persist, with internal variability necessitating large ensembles to understand trends in both rare and high-impact extreme events such as rapidly developing droughts (a term here that includes flash droughts developing on monthly scales). In this study, seven large ensemble (LE) models are employed to answer the outstanding question: How are the frequency and character of drought in the NEUS changing under a warming climate? We find that most LE models indicate the NEUS will experience a long-term wetting trend with more “extremely wet” months, but also more frequent short-term extreme droughts. These changes are associated with increasing precipitation, atmospheric water demand, and climate variability. We also conclude that discrepant trends in precipitation and evapotranspiration variability will lead to increasing anticorrelation of these variables, which is relevant to the intensification of rapidly developing drought, particularly in the spring season. These changes are associated with an increase in evapotranspiration from plants, brought by an earlier emergence of the growing season and denser vegetation. Significance Statement Droughts are extreme events with the potential to produce considerable social and economic damage. Because droughts emerge slowly and are relatively infrequent, large datasets are needed to study these features. Using multiple climate models, each producing multiple long-duration simulations, along with a novel drought index, we characterize rapidly developing droughts (those flash droughts identifiable from monthly data) and understand their drivers in the northeastern United States. We find that short-term extreme droughts are projected to be more frequent, while rapidly developing droughts will become more intense in the spring season. Intensification of rapidly developing droughts is attributed to differences in patterns of precipitation and evapotranspiration (soil evaporation and plant transpiration), which will be magnified by an extended growing season and an increase in vegetation.

Societies and ecosystems within and downstream of mountains rely on seasonal snowmelt to satisfy their water demands. Anthropogenic climate change has reduced mountain snowpacks worldwide, altering snowmelt magnitude and timing. Here the global warming level leading to widespread and persistent mountain snowpack decline, termed low-to-no snow, is estimated for the world’s most latitudinally contiguous mountain range, the American Cordillera. We show that a combination of dynamical, thermodynamical and hypsometric factors results in an asymmetric emergence of low-to-no-snow conditions within the midlatitudes of the American Cordillera. Low-to-no-snow emergence occurs approximately 20 years earlier in the southern hemisphere, at a third of the local warming level, and coincides with runoff efficiency declines (8% average) in both dry and wet years. The prevention of a low-to-no-snow future in either hemisphere requires the level of global warming to be held to, at most, +2.5 °C.

Tropical cyclones (TCs) and their associated precipitation can have devastating impacts on the areas affected, with outcomes ranging from mudslides to inland flash flooding. Previous studies have used a fixed radius around the TC to isolate storm-related precipitation. One previous study instead used a dynamic radius of 8 m s−1 winds, but the wind field of the TC can deteriorate or shift quickly after landfall or the onset of extratropical transition (ET). This study uses a dynamical radius derived from the 500-hPa geopotential height in and around the TC to define TC- and post-tropical cyclone (PTC)-related heavy precipitation, allowing for the analysis of precipitation with tropical origins after the official demise of the original TC. Climatologies are constructed, indicating a maximum in TC- and PTC-related heavy precipitation in the west North Pacific and a secondary maximum in the east North Pacific. PTC-related heavy precipitation accounts for as much as 40% of the annual heavy precipitation in the northwest portion of the west North Pacific basin and 3.13% of heavy precipitation globally. We observe that the major hurricane stage contributes on average 2.6% of the global TC- and PTC-related precipitation, while the less intense but more common tropical storm stages of the TC life cycle contribute 85.7% of this observed precipitation. This analysis framework can be further extended to assess model biases and climate projections of TC and PTC precipitation.

Cyclonic low‐pressure systems (LPS) are the dominant synoptic‐scale rain-bearing system of the South Asian summer monsoon. Traditionally categorized by intensity as monsoon lows, monsoon depressions, and more intense cyclonic storms, LPS produce intense rainfall and floods in some of the world’s most densely populated regions. Yet the contribution of the relatively weak lows vs. the stronger depressions to extreme rainfall and its trends remains unknown; this knowledge gap is particularly troubling because historical trends in LPS have been difficult to assess due to changes in the observing network. Future projections have also remained highly uncertain due to the inability of many coarse-resolution climate models to accurately simulate LPS.Here we use satellite and gauge-based precipitation estimates with atmospheric reanalyses to show that precipitation in monsoon depressions has become more intense in recent decades. This intensification has occurred as humidity over parts of India increased more rapidly than nearly anywhere else on Earth. Precipitation in depressions has risen at a relative rate larger than that of specific humidity, suggesting that upward motion in depressions has become more intense; vertical motion trends in a state-of-the-art reanalysis, which incorporates nearly all long-term climate forcings, are consistent with this hypothesis. We also examine changes in South Asian LPS precipitation simulated by an ensemble of high-resolution global models, which we find skillfully represent these storms. Future trends in total LPS precipitation, including in monsoon depressions, lie near an approximate Clausius–Clapeyron rate (7%/K) in the multi-model mean. This change in LPS rain rates contributes to a projected future increase in seasonal mean and extreme precipitation over South Asian land. Adaptation to future changes in human exposure to hydrological extremes thus requires careful monitoring, accurate multi-decadal projections, and skilful short-term forecasts of the interaction of the humidity field with the dynamics of monsoon LPS.

This paper describes and analyzes the Reed-Jablonowski (RJ) tropical cyclone (TC) test case used in the 2016 Dynamical Core Model Intercomparison Project (DCMIP2016). The intermediate complexity test case analyzes the evolution of a weak vortex into a TC in an idealized tropical environment. Simulations from 9 general circulation models (GCMs) that participated in DCMIP2016 are analyzed in this study at 50 km horizontal grid spacing, with 5 of these models also providing simulations at 25 km grid spacing for an analysis on the impact of finer grid spacing. Evolution of minimum surface pressure (MSP) and maximum 1 km azimuthally averaged wind speed (MWS), the wind-pressure relationship, radial profiles of wind speed and surface pressure, and wind composites are documented for all participating GCMs at both horizontal grid spacings. While results are generally similar between all models, some GCMs reach significantly higher storm intensities than others, ultimately impacting specific characteristics of their horizontal and vertical structure. TCs simulated at 25 km grid spacings retained these differences, but reach higher intensities and are more compact than their 50 km counterparts. These results indicate dynamical core choice is an essential factor in GCM development, and future work should be conducted to explore how specific differences within the dynamical core affect TC behavior in GCMs.

Monotonicity is an important property of remapping operators for coupled weather and climate models. However, it is often challenging to design highly accurate operators that avoid the generation of new extrema or keep a remapped field between physically prescribed bounds. To that end, this paper explores several traditional and novel approaches for both conservative and non-conservative monotone remapping on the sphere. The accuracy and effectiveness of these algorithms are evaluated in the context of several different real and idealized fields and meshes.

Despite the emerging influence of anthropogenic climate change on the global water cycle, at regional scales the combination of observational uncertainty, large internal variability, and modeling uncertainty undermine robust statements regarding the human influence on precipitation. Here, we use output from global climate models in a perfect-data sense to develop a framework for conducting regional detection and attribution (D&A) for precipitation, starting with the contiguous United States (CONUS) where observational uncertainty is lower than in other regions. Our unified approach can simultaneously detect systematic trends in mean and extreme precipitation, attribute trends to anthropogenic forcings, compute the effects of forcings as a function of time, and map the effects of individual forcings. Model output is used to conduct a set of tests that yield a parsimonious representation for characterizing seasonal precipitation over the CONUS for the historical record (1900 to present day), which ensures our D&A is insensitive to structural uncertainty. Our framework is developed using synthetic data in a Pearl-causal perspective wherein causality can be identified using intervention-based simulations. While the hypothesis-based framework and accompanying generalized D&A formula we develop should be widely applicable, we include a strong caution that the hypothesis-guided simplification of the formula for the historical climatic record of CONUS as described in this paper will likely fail to hold in other geographic regions and under future warming.

This study analyzes the quality of simulated historical precipitation across the contiguous United States (CONUS) in a 12 km Weather Research and Forecasting model version 4.2.1 (WRF v 4.2.1)-based dynamical downscaling of the fifth-generation ECMWF atmospheric reanalysis (ERA5). This work addresses the following questions. First, how well are the 3 and 24 h precipitation characteristics (diurnal and annual cycles, precipitation frequency, annual and seasonal mean and maximum precipitation, and distribution of seasonal maximum precipitation) represented in the downscaled simulation, compared to ERA5? And second, how does the performance of the simulated WRF precipitation vary across seasons, regions, and timescales? Performance is measured against the National Centers for Environmental Prediction/Environmental Modeling Center (NCEP/EMC) 4 km Stage IV and Oregon State University Parameter-Elevation Regressions on Independent Slopes Model (PRISM) data on 3 and 24 h timescales, respectively. Our analysis suggests that the 12 km WRF exhibits biases typically found in other WRF simulations, including those at convection-permitting scales. In particular, WRF simulates both the timing and magnitude of the summer diurnal precipitation peak as well as ERA5 over most of the CONUS, except for a delayed diurnal peak over the Great Plains. As compared to ERA5, both the month and the magnitude of the precipitation peak annual cycle are remarkably improved in the downscaled WRF simulation. WRF slightly overestimates 3 and 24 h precipitation maximum over the CONUS, in contrast to ERA5, which generally underestimates these quantities mainly over the eastern half of the CONUS. Notably, WRF better captures the probability density distribution (PDF) of 3 and 24 h annual and seasonal maximum precipitation. WRF exhibits seasonally dependent precipitation biases across the CONUS, while ERA5's biases are relatively consistent year round over most of the CONUS. These results suggest that dynamical downscaling to a higher resolution improves upon some precipitation metrics but is susceptible to common regional climate model biases. Consequently, if used as input data for domain-specific models, we suggest moderate bias correction be applied to the dynamically downscaled product.

The present work evaluates historical precipitation and its indices defined by the Expert Team on Climate Change Detection and Indices (ETCCDI) in suites of dynamically and statistically downscaled regional climate models (RCMs) against NOAA’s Global Historical Climatology Network Daily (GHCN-Daily) dataset over Florida. The models examined here are: (1) nested RCMs involved in the North American CORDEX (NA-CORDEX) program, (2) variable resolution Community Earth System Models (VR-CESM), (3) Coupled Model Intercomparison Project phase 5 (CMIP5) models statistically downscaled using localized constructed analogs (LOCA) technique. To quantify observational uncertainty, three in situ-based (PRISM, Livneh, CPC) and three reanalysis (ERA5, MERRA2, NARR) datasets are also evaluated against the station data. The reanalyses and dynamically downscaled RCMs generally underestimate the magnitude of the monthly precipitation and the frequency of the extreme rainfall in summer. The models forced with CanESM2 miss the phase of the seasonality of extreme precipitation. All models and reanalyses severely underestimate both the mean and interannual variability of mean wet-day precipitation (SDII), consecutive dry days (CDD), and overestimate consecutive wet days (CWD). Metric analysis suggests large uncertainty across NA-CORDEX models. Both the LOCA and VR-CESM models perform better than the majority of models. Overall, RegCM4 and WRF models perform poorer than the median model performance. The performance uncertainty across models is comparable to that in the reanalyses. Specifically, NARR performs poorer than the median model performance in simulating the mean indices and MERRA2 performs worse than the majority of models in capturing the interannual variability of the indices.

The E3SM Diagnostics Package (E3SM Diags) is a modern, Python-based Earth system model (ESM) evaluation tool (with Python module name e3sm_diags), developed to support the Department of Energy (DOE) Energy Exascale Earth System Model (E3SM). E3SM Diags provides a wide suite of tools for evaluating native E3SM output, as well as ESM data on regular latitude–longitude grids, including output from Coupled Model Intercomparison Project (CMIP) class models. E3SM Diags is modeled after the National Center for Atmospheric Research (NCAR) Atmosphere Model Working Group (AMWG, 2022) diagnostics package. In its version 1 release, E3SM Diags included a set of core essential diagnostics to evaluate the mean physical climate from model simulations. As of version 2.7, more process-oriented and phenomenon-based evaluation diagnostics have been implemented, such as analysis of the quasi-biennial oscillation (QBO), the El Niño–Southern Oscillation (ENSO), streamflow, the diurnal cycle of precipitation, tropical cyclones, ozone and aerosol properties. An in situ dataset from DOE's Atmospheric Radiation Measurement (ARM) program has been integrated into the package for evaluating the representation of simulated cloud and precipitation processes. This tool is designed with enough flexibility to allow for the addition of new observational datasets and new diagnostic algorithms. Additional features include customizable figures; streamlined installation, configuration and execution; and multiprocessing for fast computation. The package uses an up-to-date observational data repository maintained by its developers, where recent datasets are added to the repository as they become available. Finally, several applications for the E3SM Diags module were introduced to fit a diverse set of use cases from the scientific community.

As the resolution of global Earth system models increases, regional-scale evaluations are becoming ever more important. This study presents a framework for quantifying precipitation distributions at regional scales and applies it to evaluate Coupled Model Intercomparison Project (CMIP) 5 and 6 models. We employ the Intergovernmental Panel on Climate Change (IPCC) sixth assessment report (AR6) climate reference regions over land and propose refinements to the oceanic regions based on the homogeneity of precipitation distribution characteristics. The homogeneous regions are identified as heavy-, moderate-, and light-precipitating areas by K-means clustering of Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (GPM) version 6 final run product (IMERG) precipitation frequency and amount distributions. With the global domain partitioned into 62 regions, including 46 land and 16 ocean regions, we apply 10 established precipitation distribution metrics. The collection includes metrics focused on the maximum peak, lower 10th percentile, and upper 90th percentile in precipitation amount and frequency distributions; the similarity between observed and modeled frequency distributions; an unevenness measure based on cumulative amount; average total intensity on all days with precipitation; and number of precipitating days each year. We apply our framework to 25 CMIP5 and 41 CMIP6 models, as well as six observation-based products of daily precipitation. Our results indicate that many CMIP5 and 6 models substantially overestimate the observed light-precipitation amount and frequency, as well as the number of precipitating days, especially over midlatitude regions outside of some land regions in the Americas and Eurasia. Improvement from CMIP5 to 6 is shown in some regions, especially in midlatitude regions, but it is not evident globally, and over the tropics most metrics point toward degradation.

This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled Regionally Refined Model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations – a first-of-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100 km atmosphere and land, 0.125° river routing model, and 14→60 km ocean and sea ice. By design, the computational cost of NARRM is ∼3x of the uniform low-resolution (LR) model at 100 km but only ∼10–20 % of a globally uniform high-resolution model at 25 km. A novel hybrid timestep strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North polar region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better resolved topography in NARRM, whereas those over ocean are attributable to the improved air-sea interactions with finer grids for both atmosphere and ocean/sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS, and the warm-season land-atmosphere coupling at the Southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the RRM framework for a standard Earth system model release and high-resolution climate production simulations.

TempestExtremes (TE) is a multifaceted framework for feature detection, tracking, and scientific analysis of regional or global Earth system datasets on either rectilinear or unstructured/native grids. Version 2.1 of the TE framework now provides extensive support for examining both nodal (i.e., pointwise) and areal features, including tropical and extratropical cyclones, monsoonal lows and depressions, atmospheric rivers, atmospheric blocking, precipitation clusters, and heat waves. Available operations include nodal and areal thresholding, calculations of quantities related to nodal features such as accumulated cyclone energy and azimuthal wind profiles, filtering data based on the characteristics of nodal features, and stereographic compositing. This paper describes the core algorithms (kernels) that have been added to the TE framework since version 1.0, including algorithms for editing pointwise trajectory files, composition of fields around nodal features, generation of areal masks via thresholding and nodal features, and tracking of areal features in time. Several examples are provided of how these kernels can be combined to produce composite algorithms for evaluating and understanding common atmospheric features and their underlying processes. These examples include analyzing the fraction of precipitation from tropical cyclones, compositing meteorological fields around extratropical cyclones, calculating fractional contribution to poleward vapor transport from atmospheric rivers, and building a climatology of atmospheric blocks.

Intensified extreme precipitation and resulting flooding are among the most impactful consequences of climate change, especially over the northeastern US (NEUS). To project and understand the impacts of climate change (or related climate perturbations) on extreme weather events as they may occur in the future, the Pseudo-Global Warming (PGW) method has been employed with great success. However, it has never been ascertained to what degree the conclusions from PGW studies are sensitive to the design of the PGW experiment. Consequently, three key questions related to the application of the PGW method remain unanswered: At what spatial scale should climate perturbations be applied? Among the different meteorological variables available, which should be perturbed? And will PGW projections vary significantly under different experiment designs? To begin to address these questions, we examine the sensitivity and robustness of conclusions drawn from the PGW method over NEUS by conducting multiple PGW experiments. The results show that the projections of precipitation and other essential variables are consistent at the regional mean scale, with a relative difference of much less than 10\\%; however, different experimental designs nonetheless cause significant displacements among storm events. Several previously assumed advantages of modifying temperature at the regional mean scale do not appear to hold, such as the preservation of geostrophic balance. Also, we find the regional mean perturbation produces a positive precipitation bias due to overestimated warming over the ocean. Overall, PGW experiments with perturbations from temperature or the combination of temperature and wind at the gridpoint scale are both recommended, depending on the research target. The first approach can isolate the spatially-dependent thermodynamic impact, and the latter incorporates both the thermodynamic and dynamic impacts.

Traditional multimodel methods for estimating future changes in precipitation intensity, duration, and frequency (IDF) curves rely on mean or median of models’ IDF estimates. Such multimodel estimates are impaired by large estimation uncertainty, shadowing their efficacy in planning efforts. Here, assuming that each climate model is one representation of the underlying data generating process, i.e., the Earth system, we propose a novel extension of current methods through pooling model data: (i) evaluate performance of climate models in simulating the spatial and temporal variability of the observed annual maximum precipitation (AMP), (ii) bias-correct and pool historical and future AMP data of reasonably performing models, and (iii) compute IDF estimates in a nonstationary framework from pooled historical and future model data. Pooling enhances fitting of the extreme value distribution to the data and assumes that data from reasonably performing models represent samples from the “true” underlying data generating distribution. Through Monte Carlo simulations with synthetic data, we show that return periods derived from pooled data have smaller biases and lesser uncertainty than those derived from ensembles of individual model data. We apply this method to NA-CORDEX models to estimate changes in 24-h precipitation intensity–frequency (PIF) estimates over the Susquehanna watershed and Florida peninsula. Our approach identifies significant future changes at more stations compared to median-based PIF estimates. The analysis suggests that almost all stations over the Susquehanna and at least two-thirds of the stations over the Florida peninsula will observe significant increases in 24-h precipitation for 2–100-yr return periods.

Tropical cyclones (TCs) can subject an area to heavy precipitation for many hours, or even days, worsening the risk of flooding, which creates dangerous conditions for residents of the U.S. East and Gulf Coasts. To study the representation of TC-related precipitation over the eastern United States in current-generation global climate models, a novel analysis methodology is developed to track TCs and extract their associated precipitation using an estimate of their dynamical outer size. This methodology is applied to three variable-resolution (VR) configurations of the Community Atmosphere Model, version 5 (CAM5), with high-resolution domains over the North Atlantic and one low-resolution conventional configuration, as well as to a combination of reanalysis and observational precipitation data. Metrics and diagnostics such as TC counts, intensities, outer storm sizes, and annual mean total and extreme precipitation are compared between the CAM5 simulations and reanalysis/observations. The high-resolution VR configurations outperform the global low-resolution configuration for all variables in the North Atlantic. Realistic TC intensities are produced by the VR configurations. The total North Atlantic TC counts are lower than observations but better than reanalysis.