Search

Realistic simulation of the Earth's mean-state climate remains a major challenge, and yet it is crucial for predicting the climate system in transition. Deficiencies in models' process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections and incorrect attribution of the physical mechanisms governing past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and long-standing biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when it is run as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback, as well as the fidelity of the simulated base climate state, are important for constraining the climate in the past and future.

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

Marine low clouds strongly cool the planet. How this cooling effect will respond to climate change is a leading source of uncertainty in climate sensitivity, the planetary warming resulting from CO2 doubling. Here, we observationally constrain this low cloud feedback at a near-global scale. Satellite observations are used to estimate the sensitivity of low clouds to interannual meteorological perturbations. Combined with model predictions of meteorological changes under greenhouse warming, this permits quantification of spatially resolved cloud feedbacks. We predict positive feedbacks from midlatitude low clouds and eastern ocean stratocumulus, nearly unchanged trade cumulus and a near-global marine low cloud feedback of 0.19 ± 0.12 W m−2 K−1 (90% confidence). These constraints imply a moderate climate sensitivity (~3 K). Despite improved midlatitude cloud feedback simulation by several current-generation climate models, their erroneously positive trade cumulus feedbacks produce unrealistically high climate sensitivities. Conversely, models simulating erroneously weak low cloud feedbacks produce unrealistically low climate sensitivities. Marine low clouds cool the planet, but their response to warming is uncertain and dominates the spread in model-based climate sensitivities. Observational constraints suggest smaller cloud feedbacks than previously reported and imply a more moderate climate sensitivity.

The single-column model (SCM) functionality of the Energy Exascale Earth System Model version 1 (E3SMv1) is described in this paper. The E3SM SCM was adopted from the SCM used in the Community Atmosphere Model (CAM) but has evolved significantly since then. We describe changes made to the aerosol specification in the SCM, idealizations, and developments made so that the SCM uses the same dynamical core as the full general circulation model (GCM) component. Based on these changes, we describe and demonstrate the seamless capability to “replay” a GCM column using the SCM. We give an overview of the E3SM case library and briefly describe which cases may serve as useful proxies for replicating and investigate some long-standing biases in the full GCM runs while demonstrating that the E3SM SCM is an efficient tool for both model development and evaluation.

A method is proposed for combining information from several emergent constraints into a probabilistic estimate for a climate sensitivity proxy Y such as equilibrium climate sensitivity (ECS). The method is based on fitting a multivariate Gaussian PDF for Y and the emergent constraints using an ensemble of global climate models (GCMs); it can be viewed as a form of multiple linear regression of Y on the constraints. The method accounts for uncertainties in sampling this multidimensional PDF with a small number of models, for observational uncertainties in the constraints, and for overconfidence about the correlation of the constraints with the climate sensitivity. Its general form (Method C) accounts for correlations between the constraints. Method C becomes less robust when some constraints are too strongly related to each other; this can be mitigated using regularization approaches such as ridge regression. An illuminating special case, Method U, neglects any correlations between constraints except through their mutual relationship to the climate proxy; it is more robust to small GCM sample size and is appealingly interpretable. These methods are applied to ECS and the climate feedback parameter using a previously published set of 11 possible emergent constraints derived from climate models in the Coupled Model Intercomparison Project (CMIP). The ±2σ posterior range of ECS for Method C with no overconfidence adjustment is 4.3 ± 0.7 K. For Method U with a large overconfidence adjustment, it is 4.0 ± 1.3 K. This study adds confidence to past findings that most constraints predict higher climate sensitivity than the CMIP mean.

This paper presents a process-oriented evaluation of precipitating stratocumulus and its transition to cumulus in version 1 of the Energy Exascale Earth System Model (E3SMv1) using comprehensive case-study observations from a field campaign of the Atmospheric Radiation Measurement program (ARM). The E3SMv1 single-column model (SCM) of the marine boundary layer and its low clouds and precipitation are compared to observations including subcloud drizzle retrievals from a combination of Doppler radar and lidar backscatter measurements. The SCM is also compared to a large-eddy simulation (LES) of the same case. The combination of advanced remote sensing observations and LES is a powerful framework to evaluate the physical parameterizations of large-scale models. Given the observed large-scale environment, the E3SMv1 SCM realistically represents the evolution of clouds and boundary layer structure during the stratocumulus-to-cumulus transition. The model well simulates the liquid water path and its diurnal cycle in the stratocumulus period as well as the two-layer vertical thermodynamic structure and lower cloud fraction in the transition period. E3SMv1’s success in simulating the cloud in the stratocumulus period permitted examination of its precipitation processes. Here problems were identified with E3SMv1 producing an unrealistically small subcloud precipitation fraction, an unrealistic double peak in the vertical profiles of precipitation mass, and drizzle that evaporates too close to the surface. Further model diagnostics determined that these unrealistic characteristics resulted from an overly long microphysics time step and an unrealistic parameterization of the precipitation fraction. These results imply that careful consideration of these issues is needed in order to better simulate precipitation processes in marine stratocumulus.

