Search

Atmospheric rivers (ARs) are extreme weather events that can alleviate drought or cause billions of US dollars in flood damage. By transporting significant amounts of latent energy towards the poles, they are crucial to maintaining the climate system's energy balance. Since there is no first-principle definition of an AR grounded in geophysical fluid mechanics, AR identification is currently performed by a multitude of expert-defined, threshold-based algorithms. The variety of AR detection algorithms has introduced uncertainty into the study of ARs, and the thresholds of the algorithms may not generalize to new climate datasets and resolutions. We train convolutional neural networks (CNNs) to detect ARs while representing this uncertainty; we name these models ARCNNs. To detect ARs without requiring new labeled data and labor-intensive AR detection campaigns, we present a semi-supervised learning framework based on image style transfer. This framework generalizes ARCNNs across climate datasets and input fields. Using idealized and realistic numerical models, together with observations, we assess the performance of the ARCNNs. We test the ARCNNs in an idealized simulation of a shallow-water fluid in which nearly all the tracer transport can be attributed to AR-like filamentary structures. In reanalysis and a high-resolution climate model, we use ARCNNs to calculate the contribution of ARs to meridional latent heat transport, and we demonstrate that this quantity varies considerably due to AR detection uncertainty.

The deadly and destructive nature of tropical cyclones (TCs) makes understanding their response to future climate change of the utmost importance. TC genesis hinges on multiple factors, including an initial disturbance. African easterly waves (AEWs) have been shown to serve as such disturbances for TCs developing in the North Atlantic. It is therefore crucial to understand the relationship between AEWs and TCs and how this relationship may be affected by climate change. In this study, we examine the AEW-TC relationship in historical and future climates using three models from the HighResMIP PRIMAVERA simulations. The AEWs and TCs were tracked in the model data using objective tracking algorithms, and AEW and TC tracks were then matched together if they were close to each other in space and time. The strength of the AEWs was measured using the eddy kinetic energy and the curvature vorticity of the waves. TC strength and intensity were measured using potential intensity and lifetime maximum 10 m windspeed. We found that future changes in the frequency of AEWs are not a good indicator of future TC activity. However, AEW strength, as well as environmental conditions conducive to strong TCs, are good indicators of AEWs that develop into TCs in both historical and future climates.

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.

Historical trends in monsoon low pressure systems (LPS), the dominant rain-bearing weather system of South Asia, have been difficult to assess due to changes in the observing network. Future projections have also remained uncertain because prior studies concluded that many coarse-resolution climate models do not accurately simulate LPS. Here, we examine changes in South Asian monsoon LPS simulated by an ensemble of global models, including some with high spatial resolution, that we show skillfully represent LPS. In the ensemble mean, the number of strong LPS (monsoon depressions) decreased over the last 65 years (1950–2014) by about 15% while no trend was detected for weaker LPS (monsoon lows). The reduction in depression counts then moderated, yielding no trend in the periods 1980–2050 or 2015–2050. The ensemble mean projects a shift in genesis from ocean to land and an increase in LPS precipitation of at least 7% K−1, which together contribute to a projected increase in seasonal mean and extreme precipitation over central India.

Atmospheric rivers (ARs) are intensive poleward moisture transport events that are essential to the global hydrological cycle and are often linked to extreme weather events. We categorize the winter North Pacific ARs into two “flavors”: wind-dominated (windy ARs) and moisture-dominated (wet ARs) using 40 years of hourly data from fifth generation of the European Centre for Medium-Range Weather Forecasts Interim Reanalysis. We compare the differences between windy ARs and wet ARs including the lifecycle characteristics (such as genesis locations and changes of meteorological elements through the lifecycle), overall AR frequency, landfall impacts, and variability. The windy ARs are more likely to occur in the midlatitudes, while wet ARs are more active in the subtropics. Windy ARs are associated with intensive surface pressure lows, where the strong pressure gradient can support the strong wind within ARs. Due to larger size and longer lifetime, wet ARs are more likely to produce more precipitation over a lifecycle. By scaling the landfalling ARs, we show that wet ARs dominate the high-category ARs (Category 4 and 5) with higher spatial frequency and more precipitation, and windy ARs have higher contributions in the lower AR categories especially over British Columbia. Windy ARs are modulated by El Niño Southern Oscillation (ENSO) teleconnections via the anomalous geopotential height and extended subtropical jet. Wet ARs are affected by the anomalous sea surface temperature over the midlatitudes related to ENSO. Sensitivity analysis with an alternate AR detection algorithm shows consistent results on AR flavors but with disagreement on the amplitude.

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.

A central question surrounding possible human exploration of Mars is whether crewed missions can be supported by available technologies using in situ resources. Here, we show that photovoltaics-based power systems would be adequate and practical to sustain a crewed outpost for an extended period over a large fraction of the planet’s surface. Climate data were integrated into a radiative transfer model to predict spectrally-resolved solar flux across the Martian surface. This informed detailed balance calculations for solar cell devices that identified optimal bandgap combinations for maximizing production capacity over a Martian year. We then quantified power systems, manufacturing, and agricultural demands for a six-person mission, which revealed that photovoltaics-based power generation would require (Formula presented.) 10 t of carry-along mass, outperforming alternatives over ∼50% of Mars’ surface.

