Search

We propose the first unified objective framework (SyCLoPS) for detecting and classifying all types of low-pressure systems (LPSs) in a given data set. We use the state-of-the-art automated feature tracking software TempestExtremes (TE) to detect and track LPS features globally in ERA5 and compute 16 parameters from commonly found atmospheric variables for classification. A Python classifier is implemented to classify all LPSs at once. The framework assigns 16 different labels (classes) to each LPS data point and designates four different types of high-impact LPS tracks, including tracks of tropical cyclone (TC), monsoonal system, subtropical storm and polar low. The classification process involves disentangling high-altitude and drier LPSs, differentiating tropical and non-tropical LPSs using novel criteria, and optimizing for the detection of the four types of high-impact LPS. A comparison of our labels with those in the International Best Track Archive for Climate Stewardship (IBTrACS) revealed an overall accuracy of 95% in distinguishing between tropical systems, extratropical cyclones, and disturbances. SyCLoPS produces a better TC detection skill compared to the previous algorithms, highlighted by an approximately 6% reduction in the false alarm rate compared to the previous TE algorithm. The vertical cross section composite of the four types of high-impact LPS we detect each shows distinct structural characteristics. Finally, we demonstrate that SyCLoPS is valuable for investigating various aspects of LPSs in climate data, such as the evolution of a single LPS track, patterns of LPS frequencies, and precipitation or wind influence associated with a particular LPS class.

We evaluate tropical cyclones (TCs) in a set of thermodynamic global warming (TGW) simulations over the continental United States (CONUS). A 12 km simulation forced by ERA5 provides a 40-year historical (1980–2019) control. Four complimentary future scenarios are generated using thermodynamic deltas applied to lateral boundary, interior, and surface forcing. We curate a data set of 4,498 6-hourly TC snapshots in the control and find a corresponding “twin” in each counterfactual, permitting a paired comparison. Warming results in an increase in mean dynamical TC intensity and moisture-related quantities, with the latter being more pronounced. TC inner cores contract slightly but outer storm size remains unchanged. The frequency with which TCs become more intense is only moderately consistent, with snapshots having increased hazards ranging from 50% to 80% depending on warming level. The fractions of TCs undergoing rapid intensification and weakening both increase across all warming simulations, suggesting elevated short-term intensity variability.

Extreme rainfall events drive the amount and spatial distribution of rainfall in the Amazon and are a key driver of forest dynamics across the basin. This study investigates how the 3-hourly predictions in the High Resolution Model Intercomparison Project (HighResMIP, a component of the recent Coupled Model Intercomparison Project, CMIP6) represent extreme rainfall events at annual, seasonal, and sub-daily time scales. TRMM 3B42 (Tropical Rainfall Measuring Mission) 3 h data were used as observations. Our results showed that eleven out of seventeen HighResMIP models showed the observed association between rainfall and number of extreme events at the annual and seasonal scales. Two models captured the spatial pattern of number of extreme events at the seasonal and annual scales better (higher correlation) than the other models. None of the models captured the sub-daily timing of extreme rainfall, though some reproduced daily totals. Our results suggest that higher model resolution is a crucial factor for capturing extreme rainfall events in the Amazon, but it might not be the sole factor. Improving the representation of Amazon extreme rainfall events in HighResMIP models can help reduce model rainfall biases and uncertainties and enable more reliable assessments of the water cycle and forest dynamics in the Amazon.

The California-Nevada 1997 New Year’s flood was an atmospheric river (AR)-driven rain-on-snow (RoS) event and remains the costliest in their history. The joint occurrence of saturated soils, rainfall, and snowmelt generated inundation throughout northern California-Nevada. Although AR RoS events are projected to occur more frequently with climate change, the warming sensitivity of their flood drivers across scales remains understudied. We leverage the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) to recreate the 1997 New Year’s flood with horizontal grid spacings of 3.5 km across California, with forecast lead times of up to 4 days, and across six warming levels ranging from pre-industrial conditions to +3.5∘C. We describe the sensitivity of the flood drivers to warming including AR duration and intensity, precipitation phase, intensity and efficiency, snowpack mass and energy changes, and runoff efficiency. Our findings indicate current levels of climate change negligibly influence the flood drivers. At warming levels ≥1.7∘C, AR hazard potential increases, snowpack nonlinearly decreases, antecedent soil moisture decreases (except where the snowline retreats), and runoff decreases (except in the southern Sierra Nevada where antecedent snowpack persists). Storm total precipitation increases, but at rates below warming-induced increases in saturation-specific humidity. Warming intensifies short-duration, high-intensity rainfall, particularly where snowfall-to-rainfall transitions occur. This study highlights the nonlinear tradeoffs in 21st-century RoS flood hazards with warming and provides water management and infrastructure investment adaptation considerations.