This paper describes the first implementation of the Δx = 3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40‐day prescribed‐sea‐surface‐temperature simulation (January 20 through February 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non‐hydrostatic Domains (DYAMOND) Phase 2 model intercomparison. Effective resolution is found to be ∼6× the horizontal dynamics grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long‐standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g., tropical cyclones, extratropical cyclones including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.

Realistic simulation of the Earth’s mean state climate remains a major challenge and yet it is crucial for predicting the climate system in transition. Deficiencies in models’ process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections as well as incorrect attribution of the physical mechanisms governing the past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and longstanding biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when running as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback as well as the fidelity of the simulated base climate state are important for constraining the climate in the past and future.

In this study, we find significant sensitivity to the choice of time step for the Energy Exascale Earth System Model's atmospheric component, leading to large decreases in the magnitude of cloud forcing when the time step is reduced to 10 s. Reducing the time step size for the microphysics increases precipitation, leading to a drying of the atmosphere and an increase in surface evaporation. This effect is amplified when the microphysics is substepped together with other cloud physics processes. Coupling the model's dynamics and physics more frequently reduces cloud fraction at lower altitudes, while producing more cloud liquid at higher altitudes. Reducing the deep convection time step also reduces low cloud mass and cloud fraction. Together, these results suggest that cloud physics in a global circulation model can depend strongly on time step and, in particular, on the frequency with which cloud‐related processes are coupled with each other and with the model dynamics.

Climate simulations with more accurate process-level representation at finer resolutions ( km) are a pressing need in order to provide more detailed actionable information to policy makers regarding extreme events in a changing climate. Computational limitation is a major obstacle for building and running high-resolution (HR, here 0.25∘ average grid spacing at the Equator) models (HRMs). A more affordable path to HRMs is to use a global regionally refined model (RRM), which only simulates a portion of the globe at HR while the remaining is at low resolution (LR, 1∘). In this study, we compare the Energy Exascale Earth System Model (E3SM) atmosphere model version 1 (EAMv1) RRM with the HR mesh over the contiguous United States (CONUS) to its corresponding globally uniform LR and HR configurations as well as to observations and reanalysis data. The RRM has a significantly reduced computational cost (roughly proportional to the HR mesh size) relative to the globally uniform HRM. Over the CONUS, we evaluate the simulation of important dynamical and physical quantities as well as various precipitation measures. Differences between the RRM and HRM over the HR region are predominantly small, demonstrating that the RRM reproduces the precipitation metrics of the HRM over the CONUS. Further analysis based on RRM simulations with the LR vs. HR model parameters reveals that RRM performance is greatly influenced by the different parameter choices used in the LR and HR EAMv1. This is a result of the poor scale-aware behavior of physical parameterizations, especially for variables influencing sub-grid-scale physical processes. RRMs can serve as a useful framework to test physics schemes across a range of scales, leading to improved consistency in future E3SM versions. Applying nudging-to-observations techniques within the RRM framework also demonstrates significant advantages over a free-running configuration for use as a test bed and as such represents an efficient and more robust physics test bed capability. Our results provide additional confirmatory evidence that the RRM is an efficient and effective test bed for HRM development.

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean and sea ice (60km in the midlatitudes and 30km at the equator and poles), and river transport (55km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model's strong aerosol-related effective radiative forcing (ERFari+aci=−1.65W/m2) and high equilibrium climate sensitivity (ECS=5.3K).

The Single Column Atmosphere Model (SCAM) is a single column model version of the Community Atmosphere Model (CAM). Here we describe the functionality and features of SCAM6, available as part of CAM6 in the Community Earth System Model, version 2 (CESM2). SCAM6 features a wide selection of standard cases, as well as the ability to easily configure a case specified by the user based on a particular point in a CAM 3‐D simulation. This work illustrates how SCAM6 reproduces CAM6 results for physical parameterizations, mostly of moisture and clouds. We demonstrate how SCAM6 can be used for model development through different physics selections, as well as with parameter sweep experiments to highlight the sensitivity of cloud properties to the specification of the vapor deposition process in the cloud microphysics. Furthermore, we use SCAM6 to illustrate the sensitivity of CAM6 cloud radiative properties and precipitation to variable drop number (cloud microphysics properties). Finally, we illustrate how SCAM6 can be used to explore critical emergent processes such as cloud feedbacks and show that CAM6 cloud responses to surface warming in stratus and stratocumulus regimes are similar to those in CAM5. CAM6 has a larger response in the shallow cumulus regime than CAM5. CAM6 cloud feedbacks in the shallow cumulus regime are sensitive to turbulence parameters. SCAM6 is thus a valuable tool for model development, evaluation, and scientific analy sis and an important part of the model hierarchy in Community Earth System Model, version 2.

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. ISSN:1758-6798

The spread in global mean precipitation among climate models is explored in two ensembles using the complementary perspectives of surface evaporation and energy budgets. Models with higher global-m...

In this study, we analyze the realism with which tropical cyclones (TCs) are simulated in the fully coupled low‐ and high‐resolution Energy Exascale Earth System Model (E3SM) version 1, with a focus on the latter. Compared to the low‐resolution (grid spacing of ∼1°), the representation of TCs improves considerably in the high‐resolution configuration (grid spacing of ∼0.25°). Significant improvements are found in the global TC frequency, TC lifetime maximum intensities, and the relative distribution of TCs among the different basins. However, at both resolutions, spurious TC activity is found in some basins, notably in the subtropical regions. Contrasting the simulated large‐scale TC environment with observations reveals that the model environment is unrealistically conducive for TC development in those regions. Further analysis indicates that these biases are likely related to those in thermodynamic potential intensity, caused by systematic SST biases, and vertical wind shear in the coupled model. TC‐ocean interaction is also examined in the high‐resolution configuration of the model. The salient features of the ocean's response to TC‐induced mixing and the ocean's impact on TC intensification are well‐reproduced. Finally, an evaluation of the influence of El Niño Southern Oscillation (ENSO) on TCs in the high‐resolution configuration of the model reveals that the ENSO‐TC relationship in the model has the right sign and is significant for the North Atlantic and Northwest Pacific, albeit weaker than in observations. In summary, the high‐resolution configuration of the E3SM model simulates TC activity reasonably and hence could be a useful tool for TC‐related research.

The single column model (SCM) functionality of the Energy Exascale Earth System Model version 1 (E3SMv1) is described in this paper. The E3SM SCM was adopted from the SCM used in the Community Atmosphere Model (CAM), but has evolved significantly since then. We describe changes made to the aerosol specification in the SCM, idealizations, and developments made so that the SCM uses the same dynamical core as the full general circulation model (GCM) component. Based on these changes, we describe and demonstrate the seamless capability to ``replay" a GCM column using the SCM. We give an overview of the E3SM case library and briefly describe which cases may serve as useful proxies for replicating and investigate some long standing biases in the full GCM runs, while demonstrating that the E3SM SCM is an efficient tool for both model development and evaluation.

Equilibrium climate sensitivity, the global surface temperature response to CO urn:x-wiley:grl:media:grl60047:grl60047-math-0001 doubling, has been persistently uncertain. Recent consensus places it likely within 1.5–4.5 K. Global climate models (GCMs), which attempt to represent all relevant physical processes, provide the most direct means of estimating climate sensitivity via CO urn:x-wiley:grl:media:grl60047:grl60047-math-0002 quadrupling experiments. Here we show that the closely related effective climate sensitivity has increased substantially in Coupled Model Intercomparison Project phase 6 (CMIP6), with values spanning 1.8–5.6 K across 27 GCMs and exceeding 4.5 K in 10 of them. This (statistically insignificant) increase is primarily due to stronger positive cloud feedbacks from decreasing extratropical low cloud coverage and albedo. Both of these are tied to the physical representation of clouds which in CMIP6 models lead to weaker responses of extratropical low cloud cover and water content to unforced variations in surface temperature. Establishing the plausibility of these higher sensitivity models is imperative given their implied societal ramifications.

This study provides comprehensive insight into the notable differences in clouds and precipitation simulated by the Energy Exascale Earth System Model Atmosphere Model version 0 and version 1 (EAMv1). Several sensitivity experiments are conducted to isolate the impact of changes in model physics, resolution, and parameter choices on these differences. The overall improvement in EAMv1 clouds and precipitation is primarily attributed to the introduction of a simplified third‐order turbulence parameterization Cloud Layers Unified By Binormals (along with the companion changes) for a unified treatment of boundary layer turbulence, shallow convection, and cloud macrophysics, though it also leads to a reduction in subtropical coastal stratocumulus clouds. This lack of stratocumulus clouds is considerably improved by increasing vertical resolution from 30 to 72 layers, but the gain is unfortunately subsequently offset by other retuning to reach the top‐of‐atmosphere energy balance. Increasing vertical resolution also results in a considerable underestimation of high clouds over the tropical warm pool, primarily due to the selection for numerical stability of a higher air parcel launch level in the deep convection scheme. Increasing horizontal resolution from 1° to 0.25° without retuning leads to considerable degradation in cloud and precipitation fields, with much weaker tropical and subtropical short‐ and longwave cloud radiative forcing and much stronger precipitation in the intertropical convergence zone, indicating poor scale awareness of the cloud parameterizations. To avoid this degradation, significantly different parameter settings for the low‐resolution (1°) and high‐resolution (0.25°) were required to achieve optimal performance in EAMv1.

The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1∘ horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model becomes negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors in early V1 versions decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in V1 results in a very thin model layer at the Earth's surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for V1.

Emergent constraints are quantities that are observable from current measurements and have skill predicting future climate. This study explores 19 previously proposed emergent constraints related to equilibrium climate sensitivity (ECS; the global-average equilibrium surface temperature response to CO2doubling). Several constraints are shown to be closely related, emphasizing the importance for careful understanding of proposed constraints. A new method is presented for decomposing correlation between an emergent constraint and ECS into terms related to physical processes and geographical regions. Using this decomposition, one can determine whether the processes and regions explaining correlation with ECS correspond to the physical explanation offered for the constraint. Shortwave cloud feedback is generally found to be the dominant contributor to correlations with ECS because it is the largest source of intermodel spread in ECS. In all cases, correlation results from interaction between a variety of terms, reflecting the complex nature of ECS and the fact that feedback terms and forcing are themselves correlated with each other. For 4 of the 19 constraints, the originally proposed explanation for correlation is borne out by our analysis. These four constraints all predict relatively high climate sensitivity. The credibility of six other constraints is called into question owing to correlation with ECS coming mainly from unexpected sources and/or lack of robustness to changes in ensembles. Another six constraints lack a testable explanation and hence cannot be confirmed. The fact that this study casts doubt upon more constraints than it confirms highlights the need for caution when identifying emergent constraints from small ensembles.

Because weather and climate models must capture a wide variety of spatial and temporal scales, they rely heavily on parameterizations of subgrid‐scale processes. The goal of this study is to demonstrate that the assumptions used to couple these parameterizations have an important effect on the climate of version 0 of the Energy Exascale Earth System Model (E3SM) General Circulation Model (GCM), a close relative of version 1 of the Community Earth System Model (CESM1). Like most GCMs, parameterizations in E3SM are sequentially split in the sense that parameterizations are called one after another with each subsequent process feeling the effect of the preceding processes. This coupling strategy is noncommutative in the sense that the order in which processes are called impacts the solution. By examining a suite of 24 simulations with deep convection, shallow convection, macrophysics/microphysics, and radiation parameterizations reordered, process order is shown to have a big impact on predicted climate. In particular, reordering of processes induces differences in net climate feedback that are as big as the intermodel spread in phase 5 of the Coupled Model Intercomparison Project. One reason why process ordering has such a large impact is that the effect of each process is influenced by the processes preceding it. Where output is written is therefore an important control on apparent model behavior. Application of k‐means clustering demonstrates that the positioning of macro/microphysics and shallow convection plays a critical role on the model solution.

The Community Atmosphere Model (CAM) adopts Cloud Layers Unified By Binormals scheme (CLUBB) and an updated microphysics (MG2) scheme for a more unified treatment of cloud processes. This makes interactions between parameterizations tighter and more explicit. In this study, a cloudy planetary boundary layer (PBL) oscillation related to interaction between CLUBB and MG2 is identified in CAM. This highlights the need for consistency between the coupled sub-grid processes in climate model development. This oscillation occurs most often in the marine cumulus cloud regime. The oscillation occurs only if the modeled PBL is strongly decoupled and precipitation evaporates below the cloud. Two aspects of the parameterized coupling assumptions between CLUBB and MG2 schemes cause the oscillation: 1) a parameterized relationship between rain evaporation and CLUBB's sub-grid spatial variance of moisture and heat that induces an extra cooling in the lower PBL; and 2) rain evaporation which happens at a too low an altitude because of the precipitation fraction parameterization in MG2. Either one of these two conditions can overly stabilize the PBL and reduce the upward moisture transport to the cloud layer so that the PBL collapses. Global simulations prove that turning off the evaporation-variance coupling and improving the precipitation fraction parameterization effectively reduces the cloudy PBL oscillation in marine cumulus clouds. By evaluating the causes of the oscillation in CAM, we have identified the PBL processes that should be examined in models having similar oscillations. This study may draw the attention of the modeling and observational communities to the issue of coupling between parameterized physical processes.

We examine the global water cycle characteristics in the Accelerated Climate Modeling for Energy v0.3 model (a close relative to version 5.3 of the Community Atmosphere Model) in atmosphere-only simulations spanning the years 1980–2005. We evaluate the simulations using a broad range of observational and reanalysis datasets, examine how the simulations change when the horizontal resolution is increased from 1° to 0.25 $$^{\circ }$$ , and compare the simulations against models participating in the the Atmosphere Model Intercomparison Project of the 5th Coupled Model Intercomparison Project (CMIP5). Particular effort has been made to evaluate the model using the best available observational estimates and verifying model biases with additional datasets when differences are known to exist among the observations. Regardless of resolution, the model exhibits several biases: global-mean precipitation, evaporation, and precipitable water are too high, light precipitation occurs too frequently, and the atmospheric residence time of water is too short. Many of these biases are shared by the multi-model mean climate of models participating in CMIP5. The reasons behind regional biases in precipitation are discussed by examining how different fields, such as local evaporation and transport of water vapor, contribute to the bias. Although increasing the horizontal resolution does not drastically change the water cycle, it does lead to a few differences: an increase in global mean precipitation rate, an increase in the fraction of total precipitation that falls over land, more frequent heavy precipitation (>30 mm day $$^{-1}$$ ), and a decrease in precipitable water. One of the most notable changes is the shift of precipitation produced by the convective parameterization to that produced by the large-scale microphysics parameterization. We analyze how changes in moisture and circulation with resolution contribute to this shift in the precipitation partitioning. Because changing horizontal resolution requires some re-tuning, the effect of that tuning was evaluated by performing an additional simulation at 1 $$^{\circ }$$ but using the tunings from the 0.25 $$^{\circ }$$ simulation. The evaluation shows that the more frequent heavy precipitation, the decrease in precipitable water, and the shift from convective to large-scale precipitation are predominantly due to resolution changes, while tuning changes have a major influence on the global mean precipitation and the land/ocean partitioning of precipitation.

Author(s): Rasch, PJ; Xie, S; Ma, PL; Lin, W; Wang, H; Tang, Q; Burrows, SM; Caldwell, P; Zhang, K; Easter, RC; Cameron-Smith, P; Singh, B; Wan, H; Golaz, JC; Harrop, BE; Roesler, E; Bacmeister, J; Larson, VE; Evans, KJ; Qian, Y; Taylor, M; Leung, LR; Zhang, Y; Brent, L; Branstetter, M; Hannay, C; Mahajan, S; Mametjanov, A; Neale, R; Richter, JH; Yoon, JH; Zender, CS; Bader, D; Flanner, M; Foucar, JG; Jacob, R; Keen, N; Klein, SA; Liu, X; Salinger, AG; Shrivastava, M; Yang, Y | Abstract: The Energy Exascale Earth System Model Atmosphere Model version 1, the atmospheric component of the Department of Energy's Energy Exascale Earth System Model is described. The model began as a fork of the well-known Community Atmosphere Model, but it has evolved in new ways, and coding, performance, resolution, physical processes (primarily cloud and aerosols formulations), testing and development procedures now differ significantly. Vertical resolution was increased (from 30 to 72 layers), and the model top extended to 60 km (~0.1 hPa). A simple ozone photochemistry predicts stratospheric ozone, and the model now supports increased and more realistic variability in the upper troposphere and stratosphere. An optional improved treatment of light-absorbing particle deposition to snowpack and ice is available, and stronger connections with Earth system biogeochemistry can be used for some science problems. Satellite and ground-based cloud and aerosol simulators were implemented to facilitate evaluation of clouds, aerosols, and aerosol-cloud interactions. Higher horizontal and vertical resolution, increased complexity, and more predicted and transported variables have increased the model computational cost and changed the simulations considerably. These changes required development of alternate strategies for tuning and evaluation as it was not feasible to “brute force” tune the high-resolution configurations, so short-term hindcasts, perturbed parameter ensemble simulations, and regionally refined simulations provided guidance on tuning and parameterization sensitivity to higher resolution. A brief overview of the model and model climate is provided. Model fidelity has generally improved compared to its predecessors and the CMIP5 generation of climate models.

A large, long-standing, and pervasive climate model bias is excessive absorbed shortwave radiation (ASR) over the midlatitude oceans, especially the Southern Ocean. This study investigates both the underlying mechanisms for and climate impacts of this bias within the Community Earth System Model, version 1, with the Community Atmosphere Model, version 5 [CESM1(CAM5)]. Excessive Southern Ocean ASR in CESM1(CAM5) results in part because low-level clouds contain insufficient amounts of supercooled liquid. In a present-day atmosphere-only run, an observationally motivated modification to the shallow convection detrainment increases supercooled cloud liquid, brightens low-level clouds, and substantially reduces the Southern Ocean ASR bias. Tuning to maintain global energy balance enables reduction of a compensating tropical ASR bias. In the resulting preindustrial fully coupled run with a brighter Southern Ocean and dimmer tropics, the Southern Ocean cools and the tropics warm. As a result of the enhanced meridional temperature gradient, poleward heat transport increases in both hemispheres (especially the Southern Hemisphere), and the Southern Hemisphere atmospheric jet strengthens. Because northward cross-equatorial heat transport reductions occur primarily in the ocean (80%), not the atmosphere (20%), a proposed atmospheric teleconnection linking Southern Ocean ASR bias reduction and cooling with northward shifts in tropical precipitation has little impact. In summary, observationally motivated supercooled liquid water increases in shallow convective clouds enable large reductions in long-standing climate model shortwave radiation biases. Of relevance to both model bias reduction and climate dynamics, quantifying the influence of Southern Ocean cooling on tropical precipitation requires a model with dynamic ocean heat transport.

This study clarifies the causes of intermodel differences in the global-average temperature response to doubled CO2, commonly known as equilibrium climate sensitivity (ECS). The authors begin by noting several issues with the standard approach for decomposing ECS into a sum of forcing and feedback terms. This leads to a derivation of an alternative method based on linearizing the effect of the net feedback. Consistent with previous studies, the new method identifies shortwave cloud feedback as the dominant source of intermodel spread in ECS. This new approach also reveals that covariances between cloud feedback and forcing, between lapse rate and longwave cloud feedbacks, and between albedo and shortwave cloud feedbacks play an important and previously underappreciated role in determining model differences in ECS. Defining feedbacks based on fixed relative rather than specific humidity (as suggested by Held and Shell) reduces the covariances between processes and leads to more straightforward interpretations of results.

Single-column model (SCM) capability is an important tool for general circulation model development. In this study, the SCM mode of version 5 of the Community Atmosphere Model (CAM5) is shown to handle aerosol initialization and advection improperly, resulting in aerosol, cloud-droplet, and ice crystal concentrations which are typically much lower than observed or simulated by CAM5 in global mode. This deficiency has a major impact on stratiform cloud simulations but has little impact on convective case studies because aerosol is currently not used by CAM5 convective schemes and convective cases are typically longer in duration (so initialization is less important). By imposing fixed aerosol or cloud-droplet and crystal number concentrations, the aerosol issues described above can be avoided. Sensitivity studies using these idealizations suggest that the Meyers et al. (1992) ice nucleation scheme prevents mixed-phase cloud from existing by producing too many ice crystals. Microphysics is shown to strongly deplete cloud water in stratiform cases, indicating problems with sequential splitting in CAM5 and the need for careful interpretation of output from sequentially split climate models. Droplet concentration in the general circulation model (GCM) version of CAM5 is also shown to be far too low (~ 25 cm−3) at the southern Great Plains (SGP) Atmospheric Radiation Measurement (ARM) site.

Numerical weather, climate, or Earth system models involve the coupling of components. At a broad level, these components can be classified as the resolved fluid dynamics, unresolved fluid dynamical aspects (i.e., those represented by physical parameterizations such as subgrid-scale mixing), and nonfluid dynamical aspects such as radiation and microphysical processes. Typically, each component is developed, at least initially, independently. Once development is mature, the components are coupled to deliver a model of the required complexity. The implementation of the coupling can have a significant impact on the model. As the error associated with each component decreases, the errors introduced by the coupling will eventually dominate. Hence, any improvement in one of the components is unlikely to improve the performance of the overall system. The challenges associated with combining the components to create a coherent model are here termed physics–dynamics coupling. The issue goes beyond the coupling between the parameterizations and the resolved fluid dynamics. This paper highlights recent progress and some of the current challenges. It focuses on three objectives: to illustrate the phenomenology of the coupling problem with references to examples in the literature, to show how the problem can be analyzed, and to create awareness of the issue across the disciplines and specializations. The topics addressed are different ways of advancing full models in time, approaches to understanding the role of the coupling and evaluation of approaches, coupling ocean and atmosphere models, thermodynamic compatibility between model components, and emerging issues such as those that arise as model resolutions increase and/or models use variable resolutions.

Understanding the cloud response to external forcing is a major challenge for climate science. This crucial goal is complicated by intermodel differences in simulating present and future cloud cover and by observational uncertainty. This is the first formal detection and attribution study of cloud changes over the satellite era. Presented herein are CMIP5 model-derived fingerprints of externally forced changes to three cloud properties: the latitudes at which the zonally averaged total cloud fraction (CLT) is maximized or minimized, the zonal average CLT at these latitudes, and the height of high clouds at these latitudes. By considering simultaneous changes in all three properties, the authors define a coherent multivariate fingerprint of cloud response to external forcing and use models from phase 5 of CMIP (CMIP5) to calculate the average time to detect these changes. It is found that given perfect satellite cloud observations beginning in 1983, the models indicate that a detectable multivariate signal should have already emerged. A search is then made for signals of external forcing in two observational datasets: ISCCP and PATMOS-x. The datasets are both found to show a poleward migration of the zonal CLT pattern that is incompatible with forced CMIP5 models. Nevertheless, a detectable multivariate signal is predicted by models over the PATMOS-x time period and is indeed present in the dataset. Despite persistent observational uncertainties, these results present a strong case for continued efforts to improve these existing satellite observations, in addition to planning for new missions.

We examine the tropical inversion strength, measured by the estimated inversion strength (EIS), and its response to climate change in 18 models associated with phase 5 of the coupled model intercomparison project (CMIP5). While CMIP5 models generally capture the geographic distribution of observed EIS, they systematically underestimate it off the west coasts of continents, due to a warm bias in sea surface temperature. The negative EIS bias may contribute to the low bias in tropical low-cloud cover in the same models. Idealized perturbation experiments reveal that anthropogenic forcing leads directly to EIS increases, independent of “temperature-mediated” EIS increases associated with long-term oceanic warming. This fast EIS response to anthropogenic forcing is strongly impacted by nearly instantaneous continental warming. The temperature-mediated EIS change has contributions from both uniform and non-uniform oceanic warming. The substantial EIS increases in uniform oceanic warming simulations are due to warming with height exceeding the moist adiabatic lapse rate in tropical warm pools. EIS also increases in fully-coupled ocean–atmosphere simulations where $$\hbox {CO}_{2}$$ concentration is instantaneously quadrupled, due to both fast and temperature-mediated changes. The temperature-mediated EIS change varies with tropical warming in a nonlinear fashion: The EIS change per degree tropical warming is much larger in the early stage of the simulations than in the late stage, due to delayed warming in the eastern parts of the subtropical oceans. Given the importance of EIS in regulating tropical low-cloud cover, this suggests that the tropical low-cloud feedback may also be nonlinear .

Author(s): Ng, Esmond; Evans, Katherine J; Caldwell, Peter; Hoffman, Forrest M; Jackson, Charles; Kerstin, Van Dam; Leung, Ruby; Martin, Daniel F; Ostrouchov, George; Tuminaro, Raymond; Ullrich, Paul; Wild, S; Williams, Samuel | Abstract: This report presents results from the DOE-sponsored workshop titled, Advancing X-Cutting Ideas for Computational Climate Science Workshop,'' known as AXICCS, held on September 12--13, 2016 in Rockville, MD. The workshop brought together experts in climate science, computational climate science, computer science, and mathematics to discuss interesting but unsolved science questions regarding climate modeling and simulation, promoted collaboration among the diverse scientists in attendance, and brainstormed about possible tools and capabilities that could be developed to help address them. Emerged from discussions at the workshop were several research opportunities that the group felt could advance climate science significantly. These include (1) process-resolving models to provide insight into important processes and features of interest and inform the development of advanced physical parameterizations, (2) a community effort to develop and provide integrated model credibility, (3) including, organizing, and managing increasingly connected model components that increase model fidelity yet complexity, and (4) treating Earth system models as one interconnected organism without numerical or data based boundaries that limit interactions. The group also identified several cross-cutting advances in mathematics, computer science, and computational science that would be needed to enable one or more of these big ideas. It is critical to address the need for organized, verified, and optimized software, which enables the models to grow and continue to provide solutions in which the community can have confidence. Effectively utilizing the newest computer hardware enables simulation efficiency and the ability to handle output from increasingly complex and detailed models. This will be accomplished through hierarchical multiscale algorithms in tandem with new strategies for data handling, analysis, and storage. These big ideas and cross-cutting technologies for enabling breakthrough climate simulation advancements also need the "glue" of outreach and learning across the scientific domains to be successful. The workshop identified several strategies to allow productive, continuous engagement across those who have a broad knowledge of the various angles of the problem. Specific ideas to foster education and tools to make material progress were discussed. Examples include follow-on cross-cutting meetings that enable unstructured discussions of the types this workshop fostered. A concerted effort to recruit undergraduate and graduate students from all relevant domains and provide them experience, training, and networking across their immediate expertise is needed. This will broaden and expand their exposure to the future needs and solutions, and provide a pipeline of scientists with a diversity of knowledge and know-how. Providing real-world experience with subject matter experts from multiple angles may also motivate the students to attack these problems and even come up with the missing solutions.

Several recent efforts to estimate Earth's equilibrium climate sensitivity (ECS) focus on identifying quantities in the current climate which are skillful predictors of ECS yet can be constrained by observations. This study automates the search for observable predictors using data from phase 5 of the Coupled Model Intercomparison Project. The primary focus of this paper is assessing statistical significance of the resulting predictive relationships. Failure to account for dependence between models, variables, locations, and seasons is shown to yield misleading results. A new technique for testing the field significance of data-mined correlations which avoids these problems is presented. Using this new approach, all 41,741 relationships we tested were found to be explainable by chance. This leads us to conclude that data mining is best used to identify potential relationships which are then validated or discarded using physically based hypothesis testing.

In this study, a series of idealized large-eddy simulations is used to understand the relative impact of cloud-top and subcloud-layer sources of moisture on the microphysical–radiative–dynamical feedbacks in an Arctic mixed-phase stratocumulus (AMPS) cloud system. This study focuses on a case derived from observations of a persistent single-layer AMPS cloud deck on 8 April 2008 during the Indirect and Semi-Direct Aerosol Campaign near Barrow, Alaska. Moisture and moist static energy budgets are used to examine the potential impact of ice in mixed-phase clouds, specific humidity inversions coincident with temperature inversions as a source of moisture for the cloud system, and the presence of cloud liquid water above the mixed-layer top. This study demonstrates that AMPS have remarkable insensitivity to changes in moisture source. When the overlying air is dried initially, radiative cooling and turbulent entrainment increase moisture import from the surface layer. When the surface layer is dried initially, the system evolves to a state with reduced mixed-layer water vapor and increased surface-layer moisture, reducing the loss of water through precipitation and entrainment of near-surface air. Only when moisture is reduced both above and below the mixed layer does the AMPS decay without reaching a quasi-equilibrium state. A fundamental finding of this study is that, with or without cloud ice and with or without a specific humidity inversion, the cloud layer eventually extends into the temperature inversion producing a precipitation flux as a source of water into the mixed layer.