Anthropogenic climate change is decreasing seasonal snowpacks globally, with potentially catastrophic consequences on water resources, given the long-held reliance on snowpack in water management. In this Review, we examine the changes and trickle-down impacts of snow loss in the western United States (WUS). Across the WUS, snow water equivalent declines of ~25% are expected by 2050, with losses comparable with contemporary historical trends. There is less consensus on the time horizon of snow disappearance, but model projections combined with a new low-to-no snow definition suggest ~35–60 years before low-to-no snow becomes persistent if greenhouse gas emissions continue unabated. Diminished and more ephemeral snowpacks that melt earlier will alter groundwater and streamflow dynamics. The direction of these changes are difficult to constrain given competing factors such as higher evapotranspiration, altered vegetation composition and changes in wildfire behaviour in a warmer world. These changes undermine conventional WUS water management practices, but through proactive implementation of soft and hard adaptation strategies, there is potential to build resilience to extreme, episodic and, eventually, persistent low-to-no snow conditions. Federal investments offer a timely opportunity to address these vulnerabilities, but they require a concerted portfolio of activities that cross historically siloed physical and disciplinary boundaries.

This study aims to quantify the impacts of atmosphere–ocean coupling on simulated and projected tropical cyclone (TC) precipitation globally. We used global climate model (GCM) simulations over 1950–2050 from the High Resolution Model Intercomparison Project (HighResMIP) and compared its fully coupled atmosphere–ocean GCMs (AOGCMs) with atmosphere-only GCMs (AGCMs). We find that ocean coupling generally leads to decreased TC precipitation over ocean and land. Large-scale sea surface temperature (SST) biases are critical drivers of the precipitation difference, with secondary contributions from local TC–ocean feedbacks via SST cold wakes. The two driving factors, attributed to ocean coupling in the AOGCMs, influence TC precipitation in association with decreased TC intensity and specific humidity. The AOGCMs and AGCMs consistently project TC precipitation increases in 2015–2050 relative to 1950–2014 over ocean for all basins, and for landfalling TCs in the North Atlantic and western North Pacific.

Constructing an accurate, continental, in-situ-based, kilometer-scale, long-term record of the precipitation field and its spatiotemporal changes remains a significant challenge. Here, we determine the extreme-value behavior of the NEXRAD Stage IV radar-based quantitative precipitation estimate. We find that the climatology of 5-year daily return values in the contiguous United States East of the Rocky Mountains shows only slight variability on spatial scales smaller than (Formula presented.) 100 km. In light of this finding, we test whether rain-gauge-only daily precipitation data sets can produce accurate extreme-value behavior at spatial scales finer than the spacing between gauges. We find that the 5-year daily return values are accurate at locations far from rain gauges only if the interpolation between gauges is carried out appropriately for extremes. Precipitation statistics derived from in-situ rain gauge data are therefore of sufficient spatial resolution to faithfully capture daily extremes over much of the eastern United States.

A recent study found statistically significant differences in extreme precipitation distributions over the contiguous United States (CONUS) when changing the microphysics scheme in a superparameterized global climate model. Here, we repeat the analysis globally and similarly find that differences are widespread when varying the number of predicted moments in the microphysics parameterization, but not when comparing variants of the double-moment scheme. However, contrary to the previous study in which differences largely disappeared over CONUS when 5-day simulations were conducted, we found that the signal in these shorter integrations remains within the tropics, implying a direct local effect of microphysics on precipitation extremes in these regions. The effect on precipitation is traced back to changes in vertical velocity profiles changes that are then amplified in the climatological simulations compared to the 5-day ones. Finally, the superparameterized extremes, regardless of the microphysics scheme, are shown to be larger than those from the Global Precipitation Climatology Project One-Degree Daily data set and generally smaller than those from the Tropical Rainfall Measuring Mission 3B42 data set.

While various studies explore the relationship between individual sources of climate variability and extreme precipitation, there is a need for improved understanding of how these physical phenomena simultaneously influence precipitation in the observational record across the contiguous United States. In this work, we introduce a single framework for characterizing the historical signal (anthropogenic forcing) and noise (natural variability) in seasonal mean and extreme precipitation. An important aspect of our analysis is that we simultaneously isolate the individual effects of seven modes of variability while explicitly controlling for joint inter-mode relationships. Our method utilizes a spatial statistical component that uses in situ measurements to resolve relationships to their native scales; furthermore, we use a data-driven procedure to robustly determine statistical significance. In Part I of this work we focus on natural climate variability: detection is mostly limited to DJF and SON for the modes of variability considered, with the El Niño/Southern Oscillation, the Pacific–North American pattern, and the North Atlantic Oscillation exhibiting the largest influence. Across all climate indices considered, the signals are larger and can be detected more clearly for seasonal total versus extreme precipitation. We are able to detect at least some significant relationships in all seasons in spite of extremely large (> 95%) background variability in both mean and extreme precipitation. Furthermore, we specifically quantify how the spatial aspect of our analysis reduces uncertainty and increases detection of statistical significance while also discovering results that quantify the complex interconnected relationships between climate drivers and seasonal 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'.