In this work, we compare the rate of warming of summertime extreme temperatures (summer maximum value of daily maximum temperature; TXx) relative to the local mean (summer mean daily maximum temperature; TXm) over the Northern Hemisphere in observations and one set of large ensemble (LE) simulations. During the 1979–2021 historical period, observations and simulations show robust warming trends in both TXm and TXx almost everywhere in the Northern Hemisphere, except over the eastern U.S. where observations show a slight cooling trend in TXx, which may be a manifestation of internal variability. We find that the observed warming rate in TXx is significantly smaller than in TXm in North Africa, western North America, Siberia, and Eastern Asia, whereas the warming rate in TXx is significantly larger over the Eastern U.S., the U.K., and Northwestern Europe. This observed geographical pattern is successfully reproduced by the vast majority of the LE members over the historical period, and is persistent (although less intense) in future climate projections over the 2051–2100 period. We also find that these relative warming patterns are mostly driven by the local hydroclimate conditions. TXx warms slower than TXm in the hyper-arid, arid, semi-arid and moist regions, where trends in the partitioning of the turbulent surface fluxes between the latent and sensible heat flux are similar during regular and extreme hot days. In contrast, TXx warms faster than TXm in dry-subhumid regions where trends in the partitioning of the surface fluxes are significantly different between regular and extreme hot days, with a larger role of sensible heat flux during the extreme hot days.

In this study we construct and compare three different deep learning (DL) models for estimating daily snow water equivalent (SWE) from high-resolution gridded meteorological fields over the Rocky Mountain region. To train the DL models, Snow Telemetry (SNOTEL) station-based SWE observations are used as the prediction target. All DL models produce higher median Nash-Sutcliffe Efficiency (NSE) values than a conceptual SWE model and interpolated gridded data sets, although mean squared errors also tend to be higher. Sensitivity of the SWE prediction to the model's input variables is analyzed using an explainable artificial intelligence (XAI) method, yielding insight into the physical relationships learned by the models. This method reveals the dominant role precipitation and temperature play in snowpack dynamics. In applying our models to estimate SWE throughout the Rocky Mountains, an extrapolation problem arises since the statistical properties of SWE (e.g., annual maximum) and geographical properties of individual grid points (e.g., elevation) differ from the training data. This problem is solved by normalizing the SWE with its historical maximum value to alleviate extrapolation for all tested DL models. Our work shows that the DL models are promising tools for estimating SWE, and sufficiently capture relevant physical relationships to make them useful for spatial and temporal extrapolation of SWE values.

A comprehensive understanding of human-induced changes to rainfall is essential for water resource management and infrastructure design. However, at regional scales, existing detection and attribution studies are rarely able to conclusively identify human influence on precipitation. Here we show that anthropogenic aerosol and greenhouse gas (GHG) emissions are the primary drivers of precipitation change over the United States. GHG emissions increase mean and extreme precipitation from rain gauge measurements across all seasons, while the decadal-scale effect of global aerosol emissions decreases precipitation. Local aerosol emissions further offset GHG increases in the winter and spring but enhance rainfall during the summer and fall. Our results show that the conflicting literature on historical precipitation trends can be explained by offsetting aerosol and greenhouse gas signals. At the scale of the United States, individual climate models reproduce observed changes but cannot confidently determine whether a given anthropogenic agent has increased or decreased rainfall.