The social cost of methane (SC-CH4) measures the economic loss of welfare caused by emitting one tonne of methane into the atmosphere. This valuation may in turn be used in cost-benefit analyses or to inform climate policies1-3. However, current SC-CH4 estimates have not included key scientific findings and observational constraints. Here we estimate the SC-CH4 by incorporating the recent upward revision of 25 per cent to calculations of the radiative forcing of methane4, combined with calibrated reduced-form global climate models and an ensemble of integrated assessment models (IAMs). Our multi-model mean estimate for the SC-CH4 is US$933 per tonne of CH4 (5-95 per cent range, US$471-1,570 per tonne of CH4) under a high-emissions scenario (Representative Concentration Pathway (RCP) 8.5), a 22 per cent decrease compared to estimates based on the climate uncertainty framework used by the US federal government5. Our ninety-fifth percentile estimate is 51 per cent lower than the corresponding figure from the US framework. Under a low-emissions scenario (RCP 2.6), our multi-model mean decreases to US$710 per tonne of CH4. Tightened equilibrium climate sensitivity estimates paired with the effect of previously neglected relationships between uncertain parameters of the climate model lower these estimates. We also show that our SC-CH4 estimates are sensitive to model combinations; for example, within one IAM, different methane cycle sub-models can induce variations of approximately 20 per cent in the estimated SC-CH4. But switching IAMs can more than double the estimated SC-CH4. Extending our results to account for societal concerns about equity produces SC-CH4 estimates that differ by more than an order of magnitude between low- and high-income regions. Our central equity-weighted estimate for the USA increases to US$8,290 per tonne of CH4 whereas our estimate for sub-Saharan Africa decreases to US$134 per tonne of CH4.

Atmospheric rivers (ARs) account for a large portion of winter precipitation in the western US. To evaluate the sources of AR and precipitation predictability at subseasonal-to-seasonal timescales, we examine the relationships between two climate modes, the El Niño/Southern Oscillation (ENSO) and Madden-Julian oscillation (MJO), and winter hydroclimate in the western US. Our analysis uses a large ensemble of the Weather Research and Forecast (WRF) model simulations from 1981 to 2017, facilitating the assessment of uncertainty in climate mode-AR relationship due to climate variability and short length of observations. Over the North Pacific basin, we find ENSO-related latitudinal shifts of ARs, but ENSO has little effect on basin-wide averaged AR frequency. Over the western US, there is some uncertainty in the ENSO-AR connection due to the impacts of ENSO indices and datasets. However, extreme El Niño events defined by the ENSO longitude index are consistently linked to increased landfalling AR activity. The MJO can lead to significantly enhanced or suppressed landfalling AR activity in the large-ensembles, depending on the phase of MJO and time lag, while observations are too short in length to robustly show this signal. ENSO substantially modulates the MJO-AR relationship, triggering variegated responses of landfalling ARs and AR precipitation in La Niña and El Niño years. Our findings highlight the need to evaluate concurrent effects of different climate modes on ARs and precipitation, and may shed light on a path toward more accurate subseasonal-to-seasonal prediction of ARs and precipitation over the western US.

We investigate how mesoscale circulations associated with convective aggregation can modulate the sensitivity of the hydrologic cycle to warming. We quantify changes in the full distribution of rain across radiative-convective equilibrium states in a cloud-resolving model. For a given Sea Surface Temperature (SST), the shift in mean rainfall between disorganized and organized states is associated with a shift in atmospheric radiative cooling, and is roughly analogous to the effect of a 4K SST increase. With rising temperatures, the increase in mean rain rate is insensitive to the presence of organization, while extremes can intensify faster in the aggregated state, leading to a faster amplification in the sporadic nature of rain. When convection aggregates, heavy rain is enhanced by 20%-30% and nonlinear behaviors are observed as a function of SST and strength of aggregation feedbacks. First, radiative- and surface-flux aggregation feedbacks have multiplicative effects on extremes, illustrating a non-trivial sensitivity to the degree of organization. Second, alternating Clausius-Clapeyron and super-Clausius-Clapeyron regimes in extreme rainfall are found as a function of SST, corresponding to varying thermodynamic and dynamic contributions, and a large sensitivity to precipitation efficiency variations in some SST ranges. The potential for mesoscale circulations in amplifying the hydrologic cycle is established. However, these nonlinear distortions question the quantitative relevance of idealized self-aggregation. This calls for a deeper investigation of relationships which capture the coupling between global energetics, aggregation feedbacks and local convection, and for systematic tests of their sensitivity to domain configurations, surface boundary conditions, microphysics, and turbulence schemes.

Superparameterized (SP) global climate models have been shown to better simulate various features of precipitation relative to conventional models, including its diurnal cycle as well as its extremes. While various studies have focused on the effect of differing microphysics parameterizations on precipitation within limited-area cloud-resolving models, we examine here the effect on contiguous U.S. (CONUS) extremes in a global SP model. We vary the number of predicted moments for hydrometeor distributions, the character of the rimed ice species, and the representation of raindrop self-collection and breakup. Using a likelihood ratio test and accounting for the effects of multiple hypothesis testing, we find that there are some regional differences, particularly during spring and summer in the Southwest and the Midwest, in both the current climate and a warmer climate with uniformly increased sea surface temperatures. These differences are most statistically significant and widespread when the number of moments is changed. To determine whether these results are due to (fast) local effects of the different microphysics or the (slower) ensuing feedback on the large-scale atmospheric circulation, we run a series of short, 5-day simulations initialized from reanalysis data. We find that the differences largely disappear in these runs and therefore infer that the different parameterizations impact precipitation extremes indirectly via the large-scale circulation. Finally, we compare the present-day results with hourly rain gauge data and find that SP underestimates extremes relative to observations regardless of which microphysics scheme is used given a fixed model configuration and resolution.

Abstract. It has become increasingly common for researchers to utilize methods that identify weather features in climate models. There is an increasing recognition that the uncertainty associated with choice of detection method may affect our scientific understanding. For example, results from the Atmospheric River Tracking Method Intercomparison Project (ARTMIP) indicate that there are a broad range of plausible atmospheric river (AR) detectors, and that scientific results can depend on the algorithm used. There are similar examples from the literature on extratropical cyclones and tropical cyclones. It is therefore imperative to develop detection techniques that explicitly quantify the uncertainty associated with the detection of events. We seek to answer the question: given a ‘plausible’ AR detector, how does uncertainty in the detector quantitatively impact scientific results? We develop a large dataset of global AR counts, manually identified by a set of 8 researchers with expertise in atmospheric science, which we use to constrain parameters in a novel AR detection method. We use a Bayesian framework to sample from the set of AR detector parameters that yield AR counts similar to the expert database of AR counts; this yields a set of 'plausible' AR detectors from which we can assess quantitative uncertainty. This probabilistic AR detector has been implemented in the Toolkit for Extreme Climate Analysis (TECA), which allows for efficient processing of petabyte-scale datasets. We apply the TECA Bayesian AR Detector, TECA-BARD v1.0, to the MERRA2 reanalysis and show that the sign of the correlation between global AR count and El Nino Southern Oscillation depends on the set of parameters used.

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.

The gridding of daily accumulated precipitation—especially extremes—from ground-based station observations is problematic due to the fractal nature of precipitation, and therefore estimates of long period return values and their changes based on such gridded daily datasets are generally underestimated. In this paper, we characterize high-resolution changes in observed extreme precipitation from 1950 to 2017 for the contiguous United States (CONUS) based on in situ measurements only. Our analysis utilizes spatial statistical methods that allow us to derive gridded estimates that do not smooth extreme daily measurements and are consistent with statistics from the original station data while increasing the resulting signal-to-noise ratio. Furthermore, we use a robust statistical technique to identify significant pointwise changes in the climatology of extreme precipitation while carefully controlling the rate of false positives. We present and discuss seasonal changes in the statistics of extreme precipitation: the largest and most spatially coherent pointwise changes are in fall (SON), with approximately 33% of CONUS exhibiting significant changes (in an absolute sense). Other seasons display very few meaningful pointwise changes (in either a relative or absolute sense), illustrating the difficulty in detecting pointwise changes in extreme precipitation based on in situ measurements. While our main result involves seasonal changes, we also present and discuss annual changes in the statistics of extreme precipitation. In this paper we only seek to detect changes over time and leave attribution of the underlying causes of these changes for future work.