The water cycle is an important component of the earth system and it plays a key role in many facets of society, including energy production, agriculture, and human health and safety. In this study, the Energy Exascale Earth System Model version 1 (E3SMv1) is run with low-resolution (roughly 110 km) and high-resolution (roughly 25 km) configurations—as established by the High Resolution Model Intercomparison Project protocol—to evaluate the atmospheric and terrestrial water budgets over the conterminous United States (CONUS) at the large watershed scale. The warm season water cycle slows down in the HR experiment relative to the LR, with decreasing fluxes of precipitation, evapotranspiration, atmospheric moisture convergence, and runoff. The reductions in these terms exacerbate biases for some watersheds, while reducing them in others. For example, precipitation biases are exacerbated at HR over the Eastern and Central CONUS watersheds, while precipitation biases are reduced at HR over the Western CONUS watersheds. The most pronounced changes with resolution to the water cycle come from reductions in precipitation and evapotranspiration. The reduction in evapotranspiration reduces the biases across nearly all of the CONUS. Additional exploratory metrics show improvements to water cycle extremes (both in precipitation and streamflow), fractional contributions of different storm types to total precipitation, and mountain snowpack.

The 1997 New Year's flood event was the most costly in California's history. This compound extreme event was driven by a category 5 atmospheric river that led to widespread snowmelt. Extreme precipitation, snowmelt, and saturated soils produced heavy runoff causing widespread inundation in the Sacramento Valley. This study recreates the 1997 flood using the Regionally Refined Mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) under prescribed ocean conditions. Understanding the processes causing extreme events informs practical efforts to anticipate and prepare for such events in the future, and also provides a rich context to evaluate model skill in representing extremes. Three California-focused RRM grids, with horizontal resolution refinement of 14 km down to 3.5 km, and six forecast lead times, 28 December 1996 at 00Z through 30 December 1996 at 12Z, are assessed for their ability to recreate the 1997 flood. Planetary to synoptic scale atmospheric circulations and integrated vapor transport are weakly influenced by horizontal resolution refinement over California. Topography and mesoscale circulations, such as the Sierra barrier jet, are better represented at finer horizontal resolutions resulting in better estimates of storm total precipitation and storm duration snowpack changes. Traditional time-series and causal analysis frameworks are used to examine runoff sensitivities state-wide and above major reservoirs. These frameworks show that horizontal resolution plays a more prominent role in shaping reservoir inflows, namely the magnitude and time-series shape, than forecast lead time, 2-to-4 days prior to the 1997 flood onset.

Dynamical downscaling is a crucial process for providing regional climate information for broad uses, using coarser-resolution global models to drive higher-resolution regional climate simulations. The pool of global climate models (GCMs) providing the fields needed for dynamical downscaling has increased from the previous generations of the Coupled Model Intercomparison Project (CMIP). However, with limited computational resources, the need for prioritizing the GCMs for subsequent downscaling studies remains. GCM selection for dynamical downscaling should focus on evaluating processes relevant for providing boundary conditions to the regional models and be inspired by regional uses such as the response of extremes to changes in the boundary conditions. This leads to the need for metrics representing processes of relevance to diverse stakeholders and subregions of a domain. Procedures to account for metric redundancy and the statistical distinguishability of GCM rankings are required. Further, procedures for selecting realizations from ensembles of top-performing GCM simulations can be used to span the range of climate change signals in multiple ways. As a result, distinct weighting of metrics and prioritization of particular realizations may depend on user needs. We provide high-level guidelines for such region-specific evaluations and address how CMIP7 might enable dynamical downscaling of a representative sample of high-quality models across representative shared socioeconomic pathways (SSPs).

Developing actionable climate information and integrating it into decision-making are two crucial elements for promoting effective societal responses to climate change. However, what constitutes actionable climate information, and how it is used, varies based on the actors, systems, and scales that are relevant to specific decisions. Yet, the terms ‘actionable climate information’ or ‘use of climate information’ are used abstractly. There is a lack of holistic understanding of the various types of information that can be deemed as usable by different users, and the different ways in which they may be used in decision-making. Typologies or generalizable categorizations can help both knowledge producers and users to better envision the entire landscape of climate information and its uses and can help to reduce the time and cost of actionable knowledge production. Through systematic coding and analysis of ∼ 4 years of co-production engagements between climate scientists and resource managers, this paper presents empirically derived typologies of actionable climate information and its use, and explores whether certain uses are better informed by specific types of climate information. These typologies provide a valuable starting point for climate information producers, users, and boundary spanners working on climate-informed resource management, to reduce some of the time-intensive elements of the process.