Gridded data products, for example interpolated daily measurements of precipitation from weather stations, are commonly used as a convenient substitute for direct observations because these products provide a spatially and temporally continuous and complete source of data. However, when the goal is to characterize climatological features of extreme precipitation over a spatial domain (e.g., a map of return values) at the native spatial scales of these phenomena, then gridded products may lead to incorrect conclusions because daily precipitation is a fractal field and hence any smoothing technique will dampen local extremes. To address this issue, we create a new “probabilistic” gridded product specifically designed to characterize the climatological properties of extreme precipitation by applying spatial statistical analysis to daily measurements of precipitation from the Global Historical Climatology Network over the contiguous United States. The essence of our method is to first estimate the climatology of extreme precipitation based on station data and then use a data-driven statistical approach to interpolate these estimates to a fine grid. We argue that our method yields an improved characterization of the climatology within a grid cell because the probabilistic behavior of extreme precipitation is much better behaved (i.e., smoother) than daily weather. Furthermore, the spatial smoothing innate to our approach significantly increases the signal-to-noise ratio in the estimated extreme statistics relative to an analysis without smoothing. Finally, by deriving a data-driven approach for translating extreme statistics to a spatially complete grid, the methodology outlined in this paper resolves the issue of how to properly compare station data with output from earth system models. We conclude the paper by comparing our probabilistic gridded product with a standard extreme value analysis of the Livneh gridded daily precipitation product. Our new data product is freely available on the Harvard Dataverse (https://bit.ly/2CXdnuj).

Earth system models are complex and represent a large number of processes, resulting in a persistent spread across climate projections for a given future scenario. Owing to different model performances against observations and the lack of independence among models, there is now evidence that giving equal weight to each available model projection is suboptimal. This Perspective discusses newly developed tools that facilitate a more rapid and comprehensive evaluation of model simulations with observations, process-based emergent constraints that are a promising way to focus evaluation on the observations most relevant to climate projections, and advanced methods for model weighting. These approaches are needed to distil the most credible information on regional climate changes, impacts, and risks for stakeholders and policy-makers.

Changes in land and ocean carbon storage in response to elevated atmospheric carbon dioxide concentrations and associated climate change, known as the concentration-carbon and climate-carbon feedbacks, are principal controls on the response of the climate system to anthropogenic greenhouse gas emissions. Such feedbacks have typically been quantified in the context of natural ecosystems, but land management activities are also responsive to future atmospheric carbon and climate changes. Here we show that inclusion of such human-driven responses within an Earth system model shifts both the terrestrial concentration-carbon and climate-carbon feedbacks toward increased carbon storage. We introduce a conceptual framework for decomposing these changes into separate concentration-land cover, climate-land cover, and land cover-carbon effects, providing a parsimonious means to diagnose sources of variation across numerical models capable of estimating such feedbacks.

Recently, it was recognized that widely used calculations of methane radiative forcing systematically underestimated its global value by 15% by omitting its shortwave effects. We show that shortwave forcing by methane can be accurately calculated despite considerable uncertainty and large gaps in its shortwave spectroscopy. We demonstrate that the forcing is insensitive, even when confronted with much more complete methane absorption spectra extending to violet light wavelengths derived from observations of methane-rich Jovian planets. We undertake the first spatially resolved global calculations of this forcing and find that it is dependent on bright surface features and clouds. Localized annual mean forcing from preindustrial to present-day methane increases approaches +0.25 W/m2, 10 times the global annualized shortwave forcing and 43% of the total direct CH4 forcing. Shortwave forcing by anthropogenic methane is sufficiently large and accurate to warrant its inclusion in historical analyses, projections, and mitigation strategies for climate change.

Dynamical downscaling is a widely used technique to properly capture regional surface heterogeneities that shape the local hydroclimatology. However, in the context of dynamical downscaling, the impacts on simulation fidelity have not been comprehensively evaluated across many user-specified factors, including the refinements of model horizontal resolution, large-scale forcing datasets, and dynamical cores. Two global-to-regional downscaling methods are used to assess these: specifically, the variable-resolution Community Earth System Model (VR-CESM) and the Weather Research and Forecasting (WRF) Model with horizontal resolutions of 28, 14, and 7 km. The modeling strategies are assessed by comparing the VR-CESM and WRF simulations with consistent physical parameterizations and grid domains. Two groups of WRF Models are driven by either the NCEP reanalysis dataset (WRF_NCEP) or VR-CESM7 results (WRF_VRCESM) to evaluate the effects of large-scale forcing datasets. The simulated hydroclimatologies are compared with reference datasets for key properties including total precipitation, snow cover, snow water equivalent (SWE), and surface temperature. The large-scale forcing datasets are critical to the WRF simulations of total precipitation but not surface temperature, controlled by the wind field and atmospheric moisture transport at the ocean boundary. No significant benefit is found in the regional average simulated hydroclimatology by increasing horizontal resolution refinement from 28 to 7 km, probably due to the systematic biases from the diagnostic treatment of rainfall and snowfall in the microphysics scheme. The choice of dynamical core has little impact on total precipitation but significantly determines simulated surface temperature, which is affected by the snow-albedo feedback in winter and soil moisture estimations in summer.

Implementing accurate calculations of radiative forcing can improve climate projections

Mountains are natural dams that impede atmospheric moisture transport and water towers that cool, condense, and store precipitation. They are essential in the western United States where precipitation is seasonal, and snowpack is needed to meet water demand. With anthropogenic climate change increasingly threatening mountain snowpack, there is a pressing need to better understand the driving climatological processes. However, the coarse resolution typical of modern global climate models renders them largely insufficient for this task, and signals a need for an advanced strategy. This paper continues the assessment of variable-resolution in the Community Earth System Model (VR-CESM) in modeling mountain hydroclimatology to understand the role of grid-spacing at 55, 28, 14, and 7 km and microphysics, specifically the Morrison and Gettelman (, MG1, https://doi.org/10.1175/2008JCLI2105.1) scheme versus the Gettelman and Morrison (, MG2, https://doi.org/10.1175/JCLI-D-14-00102.1) scheme. Eight VR-CESM simulations were performed from 1999 to 2015 with the F_AMIP_CAM5 component set, which couples the atmosphere-land models and prescribes ocean data. Refining horizontal grid-spacing from 28 to 7 km with the MG1 scheme did not improve the simulated mountain hydroclimatology. Substantial improvements occurred with the use of MG2 at grid-spacings ≤28 km compared to MG1 as shown with subsequent statistics. Average SWE bias diminished by 9.4X, 4.9X, and 3.5X from 55 to 7 km. The range in minimum (maximum) DJF spatial correlations increased by 0.1–0.2 in both precipitation and SWE. Mountain windward/leeward distributions and elevation profiles improved across hydroclimate variables, however not always with model resolution alone. Disconcertingly, all VR-CESM simulations exhibited a systemic mountain cold bias that worsened with elevation and will require further examination.

This paper presents two contributions for research into better understanding the role of anthropogenic warming in extreme weather. The first contribution is the generation of a large number of multi-decadal simulations using a medium-resolution atmospheric climate model, CAM5.1-1degree, under two scenarios of historical climate following the protocols of the C20C+ Detection and Attribution project: the one we have experienced (All-Hist), and one that might have been experienced in the absence of human interference with the climate system (Nat-Hist). These simulations are specifically designed for understanding extreme weather and atmospheric variability in the context of anthropogenic climate change. The second contribution takes advantage of the duration and size of these simulations in order to identify features of variability in the prescribed ocean conditions that may strongly influence calculated estimates of the role of anthropogenic emissions on extreme weather frequency (event attribution). There is a large amount of uncertainty in how much anthropogenic emissions should warm regional ocean surface temperatures, yet contributions to the C20C+ Detection and Attribution project and similar efforts so far use only one or a limited number of possible estimates of the ocean warming attributable to anthropogenic emissions when generating their Nat-Hist simulations. Thus, the importance of the uncertainty in regional attributable warming estimates to the results of event attribution studies is poorly understood. The identification of features of the anomalous ocean state that seem to strongly influence event attribution estimates should therefore be able to serve as a basis set for effective sampling of other plausible attributable warming patterns. The identification performed in this paper examines monthly temperature and precipitation output from the CAM5.1-1degree simulations averaged over 237 land regions, and compares interannual anomalous variations in the ratio between the frequencies of extremes in the All-Hist and Nat-Hist simulations against variations in ocean temperatures.

The Colorado floods of September 2013 caused severe damage and fatalities, and resulted from prolonged heavy rainfall unusual for that time of year – both in its record-breaking amounts and associated weather systems. We investigate the possible role of anthropogenic climate change in this extreme event. The unusual hydrometeorology of the event, however, challenges standard frameworks for attributing extreme events to anthropogenic climate change, because they typically struggle to simulate and connect the large-scale meteorology associated with local weather processes. Therefore we instead employ a part dynamical modelling- part observational- based event attribution approach, which simulates regional Colorado rainfall conditional on boundary conditions prescribed from the observed synoptic-scale meteorology in September 2013 – and assumes these conditions would have been similar in the absence of anthropogenic forcing. Using this ‘conditional event attribution’ approach we find that our regional climate model simulations indicate that anthropogenic drivers increased the magnitude of heavy northeast Colorado rainfall for the wet week in September 2013 by 30%, with the occurrence probability of a week at least that wet increasing by at least a factor of 1.3. By comparing the convective and large-scale components of rainfall, we find that this increase resulted in part from the additional moisture-carrying capacity of a warmer atmosphere – allowing more intense local convective rainfall that induced a dynamical positive feedback in the existing larger scale moisture flow – and also in part from additional moisture transport associated with larger scale circulation change. Our approach precludes assessment of changes in the frequency of the observed synoptic meteorological conditions themselves, and thus does not assess the effect of anthropogenic climate drivers on the statistics of heavy Colorado rainfall events. However, tailoring analysis tools to diagnose particular aspects of localized extreme weather events, conditional on the observed large-scale meteorology, can prove useful for diagnosing the physical effects of anthropogenic climate change on severe weather events – especially given large uncertainties in assessments of anthropogenic driven changes in atmospheric circulation.

Fossil fuel combustion and land-use change are the two largest contributors to industrial-era increases in atmospheric CO2 concentration. Projections of these are thus fundamental inputs for coupled Earth system models (ESMs) used to estimate the physical and biological consequences of future climate system forcing. While historical data sets are available to inform past and current climate analyses, assessments of future climate change have relied on projections of energy and land use from energy-economic models, constrained by assumptions about future policy, land-use patterns and socio-economic development trajectories. Here we show that the climatic impacts on land ecosystems drive significant feedbacks in energy, agriculture, land use and carbon cycle projections for the twenty-first century. We find that exposure of human-appropriated land ecosystem productivity to biospheric change results in reductions of land area used for crops; increases in managed forest area and carbon stocks; decreases in global crop prices; and reduction in fossil fuel emissions for a low-mid-range forcing scenario. The feedbacks between climate-induced biospheric change and human system forcings to the climate system - demonstrated here - are handled inconsistently, or excluded altogether, in the one-way asynchronous coupling of energy-economic models to ESMs used to date.

This study explores the use of nonuniform fast spherical Fourier transforms on meteorological data that are arbitrarily distributed on the sphere. The applicability of this methodology in the atmospheric sciences is demonstrated by estimating spectral coefficients for nontrivial subsets of reanalysis data on a uniformly spaced latitude-longitude grid, a global cloud resolving model on an icosahedral mesh with 3-km horizontal grid spacing, and for temperature anomalies from arbitrarily distributed weather stations over the United States. A spectral correction technique is developed that can be used in conjunction with the inverse transform to yield data interpolated onto a uniformly spaced grid, with optional triangular truncation, at reduced computational cost compared to other variance conserving interpolation methods, such as kriging or natural spline interpolation. The spectral correction yields information that can be used to deduce gridded observational biases not directly available from other methods.

Understanding potential impacts of climate change is complicated by spatially mismatched land representations between gridded datasets and models, and land use models with larger regions defined by geopolitical and/or biophysical criteria. Here we quantify the sensitivity of Global Change Assessment Model (GCAM) outputs to the delineation of Agro-Ecological Zones (AEZs), which are normally based on historical (1961–1990) climate. We reconstruct GCAM's land regions using projected (2071–2100) climate, and find large differences in estimated future land use that correspond with differences in agricultural commodity prices and production volumes. Importantly, historically delineated AEZs experience spatially heterogeneous climate impacts over time, and do not necessarily provide more homogenous initial land productivity than projected AEZs. We conclude that non-climatic criteria for land use region delineation are likely preferable for modeling land use change in the context of climate change, and that uncertainty associated with land delineation needs to be quantified.

An ensemble of six pairs of RCM experiments performed at 25 and 50 km for the period 1961–2000 over a large European domain is examined in order to evaluate the effects of resolution on the simulation of daily precipitation statistics. Application of the non-parametric two-sample Kolmorgorov–Smirnov test, which tests for differences in the location and shape of the probability distributions of two samples, shows that the distribution of daily precipitation differs between the pairs of simulations over most land areas in both summer and winter, with the strongest signal over southern Europe. Two-dimensional histograms reveal that precipitation intensity increases with resolution over almost the entire domain in both winter and summer. In addition, the 25 km simulations have more dry days than the 50 km simulations. The increase in dry days with resolution is indicative of an improvement in model performance at higher resolution, while the more intense precipitation exceeds observed values. The systematic increase in precipitation extremes with resolution across all models suggests that this response is fundamental to model formulation. Simple theoretical arguments suggest that fluid continuity, combined with the emergent scaling properties of the horizontal wind field, results in an increase in resolved vertical transport as grid spacing decreases. This increase in resolution-dependent vertical mass flux then drives an intensification of convergence and resolvable-scale precipitation as grid spacing decreases. This theoretical result could help explain the increasingly, and often anomalously, large stratiform contribution to total rainfall observed with increasing resolution in many regional and global models.

Numerous facets of scientific research implicitly or explicitly call for the estimation of probability densities. Histograms and kernel density estimates (KDEs) are two commonly used techniques for estimating such information, with the KDE generally providing a higher fidelity representation of the probability density function (PDF). Both methods require specification of either a bin width or a kernel bandwidth. While techniques exist for choosing the kernel bandwidth optimally and objectively, they are computationally intensive, since they require repeated calculation of the KDE. A solution for objectively and optimally choosing both the kernel shape and width has recently been developed by Bernacchia and Pigolotti (2011). While this solution theoretically applies to multidimensional KDEs, it has not been clear how to practically do so. A method for practically extending the Bernacchia-Pigolotti KDE to multidimensions is introduced. This multidimensional extension is combined with a recently-developed computational improvement to their method that makes it computationally efficient: a 2D KDE on 105 samples only takes 1 s on a modern workstation. This fast and objective KDE method, called the fastKDE method, retains the excellent statistical convergence properties that have been demonstrated for univariate samples. The fastKDE method exhibits statistical accuracy that is comparable to state-of-the-science KDE methods publicly available in R, and it produces kernel density estimates several orders of magnitude faster. The fastKDE method does an excellent job of encoding covariance information for bivariate samples. This property allows for direct calculation of conditional PDFs with fastKDE. It is demonstrated how this capability might be leveraged for detecting non-trivial relationships between quantities in physical systems, such as transitional behavior.

NASA Aqua-derived thermodynamic profiles, calculated spectral clear-sky outgoing longwave radiation (OLR), and vertical velocity fields from meteorological reanalyses are combined to determine the relative proportion of the far-infrared (FIR) and mid-infrared (MIR) spectral contributions to the total clear-sky OLR during different phases of El Niño–Southern Oscillation (ENSO). In the ascending branch of the tropical circulation, the spatial variance of upper tropospheric water vapor is shown to be larger during La Niña than El Niño and is consistent with zonal symmetry changes in the tropical waveguide and associated tropical-extratropical mixing. In the descending branch, upper tropospheric water vapor shows weaker coupling to lower layers that is evidenced by changes in the ratio of FIR to MIR in the clear-sky OLR. Diagnostics from the Geophysical Fluid Dynamics Laboratory AM3 model simulation are generally similar to satellite data, but the ratio of FIR to MIR is 5–10% larger with respect to dynamic regime.

Numerous studies have shown that atmospheric models with high horizontal resolution better represent the physics and statistics of precipitation in climate models. While it is abundantly clear from these studies that high-resolution increases the rate of extreme precipitation, it is not clear whether these added extreme events are “realistic”; whether they occur in simulations in response to the same forcings that drive similar events in reality. In order to understand whether increasing horizontal resolution results in improved model fidelity, a hindcast-based, multiresolution experimental design has been conceived and implemented: the InitiaLIzed-ensemble, Analyze, and Develop (ILIAD) framework. The ILIAD framework allows direct comparison between observed and simulated weather events across multiple resolutions and assessment of the degree to which increased resolution improves the fidelity of extremes. Analysis of 5 years of daily 5 day hindcasts with the Community Earth System Model at horizontal resolutions of 220, 110, and 28 km shows that: (1) these hindcasts reproduce the resolution-dependent increase of extreme precipitation that has been identified in longer-duration simulations, (2) the correspondence between simulated and observed extreme precipitation improves as resolution increases; and (3) this increase in extremes and precipitation fidelity comes entirely from resolved-scale precipitation. Evidence is presented that this resolution-dependent increase in precipitation intensity can be explained by the theory of Rauscher et al. (), which states that precipitation intensifies at high resolution due to an interaction between the emergent scaling (spectral) properties of the wind field and the constraint of fluid continuity.

Runoff and soil moisture are two key components of the global hydrologic cycle that should be validated at local to global scales in Earth System Models (ESMs) used for climate projection. We have evaluated the runoff and surface soil moisture output by the Community Climate System Model (CCSM) along with 8 other models from the Coupled Model Intercomparison Project (CMIP5) repository using satellite soil moisture observations and stream gauge corrected runoff products. A series of Community Land Model (CLM) runs forced by reanalysis and coupled model outputs was also performed to identify atmospheric drivers of biases and uncertainties in the CCSM. Results indicate that surface soil moisture simulations tend to be positively biased in high latitude areas by most selected CMIP5 models except CCSM, FGOALS, and BCC, which share similar land surface model code. With the exception of GISS, runoff simulations by all selected CMIP5 models were overestimated in mountain ranges and in most of the Arctic region. In general, positive biases in CCSM soil moisture and runoff due to precipitation input error were offset by negative biases induced by temperature input error. Excluding the impact from atmosphere modeling, the global mean of seasonal surface moisture oscillation was out of phase compared to observations in many years during 1985–2004. The CLM also underestimated runoff in the Amazon, central Africa, and south Asia, where soils all have high clay content. We hypothesize that lack of a macropore flow mechanism is partially responsible for this underestimation. However, runoff was overestimated in the areas covered by volcanic ash soils (i.e., Andisols), which might be associated with poor soil porosity representation in CLM. Our results indicate that CCSM predictability of hydrology could be improved by addressing the compensating errors associated with precipitation and temperature and updating the CLM soil representation.

The radiative forcing (RF) of carbon dioxide (CO2) is the leading contribution to climate change from anthropogenic activities. Calculating CO2 RF requires detailed knowledge of spectral line parameters for thousands of infrared absorption lines. A reliable spectroscopic characterization of CO2 forcing is critical to scientific and policy assessments of present climate and climate change. Our results show that CO2 RF in a variety of atmospheres is remarkably insensitive to known uncertainties in the three main CO2 spectroscopic parameters: the line shapes, line strengths, and half widths. We specifically examine uncertainty in RF due to line mixing as this process is critical in determining line shapes in the far wings of CO2 absorption lines. RF computed with a Voigt line shape is also examined. Overall, the spectroscopic uncertainty in present-day CO2 RF is less than 1%, indicating a robust foundation in our understanding of how rising CO2 warms the climate system.

Greenhouse gas emissions from fossil fuel combustion and land use are changing the radiative budget of the Earth and changing its climate. The negative impacts of this climate change on natural and human systems are already emergent. The solution is to eliminate greenhouse gas emissions altogether as soon as possible, but the rate at which these emissions can decrease is limited by human reliance on fossil fuels for energy and the infrastructural, socio-economic, and behavioral inertia of current systems around the world. In this chapter, we discuss the physical impacts as well as the many challenges and obstacles to ‘bending the curve’, and provide a framework of possible solutions.

The superparameterized Community Atmosphere Model (SPCAM) is used to investigate the impact and geographic sensitivity of positive Indian Ocean Dipole (+IOD) sea-surface temperatures (SSTs) on Madden-Julian oscillation (MJO) propagation. The goal is to clarify potentially appreciable +IOD effects on MJO dynamics detected in prior studies by using a global model with explicit convection representation. Prescribed climatological October SSTs and variants of the SST distribution from October 2006, a +IOD event, force the model. Modest MJO convection weakening over the Maritime Continent occurs when either climatological SSTs, or +IOD SST anomalies restricted to the Indian Ocean, are applied. However, severe MJO weakening occurs when either +IOD SST anomalies are applied globally or restricted to the equatorial Pacific. MJO disruption is associated with time-mean changes in the zonal wind profile and lower moist static energy (MSE) in subsiding air masses imported from the Subtropics by Rossby-like gyres. On intraseasonal scales, MJO disruption arises from significantly smaller MSE accumulation, weaker meridional advective moistening, and overactive submonthly eddies that mix drier subtropical air into the path of MJO convection. These results (1) demonstrate that SPCAM reproduces observed time-mean and intraseasonal changes during +IOD episodes, (2) reaffirm the role that submonthly eddies play in MJO propagation and show that such multiscale interactions are sensitive to interannual SST states, and (3) suggest that boreal fall +IOD SSTs local to the Indian Ocean have a significantly smaller impact on Maritime Continent MJO propagation compared to contemporaneous Pacific SST anomalies which, for October 2006, resemble El Ninõ-like conditions.

Projected increases in mean precipitation are constrained by the atmospheric energy budget through radiative-convective equilibrium. However, significant differences persist between climate models on the rate of increase in precipitation per unit warming, mostly arising from the clear-sky radiative response. While the intermodel spread in clear-sky longwave cooling has been explained by climate feedbacks, the sources of spread in clear-sky shortwave heating are still unclear. This article focuses on the latter. Since water vapor contributes most of the atmospheric shortwave absorption, both intermodel differences in its spatial distribution and in radiative transfer parameterizations are plausible hypotheses for the spread. This work reestablishes the primary contribution from water vapor relative to other shortwave-absorbing species and evaluates the validity of both hypotheses. It is found that the intermodel spread in shortwave absorption change most likely originates from the radiation schemes, possibly because of simplifications induced by their low spectral resolutions. Key Points Water vapor appears as the main cause of intermodel differences in atmospheric shortwave absorption Such discrepancies seem to come from the shortwave radiative transfer formulation for water vapor Errors arising from the coarse spectral resolution in radiation schemes could explain this spread

This paper is among the first to report on the full integration of basin-scale models that include projections of the demand and supply of water and energy for residential, commercial, industrial, and agricultural sector users. We link two widely used regional planning models that allow one to study the impact of rising climate variability on water and electricity use in Sacramento, California. Historic data combined with the current energy and water system configuration was used to assess the implications of changes in temperature and precipitation. Climate simulations suggest that electricity imports to the region would increase during hot dry spells, when regional power production is most constrained. In particular, regional imports of electricity would increase over 35 % in hot dry years, assuming a 4 °C increase in average temperature and a 25 % decrease in average precipitation.

We demonstrate the effectiveness of a new method for quantifying radiative forcing from land use and land cover change (LULCC) within an integrated assessment model, the Global Change Assessment Model (GCAM). The method relies on geographically differentiated estimates of radiative forcing from albedo change associated with major land cover transitions derived from the Community Earth System Model. We find that conversion of 1 km2 of woody vegetation (forest and shrublands) to non-woody vegetation (crops and grassland) yields between 0 and −0.71 nW/m2 of globally averaged radiative forcing determined by the vegetation characteristics, snow dynamics, and atmospheric radiation environment characteristic within each of 151 regions we consider globally. Across a set of scenarios designed to span a range of potential future LULCC, we find LULCC forcing ranging from −0.06 to −0.29 W/m2 by 2070 depending on assumptions regarding future crop yield growth and whether climate policy favors afforestation or bioenergy crops. Inclusion of this previously uncounted forcing in the policy targets driving future climate mitigation efforts leads to changes in fossil fuel emissions on the order of 1.5 PgC/yr by 2070 for a climate forcing limit of 4.5 Wm−2, corresponding to a 12–67 % change in fossil fuel emissions depending on the scenario. Scenarios with significant afforestation must compensate for albedo-induced warming through additional emissions reductions, and scenarios with significant deforestation need not mitigate as aggressively due to albedo-induced cooling. In all scenarios considered, inclusion of albedo forcing in policy targets increases forest and shrub cover globally.