Regional climate models can be used to examine how past weather events might unfold under different climate conditions by simulating analogue versions of those events with modified thermodynamic conditions (i.e., warming signals). Here, we apply this approach by dynamically downscaling a 40-year sequence of past weather from 1980-2019 driven by atmospheric re-analysis, and then repeating this 40-year sequence a total of 8 times using a range of time-evolving thermodynamic warming signals that follow 4 80-year future warming trajectories from 2020-2099. Warming signals follow two emission scenarios (SSP585 and SSP245) and are derived from two groups of global climate models based on whether they exhibit relatively high or low climate sensitivity. The resulting dataset, which contains 25 hourly and over 200 3-hourly variables at 12 km spatial resolution, can be used to examine a plausible range of future climate conditions in direct reference to previously observed weather and enables a systematic exploration of the ways in which thermodynamic change influences the characteristics of historical extreme events.

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 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-its-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, a 0.125 river-routing model, and 14→60¯km ocean and sea ice. By design, the computational cost of NARRM is 1/43× of the uniform low-resolution (LR) model at 100¯km but only 1/4¯10¯%-20¯% of a globally uniform high-resolution model at 25¯km. A novel hybrid time step 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 Pole 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 and 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 fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations.

Most extreme precipitation in the densely populated region of central India is produced by atmospheric vortices called monsoon lows and monsoon depressions. Here we use satellite and gauge-based precipitation estimates with atmospheric reanalyses to assess 40-year trends in the rain rates of these storms, which have remained unknown. We show that rain rates increased in the rainiest quadrant of monsoon depressions, southwest of the vortex center; precipitation decreased in eastern quadrants, yielding no clear trend in precipitation averaged over the entire storm diameter. In an atmospheric reanalysis, ascent increased in the region of amplifying precipitation, but we could not detect trends in the intensity of rotational winds around the storm center. These storm changes occurred in a background environment where humidity increased rapidly over land while warming was more muted. Monsoon lows, which we show produce less precipitation than depressions, exhibit weaker trends that are less statistically robust.

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.

In climate science and applications, the term “metric” is used to describe the distillation of complex, multifaceted evaluations to summarize the overall quality of a model simulation, or other data product, and/or as a means to quantify some response to climate change. Metrics provide insights into the fidelity of processes and outcomes from climate models and can assist with both differentiating models' representation of variables or processes and informing whether models are “fit for purpose.” Metrics can also provide a valuable reference point for co-production of knowledge between climate scientists and climate impact practitioners. Although continued metric developments enable model developers to better understand the impacts of decisions made in the model design process, metrics also have implications for the characterization of uncertainty and facilitating analyses of underlying physical processes. As a result, comprehensive evaluation with multiple metrics enhances usability of climate information by both scientific and stakeholder communities. This paper presents examples of insights gained from the development and appropriate use of metrics, and provides examples of how metrics can be used to engage with stakeholders and inform decision-making. This article is categorized under: Climate Models and Modeling > Knowledge Generation with Models The Social Status of Climate Change Knowledge > Climate Science and Decision Making Assessing Impacts of Climate Change > Evaluating Future Impacts of Climate Change.

We introduce variable-resolution enabled Community Earth System Model (VR-CESM) results simulating historical and future climate conditions at 28 km over South America and 14 km over the Andes. Three 30-year simulations are performed: a historic (1985–2014), a near future (2030–2059), and an end-century (2070–2099) simulation under the RCP8.5 scenario. Historic results compare favourably to several temperature and precipitation reanalysis products, though local biases are present, particularly during austral summer. Future simulations highlight broad warming patterns (+3–6°C by end-century) and heterogeneous precipitation responses across South America that qualitatively agree with prior modelling efforts. Our results reveal that the interaction between temperature and precipitation changes produce shifts in several Köppen–Geiger climates. Notable changes include the near-elimination of the Andean Tundra or Alpine climates, a 15% decrease in Tropical Rainforests and a Tropical Savannah expansion of 20%. To provide a regionally focused analysis of projected climate change and to illustrate the benefits of variable resolution modelling, we analyse changes in the magnitude and trend in seasonal and daily temperature and precipitation in Chile. We also examined several metrics [e.g., snow water equivalent (SWE), temperatures on wet days, and days below 0°C] to evaluate potential impacts of climate change on the Chilean cryosphere between the end-of-century and historic periods, finding wide-ranging indications of cryospheric decline. These changes are interpreted through reductions in the timing (1–2.5 months earlier peak SWE) and magnitude (200–1,000 mm SWE decreases) of water stored as snow in the Andes, a 10–30% decrease in number of cool season wet days with temperatures below 1°C, and 50–200 fewer days (annually) with minimum temperatures below 0°C. Our aim in producing a high-resolution dataset of climate projections from VR-CESM is to support analyses of climate change throughout South America but especially in vulnerable montane regions and to provide additional results for comparison with previous, ongoing, and upcoming modelling efforts.

The North American monsoon is frequently affected by transient, propagating upper tropospheric vorticity anomalies. Sometimes called Tropical Upper-Tropospheric Troughs (TUTTs), these features have been claimed to episodically enhance monsoon rainfall. Here, we track long-lived TUTTs in 40 years of reanalysis data, producing composites and case studies from 340 TUTTs which last, on average, 7 days as they move westward across the North American monsoon region. TUTTs are thought to form from midlatitude Rossby wave breaking; case studies from our dataset support this theory. TUTTs move westward within the easterly upper-level flow in which they are embedded. In vortex-centered composites along the full tracks of long-lived TUTTs, we find no detectable increase in rainfall within the main TUTT circulation. Instead, negative precipitation anomalies lie within about 500 km of the TUTT center. Quasi-geostrophic ascent occurs in the southeast quadrant of TUTTs but is confined to the upper troposphere and does not appear to interact with precipitation. Positive anomalies of ascent and rainfall occur south and southeast of TUTTs but lie outside the main TUTT vortex, perhaps indicating concurrent variations in nearby climatological precipitation maxima. In contrast with previous case studies and subjective analyses that showed TUTTs enhance precipitation in parts of northwestern Mexico, our composites along the tracks of long-lived TUTTs portray these systems, to first order, as strong vorticity anomalies trapped in the upper troposphere that interact only weakly and indirectly with precipitation.

Atmospheric rivers (ARs) are long and narrow filaments of vapor transport that are responsible for most poleward moisture transport outside of the tropics. Many AR detection algorithms have been developed to automatically identify ARs in climate data. The diversity of these algorithms has introduced appreciable uncertainties in quantitative measures of AR properties and thereby impedes the construction of a unified and internally consistent climatology of ARs. This paper compares nine global AR detection algorithms from the perspective of AR lifecycles following the propagation of ARs from origin to termination in the MERRA2 reanalysis over the period 1980–2017. Uncertainties in AR lifecycle characteristics, including event number, lifetime, intensity, and frequency distribution are discussed. Notably, the number of AR events per year in the Northern Hemisphere can vary by a factor of 5 with different algorithms. Although all algorithms show that the maximum AR origin (termination) frequency is located over the western (eastern) portion of ocean basins, significant disagreements appear in regional distribution. Spreads are large in AR lifetime and intensity. The number of landfalling AR events produced by the algorithms can vary from 16 to 80 events per year, although the agreement improves for stronger ARs. By examining the ARs' connections with the Madden-Julian Oscillation and El Niño Southern Oscillation, we find that the overall responses of ARs (such as changes in AR frequency, origin, and landfall activity) to climate variability are consistent among algorithms.

We examine the resolution dependence of errors in extreme sub-daily precipitation in available high-resolution climate models. We find that simulated extreme precipitation increases as horizontal resolution increases but that appropriately constructed model skill metrics do not significantly change. We find little evidence that simulated extreme winter or summer storm processes significantly improve with the resolution because the model performance changes identified are consistent with expectations from scale dependence arguments alone. We also discuss the implications of these scale-dependent limitations on the interpretation of simulated extreme precipitation. This article is part of a discussion meeting issue 'Intensification of short-duration rainfall extremes and implications for flash flood risks'.

Filaments of intense vapor transport called atmospheric rivers (ARs) are responsible for the majority of poleward vapor transport in the midlatitudes. Despite their importance to the hydrologic cycle, there remain many unanswered questions about changes to ARs in a warming climate. In this study we perform a series of escalating uniform SST increases (+2, +4, and +6K, respectively) in the Community Atmosphere Model version 5 in an aquaplanet configuration to evaluate the thermodynamic and dynamical response of AR vapor content, transport, and precipitation to warming SSTs. We find that AR column integrated water vapor (IWV) is especially sensitive to SST and increases by 6.3-9.7% per degree warming despite decreasing relative humidity through much of the column. Further analysis provides a more nuanced view of AR IWV changes: Since SST warming is modest compared to that in the midtroposphere, computing fractional changes in IWV with respect to SST results in finding spuriously large increases. Meanwhile, results here show that AR IWV transport increases relatively uniformly with temperature and at consistently lower rates than IWV, as modulated by systematically decreasing low-level wind speeds. Similarly, changes in AR precipitation are related to a compensatory relationship between enhanced near-surface moisture and damped vertical motions.

Atmospheric rivers (ARs) can be a boon and bane to water resource managers as they have the ability to replenish water reserves, but they can also generate million-to-billion-dollar flood damages. To investigate how anthropogenic climate change may influence AR characteristics in the coastal western United States by end century, we employ a suite of novel tools such as variable resolution in the Community Earth System Model (VR-CESM), the TempestExtremes AR detection algorithm, and the Ralph, Rutz, et al. (2019, https://doi.org/10.1175/BAMS-D-18-0023.1) AR category scale. We show that end-century ARs primarily shift from being “mostly or primarily beneficial” to “mostly or primarily hazardous” with a concomitant sharpening and intensification of winter season precipitation totals. Changes in precipitation totals are due to a significant increase in AR (+260%) rather than non-AR (+7%) precipitation, largely through increases in the most intense category of AR events and a decrease in the interval between landfalling ARs.

Future changes in tropical cyclone properties are an important component of climate change impacts and risk for many tropical and midlatitude countries. In this study we assess the performance of a multimodel ensemble of climate models, at resolutions ranging from 250 to 25 km. We use a common experimental design including both atmosphere-only and coupled simulations run over the period 1950-2050, with two tracking algorithms applied uniformly across the models. There are overall improvements in tropical cyclone frequency, spatial distribution, and intensity in models at 25 km resolution, with several of them able to represent very intense storms. Projected tropical cyclone activity by 2050 generally declines in the South Indian Ocean, while changes in other ocean basins are more uncertain and sensitive to both tracking algorithm and imposed forcings. Coupled models with smaller biases suggest a slight increase in average TC 10 m wind speeds by 2050.

Variable-resolution global climate models (VRGCMs) are a dynamical downscaling method that can reach spatiotemporal scales needed for regional climate assessments. Over the years, several users of VRGCMs have assumed where the location and extent of the refinement domain should be based on knowledge of the prevailing storm tracks and resolution dependence of important regional climate processes (e.g., atmospheric rivers [ARs] and orographic uplift), but the effect of high-resolution domain size and extent on the simulation of downstream hydroclimatic phenomena has not been systematically evaluated. Here, we use variable resolution in the Community Earth System Model (VR-CESM) to perform such a test. To do this, three VR-CESM grids were generated that span the entire, two thirds, and one third of the North Pacific and evaluated for a 30-year climatology using Atmospheric Model Intercomparison Project protocols. Simulations are compared with reanalysis products offshore (fifth-generation of the European Centre for Medium-Range Weather Forecasts atmospheric reanalysis [ERA5]) and onshore (Livneh, 2015, https://doi.org/10.1038/sdata.2015.42, and Parameter-elevation Regressions on Independent Slopes Model [PRISM]) of the western United States. The westward expansion of refinement domain influenced integrated vapor transport (IVT), which was generally high biased but minimally impacted AR characteristics. Due to slight differences in landfalling AR counts in the western United States, California winter precipitation generally improved with westward expansion of the refinement domain. Western U.S. mountain snowpack and surface temperatures were insensitive to refinement domain size and were more influenced by changes in topographic resolution and/or land surface model version. Given minimal dependence of simulated western U.S. hydroclimate on refinement domain size over the North Pacific, we advise future VR-CESM studies to focus grid resolution on better resolving land surface heterogeneity.

Wind energy production is expected to be affected by shifts in wind patterns that will accompany climate change. However, many questions remain on the magnitude and character of this impact, especially on regional scales. In this study, clustering is used to group and analyze large-scale wind patterns in California using model simulations from the variable-resolution Community Earth System Model (VR-CESM). Specifically, simulations have been produced that cover historical (1980–2000), mid-century (2030–2050), and end-of-century (2080–2100) time periods. Once clustered, observed changes to wind patterns can be analyzed in terms of both the change in frequency of those clusters and changes to winds within-clusters. Statistically significant capacity factors changes have been found at all five wind plant sites. Decomposition of the capacity factor changes into frequency changes and within-cluster changes enables a better understanding of their drivers. A further examination of the synoptic-scale fields associated with each cluster then provides a better understanding of how changes to large-scale meteorological fields are important for driving changes in localized wind speeds.

Until recently, the El Niño–Southern Oscillation (ENSO) was considered a reliable source of winter precipitation predictability in the western US, with a historically strong link between extreme El Niño events and extremely wet seasons. However, the 2015–2016 El Niño challenged our understanding of the ENSO-precipitation relationship. California precipitation was near-average during the 2015–2016 El Niño, which was characterized by warm sea surface temperature (SST) anomalies of similar magnitude compared to the extreme 1997–1998 and 1982–1983 El Niño events. We demonstrate that this precipitation response can be explained by El Niño’s spatial pattern, rather than internal atmospheric variability. In addition, observations and large-ensembles of regional and global climate model simulations indicate that extremes in seasonal and daily precipitation during strong El Niño events are better explained using the ENSO Longitude Index (ELI), which captures the diversity of ENSO’s spatial patterns in a single metric, compared to the traditional Niño3.4 index, which measures SST anomalies in a fixed region and therefore fails to capture ENSO diversity. The physically-based ELI better explains western US precipitation variability because it tracks the zonal shifts in tropical Pacific deep convection that drive teleconnections through the response in the extratropical wave-train, integrated vapor transport, and atmospheric rivers. This research provides evidence that ELI improves the value of ENSO as a predictor of California’s seasonal hydroclimate extremes compared to traditional ENSO indices, especially during strong El Niño events.

Atmospheric rivers (ARs) are now widely known for their association with high-impact weather events and long-term water supply in many regions. Researchers within the scientific community have developed numerous methods to identify and track of ARs—a necessary step for analyses on gridded data sets, and objective attribution of impacts to ARs. These different methods have been developed to answer specific research questions and hence use different criteria (e.g., geometry, threshold values of key variables, and time dependence). Furthermore, these methods are often employed using different reanalysis data sets, time periods, and regions of interest. The goal of the Atmospheric River Tracking Method Intercomparison Project (ARTMIP) is to understand and quantify uncertainties in AR science that arise due to differences in these methods. This paper presents results for key AR-related metrics based on 20+ different AR identification and tracking methods applied to Modern-Era Retrospective Analysis for Research and Applications Version 2 reanalysis data from January 1980 through June 2017. We show that AR frequency, duration, and seasonality exhibit a wide range of results, while the meridional distribution of these metrics along selected coastal (but not interior) transects are quite similar across methods. Furthermore, methods are grouped into criteria-based clusters, within which the range of results is reduced. AR case studies and an evaluation of individual method deviation from an all-method mean highlight advantages/disadvantages of certain approaches. For example, methods with less (more) restrictive criteria identify more (less) ARs and AR-related impacts. Finally, this paper concludes with a discussion and recommendations for those conducting AR-related research to consider.

This study provides an overview of the coupled high-resolution Version 1 of the Energy Exascale Earth System Model (E3SMv1) and documents the characteristics of a 50-year-long high-resolution control simulation with time-invariant 1950 forcings following the HighResMIP protocol. In terms of global root-mean-squared error metrics, this high-resolution simulation is generally superior to results from the low-resolution configuration of E3SMv1 (due to resolution, tuning changes, and possibly initialization procedure) and compares favorably to models in the CMIP5 ensemble. Ocean and sea ice simulation is particularly improved, due to better resolution of bathymetry, the ability to capture more variability and extremes in winds and currents, and the ability to resolve mesoscale ocean eddies. The largest improvement in this regard is an ice-free Labrador Sea, which is a major problem at low resolution. Interestingly, several features found to improve with resolution in previous studies are insensitive to resolution or even degrade in E3SMv1. Most notable in this regard are warm bias and associated stratocumulus deficiency in eastern subtropical oceans and lack of improvement in El Niño. Another major finding of this study is that resolution increase had negligible impact on climate sensitivity (measured by net feedback determined through uniform +4K prescribed sea surface temperature increase) and aerosol sensitivity. Cloud response to resolution increase consisted of very minor decrease at all levels. Large-scale patterns of precipitation bias were also relatively unaffected by grid spacing.

The accurate characterization of seasonal and inter-annual site-level wind energy variability is essential during wind project development. Understanding the features and probability of low-wind years is of particular interest to developers and financers. However, a dearth of long-term, hub-height wind observations makes these characterizations challenging, and thus techniques to improve these characterizations are of great value. To improve resource characterization, we explicitly link wind resource variability (at hub-height, and at specific sites) to regional and synoptic scale wind regimes. Our approach involves statistical clustering of high-resolution modeled wind data, and is applied to California for a period covering 1980–2015. With this approach, we investigate the unique meteorological patterns driving low and high wind years at five separate wind project sites. We also find wind regime changes over the 36-year period consistent with global warming: wind regimes associated with anomalously hot summer days increased at half a day per year and stagnant conditions increased at one-third days per year. Despite these changes, the average annual resource potential remained constant at all project sites. Additionally, we identify correlations between climate modes and wind regime frequency, a linkage valuable for resource characterization and forecasting. Our general approach can be applied in any location and may benefit many aspects of wind energy resource evaluation and forecasting.

The mountains of the Western United States provide a vital natural service through the storage and release of mountain snowpack, lessening impacts of seasonal aridity and satiating summer water demand. However, climate change continues to undermine these important processes. To understand how snowpack may change in the headwaters of California's major reservoirs, the North American Coordinated Regional Climate Downscaling Experiment is analyzed to assess peak water volume, peak timing, accumulation rate, melt rate, and snow season length across both latitudinal and elevational gradients. Under a high-emissions scenario, end-of-century peak snowpack timing occurs 4 weeks earlier and peak water volume is 79.3% lower. The largest reductions are above Shasta, Oroville, and Folsom and between 0- and 2,000-m elevations. Regional climate model and global forcing data set choice is important in determining historical snowpack character, yet by end century all models show a significant and similar decline in mountain snowpack.

Shifting wind patterns are an expected consequence of global climate change, with direct implications for wind energy production. However, wind is notoriously difficult to predict, and significant uncertainty remains in our understanding of climate change impacts on existing wind generation capacity. In this study, historical and future wind climatology and associated capacity factors at five wind turbine sites in California are examined. Historical (1980-2000) and mid-century (2030-2050) simulations were produced using the Variable-Resolution Community Earth System Model (VR-CESM) to understand how these wind generation sites are expected to be impacted by climate change. A high-resolution statistically downscaled WRF product provided by DNV GL, reanalysis datasets MERRA-2, CFSR, NARR, and observational data were used for model validation and comparison. These projections suggest that wind power generation capacity throughout the state is expected to increase during the summer, and decrease during fall and winter, based on significant changes at several wind farm sites. This study improves the characterization of uncertainty around the magnitude and variability in space and time of California's wind resources in the near future, and also enhances our understanding of the physical mechanisms related to the trends in wind resource variability.