Search

The atmospheric Green's function method is a technique for modeling the response of the atmosphere to changes in the spatial field of surface temperature. While early studies applied this method to changes in atmospheric circulation, it has also become an important tool to understand changes in radiative feedbacks due to evolving patterns of warming, a phenomenon called the "pattern effect." To better study this method, this paper presents a protocol for creating atmospheric Green's functions to serve as the basis for a model intercomparison project, GFMIP. The protocol has been developed using a series of sensitivity tests performed with the HadAM3 atmosphere‐only general circulation model, along with existing and new simulations from other models. Our preliminary results have uncovered nonlinearities in the response of the atmosphere to surface temperature changes, including an asymmetrical response to warming versus cooling patch perturbations, and a change in the dependence of the response on the magnitude and size of the patches. These nonlinearities suggest that the pattern effect may depend on the heterogeneity of warming as well as its location. These experiments have also revealed tradeoffs in experimental design between patch size, perturbation strength, and the length of control and patch simulations. The protocol chosen on the basis of these experiments balances scientific utility with the simulation time and setup required by the Green's function approach. Running these simulations will further our understanding of many aspects of atmospheric response, from the pattern effect and radiative feedbacks to changes in circulation, cloudiness, and precipitation. Plain Language Summary: Many properties of the atmosphere are affected by the temperature of the ocean surface. Knowing how strong these effects are would help us to better predict global warming. The response to a given surface warming depends on where the warming occurs. To account for this, researchers sometimes simulate the response to individual patches of warming and then assume the response to an arbitrary warming pattern can be summed together from these patch responses. This is sometimes called the atmospheric Green's function method, and it works well at recreating the atmospheric response to historical temperature changes. We are organizing a Green's Function Model Intercomparison Project (GFMIP), in which participants will apply the method consistently for many climate models. This paper presents the GFMIP protocol. In the course of developing this protocol, we found that the atmospheric response to warming is not proportional in all cases: the response to surface warming is not the opposite of the response to surface cooling; warming twice as much doesn't cause twice as much of a response; and making a patch of warming twice as large doesn't cause twice as large a response. GFMIP will help us figure out how to account for this nonlinearity. Key Points: The Green's Function Model Intercomparison Project (GFMIP) explores the atmospheric response to surface temperature patch perturbationsThis paper presents the GFMIP protocol, which was generated using insights from past studies and new sensitivity testsGreen's functions reconstruct the response to historical temperatures, but nonlinearities can affect responses to other warming patterns [ABSTRACT FROM AUTHOR]

When the negatively buoyant air in the cloud downdrafts reaches the surface, it spreads out horizontally, producing cold pools. A cold pool can trigger new convective cells. However, when combined with the ambient vertical wind shear, it can also connect and upscale them into large mesoscale convective systems (MCS). Given the broad spectrum of scales of the atmospheric phenomenon involving the interaction between cold pools and the MCS, a parameterization was designed here. Then, it is coupled with a classical convection parameterization to be applied in an atmospheric model with an insufficient spatial resolution to explicitly resolve convection and the sub‐cloud layer. A new scalar quantity related to the deficit of moist static energy detrained by the downdrafts mass flux is proposed. This quantity is subject to grid‐scale advection, mixing, and a sink term representing dissipation processes. The model is then applied to simulate moist convection development over a large portion of tropical land in the Amazon Basin in a wet and dry‐to‐wet 10‐days period. Our results show that the cold pool edge parameterization improves the organization, longevity, propagation, and severity of simulated MCS over the Amazon and other different continental areas. Plain Language Summary: In nature, cold pools are formed by cold air masses descending from the low to mid‐troposphere in thunderstorms. When these drafts reach the surface, they spread out horizontally. A manifestation of cold pools is the relatively high speed at the gust front, which can lift environmental air producing new convective cells. Moreover, depending on ambient conditions, the cold pools may help organize the new convective cells, increasing their aggregation and forming the so‐called mesoscale convective systems (MCSs). MCSs, which cover hundreds to thousands of km2, significantly impact the global scale circulation, energy budget, hydrological cycle, and population safety. Forecasting MCSs is challenging for global circulation models (GCM) due to the broad spectrum of scales of the involved atmospheric phenomenon. The computational limitations, at present and for some time to come, do not allow running in real‐time GCMs, which explicitly solves all relevant scales of motion. This paper describes a methodology to account for essential interplays between cold pools edges and moist convection to be applied in the GCMs of weather and climate forecasting. We show that the method improves the model simulation of the main types of MCSs over the Amazon Basin and other continental areas. Key Points: A model scheme for including effect of cold pools edges in triggering new convective cells and storm propagation is presentedThe scheme is coupled with a convection parameterization and applied in the modeling of moist convective systemsThe method improves the organization, longevity, propagation, and severity of the simulated mesoscale convective systems [ABSTRACT FROM AUTHOR]

We address the question of how to use a machine learned (ML) parameterization in a general circulation model (GCM), and assess its performance both computationally and physically. We take one particular ML parameterization (Guillaumin & Zanna, 2021, https://doi.org/10.1002/essoar.10506419.1) and evaluate the online performance in a different model from which it was previously tested. This parameterization is a deep convolutional network that predicts parameters for a stochastic model of subgrid momentum forcing by mesoscale eddies. We treat the parameterization as we would a conventional parameterization once implemented in the numerical model. This includes trying the parameterization in a different flow regime from that in which it was trained, at different spatial resolutions, and with other differences, all to test generalization. We assess whether tuning is possible, which is a common practice in GCM development. We find the parameterization, without modification or special treatment, to be stable and that the action of the parameterization to be diminishing as spatial resolution is refined. We also find some limitations of the machine learning model in implementation: (a) tuning of the outputs from the parameterization at various depths is necessary; (b) the forcing near boundaries is not predicted as well as in the open ocean; (c) the cost of the parameterization is prohibitively high on central processing units. We discuss these limitations, present some solutions to problems, and conclude that this particular ML parameterization does inject energy, and improve backscatter, as intended but it might need further refinement before we can use it in production mode in contemporary climate models. Plain Language Summary: This paper discusses how machine learning can be used to make climate models more accurate. Specifically, we import an existing machine learning model that predicts how small eddies (in the order of 10–100 km) in the ocean affect larger currents. We test this machine learning model in a different ocean circulation model than the one it was originally designed for, and found that it worked well. However, we also found some limitations: the model works differently at different depths in the ocean, and it does not work as well near the coasts of the ocean. We also found that the model takes a long time to run on normal computers. Overall, we concluded that the model is promising, but more work is needed to make it work well in realistic situations. Key Points: A stochastic‐deep learning model is implemented in an ocean circulation model, MOM6We evaluate the online performance of the stochastic‐deep learning model as a subgrid parameterizationWe identify certain limitations of the machine learned parameterization which otherwise has the potential to improve specific metrics [ABSTRACT FROM AUTHOR]

Precipitation is highly variable on spatial scales smaller than a typical general circulation model (GCM) grid box. Neglecting this subgrid‐scale variability can impact the simulation of the underlying land surface and its coupling with the atmosphere. In this study we present a scheme to stochastically disaggregate precipitation from an atmospheric grid to an underlying mosaic of surface tiles, using an assumed probability distribution of subgrid precipitation derived from observations. Unlike previous GCM‐based schemes, our approach includes a treatment of memory to allow persistence of subgrid features. In single column model experiments, we demonstrate the ability of the scheme to reproduce observed precipitation statistics at the spatial scale of the surface tiles. The root mean square error of hourly spatial standard deviation of precipitation is reduced from 1.90 to 0.96 mm hr−1 when disaggregation is applied. The scheme increases the spatial standard deviations of surface heat fluxes, soil moisture and temperature, and we show that incorporating memory amplifies these increases by a factor of 2–4. We also document increases in mean precipitation runoff and the Bowen ratio. Plain Language Summary: Global climate and weather models (GCMs), used for global warming projections and near‐term weather prediction, represent the atmosphere on a grid, while the underlying land surface is often composed of smaller elements. In most models, all surface elements under a given atmospheric grid point receive the same precipitation rate. However, real‐world precipitation rates vary significantly over the scale of a GCM grid cell. This paper presents an approach to represent subgrid precipitation variability in a GCM. Grid mean precipitation is randomly distributed across surface elements, such that surface precipitation reproduces observed statistics. The approach is found to increase grid cell‐averaged surface runoff as well as the spatial variances of surface temperature, soil moisture, and surface sensible heat and moisture fluxes. Key Points: A parameterization is developed to stochastically distribute precipitation across subgrid surface tiles within a general circulation model (GCM)In single column experiments, the scheme reproduces observed precipitation statistics at the scale of land surface tilesThe scheme increases surface runoff, Bowen ratio, and the spatial variance of surface properties at the ARM Southern Great Plains site [ABSTRACT FROM AUTHOR]

This paper presents the response to anthropogenic forcing in the GISS‐E2.1 climate models for the 21st century Shared Socioeconomic Pathways emission scenarios within the Coupled Model Intercomparison Project Phase 6 (CMIP6). The experiments were performed using an updated and improved version of the NASA Goddard Institute for Space Studies (GISS) coupled general circulation model that includes two different versions for atmospheric composition: A non‐interactive version (NINT) with prescribed composition and a tuned aerosol indirect effect and the One‐Moment Aerosol model (OMA) version with fully interactive aerosols which includes a parameterized first indirect aerosol effect on clouds. The effective climate sensitivities are 3.0°C and 2.9°C for the NINT and OMA models, respectively. Each atmospheric version is coupled to two different ocean general circulation models: The GISS ocean model (E2.1‐G) and HYCOM (E2.1‐H). We describe the global mean responses for all future scenarios and spatial patterns of change for surface air temperature and precipitation for four of the marker scenarios: SSP1‐2.6, SSP2‐4.5, SSP4‐6.0, and SSP5‐8.5. By 2100, global mean warming ranges from 1.5°C to 5.2°C relative to 1,850–1,880 mean temperature. Two high‐mitigation scenarios SSP1‐1.9 and SSP1‐2.6 limit the surface warming to below 2°C by the end of the 21st century, except for the NINT E2.1‐H model that simulates 2.2°C of surface warming. For the high emission scenario SSP5‐8.5, the range is 4.6–5.2°C at 2100. Due to about 15% larger effective climate sensitivity and stronger transient climate response in both NINT and OMA CMIP6 models compared to CMIP5 versions, there is a stronger warming by 2100 in the SSP emission scenarios than in the comparable Representative Concentration Pathway (RCP) scenarios in CMIP5. Changes in sea ice area are highly correlated to global mean surface air temperature anomalies and show steep declines in both hemispheres, with the largest sea ice area decreases occurring during September in the Northern Hemisphere in both E2.1‐G (−1.21 × 106 km2/°C) and E2.1‐H models (−0.94 × 106 km2/°C). Both coupled models project decreases in the Atlantic overturning stream function by 2100. The largest decrease of 56%–65% in the 21st century overturning stream function is produced in the warmest scenario SSP5‐8.5 in the E2.1‐G model, comparable to the reduction in the corresponding CMIP5 GISS‐E2 RCP8.5 simulation. Both low‐end scenarios SSP1‐1.9 and SSP1‐2.6 also simulate substantial reductions of the overturning (9%–37%) with slow recovery of about 10% by the end of the 21st century (relative to the maximum decrease at the middle of the 21st century). Plain Language Summary: The projections of future climate change are uncertain because they are dependent on different possible scenarios of human‐caused emissions and their interaction with natural forcings, internal climate variability, and inter‐model differences. This paper presents the results of the climate model of the NASA Goddard Institute for Space Studies, GISS‐E2.1, for the anthropogenically forced climate response for the twenty first century Shared Socioeconomic Pathways emission scenarios within the Coupled Model Intercomparison Project Phase 6 (CMIP6). The sensitivity of the model response to different magnitudes of anthropogenic forcings in the twenty‐first‐century scenarios were performed using two different versions for the atmospheric composition and two different ocean general circulation models. Compared to CMIP5 GISS‐E2 versions, the CMIP6 GISS‐E2.1 climate model shows a stronger warming by 2100 in comparable scenarios due to larger effective climate sensitivity and transient climate response. Both climate models with two different ocean components project decreases in the Atlantic overturning stream function by 2100. Key Points: GISS E2.1 model with different configurations is used to carry out 134 Shared Socioeconomic Pathway (SSP) experimentsGISS‐E2.1 climate model shows a stronger warming by 2,100 in comparable Representative Concentration Pathway scenarios in CMIP5 due to larger effective climate sensitivity and stronger transient climate responseBoth coupled models, E2.1‐G and E2.1‐H, project decreases in the Atlantic overturning stream function by 2100 with the largest decrease in the warmest scenario SSP5‐8.5 in the E2.1‐G model [ABSTRACT FROM AUTHOR]

In geoscience and other fields, researchers use models as a simplified representation of reality. The models include processes that often rely on uncertain parameters that reduce model performance in reflecting real‐world processes. The problem is commonly addressed by adapting parameter values to reach a good match between model simulations and corresponding observations. Different optimization tools have been successfully applied to address this task of model calibration. However, seeking one best value for every single model parameter might not always be optimal. For example, if model equations integrate over multiple real‐world processes which cannot be fully resolved, it might be preferable to consider associated model parameters as random parameters. In this paper, a random parameter is drawn from a wide probability distribution for every singe model simulation. We developed an optimization approach that allows us to declare certain parameters random while optimizing those that are assumed to take fixed values. We designed a corresponding variant of the well known Covariance Matrix Adaption Evolution Strategy (CMA‐ES). The new algorithm was applied to a global biogeochemical circulation model to quantify the impact of zooplankton mortality on the underlying biogeochemistry. Compared to the deterministic CMA‐ES, our new method converges to a solution that better suits the credible range of the corresponding random parameter with less computational effort. Plain Language Summary: In geoscience and other fields, researchers use models as a simplified representation of reality. The models include processes that often rely on uncertain parameters that reduce model performance in reflecting real‐world processes. The problem is commonly addressed by adapting parameter values to reach a good match between model simulations and corresponding observations. However, seeking one best value for every single model parameter might not always be optimal. For example, if model equations integrate over multiple real‐world processes it might be preferable to consider associated model parameters as random parameters that are drawn from a wide probability distribution for every single model simulation. To explore this further, we modified a popular variant of the Covariance Matrix Adaption Evolution Strategy (CMA‐ES) in order to tune the other parameters. The new algorithm was applied to a global biogeochemical circulation model to quantify the impact of zooplankton mortality on the underlying biogeochemistry. Compared to the deterministic CMA‐ES our new method converges to a solution that better suits the credible range of the corresponding random parameter. Key Points: Model calibration is an important task in the view of parametric uncertaintiesWe allow specific uncertain model parameters to be random while optimizing other model parameters, providing an efficient tool for the taskWe apply it to a global biogeochemical circulation model to quantify the impact of zooplankton mortality on the underlying biogeochemistry [ABSTRACT FROM AUTHOR]

We develop a doubly periodic version of the Simple Convection‐Permitting E3SM Atmosphere Model (SCREAM) to provide an efficient configuration for this global convection permitting model (GCPM), akin to a single column model often found in conventional general circulation models. The design details are explained, in addition to the extensive case library associated with the doubly periodic SCREAM (DP‐SCREAM) configuration. We demonstrate that doubly periodic cloud resolving models are useful tools to explore the horizontal resolution sensitivity of GCPMs, in addition to replicating biases seen in the global models. Using DP‐SCREAM, we show that SCREAM exhibits behaviors of a scale aware model as it is able to naturally partition between sub‐grid scale (SGS) and resolved vertical transport across the gray zone of turbulence. We show that SCREAM is reasonably scale insensitive when run at resolutions from 1 to 5 km, but can exhibit sensitivity, particularly for the shallow convective regime, when run at resolutions approaching that of large eddy simulations. We conclude that SGS parameterization improvements are likely needed to reduce this scale sensitivity. Plain Language Summary: Advances in computational resources have allowed climate simulations to be performed with very high resolution, which provides higher quality results. However, these simulations require a lot of time and computer resources to perform, which makes these models difficult to use for the common scientist. In this paper we develop a high‐resolution configuration which focuses on a specific point on the globe, enabling it to run fast and to use minimal computational resources. This allows users and developers to gauge how the model may perform before doing a computationally intensive global simulation. We show that this faster configuration is a useful tool to replicate problems that are found in the global model and is a valuable way to assess sensitivities of the model, particularly pertaining to choices made in its resolution. Key Points: Doubly periodic configurations represent an efficient framework to supplement global cloud resolving modelsShallow convection is sensitive to horizontal resolution in Simple Convection‐Permitting E3SM Atmosphere Model (SCREAM) but stratiform clouds are more robustSCREAM exhibits the behaviors of a scale aware model and can represent a range of cloud regimes [ABSTRACT FROM AUTHOR]

The mixing of tracers by mesoscale eddies, parameterized in many ocean general circulation models (OGCMs) as a diffusive‐advective process, contributes significantly to the distribution of tracers in the ocean. In the ocean interior, diffusive contribution occurs mostly along the direction parallel to local neutral density surfaces. However, near the surface of the ocean, small‐scale turbulence and the presence of the boundary itself break this constraint and the mesoscale transport occurs mostly along a plane parallel to the ocean surface (horizontal). Although this process is easily represented in OGCMs with geopotential vertical coordinates, the representation is more challenging in OGCMs that use a general vertical coordinate, where surfaces can be tilted with respect to the horizontal. We propose a method for representing the diffusive horizontal mesoscale fluxes within the surface boundary layer of general vertical coordinate OGCMs. The method relies on regridding/remapping techniques to represent tracers in a geopotential grid. Horizontal fluxes are calculated on this grid and then remapped back to the native grid, where fluxes are applied. The algorithm is implemented in an ocean model and tested in idealized and realistic settings. Horizontal diffusion can account for up to 10% of the total northward heat transport in the Southern Ocean and Western boundary current regions of the Northern Hemisphere. It also reduces the vertical stratification of the upper ocean, which results in an overall deepening of the surface boundary layer depth. Finally, enabling horizontal diffusion leads to meaningful reductions in the near‐surface global bias of potential temperature and salinity. Plain Language Summary: Mesoscale ocean eddies, which are analogous to the weather systems in the atmosphere, are crucial to the distribution of heat, salt, carbon, and nutrients throughout the global ocean. Most of the ocean models used in climate simulations do not have enough horizontal resolution to resolve these eddies and, therefore, their effects must be parameterized. In the ocean interior, where no appreciable heat or mass is exchanged across density surfaces, the mixing of tracers due to mesoscale eddies occurs along surfaces of constant density. However, as ocean boundaries are approached, mixing occurs on a plane parallel to the boundary. For example, near the surface of the ocean, which is where the scheme presented here is designed to work, this plane is mostly horizontal. There is a class of ocean models that use a vertical coordinate system whose layer thicknesses can vary horizontally, thus complicating the implementation of horizontal diffusive parameterizations. This paper presents and evaluates a method that allows horizontal fluxes to be calculated and applied within the surface layer of this class of ocean models. Key Points: General vertical coordinate ocean models lack the option to shift diffusive fluxes from neutral surfaces to a purely horizontal directionRegridding/remapping techniques are used to represent tracers in z*‐coordinates, where horizontal fluxes are easily appliedThe method reduces global tracer biases in forced simulations, regardless of the coordinate system employed [ABSTRACT FROM AUTHOR]

This paper describes an implementation of the combined hybrid‐parallel prediction (CHyPP) approach of Wikner et al. (2020), https://doi.org/10.1063/5.0005541 on a low‐resolution atmospheric global circulation model (AGCM). The CHyPP approach combines a physics‐based numerical model of a dynamical system (e.g., the atmosphere) with a computationally efficient type of machine learning (ML) called reservoir computing to construct a hybrid model. This hybrid atmospheric model produces more accurate forecasts of most atmospheric state variables than the host AGCM for the first 7–8 forecast days, and for even longer times for the temperature and humidity near the earth's surface. It also produces more accurate forecasts than a model based only on ML, or a model that combines linear regression, rather than ML, with the AGCM. The potential of the CHyPP approach for climate research is demonstrated by a 10‐year long hybrid model simulation of the atmospheric general circulation, which shows that the hybrid model can simulate the general circulation with substantially smaller systematic errors and more realistic variability than the host AGCM. Plain Language Summary: This paper presents a computationally efficient novel approach to construct a hybrid model of the atmosphere by combining a physics‐based model of the global atmospheric circulation with a machine learning component. The primary purpose of the hybrid model is to produce quantitative weather forecasts on the same grid as the physics‐based model. It is found that the hybrid model produces more accurate forecasts than the host physics‐based model for the first 7–8 forecast days for most forecast variables, and for even longer times for the temperature and humidity near the Earth's surface. Furthermore, the hybrid model is found to simulate the climate with substantially smaller systematic errors and more realistic temporal variability than the host model. Key Points: A hybrid model incorporating machine learning produces more accurate forecasts and more realistic climate than the host physics‐based modelThe hybrid model states are more realistically balanced and have substantially lower biases than the host modelThe hybrid model produces more realistic atmospheric variability than the host model at time scales shorter than about a week [ABSTRACT FROM AUTHOR]

Radiative transfer models are widely applied in climate models to simulate vertical temperature perturbations caused by external radiative forcings. A large part of radiative transfer models is the infrared gaseous spectral transmittance scheme, which quantifies the longwave atmospheric absorption. A rapid infrared gaseous spectral transmittance scheme, called the Optimized alternate Mapping Correlated K‐Distribution model (OMCKD), is introduced in this paper. To improve the accuracy of our scheme without increasing pseudo‐monochromatic calculations, we introduce the optimal iteration method to automatically tune the equivalent absorption coefficients in the cumulative probability subspace. In addition, a new expression weighted by black‐body radiation is introduced to calculate the equivalent absorption coefficient. The OMCKD simulates heating rate and radiation flux with errors of less than 0.12 K d−1 and 0.35 W m−2, respectively, below stratopause for standard atmospheric profiles. The OMCKD is also evaluated and compared with the rapid radiative transfer model for general circulation models (RRTMG) in realistic atmospheric profiles. We found that OMCKD can accurately produce heating rates and generally captures radiative forcings associated with large perturbations to the concentrations of main greenhouse gases. Furthermore, the number of pseudo‐monochromatic calculations in OMCKD is 11.4% less than that in RRTMG, which indicates less computational cost. Plain Language Summary: Radiative models are practical tools for studying the greenhouse effect. The thermal infrared radiation absorbed by greenhouse gases is quantified by a specially designed program called the radiation scheme. The scheme divides the thermal infrared band into hundreds of intervals and calculates each interval's gaseous radiation absorption ability separately. Generally, increasing the number of intervals improves simulation accuracy, but the computational cost is proportional to the number of intervals. We have developed a rapid infrared gaseous spectral transmittance scheme called the Optimized alternate Mapping Correlated K‐Distribution model (OMCKD). Our scheme introduces an optimal method to adjust the radiation absorption ability of each interval with the line‐by‐line integration scheme (accurate but extremely time‐consuming) as the standard. We succeeded in improving the accuracy of our scheme without increasing the number of intervals. The OMCKD simulates heating rate and radiation flux with errors of less than 0.12 K d−1 and 0.35 W m−2, respectively, below stratopause for standard atmospheric conditions. The OMCKD generally captures the changes in the net fluxes and heating rates associated with perturbations to the concentrations of five main greenhouse gases. The stability of OMCKD implies its potential to improve the accuracy of climate predictions under higher CO2 concentration scenarios. Key Points: A rapid infrared gaseous spectral transmittance scheme based on the correlated K‐distribution method and the optimal method is developedThe optimal iteration method can significantly improve the accuracy of heating rates simulated with the correlated K‐distribution modelOur infrared gaseous spectral transmittance scheme maintains high accuracy across various carbon dioxide concentration scenarios [ABSTRACT FROM AUTHOR]

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. Plain Language Summary: This paper describes the new global 3.25 km version of the Energy Exascale Earth System Model (E3SM) atmosphere model and its behavior in a 40‐day northern‐hemisphere wintertime simulation. In exchange for huge computational expense, this high‐resolution model avoids many but not all biases common in lower‐resolution models. It also captures several types of extreme weather that would simply not be resolved in lower‐resolution models. Several opportunities for further development are identified. Key Points: Describes the Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM)SCREAM performs well in a 40‐day boreal winter simulation at 3.25 km ΔxResolving deep convection solves many long‐standing climate model biases [ABSTRACT FROM AUTHOR]

This study further evaluates the modeling approach of Jia et al. (2019), https://doi.org/10.1029/2019gl081960 (JEA19) to investigate the potential effects of mesoscale sea surface temperature (SST) variability on the atmospheric circulation. The approach employs a global atmospheric circulation model coupled to a slab ocean model to produce two ensembles of simulations: one in which the SST exhibits realistic mesoscale variability, and another in which the mesoscale SST variability is suppressed. The results of the present study suggest that the modeling approach yields the desired SST differences between the two ensembles at the mesoscales. They also show, however, that the approach can lead to undesirable SST differences at the large scales, if the prescribed pair of oceanic heat flux fields have large scale differences. In the experiments of JEA19, such forced, large scale SST differences dominate over the large scale differences that may develop in response to changes in the large scale atmospheric circulation by nonlinear interactions with the SST. This result suggests that finding a proper pair of estimates of the oceanic heat flux fields is necessary for the investigation of the potential upscale impact of mesoscale SST variability on the atmosphere by the approach of JEA19. The paper concludes with proposing a potential improved strategy to obtain such a pair of estimates. Plain Language Summary: This study evaluates a modeling approach to investigate the potential effects of mesoscale (scales smaller than 100 km) ocean sea surface temperature (SST) variability on the atmospheric circulation. The approach employs a global atmospheric circulation model coupled to a thermodynamical model of the ocean. Two ensembles of model simulations are prepared: one in which the ocean exhibits realistic mesoscale SST variability, and another in which that variability is suppressed. The results suggest that the modeling approach yields the desired mesoscale SST differences. They also show, however, that the approach can lead to undesirable large scale SST differences, if there are large scale differences between the prescribed oceanic heat transport fields of the two ensembles. In the experiments of JEA19, such forced, large scale SST differences dominate over the large scale differences that may develop in response to changes in the large scale atmospheric circulation by two‐way interactions with the SST. This result suggests that finding a proper pair of estimates of the oceanic heat transport fields is necessary for the investigation of the potential upscale impact of mesoscale SST variability on the atmosphere by the approach of JEA19. The paper concludes with proposing a potential improved strategy to obtain such a pair of estimates. Key Points: The large scale response in an uncoupled atmospheric simulation is the response of the model rather than that of the atmosphereA coupled modeling approach is necessary to investigate the large scale atmospheric response to mesoscale sea surface temperature variabilityA pair of oceanic heat fluxes that do not have large scale differences is critical to produce a realistic large scale atmospheric response [ABSTRACT FROM AUTHOR]

To assess the effect of ocean‐atmosphere coupling in the climate response to forced sea ice loss, the Polar Amplification Model Intercomparison Project protocol includes centennial coupled atmosphere‐ocean general circulation model simulations with imposed sea ice loss. The protocol, which specifies sea ice concentration and thickness distribution targets, does not prescribe a method for achieving them. Although different methods for imposing sea ice loss (or growth) in models have been documented, testing of the method‐dependence of the resulting climate responses has been limited. Achieving the targeted sea ice state has proven to be challenging using the so‐called ghost‐flux nudging method, which induces ice melt from below, as this method does not constrain the partitioning between thickness and concentration. We propose, describe and test a simple method that combines the advantages of direct sea ice nudging and ghost‐flux nudging. The hybrid nudging method better captures the partitioning between thickness and concentration while conserving total water content. We document how this novel sea ice constraining method reaches specific targets, enhances surface turbulent heat flux responses to sea ice loss, and induces tropospheric warming for both polar regions. Plain Language Summary: The Arctic is warming faster than the global average due to several processes that, once combined, lead to so‐called Arctic Amplification. Part of this anomalous polar warming comes from an intense reduction in ice cover allowing heat into the ocean, warming the Arctic ocean near the surface, hence melting more ice. A joint effort by several climate modeling groups called the Polar Amplification Model Intercomparison Project (PAMIP) aims at better understanding Arctic Amplification through a coordinated set of climate simulations. Among this ensemble of simulations is a set of centennial simulations performed with fully coupled state‐of‐the‐art climate models. In these experiments, Arctic (and Antarctic) sea ice are forced to reach specific states in order to better isolate Arctic Amplification and sea ice loss from the rest of anthropogenic global warming. In this paper, we propose a simple technique to nudge sea ice models to specific states such as prescribed by PAMIP. This new method combines advantages from existing techniques to improve the control over the extent and the thickness of the ice. We document how our novel method produces atmospheric warming and circulation responses that are more consistent with previously published results than the existing techniques. Key Points: Already established sea ice nudging method leads to issues in the partitioning between sea ice concentration (SIC) and sea ice thickness (SIT)New hybrid nudging method significantly improves capture of both SIC and SIT targets for Polar Amplification Model Intercomparison Project, in particular for Arctic sea iceHybrid scheme increases Arctic turbulent heat flux and warming [ABSTRACT FROM AUTHOR]

Abstract: The Global Coupled 3 (GC3) configuration of the Met Office Unified Model is presented. Among other applications, GC3 is the basis of the United Kingdom's submission to the Coupled Model Intercomparison Project 6 (CMIP6). This paper documents the model components that make up the configuration (although the scientific descriptions of these components are in companion papers) and details the coupling between them. The performance of GC3 is assessed in terms of mean biases and variability in long climate simulations using present‐day forcing. The suitability of the configuration for predictability on shorter time scales (weather and seasonal forecasting) is also briefly discussed. The performance of GC3 is compared against GC2, the previous Met Office coupled model configuration, and against an older configuration (HadGEM2‐AO) which was the submission to CMIP5. In many respects, the performance of GC3 is comparable with GC2, however, there is a notable improvement in the Southern Ocean warm sea surface temperature bias which has been reduced by 75%, and there are improvements in cloud amount and some aspects of tropical variability. Relative to HadGEM2‐AO, many aspects of the present‐day climate are improved in GC3 including tropospheric and stratospheric temperature structure, most aspects of tropical and extratropical variability and top‐of‐atmosphere and surface fluxes. A number of outstanding errors are identified including a residual asymmetric sea surface temperature bias (cool northern hemisphere, warm Southern Ocean), an overly strong global hydrological cycle and insufficient European blocking. [ABSTRACT FROM AUTHOR]

We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed‐layer model (MLM) with entrainment and detrainment, an Ekman MLM, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting three‐dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed‐layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid‐latitude westerly jet and of the Inter‐tropical Convergence Zone (ITCZ) in response to sea‐ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed‐layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization. Plain Language Summary: Coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the climate system. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this paper, we manage to overcome this problem and develop a globally coupled ocean model hierarchy that can turn on and off the processes of mixed‐layer entrainment and Ekman flow. We use the hierarchy to study the impact of Arctic sea‐ice loss on the climate system. We find that the effect of mixed‐layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid‐latitude westerly jet and the Inter‐tropical Convergence Zone. Key Points: We develop a flux‐corrected globally coupled ocean‐model hierarchy that explicitly includes mixed‐layer entrainment and Ekman flowFlux corrections need to contain components due to interannual variability to reproduce vertical gradient of ocean temperature and salinityThe hierarchy applied to Arctic sea‐ice loss reveals that the Inter‐tropical Convergence Zone is sensitive to the entrainment and Ekman flow parameterization [ABSTRACT FROM AUTHOR]

This paper describes the GISS‐E2.1 contribution to the Coupled Model Intercomparison Project, Phase 6 (CMIP6). This model version differs from the predecessor model (GISS‐E2) chiefly due to parameterization improvements to the atmospheric and ocean model components, while keeping atmospheric resolution the same. Model skill when compared to modern era climatologies is significantly higher than in previous versions. Additionally, updates in forcings have a material impact on the results. In particular, there have been specific improvements in representations of modes of variability (such as the Madden‐Julian Oscillation and other modes in the Pacific) and significant improvements in the simulation of the climate of the Southern Oceans, including sea ice. The effective climate sensitivity to 2 × CO2 is slightly higher than previously at 2.7–3.1°C (depending on version) and is a result of lower CO2 radiative forcing and stronger positive feedbacks. Plain Language Summary: This paper describes the latest iteration of the National Aeronautics and Space Administration (NASA) Goddard Institute for Space Studies (GISS) climate model, which will be used for understanding historical climate change and to make projections for the future. We compare the model output to a wide range of observations over the recent era (1979–2014) and show that there has been a significant increase in how well the model performs compared to the previous version from 2014, particularly in the Southern Ocean, though some persistent biases remain. The model has a temperature response to the increase of carbon dioxide that is slightly higher than previous versions but is well within the range expected from observational and past climate constraints. Key Points: GISS‐E2.1 is an updated climate model version for use within the CMIP6 projectAtmospheric composition is calculated consistently in all model versionsResults demonstrate a significant improvement in skill in a climate model without changes to atmospheric resolution [ABSTRACT FROM AUTHOR]

Data‐driven parameterizations offer considerable potential for improving the fidelity of General Circulation Models. However, ensuring that these remain consistent with the governing equations while still producing stable simulations remains a challenge. In this paper, we propose a combined Variational‐Multiscale (VMS) Artificial Neural Network (ANN) discretization which makes no a priori assumptions on the model form, and is only restricted in its accuracy by the precision of the ANN. Using a simplified problem, we demonstrate that good predictions of the required closure terms can be obtained with relatively compact ANN architectures. We then turn our attention to the stability of the VMS‐ANN discretization in the context of a single implicit time step. It is demonstrated that the ANN parameterization introduces nonphysical solutions to the governing equations that can significantly affect or prevent convergence. We show that enriching the training data with nonphysical states from intra‐time step iterations is an effective remedy. This indicates that the lack of representative ANN‐induced errors in our original, exact training inputs underpin the observed instabilities. In turn, this suggests that data set enrichment might aid in resolving instabilities that develop over several time steps. Plain Language Summary: Computer models of weather and climate have coarse resolutions, which prevent them from accurately predicting the effect of small atmospheric motions, for instance in low‐lying clouds, on the global climate. Recent studies indicate that improvements might be had by using models for the small motions that are informed purely by data, for example Artificial Neural Network (ANN) models. However, capitalizing on this potential is challenging in practice, since ANNs can introduce instabilities to the numerical model for larger‐scale motions. In this study, we zero in on these instabilities by introducing a model framework in which prediction accuracy is only limited by an ANN's ability to predict the net influence of small atmospheric motions. We contextualize this framework by contrasting it to more commonly used approaches. Using a simple test case based on the motions underpinning the development of low clouds, we illustrate that standard training procedures do not prepare ANN models for their interaction with the rest of the numerical model, and that this can contribute to numerical instability. We demonstrate an effective remedy for a specific type of instability, and suggest how this technique might also be used to treat other types of instabilities prompted by ANN models. Key Points: Artificial Neural Networks (ANNs) can parameterize exact closure terms of a Variational Multiscale method without a priori model form assumptionsThe ANN parameterization outperforms state‐of‐the‐art schemes in offline testsInstabilities in the numerical model comprising the ANN parameterization are amendable through data set enrichment [ABSTRACT FROM AUTHOR]

This paper describes the main characteristics of CNRM‐CM6‐1, the fully coupled atmosphere‐ocean general circulation model of sixth generation jointly developed by Centre National de Recherches Météorologiques (CNRM) and Cerfacs for the sixth phase of the Coupled Model Intercomparison Project 6 (CMIP6). The paper provides a description of each component of CNRM‐CM6‐1, including the coupling method and the new online output software. We emphasize where model's components have been updated with respect to the former model version, CNRM‐CM5.1. In particular, we highlight major improvements in the representation of atmospheric and land processes. A particular attention has also been devoted to mass and energy conservation in the simulated climate system to limit long‐term drifts. The climate simulated by CNRM‐CM6‐1 is then evaluated using CMIP6 historical and Diagnostic, Evaluation and Characterization of Klima (DECK) experiments in comparison with CMIP5 CNRM‐CM5.1 equivalent experiments. Overall, the mean surface biases are of similar magnitude but with different spatial patterns. Deep ocean biases are generally reduced, whereas sea ice is too thin in the Arctic. Although the simulated climate variability remains roughly consistent with CNRM‐CM5.1, its sensitivity to rising CO2 has increased: the equilibrium climate sensitivity is 4.9 K, which is now close to the upper bound of the range estimated from CMIP5 models. Key Points: Description of CNRM‐CM6‐1 model components, their coupling, and tuning procedures are describedHistorical simulations and DECK experiments are assessedPreindustrial simulation is stable and mean climate and variability in historical runs is realistic [ABSTRACT FROM AUTHOR]

Parameters in climate models are usually calibrated manually, exploiting only small subsets of the available data. This precludes both optimal calibration and quantification of uncertainties. Traditional Bayesian calibration methods that allow uncertainty quantification are too expensive for climate models; they are also not robust in the presence of internal climate variability. For example, Markov chain Monte Carlo (MCMC) methods typically require O(105) model runs and are sensitive to internal variability noise, rendering them infeasible for climate models. Here we demonstrate an approach to model calibration and uncertainty quantification that requires only O(102) model runs and can accommodate internal climate variability. The approach consists of three stages: (a) a calibration stage uses variants of ensemble Kalman inversion to calibrate a model by minimizing mismatches between model and data statistics; (b) an emulation stage emulates the parameter‐to‐data map with Gaussian processes (GP), using the model runs in the calibration stage for training; (c) a sampling stage approximates the Bayesian posterior distributions by sampling the GP emulator with MCMC. We demonstrate the feasibility and computational efficiency of this calibrate‐emulate‐sample (CES) approach in a perfect‐model setting. Using an idealized general circulation model, we estimate parameters in a simple convection scheme from synthetic data generated with the model. The CES approach generates probability distributions of the parameters that are good approximations of the Bayesian posteriors, at a fraction of the computational cost usually required to obtain them. Sampling from this approximate posterior allows the generation of climate predictions with quantified parametric uncertainties. Plain Language Summary: Calibrating climate models with available data and quantifying their uncertainties are essential to make climate predictions accurate and actionable. A primary source of uncertainties in climate models comes from representation of small‐scale processes such as moist convection. Parameters in convection schemes and other parameterizations are usually calibrated by hand, using only a small fraction of data that are available. As a result, the calibration process may miss information about the small‐scale processes in question. This paper presents a proof‐of‐concept, in an idealized setting, of how parameters in climate models can be calibrated using a substantial fraction of the available data, and how uncertainties in the parameters can be quantified. We employ a new algorithm, called calibrate‐emulate‐sample (CES), which makes such calibration and uncertainty quantification feasible for computationally expensive climate models. CES reduces the hundreds of thousands of model runs usually required to quantify uncertainties in computer models to hundreds, thereby achieving about a factor 1,000 speedup. It leads to more robust calibration and uncertainty quantification in the presence of noise arising from chaotic variability of the climate system. We show how uncertainties in climate model parameters can be translated into quantified uncertainties of climate predictions through ensemble integrations. Key Points: We use time averaged climate statistics to calibrate convective parameters and quantify their uncertaintiesWe demonstrate use of the calibrate‐emulate‐sample algorithm to provide efficient calibration and uncertainty quantificationParametric uncertainty in climate predictions is quantified by sampling from the learnt convective parameter distribution [ABSTRACT FROM AUTHOR]

We present a neutral diffusion operator appropriate for an ocean model making use of general vertical coordinates. The diffusion scheme uses polynomial reconstructions in the vertical, along with a horizontally local but vertically nonlocal stencil for estimates of tracer fluxes. These fluxes are calculated on a vertical grid that is the superset of model columns in a neutral density space. Using flux‐limiters, the algorithm dissipates tracer extrema locally, and no new extrema are created. A demonstration using a linear equation of state in an idealized configuration shows that the algorithm is perfectly neutral. When using the nonlinear TEOS‐10 equation of state with a constant reference pressure, the algorithm compares nearly exactly to a case discretized onto isopycnal surfaces and using along‐layer diffusion. The algorithm's cost is comparable to that of tracer advection and can be readily implemented into ocean general circulation models. Plain Language Summary: One prominent effect of ocean turbulence is the mixing of quantities such as temperature, salinity, and other tracers carried by the seawater. For turbulence associated with geostrophic motions, this mixing occurs along neutral directions, which are directions where buoyancy does not change. In computer models that do not explicitly represent such geostrophic eddy motions, this type of turbulent mixing is parameterized by diffusion oriented along the neutral directions. A common means to represent neutral diffusion is via rotating the diffusion fluxes to align with neutral directions (which are typically not aligned with the model grid lines). Unfortunately, this rotation can result in spurious numerical artifacts such as the creation of new minimum or maximum values of tracer concentration (referred to as extrema). In this paper, we propose a new algorithm for calculating and applying neutral diffusion in sublayers of a numerical model. We show that this new method does not create extrema, and it provides a physically accurate representation of turbulent diffusive mixing along neutral directions. Key Points: A new discretization of a neutral diffusion operator is detailed, with the new method not based on a rotated tensorThe new algorithm overcomes many numerical deficiencies of previous approachesIdealized test cases show minimal spurious dianeutral diffusion for linear and nonlinear equations of state [ABSTRACT FROM AUTHOR]

This study presents the version of the LMDZ global atmospheric model used as the atmospheric component of the Institut Pierre Simon Laplace coupled model (IPSL‐CM6A‐LR) to contribute to the 6th phase of the international Coupled Model Intercomparison Project (CMIP6). This LMDZ6A version includes original convective parameterizations that define the LMDZ "New Physics": a mass flux parameterization of the organized structures of the convective boundary layer, the "thermal plume model," and a parameterization of the cold pools created by reevaporation of convective rainfall. The vertical velocity associated with thermal plumes and gust fronts of cold pools are used to control the triggering and intensity of deep convection. Because of several shortcomings, the early version 5B of this New Physics was worse than the previous "Standard Physics" version 5A regarding several classical climate metrics. To overcome these deficiencies, version 6A includes new developments: a stochastic triggering of deep convection, a modification of the thermal plume model that allows the representation of stratocumulus and cumulus clouds in a unified framework, an improved parameterization of very stable boundary layers, and the modification of the gravity waves scheme targeting the quasi‐biennal oscillation in the stratosphere. These improvements to the physical content and a more well‐defined tuning strategy led to major improvements in the LMDZ6A version model climatology. Beyond the presentation of this particular model version and documentation of its climatology, the present paper underlines possible methodological pathways toward model improvement that can be shared across modeling groups. Plain Language Summary: The improvement of global numerical models is essential for the anticipation of future climate changes. We present significant advances in the physical content of a particular atmospheric model which contributes to the simulations of the Coupled Model Intercomparison project CMIP that feed reports from the IPCC. We document in particular the improvements of the representation through "parameterizations" of convective and cloudy processes. The article emphasizes the importance of strengthening the formalization of the methodology of development and tuning of models, so that new physical ideas can be translated into effective improvement of the climate representation. Key Points: The development strategy of the LMDZ6A global atmospheric circulation model is presentedImprovements with respect to previous versions are documented in the context of the Coupled Model Intercomparison Project, CMIPThe improvements are based on significant changes of the physics content as well as on a better controlled tuning strategy [ABSTRACT FROM AUTHOR]

Climate models rely on parameterizations of a variety of processes in the atmospheric physics, but a common concern is that the temporal resolution is too coarse to consistently resolve the behavior that individual parameterizations are designed to capture. This study examines timescales numerically derived from the Morrison‐Gettelman (MG2) microphysics as implemented within the Energy Exascale Earth System Model, version 1 (E3SMv1). Numerically relevant timescales in MG2 are derived by computing the eigenspectrum of its Jacobian. These timescales are found to often be smaller than the default 5 min time step used for MG2. The fast timescales are then heuristically connected to individual microphysics processes. By substepping a few particular rain processes within MG2, the time discretization error for those processes was considerably reduced with minimal additional expense to the overall microphysics. While this improvement has a substantial effect on the target processes and on the vertical distribution of stratiform‐derived rain within E3SMv1, the overall model climate is found to not be sensitive to MG2 time step. We hypothesize that this is because the surface climate does not depend strongly on certain process rates, especially MG2's rain evaporation rate. Plain Language Summary: The atmospheric components of climate models contain a number of physics parameterizations, subcomponents that are designed to capture particular aspects of the atmospheric physics. Cloud microphysics models are parameterizations designed to represent very small scale cloud processes, including phase changes and the formation of precipitation. The accuracy of these parameterizations depends on the model time step: A shorter time step typically requires more computational resources, but also improves the model's accuracy. This paper examines a particular microphysics model, MG2, used with a time step of 5 min in the E3SMv1 climate model. By linearizing MG2, we can find characteristic timescales associated with this model, which are often much shorter than 5 min. This suggests that the usual time step is too large to fully capture the physics that MG2 represents. We also experiment with using a shorter time step for parts of the rain physics and find that in many cases the rain mass is strongly affected. However, reducing the MG2 time step does not have much effect on the overall climate of E3SMv1. Key Points: MG2 contains several processes that are poorly resolved at a time step of 5 minSubstepping MG2's rain‐related processes changes their rates considerablyHowever, reducing MG2's time step within E3SM has little to no impact on the climate [ABSTRACT FROM AUTHOR]

Canonical understanding based on general circulation models (GCMs) is that the atmospheric circulation response to midlatitude sea‐surface temperature (SST) anomalies is weak compared to the larger influence of tropical SST anomalies. However, the ∼100‐km horizontal resolution of modern GCMs is too coarse to resolve strong updrafts within weather fronts, which could provide a pathway for surface anomalies to be communicated aloft. Here, we investigate the large‐scale atmospheric circulation response to idealized Gulf Stream SST anomalies in Community Atmosphere Model (CAM6) simulations with 14‐km regional grid refinement over the North Atlantic, and compare it to the responses in simulations with 28‐km regional refinement and uniform 111‐km resolution. The highest resolution simulations show a large positive response of the wintertime North Atlantic Oscillation (NAO) to positive SST anomalies in the Gulf Stream, a 0.4‐standard‐deviation anomaly in the seasonal‐mean NAO for 2°C SST anomalies. The lower‐resolution simulations show a weaker response with a different spatial structure. The enhanced large‐scale circulation response results from an increase in resolved vertical motions with resolution and an associated increase in the influence of SST anomalies on transient‐eddy heat and momentum fluxes in the free troposphere. In response to positive SST anomalies, these processes lead to a stronger and less variable North Atlantic jet, as is characteristic of positive NAO anomalies. Our results suggest that the atmosphere responds differently to midlatitude SST anomalies in higher‐resolution models and that regional refinement in key regions offers a potential pathway to improve multi‐year regional climate predictions based on midlatitude SSTs. Plain Language Summary: Variations in the ocean surface temperature (SST) influence the atmospheric circulation and thus climate over land. Canonical understanding is that tropical SSTs are more important than SSTs in midlatitudes. However, this understanding is based on climate models that don't resolve processes at scales less than 100 km. Here, we show that by increasing the atmospheric model resolution to resolve features on smaller scales, such as weather fronts, we find a larger atmospheric circulation response to midlatitude SST anomalies in the North Atlantic. North Atlantic SST anomalies can be predicted multiple years in advance, and a larger atmospheric circulation response to these predictable SST anomalies therefore implies increased predictability of climate over the surrounding land regions. Key Points: There is a large circulation response to idealized Gulf Stream sea‐surface temperature (SST) anomalies in an atmospheric model with 14‐km regional grid refinementThis response is weaker or absent in simulations with 28‐km or coarser resolution, which do not fully resolve mesoscale frontal processesTransient‐eddy fluxes of heat and momentum are modified as fronts pass over warm SSTs, leading to a large‐scale circulation response [ABSTRACT FROM AUTHOR]

Current eddy‐permitting and eddy‐resolving ocean models require dissipation to prevent a spurious accumulation of enstrophy at the grid scale. We introduce a new numerical scheme for momentum advection in large‐scale ocean models that involves upwinding through a weighted essentially non‐oscillatory (WENO) reconstruction. The new scheme provides implicit dissipation and thereby avoids the need for an additional explicit dissipation that may require calibration of unknown parameters. This approach uses the rotational, "vector invariant" formulation of the momentum advection operator that is widely employed by global general circulation models. A novel formulation of the WENO "smoothness indicators" is key for avoiding excessive numerical dissipation of kinetic energy and enstrophy at grid‐resolved scales. We test the new advection scheme against a standard approach that combines explicit dissipation with a dispersive discretization of the rotational advection operator in two scenarios: (a) two‐dimensional turbulence and (b) three‐dimensional baroclinic equilibration. In both cases, the solutions are stable, free from dispersive artifacts, and achieve increased "effective" resolution compared to other approaches commonly used in ocean models. Plain Language Summary: High‐resolution climate models that resolve the cyclones and anticyclones in the ocean, often called "eddies," must prevent an artificial build‐up of whirl‐like movements, or "enstrophy," at the model's grid‐scale. Even though methods that prevent artificial accumulation of enstrophy are included only to ensure numerical stability, they unfortunately also negatively impact the accuracy of model predictions even at scales larger than the grid‐scale. Here, we devise a novel numerical method to overcome this deficiency. Our method has the best of both worlds: it removes just enough enstrophy so that the flow is as close to reality as possible and it achieves this without accumulating enstrophy at grid‐scale. Key Points: We describe a new momentum advection scheme based on upwind‐biased, weighted essentially non‐oscillatory reconstructionsThe new scheme automatically adapts to horizontal resolution without generating grid‐scale noise typical of low order oscillatory schemesWe demonstrate a higher "effective" resolution compared to other widely used approaches [ABSTRACT FROM AUTHOR]

The non‐local model of mixing based on internal wave breaking, IDEMIX, is implemented as an enhancement of a turbulent kinetic energy closure model in three non‐eddy resolving general circulation ocean models that differ in the discretization and choice of computational grids. In IDEMIX internal wave energy is generated by an energy flux resulting from near‐inertial waves induced by wind forcing at the surface, and at the bottom, by an energy flux that parameterizes the transfer of energy between baroclinic and barotropic tides. In all model simulations with IDEMIX, the mixing work is increased compared to the reference solutions without IDEMIX, reaching values in better agreement with finestructure observations. Furthermore, the horizontal structure of the mixing work is more realistic as a consequence of the heterogeneous forcing functions. All models with IDEMIX simulate deeper thermocline depths related to stronger shallow overturning cells in the Indo‐Pacific. In the North Atlantic, deeper mixed layers in simulations with IDEMIX are associated with an increased Atlantic overturning circulation and an increase of northward heat transports toward more realistic values. The response of the deep Indo‐Pacific overturning circulation and the weak bottom cell of the Atlantic to the inclusion of IDEMIX is incoherent between the models, suggesting that additional unidentified processes and numerical mixing may confound the analysis. Applying different tidal forcing functions leads to simulation differences that are small compared to differences between the different models or between simulations with IDEMIX and without IDEMIX. Plain Language Summary: Waves in the ocean interior play a fundamental role for ocean dynamics since they can carry energy over long distances and, once they break, lead to turbulent mixing. This turbulent mixing can cause dense water masses to rise from the deep ocean with a direct impact on large‐scale currents. The wave dynamics occur on spatial scales that cannot be resolved in global ocean or climate models. To account for these processes, we apply the new parameterization IDEMIX that describes internal wave generation, propagation, and mixing. Using three different ocean models with and without IDEMIX ensures that we can identify model‐specific effects of the parameterization and discriminate them from those independent of the model. We find that the simulated mixing patterns agree better with observations once IDEMIX is applied. Large‐scale currents and the vertical temperature distribution are substantially affected by the internal wave parameterization. Whether this leads to an improved agreement with observed currents and water mass properties depends on the specific model and on numerical effects. In most cases, simulations with IDEMIX are not very sensitive to details of how the internal wave model is driven by tidal energy input. Key Points: the IDEMIX closure for a consistent representation of internal wave‐induced mixing is evaluated in three state‐of‐the‐art ocean modelsonly in simulations with IDEMIX can the models reproduce the magnitude and spatial variability of the observed mixing workmost changes with IDEMIX can be attributed to stronger mixing, but some effects are confounded by other processes and numerical mixing [ABSTRACT FROM AUTHOR]

The internal variability pertains to fluctuations originating from processes inherent to the climate component and their mutual interactions. On the other hand, forced variability delineates the influence of external boundary conditions on the physical climate system. A methodology is formulated to distinguish between internal and forced variability within the surface air temperature. The noise‐to‐noise approach is employed for training a neural network, drawing an analogy between internal variability and image noise. A large training data set is compiled using surface air temperature data spanning from 1901 to 2020, obtained from an ensemble of Atmosphere‐Ocean General Circulation Model simulations. The neural network utilized for training is a U‐Net, a widely adopted convolutional network primarily designed for image segmentation. To assess performance, comparisons are made between outputs from two single‐model initial‐condition large ensembles, the ensemble mean, and the U‐Net's predictions. The U‐Net reduces internal variability by a factor of four, although notable discrepancies are observed at the regional scale. While demonstrating effective filtering of the El Niño Southern Oscillation, the U‐Net encounters challenges in capturing the changes in the North Atlantic Ocean. This methodology holds potential for extension to other physical variables, facilitating insights into the climate change triggered by external forcings over the long term. Plain Language Summary: To anticipate future climate change, it is crucial to detect and understand the impacts of human activities. However, distinguishing the effects of anthropogenic forcing from natural climate variations in observational data is challenging. Natural climate variability, known as internal variability, arises from the chaotic nature of atmospheric and oceanic circulation, and from the interactions among the ocean, atmosphere, and land. Here, a novel approach is introduced to distinguish the changes caused by human activities from internal variability. It is applied to the surface air temperature evolution from 1901 to 2020. This method uses an artificial neural network designed to separate the internal from the human‐induced variability. An unprecedented number of climate model simulations are used, enabling precise estimation of human‐forced variability in these climate models. The spatio‐temporal variations are distinguished by applying a well‐known methodology previously used to remove noise from images. The method's performance is evaluated, revealing errors regarding the internal variability that are typically one‐fourth of the actual variations. Regions with important internal variability or with low agreement among models exhibit the largest errors. Overall, the skills are comparable to other existing approaches, but improvements are anticipated. Key Points: We present a new method to separate the forced and internal variability of the surface air temperatureWe utilize a U‐Net trained with global climate models outputs and implement a noise to noise methodology to eliminate internal variabilityThe results are assessed through the utilization of very large ensemble simulations of two distinct climate models [ABSTRACT FROM AUTHOR]

A single column model with parameterized large‐scale (LS) dynamics is used to better understand the response of steady‐state tropical precipitation to relative sea surface temperature under various representations of radiation, convection, and circulation. The large‐scale dynamics are parametrized via the weak temperature gradient (WTG), damped gravity wave (DGW), and spectral weak temperature gradient (Spectral WTG) method in NCAR's Single Column Atmosphere Model (SCAM6). Radiative cooling is either specified or interactive, and the convective parameterization is run using two different values of a parameter that controls the degree of convective inhibition. Results are interpreted in the context of the Global Atmospheric System Studies ‐Weak Temperature Gradient (GASS‐WTG) Intercomparison project. Using the same parameter settings and simulation configuration as in the GASS‐WTG Intercomparison project, SCAM6 under the WTG and DGW methods produces erratic results, suggestive of numerical instability. However, when key parameters are changed to weaken the large‐scale circulation's damping of tropospheric temperature variations, SCAM6 performs comparably to single column models in the GASS‐WTG Intercomparison project. The Spectral WTG method is less sensitive to changes in convection and radiation than are the other two methods, performing qualitatively similarly across all configurations considered. Under all three methods, circulation strength, represented in 1D by grid‐scale vertical velocity, is decreased when barriers to convection are reduced. This effect is most extreme under specified radiative cooling, and is shown to come from increased static stability in the column's reference radiative‐convective equilibrium profile. This argument can be extended to interactive radiation cases as well, though perhaps less conclusively. Plain Language Summary: Single column models, as the name suggests, only model the vertical dimension of the atmosphere. They are simpler than full‐scale 3D global circulation models, but nonetheless play an important role in model development and in better understanding physical phenomena. We use NCAR's Single Column Atmosphere Model (SCAM6) to better understand tropical rainfall. In a single column model, the atmospheric wind coming from other locations (the large‐scale circulation) must be either specified or approximated using a parameterization. We implement three different parameterizations of the large‐scale circulation into SCAM6 and assess how SCAM6 responds to various changes while using these parameterizations. The Global Atmospheric System Studies‐Weak Temperature Gradient (GASS‐WTG) Intercomparison project assessed the performance of various other single column models using two parameterizations of the large‐scale circulation. We fold SCAM6 into this Intercomparison, and we find that SCAM6 performs substantially differently to its peers under the settings used in the GASS‐WTG Intercomparison project, but more comparably to its peers when key circulation parameters are relaxed. Another notable finding is that, across all three parameterizations of large‐scale dynamics, circulation strength decreases when we reduce SCAM6's barrier to convection. Key Points: We use three different methods to parameterize the large‐scale (LS) dynamics in NCAR's single column atmospheric model (SCAM6)As in the Global Atmospheric System Studies Intercomparison, we test SCAM6's response to various boundary conditions and model parametersUnder all three methods, circulation strength is decreased when barriers to convection are reduced [ABSTRACT FROM AUTHOR]

Many general circulation models (GCMs) assume some heterogeneity of water amounts in their grid boxes and use probability density functions to parameterize cloud fractions CF and amounts of condensed water qc. Most GCM cloud schemes calculate the CF as the volume of the grid box that contains clouds (CFvol), whereas radiative fluxes primarily depend on the CF by surface (CFsurf), that is, the surface of the grid box covered by clouds when looking from above. This discrepancy matters as previous findings suggest that CFsurf is typically greater than CFvol by about 30%. In this paper we modify the single column model version of the LMDz GCM cloud scheme by introducing the vertical subgrid‐scale heterogeneity of water content. This allows to distinctly compute the two fractions, CFvol and CFsurf, as well as the amount of condensed water qc. This study is one of the first to take into account such vertical subgrid‐scale heterogeneity in a GCM cloud scheme. Three large eddy simulation cases of cumuliform boundary layer clouds are used to test and calibrate two different parameterizations. These new developments increase cloud cover by about 10% for the oceanic cases RICO and Barbados Oceanographic Meteorological Experiment and by up to 50% for the continental case ARM. The change in condensed water reduces the liquid water path by 10–20% and therefore the cloud opacity by 5–50%. These results show the potential of the new scheme to reduce the too few, too bright bias by increasing low‐level CF and decreasing cloud reflectance. Key Points: Most cloud schemes assume vertically homogeneous cloud fractions in GCM grid boxesThis approximation may induce an underestimation of cloud cover and overestimation of reflectanceWe parameterize 3‐D subgrid‐scale heterogeneities for low‐level clouds and reduce these biases [ABSTRACT FROM AUTHOR]

Abstract: NASA's Goddard Earth Observing System (GEOS) Earth System Model (ESM) is a modular, general circulation model (GCM), and data assimilation system (DAS) that is used to simulate and study the coupled dynamics, physics, chemistry, and biology of our planet. GEOS is developed by the Global Modeling and Assimilation Office (GMAO) at NASA Goddard Space Flight Center. It generates near‐real‐time analyzed data products, reanalyses, and weather and seasonal forecasts to support research targeted to understanding interactions among Earth System processes. For chemistry, our efforts are focused on ozone and its influence on the state of the atmosphere and oceans, and on trace gas data assimilation and global forecasting at mesoscale discretization. Several chemistry and aerosol modules are coupled to the GCM, which enables GEOS to address topics pertinent to NASA's Earth Science Mission. This paper describes the atmospheric chemistry components of GEOS and provides an overview of its Earth System Modeling Framework (ESMF)‐based software infrastructure, which promotes a rich spectrum of feedbacks that influence circulation and climate, and impact human and ecosystem health. We detail how GEOS allows model users to select chemical mechanisms and emission scenarios at run time, establish the extent to which the aerosol and chemical components communicate, and decide whether either or both influence the radiative transfer calculations. A variety of resolutions facilitates research on spatial and temporal scales relevant to problems ranging from hourly changes in air quality to trace gas trends in a changing climate. Samples of recent GEOS chemistry applications are provided. [ABSTRACT FROM AUTHOR]

Tropical high cloud cover decreases with surface warming in most general circulation models. This reduction, according to the "stability‐iris" hypothesis, is thermodynamically controlled and linked to a decrease in the radiatively‐driven clear‐sky convergence, when the peak anvil clouds rise because of the rising isotherms. The influence of the large‐scale dynamical changes on the tropical high cloud fraction remains difficult to disentangle from the local thermodynamic influence, given that the mean meridional circulation remains inextricably tied to the local thermodynamic structure of the atmosphere. However, using idealized general circulation model simulations, we propose a novel method to segregate the dynamical impact from the thermodynamic impact on the tropical high cloud fraction. To this end, our investigation primarily focuses on the mechanisms underpinning changes in the high cloud cover in the deep tropics in response to extratropical surface warming, when the tropical sea surface temperatures remain invariant. Net convective detrainment of ice cloud condensates decreases at the peak detrainment region, without a rise in its altitude. We also find that the importance of depositional growth of ice cloud condensates in controlling the high cloud fraction response in the deep tropics varies with altitude. Plain Language Summary: The cloud feedback associated with changes in the tropical high cloud cover is one of the major uncertainties in calculating the current estimates of climate sensitivity, which is a measure of how much the Earth's average surface temperature would increase if we double the amount of atmospheric carbon dioxide. When the surface becomes warmer in the tropics, the tropical high cloud cover decreases. However, this raises an important question: how do circulation changes independent of the temperature changes within the tropics impact the tropical high cloud cover? Using idealized general circulation model simulations, we found that the tropical high cloud fraction decreases as a result of circulation changes induced by extratropical warming, even when the tropical sea surface temperatures are held constant. Both convective and cloud microphysical processes play significant roles in controlling the changes in tropical high cloud fraction, and their relative importance varies with altitude. Understanding the different factors responsible for the changes in high cloud cover is important, as the area covered by these tropical high clouds can affect how much the Earth warms under climate change. Key Points: An indirect weakening of the Hadley circulation decreases the convectively‐detrained ice cloud condensates in the deep tropicsA concurrent reduction in the net vertical transport of water vapor limits the depositional growth of ice cloud condensatesThe relative influence of depositional growth to net convective detrainment on the tropical high cloud response is altitude‐dependent [ABSTRACT FROM AUTHOR]

To improve cloud simulation in a general circulation model, we develop a new cloud macrophysics scheme that treats detrained cumulus generated by convective detrainment process separately from pure stratus; prognoses two symmetric triangular probability density functions of total specific humidity for each detrained cumulus and pure stratus; and diagnoses both cloud fraction and cloud condensate in all liquid, ice, and mixed‐phases in a consistent and unified way without any adjustment to remove empty or very dense cloud. Supersaturation is allowed within ice cloud. The new scheme ("NEW") is compared with the previous model ("OLD") using single‐column simulations for subtropical marine stratocumulus (DYCOMS2) and continental deep convection (ARM97) cases. In DYCOMS2, both NEW and OLD produce vertical profiles of grid‐mean cloud condensate similar to large‐eddy simulation. In ARM97, compared with OLD, NEW simulates less sporadic vertical profiles of in‐cloud condensates, due to consistent diagnosis of cloud fraction and cloud condensate; more continuously‐varying detrained cumulus with time, due to prognostic treatment of detrained cumulus; and ice cloud fraction and ice condensate similar to those of OLD, in spite of completely different treatment of ice cloud processes. The global performance of NEW is similar to OLD with improved relative humidity. Compared to OLD, NEW simulates more and improved cloud condensate, but less and degraded cloud fraction, particularly, in the lower troposphere. Detrained cumulus is moister and colder and has a larger moisture variance than pure stratus. Overall, NEW simulates stronger condensation‐deposition rates than OLD, due in part to the separate treatment of detrained cumulus and pure stratus. Plain Language Summary: This study suggests a new cloud parameterization scheme in a general circulation model that can represent clouds in nature in a realistic way. Our new cloud scheme separates two types of stratus, called detrained cumulus and pure stratus, by using a set of prognostic equations. Moreover, both liquid, ice, and mixed‐phase clouds are treated in a consistent way. We found that the new scheme successfully operates for two single‐column simulations for subtropical marine stratocumulus and continental deep convection cases. In addition, we performed the global simulation and found that the new scheme simulated more cloud condensate but less cloud fraction than the old scheme with some improvements. Key Points: We develop a new cloud macrophysics scheme by prognosing two triangular probability density functions for each detrained cumulus and pure stratusThe new scheme computes both cloud fraction and cloud condensate in all liquid, ice, and mixed‐phases in a consistent and unified wayThe new scheme simulates more cloud condensate but less cloud fraction than the old scheme with improved relative humidity and global performance [ABSTRACT FROM AUTHOR]

In recent years, efforts have been made to include tides in both operational ocean models as well as climate and earth system models. The accuracy of the barotropic tides is often limited by the model topography, which is in turn limited by model horizontal resolution. In this work, we explore the reduction of barotropic tidal errors in an ocean general circulation model (Modular Ocean Model version 6; MOM6) using sub‐grid scale topography representation. We follow the methodology from Adcroft (2013, https://doi.org/10.1016/j.ocemod.2013.03.002), which utilizes statistics from finer resolution topographic data sets to represent sub‐grid scale features with a light computational cost in a structured finite volume formulation. The geometric effect from sub‐grid scale topography can be introduced to the model with only a few parameters at each grid cell. The porous barriers, which are implemented at the walls of the grid cells, are used to modify transport between grid cells. Our results show that the globally averaged tidal error in lower‐resolution simulations is significantly reduced with the use of porous barriers. We argue this method is a potentially useful tool to improve simulations of tides (and other flows) in low‐resolution simulations. Plain Language Summary: In recent years, significant attention has been given to simulating tides in diverse ocean and climate models, because of the roles of tides in impacting sea levels. In order to better predict future sea level changes, it is therefore essential to improve the accuracy of surface tides simulated in numerical ocean models. One known source of errors that could impact the accuracy of modeled tides stems from the accuracy of ocean topography in the models, which relies on how many grid boxes (i.e., model resolution) are used to represent the complicated seafloor elevations. The model's resolution, however, is often constrained by the efficiency of computations. Our objective in this study is to enhance the portrayal of ocean topography at a given resolution to minimize tidal errors. To achieve this, we employ a technique known as "porous barriers," which uses approximations to mimic the effects of ocean topography that would be omitted in the models, without imposing significant computational burdens. Our findings demonstrate that this technique of porous barriers can substantially reduce tidal errors in our ocean model. We also expect a wider application of this technique to numerous other physical processes within the ocean. Key Points: The accuracy of tides in ocean models is affected by the representation of topography, which is constrained by model horizontal resolutionThe geometric effect of the sub‐grid scale topographic features can be approximated in the model with the porous barriersThe porous barriers can significantly reduce globally averaged errors for barotropic tides [ABSTRACT FROM AUTHOR]

Gravity wave parameterizations contribute to uncertainties in middle atmosphere modeling. To investigate the potential for using machine learning to represent atmospheric gravity waves and the impact of implementing such schemes in a general circulation model (GCM), we train a random forest (RF) emulator on outputs from an existing complex parameterization scheme for orographic gravity wave drag (GWD). The performance of the RF emulator is then evaluated with a focus on stratospheric climatology and variability in climate simulations from the middle atmosphere resolving Beijing Climate Center Atmospheric Circulation Model. In offline tests, the predicted orographic GWD by the RF agrees remarkably well with the target GWD throughout the troposphere and the middle atmosphere. The RF emulator can reproduce the observed climatology of zonal‐mean wind and air temperature in the GCM simulation, as well as its target scheme. Compared to the target orographic GWD parameterization scheme, the RF emulator can reproduce the breakdown of the polar vortex in the Southern Hemisphere stratosphere. This study demonstrates the feasibility of using machine learning to emulate parameterized orographic GWD for modeling the stratosphere with a GCM. Plain Language Summary: Machine learning has been utilized to learn atmospheric physical parameterizations such as moist convection, cloud microphysics, radiation, nonorographic gravity waves, etc. However, it remains unknown whether it is feasible to apply machine learning algorithms to parameterize orographic gravity waves and how such schemes perform when interactively coupled in atmospheric general circulation models. Here we employ the widely used random forest algorithm to learn from an existing complex parameterization scheme for orographic gravity wave drag. Then, we implement the random forest scheme in a high‐top atmospheric general circulation model to test its performance. We find that the random forest scheme acts as well as the target parameterization scheme. This study demonstrates the potential of using machine learning to emulate parameterized orographic gravity waves for modeling the middle atmosphere. In future work, it is recommended that machine learning may be used to develop new parameterization schemes directly from observations or high‐resolution model outputs to overcome critical deficiencies in current orographic gravity wave parameterizations, helping to reduce the common biases in simulating the stratosphere with general circulation models. Key Points: Parameterized gravity wave drag from unresolved orography can be emulated with a random forestThe random forest emulator simulates vertical distributions of the zonal wind and air temperature quite well in a general circulation model (GCM)The GCM configured with the random forest scheme reproduces a reasonable breakdown of the austral polar vortex in the stratosphere [ABSTRACT FROM AUTHOR]

To improve the simulation of nocturnal precipitation, we develop a parameterization for the elevated convection that is launched from the level of maximum grid‐mean moist static energy, when grid‐mean vertical flow is upward at the launching interface. The parameterized elevated convection is forced by both cold pool‐driven and partially resolved external mesoscale organized flows. Properties of the external mesoscale flow are estimated from grid‐mean vertical velocity and three dimensional advection tendencies of temperature and moisture. The new parameterization is implemented into a unified convection scheme (UNICON) and tested for both the single‐column case at the Southern Great Plain (SGP) site in US and in global simulations. At the SGP site, the parameterized elevated convection strengthens nocturnal convection, increases nocturnal precipitation, and better simulates the observed diurnal cycle of precipitation. It appears that the elevated and surface‐based convections interact with each other by stabilizing the atmospheric column, which affects subsequent convection. Global simulation shows that the elevated convection mostly occurs over the continents during the night, and also over the oceanic mid‐latitude warm air advection and storm track regions during summer. Without degrading global mean climate, the elevated convection improves the simulation of nocturnal precipitation in the mid‐US but simulates somewhat strong nocturnal precipitation in northeastern Asia. It seems that a key to simulate observed nocturnal precipitation is to appropriately parameterize the impacts of external organized flow on the elevated convection. Plain Language Summary: Precipitation has a strong diurnal cycle as a result of complex interactions between the atmosphere and the surface through moist convection. The diurnal cycle of precipitation has a significant impact on surface energy and water budgets, and its accumulated effect may appear even on the climate timescale. However, most general circulation models (GCM) poorly simulate the observed diurnal cycle of precipitation. Particularly, it has been difficult to simulate nocturnal precipitation over the central US during summer. To improve the simulation of nocturnal precipitation, we develop a new parameterization for the elevated convection. Both single‐column and global test simulations show that the new parameterization increases nocturnal precipitation and better simulates the observed diurnal cycle of precipitation. Our research will contribute to improving the performance of GCM and so better simulating current and future climates. Key Points: To improve the simulation of nocturnal precipitation, we develop a parameterization for the elevated convectionIn the mid‐US, the new parameterization increases nocturnal precipitation and better simulates the observed diurnal cycle of precipitationA key to simulate observed nocturnal precipitation is to parameterize the impacts of external mesoscale flow on the elevated convection [ABSTRACT FROM AUTHOR]

The accuracy of Earth system models is compromised by unknown and/or unresolved dynamics, making the quantification of systematic model errors essential. While a model parameter estimation, which allows parameters to change spatio‐temporally, shows promise in quantifying and mitigating systematic model errors, the estimation of the spatio‐temporally distributed model parameters has been practically challenging. Here we present an efficient and practical method to estimate time‐varying parameters in high‐dimensional spaces. In our proposed method, Hybrid Offline and Online Parameter Estimation with ensemble Kalman filtering (HOOPE‐EnKF), model parameters estimated by EnKF are constrained by results of offline batch optimization, in which the posterior distribution of model parameters is obtained by comparing simulated and observed climatological variables. HOOPE‐EnKF outperforms the original EnKF in synthetic experiments using a two‐scale Lorenz96 model and a simple global general circulation model. One advantage of HOOPE‐EnKF over traditional EnKFs is that its performance is not greatly affected by inflation factors for model parameters, thus eliminating the need for extensive tuning of inflation factors. We thoroughly discuss the potential of HOOPE‐EnKF as a practical method for improving parameterizations of process‐based models and prediction in real‐world applications such as numerical weather prediction. Plain Language Summary: Earth system models help us understand and predict the behavior of our planet, but their accuracy is limited due to unknown or unresolved factors. Adjusting the parameters of these models based on the changing patterns in time and space can help improve their accuracy. However, this has been a challenging task. In this study, we introduce a new method called Hybrid Offline and Online Parameter Estimation with ensemble Kalman filtering (HOOPE‐EnKF) that efficiently estimates these changing parameters. This method combines an existing EnKF with offline batch optimization, which compares long‐term simulation with observations to fine‐tune the model. Our experiments show that HOOPE‐EnKF performs better than the original EnKF and is less sensitive to certain hyperparameters, reducing the need for time‐consuming calibrations. HOOPE‐EnKF has a great potential for improving the accuracy of models and predictions in real‐world applications, like weather forecasting, and can help scientists better understand our planet's complex systems. Key Points: Ensemble Kalman filter (EnKF) is extended to estimate high‐dimensional time‐varying model parametersThe proposed method is insensitive to the choice of covariance inflation factorsThe proposed method successfully quantifies systematic model errors from observations [ABSTRACT FROM AUTHOR]

This paper describes the initial implementation of a new toolbox that seeks to balance accuracy, efficiency, and flexibility in radiation calculations for dynamical models. The toolbox consists of two related code bases: Radiative Transfer for Energetics (RTE), which computes fluxes given a radiative transfer problem defined in terms of optical properties, boundary conditions, and source functions; and RRTM for General circulation model applications--Parallel (RRTMGP), which combines data and algorithms to map a physical description of the gaseous atmosphere into such a radiative transfer problem. The toolbox is an implementation of well-established ideas, including the use of a k-distribution to represent the spectral variation of absorption by gases and the use of two-stream, plane-parallel methods for solving the radiative transfer equation. The focus is instead on accuracy, by basing the k-distribution on state-of-the-art spectroscopy and on the sometimes-conflicting goals of flexibility and efficiency. Flexibility is facilitated by making extensive use of computational objects encompassing code and data, the latter provisioned at runtime and potentially tailored to specific problems. The computational objects provide robust access to a set of high-efficiency computational kernels that can be adapted to new computational environments. Accuracy is obtained by careful choice of algorithms and through tuning and validation of the k-distribution against benchmark calculations. Flexibility with respect to the host model implies user responsibility for maps between clouds and aerosols and the radiative transfer problem, although comprehensive examples are provided for clouds. [ABSTRACT FROM AUTHOR]

Modeling deep convection accurately in tropical regions is important. However, biases remain in current trigger functions. To alleviate the overestimation of frequency and wrong depiction of the diurnal cycle, we propose a deep convection trigger function, ResU‐Deep, based on the framework of U‐net with three modifications to better suit the problem of deep convection identification: (a) adding the upsampling process into the encoder part, (b) replacing the double convolution block with a residual‐convolutional block, and (c) adding a dynamic weight into the loss function. Thirty‐three environmental variables within tropical regions are used in ResU‐Deep, including 31 features from ECMWF atmospheric reanalysis (ERA5) data set, and two historical convection fields. Tropical Rainfall Measuring Mission 3B42 data set is used as the precipitation observation. Central America, North Africa, South and East Asia, and West Pacific Ocean within 0°∼30°N are selected as the study regions for the high frequency of deep convection activities. ResU‐Deep, incorporating the surrounding information, is separately trained and evaluated in four regions and has the F1‐scores of 58%, 53%, 60%, and 63% for the occurrence, outperforming the single‐column‐based machine learning methods. Also, a unified model has similar performance in four regions. Further comparisons are made with convective available potential energy‐based trigger functions in Southern Great Plains. Results show that ResU‐Deep can capture the trends and peaks of diurnal cycles on complex terrains in large regions. According to feature importance test, the contribution levels of environmental features are different in four regions, indicating the model can learn the mechanisms of deep convection in specific region, thus improving the prediction accuracy. Plain Language Summary: The deep convection trigger function is a widely‐used method to determine the occurrence of deep convection in general circulation models. However, many trigger functions face the overestimation of deep convections and cannot accurately capture the variations of occurrence frequency. We propose a convection trigger function based on deep learning, called ResU‐Deep, to improve the trigger function in tropical regions. Three modifications are made based on the framework of U‐net. Both the architecture and the loss function are modified to improve the model's performance and increase the training speed. Thirty‐one environmental variables from ERA5 data set and two historical deep convection occurrence information are used as the features. Central America, North Africa, South and East Asia, and West Pacific Ocean within the tropics have high deep convection frequency. The model is trained and tested independently on each site, reflecting the non‐uniform physical mechanisms for deep convection. Also, a unified model has similar performance. Besides, incorporating surrounding information in ResU‐Deep can help improve prediction accuracy. To compare with the convective available potential energy‐based trigger functions, we test the model on Southern Great Plains in central United States. Results show that deep learning is a promising solution to identify deep convection. Key Points: A location‐aware and deep‐learning‐based deep convection trigger function is proposed to improve the diurnal cycle simulation in tropicsResults show that terrain information, temporal dependence of convection, and water vapor content are essential for predicting convectionIncluding information from neighboring atmospheric columns can improve the performance of the deep convection trigger function [ABSTRACT FROM AUTHOR]

Deep learning (DL) has emerged as a promising tool to downscale climate projections at regional‐to‐local scales from large‐scale atmospheric fields following the perfect‐prognosis approach. Given their complexity, it is crucial to properly evaluate these methods, especially when applied to changing climatic conditions where the ability to extrapolate/generalize is key. In this work, we intercompare several DL models extracted from the literature for the same challenging use‐case (downscaling temperature in the CORDEX North America domain) and expand standard evaluation methods building on eXplainable Artificial Intelligence (XAI) techniques. Specifically, we introduce two novel XAI‐based diagnostics—Aggregated Saliency Map and Saliency Dispersion Maps—and show how they can be used to unravel the internal behavior of these models, aiding in their design and evaluation. This work advocates for the introduction of XAI techniques into deep downscaling evaluation frameworks, especially when working with large regions and/or under climate change conditions. Plain Language Summary: Due to limitations in the computational resources available, General Circulation Models (GCMs) are often used to simulate the climate system over coarse resolution grids. This hampers the applicability of GCM products in the regional‐to‐local scale, highly demanded by different socio‐economic sectors. Statistical downscaling aims to solve this problem by generating high‐resolution climate fields. Recently, machine learning techniques—particularly deep learning (DL) models—have shown promising results in this task. These models are first trained in a historical period through observational data sets, and then applied to the GCM outputs of plausible future scenarios, thus generating high‐resolution climate change products. To assess the performance of these methods, a number of evaluation metrics have been proposed considering both the skill to reproduce present climate conditions and the ability to generalize changing conditions. Here, we illustrate the possibilities of eXplainable Artificial Intelligence (XAI) techniques to expand the evaluation framework for deep downscaling methods, introducing new XAI‐derived diagnostics to unravel their internal behavior. The results show the usefulness of incorporating XAI techniques into statistical downscaling evaluation frameworks, especially when working with large regions and/or under climate change conditions. Key Points: Explainable artificial intelligence (XAI) facilitates the evaluation of deep downscaling models by unraveling their internal behaviorXAI techniques can detect structural problems not revealed by standard evaluation, which may be relevant for understanding model differencesXAI techniques can assess the relevance and locality of the predictors of deep learning models, helping to assess their physical consistency [ABSTRACT FROM AUTHOR]

The ECHAM6 atmospheric general circulation model is the atmosphere component of the Max Planck Institute Earth System Model (MPI-ESM) that is used in the Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations. As ECHAM6 has its uppermost layer centered at 0.01 hPa in the upper mesosphere, these simulations offer the opportunity to study the middle atmosphere climate change and its relation to the troposphere on the basis of a very comprehensive set of state-of-the-art model simulations. The goals of this paper are (a) to introduce those new features of ECHAM6 particularly relevant for the middle atmosphere, including external forcing data, and (b) to evaluate the simulated middle atmosphere and describe the simulated response to natural and anthropogenic forcings. New features in ECHAM6 with respect to ECHAM5 include a new short-wave radiation scheme, the option to vary spectral irradiance independent of total solar irradiance, and a latitude-dependent gravity-wave source strength. The description of external forcing data focuses on solar irradiance and ozone. Stratospheric temperature trends simulated with the MPI-ESM for the last decades of the 20th century agree well with observations. The future projections depend strongly on the scenario. Under the high emission scenario RCP8.5, simulated temperatures are locally lower by more than 20 K than preindustrial values. Many of the simulated patterns of the responses to natural forcings as provided by solar variability, volcanic aerosols, and El Niño-Southern Oscillation, largely agree with the observations. [ABSTRACT FROM AUTHOR]

An investigation of the impact of the subgrid-scale variability of cloud liquid water on the autoconversion process as parameterized in a general circulation model is presented in this paper. For this purpose, a prognostic statistical probability density distribution (PDF) of the subgrid scale variability of cloud water is incorporated in a continuous autoconversion parameterization. Thus, the revised autoconversion rate is calculated by an integral of the autoconversion equation over the PDF of total water mixing ratio from the saturation vapor mixing ratio to the maximum of total water mixing ratio. An evaluation of the new autoconversion parameterization is carried out by means of one year simulations with the ECHAM5 climate model. The results indicate that the new autoconversion scheme causes an increase of the frequency of occurrence of high autoconversion rates and a decrease of low ones compared to the original scheme. This expected result is due to the emphasis on areas of high cloud liquid water in the new approach, and the non-linearity of the autoconversion with respect to liquid water mixing ratio. A similar trend as in the autoconversion is observed in the accretion process resulting from the coupling of both processes. As a consequence of the altered autoconversion, large-scale surface precipitation also shows a shift of occurrence from lower to higher rates. The vertically integrated cloud liquid water estimated by the model shows slight improvements compared to satellite data. Most importantly, the artificial tuning factor for autoconversion in the continuous parameterization could be reduced by almost an order of magnitude using the revised parameterization. [ABSTRACT FROM AUTHOR]

Soil water‐groundwater interactions are important in determining the moisture profile in a soil column. However, the Noah Land Surface Model (Noah LSM), the land surface component of several general circulation models, does not consider groundwater effects. The present study investigates the impact of integrating a groundwater model into the Noah LSM. By performing experiments forced by observations at the surface, insights into the characteristics of soil moisture evolution with and without the aquifer are investigated. In the absence of surface precipitation, soil moisture is shown to exhibit an exponential decay rate without the aquifer (termed as Noah‐Cntl). This is because evapotranspiration dominates other processes governing the soil moisture evolution. Based on scale analysis, we derived an analytical equation that could well represent the aforementioned exponential variation in soil moisture. Presence of the aquifer (termed as Noah‐GW) makes the capillary processes important as well. The result is a slower decay of soil moisture as compared to Noah‐Cntl without surface precipitation. Soil moisture evolution in both the models is a function of vegetation and soil types. As a result, the aforementioned decay timescale displays significant heterogeneity over the Indian region, with Central India exhibiting relatively faster decay rates. We also performed 24‐year long experiments with observed interannually varying forcing that indicates an enhancement in the annual mean soil moisture at all levels due to groundwater effects. This enhancement is particularly prominent during the post‐monsoon season, reaffirming the results of the drain‐out experiments mentioned above. Plain Language Summary: Land surface processes are integral parts of weather and climate. Hence, an accurate representation of these processes is crucial from prediction and understanding points of view. However, state‐of‐the‐art land surface models do not incorporate some essential components that make these models erroneous. This study integrates a groundwater model into an existing land surface model and characterizes its impact on the evolution of soil moisture. We discuss factors affecting soil moisture decay rate over the Indian subcontinent without surface forcing. We show that the post‐monsoon and annual mean soil moisture increases in the presence of a groundwater model. However, spatial heterogeneity of groundwater effects is dependent on parameters like soil type and surface rain rate. Key Points: A prognostic model for groundwater is developed and integrated inside the Noah Land Surface ModelAnalytical solutions for the rate of change of soil moisture in the absence of surface water forcing are derived to explain model resultsPresence of an aquifer increases the soil moisture retention capability in the Noah Land surface model [ABSTRACT FROM AUTHOR]

Cloud fraction (CF) significantly affects the short‐ and long‐wave radiation. Its realistic representation in general circulation models (GCMs) still poses great challenges in modeling the atmosphere. Here, we present a neural network‐based (NN‐based) diagnostic scheme that uses the grid‐mean temperature, pressure, liquid and ice water mixing ratios, and relative humidity to simulate the sub‐grid CF. The scheme, trained using CloudSat data with explicit consideration of grid sizes, realistically simulates the observed CF with a correlation coefficient >0.9 for liquid‐, mixed‐, and ice‐phase clouds. The scheme also captures the observed non‐monotonic relationship between CF and relative humidity and is computationally efficient, and robust for GCMs with a variety of horizontal and vertical resolutions. For illustrative purposes, we conducted comparative analyses of the 2006–2019 climatological‐mean cloud fractions among CloudSat, and simulations from the NN‐based scheme and the Xu‐Randall scheme (optimized the same way as the NN‐based scheme). The NN‐based scheme improves not only the spatial distribution of the total CF but also the cloud vertical structure. For example, the biases of too‐many high‐level clouds over the tropics and too‐many low‐level clouds over regions around 60°S and 60°N in the Xu‐Randall scheme are significantly reduced. These improvements are also found to be insensitive to the spatio‐temporal variability of large‐scale meteorology conditions, implying that the scheme can be used in different climate regimes. Plain Language Summary: Cloud fraction (CF) is crucial to cloud climate effects of both reflecting shortwave radiation and absorbing/emitting longwave radiation. However, the simulation of CF in general circulation models (GCMs) has been difficult, because most clouds are smaller than the typical scales of GCM grids and cannot be resolved by the grid‐scale physics, while the physical understanding of sub‐grid processes is still inadequate. Thus, this study uses a data‐driven approach, that is, a neural network, to parameterize the sub‐grid CF in climate models. The database for training and evaluating this NN‐based scheme is obtained by upscaling the CloudSat (quasi‐) observational data and emulating the GCM "grid‐mean" properties required for cloud‐fraction parameterization, minimizing the data‐oriented biases. Moreover, the effects of the GCM horizontal and vertical grid sizes are both considered in the network, increasing the scheme adaptivity for use in GCMs with different resolutions. Results show that the NN‐based scheme correctly predicts observed features of cloud‐fraction variation with cloud condensate content and relative humidity for clouds of different phases and better predicts total cloud‐fraction spatial distribution and cloud vertical structure than the conventional Xu‐Randall scheme. This suggests that the NN‐based scheme has the potential to reduce the biases of cloud radiative effects existing in current GCMs. Key Points: A neural network‐based scheme for parameterizing sub‐grid cloud fraction in climate models is developed using the CloudSat dataThe scheme considers the effects of both horizontal and vertical grid sizes on cloud‐fraction parameterizationThe scheme better predicts total cloud‐fraction spatial distribution and cloud vertical structure than the Xu‐Randall scheme [ABSTRACT FROM AUTHOR]

We conduct a series of eight 45‐day experiments with a global storm‐resolving model (GSRM) to test the sensitivity of relative humidity R $\mathcal{R}$ in the tropics to changes in model resolution and parameterizations. These changes include changes in horizontal and vertical grid spacing as well as in the parameterizations of microphysics and turbulence, and are chosen to capture currently existing differences among GSRMs. To link the R $\mathcal{R}$ distribution in the tropical free troposphere with processes in the deep convective regions, we adopt a trajectory‐based assessment of the last‐saturation paradigm. The perturbations we apply to the model result in tropical mean R $\mathcal{R}$ changes ranging from 0.5% to 8% (absolute) in the mid troposphere. The generated R $\mathcal{R}$ spread is similar to that in a multi‐model ensemble of GSRMs and smaller than the spread across conventional general circulation models, supporting that an explicit representation of deep convection reduces the uncertainty in tropical R $\mathcal{R}$. The largest R $\mathcal{R}$ changes result from changes in parameterizations, suggesting that model physics represent a major source of humidity spread across GSRMs. The R $\mathcal{R}$ in the moist tropical regions is particularly sensitive to vertical mixing processes within the tropics, which impact R $\mathcal{R}$ through their effect on the last‐saturation temperature rather than their effect on the evolution of the humidity since last‐saturation. In our analysis the R $\mathcal{R}$ of the dry tropical regions strongly depends on the exchange with the extratropics. The interaction between tropics and extratropics could change with warming and presage changes in the radiatively sensitive dry regions. Plain Language Summary: Water vapor is the most important greenhouse gas in the atmosphere. Therefore, for the prediction of future warming it is important that climate models capture the distribution of atmospheric humidity and its change under warming. However, climate models currently strongly disagree in their representation of humidity, causing uncertainty in climate predictions. A recent study has shown that, while there is better agreement among the newest generation of climate models, so called global storm‐resolving models, the remaining inter‐model differences are still relevant and therefore need to be better understood. To narrow down the causes of these differences, in this study we examine how much the humidity in a storm‐resolving model changes in response to changes in different model components, which are chosen to reflect the differences that currently exist between models. We find the largest humidity changes in response to changes in the model's representation of sub‐grid scale processes. In storm‐resolving models these are turbulent motions and cloud microphysics. Our results suggest that differences in the representation of these processes cause a major part of the humidity differences between storm‐resolving models. Key Points: Sensitivity experiments suggest that parameterizations are the major source of relative humidity spread across global storm‐resolving modelsVertical mixing processes strongly impact the humidity of the moist tropics by affecting last‐saturation statistics within the tropicsCompared to the rest of the tropics, the humidity of the dry tropics is more sensitive to the pathways of exchange with the extratropics [ABSTRACT FROM AUTHOR]

Atmospheric gravity waves (GWs) span a broad range of length scales. As a result, the un‐resolved and under‐resolved GWs have to be represented using a sub‐grid scale (SGS) parameterization in general circulation models (GCMs). In recent years, machine learning (ML) techniques have emerged as novel methods for SGS modeling of climate processes. In the widely used approach of supervised (offline) learning, the true representation of the SGS terms have to be properly extracted from high‐fidelity data (e.g., GW‐resolving simulations). However, this is a non‐trivial task, and the quality of the ML‐based parameterization significantly hinges on the quality of these SGS terms. Here, we compare three methods to extract 3D GW fluxes and the resulting drag (Gravity Wave Drag [GWD]) from high‐resolution simulations: Helmholtz decomposition, and spatial filtering to compute the Reynolds stress and the full SGS stress. In addition to previous studies that focused only on vertical fluxes by GWs, we also quantify the SGS GWD due to lateral momentum fluxes. We build and utilize a library of tropical high‐resolution (Δx = 3 km) simulations using weather research and forecasting model. Results show that the SGS lateral momentum fluxes could have a significant contribution to the total GWD. Moreover, when estimating GWD due to lateral effects, interactions between the SGS and the resolved large‐scale flow need to be considered. The sensitivity of the results to different filter type and length scale (dependent on GCM resolution) is also explored to inform the scale‐awareness in the development of data‐driven parameterizations. Plain Language Summary: Gravity waves (GWs) present a challenge to climate prediction: waves on scales of O(1)–O(100) km can neither be systematically measured with conventional observational systems, nor properly represented (resolved) in operational climate models, which have a typical grid spacing on the order of 100 km. Therefore, in these climate models, small‐scale GWs must be parameterized, or estimated, based on the resolved (large‐scale) flow. The primary effects of these small‐scale waves on the resolved flow is the so‐called sub‐grid scale drag (Gravity Wave Drag [GWD]), resulting from the propagation and breaking of these waves. Existing GW parameterizations in general circulation models are all highly simplified; for example, they only account for vertical propagation of GWs. With growing computing power, a promising alternative approach is to use machine learning to develop data‐driven parameterizations. However, this requires to first generate reliable high‐resolution computer simulations and then extract GWD from these simulations. This study follows these steps, compares different extraction methods, and describes some challenges and pathways to make advances. Furthermore, our results suggest that the horizontal propagation of GWs should be included in parameterizations too, however, extra care is needed in order to extract the resulting GWD from high‐resolution data. Key Points: In a library of weather research and forecasting model simulations, we compare methods for estimating 3D gravity wave drag force that are un‐ and under‐resolved by general circulation modelsFor drag associated with vertical fluxes, different methods agree on time‐ and zonal‐mean but not on instantaneous spatiotemporal patternsDrag associated with horizontal fluxes is significant but is very sensitive to the estimation methodology [ABSTRACT FROM AUTHOR]

The ocean circulation dynamics can be represented as a high‐dimensional multi‐scale nonlinear system with inhomogenous meanfields and topography. Numerical simulations of the ocean dynamics resolve the large scales of motion on a computational grid, with the unresolved subgrid interactions parameterized. These simulations are highly dependent upon the grid resolution, unless the subgrid terms are appropriately parameterized. There are five fundamental classes of subgrid interactions: eddy‐eddy; eddy‐topographic; eddy‐meanfield; meanfield‐meanfield; and meanfield‐topographic. Scale dependent parameterizations representing each of these interaction classes are presented here in oceanic flows for the first time. Subgrid parameterizations are calculated and validated in baroclinic quasi‐geostrophic simulations of idealized Antarctic Circumpolar Current flows with representative mean currents and ocean floor topography. The parameterization coefficients are derived from the coarse grained statistics of high resolution reference simulations in spectral space. Stochastic and deterministic parameterizations are developed for the eddy‐eddy interactions, and deterministic forms for the remaining classes. The kinetic energy spectra and meanfield resulting from large eddy simulations (LES) adopting these coefficients accurately replicate those of the reference simulation. The eddy‐eddy interactions are dominant, but all classes need to be parameterized for best results. For LES where baroclinic instability is explicitly resolved the stochastic variants out‐perform the deterministic ones across all scales. When baroclinic instability is not explicitly resolved, the stochastic variants out‐perform the deterministic cases at the large scales, but introduce some distortions at the smallest resolved scales. This study provides an assessment of the relevant strengths of the subgrid interaction classes in an idealized yet representative ocean. Plain Language Summary: Numerical simulation of oceanic circulations is incredibly challenging due in part to the large range of scales. There are oscillations on the size of entire ocean basins, through to fine scale centimeter sized turbulence, and everything in between. Global simulations of the ocean typically capture the large scale features on a discrete grid. To ensure these simulations are physically representative, the influence of the missing smaller scales needs to be captured accurately. Here we present a first calculation of parameterizations governing the interactions between the resolved dynamics and the missing small scale eddies, meanfield and topographic features. Reduced resolution simulations incorporating these parameterizations without tuning have statistics in very good agreement with the high resolution simulations from which they were derived. Key Points: Subgrid models of all interaction classes between eddies, meanfield and topography are calculated from reference global ocean simulationsInter‐eddy interactions are modeled both stochastically and deterministically, with deterministic models used for the other interactionsCoarse simulations using these subgrid coefficients without tuning are in very good agreement with those of the high‐resolution simulation [ABSTRACT FROM AUTHOR]

Previous work has found that as the surface warms the large‐scale tropical circulations weaken, convective anvil cloud fraction decreases, and atmospheric static stability increases. Circulation changes inevitably lead to changes in the humidity and cloud fields which influence the surface energetics. The exchange of mass between the boundary layer (BL) and the midtroposphere has also been shown to weaken in global climate models. What has remained less clear is how robust these changes in the circulation are to different representations of convection, clouds, and microphysics in numerical models. We use simulations from the Radiative‐Convective Equilibrium Model Intercomparison Project to investigate the interaction between overturning circulations, surface temperature, and atmospheric moisture. We analyze the underlying mechanisms of these relationships using a 21‐member model ensemble that includes both General Circulation Models and Cloud‐system Resolving Models. We find a large spread in the change of intensity of the overturning circulation. Both the range of the circulation intensity, and its change with warming can be explained by the range of the mean upward vertical velocity. There is also a consistent decrease in the exchange of mass between the BL and the midtroposphere. However, the magnitude of the decrease varies substantially due to the range of responses in both mean precipitation and mean precipitable water. We hypothesize based on these results that despite well understood thermodynamic constraints, there is still a considerable ability for the cloud fields and the precipitation efficiency to drive a substantial range of tropical convective responses to warming. Plain Language Summary: Tropical large‐scale overturning circulations are expected to weaken with warming. This weakening is the result of precipitation increasing at a slower rate than the atmospheric water vapor. Because precipitation and water vapor are important measures of how energy flows through the atmosphere it is important to understand how they will respond to a warming climate. We use two methods to calculate the change of the overturning circulation in 21 different simulations of the tropical atmosphere. This group of 21 models includes high resolution models that resolve cloud systems, and global models with grid‐spacing of about 100 km. We show that a weakening circulation that results from increasing atmospheric stability and increased water vapor is a robust result across most models. But across the group of models there is a large range of magnitudes in the response of the circulation to warming. This variability is well explained by the magnitude of the mean upward vertical velocity. Higher resolution models do not have a narrower range of responses. Narrowing this range of responses will depend on developing a better understanding of what drives the variations in atmospheric stability, surface fluxes of latent energy, and relative humidity. Key Points: The overturning tropical circulation weakens as the surface warms in the majority of radiative convective equilibrium models examinedInter‐model spread of circulation intensity and change of intensity are highly correlated with ascending velocity spread at 500 hPaVariability of both the clear‐sky heating and static stability result in large variations of the subsidence velocity [ABSTRACT FROM AUTHOR]

Proxy records that document the last 2000 years of climate provide evidence for the wide range of the natural climate variability from inter‐annual to secular timescales not captured by the short window of recent direct observations. Assessing climate models ability to reproduce such natural variations is crucial to understand climate sensitivity and impacts of future climate change. Paleoclimate data assimilation (PDA) offers a powerful way to extend the short instrumental period by optimally combining the physics described by General Circulation Climate Models (GCMs) with information from available proxy records while taking into account their uncertainties. Here we present a new PDA approach based on a sequential importance resampling (SIR) Particle filter (PF) that uses Linear Inverse Modeling (LIM) as an emulator of several CMIP‐class GCMs. We examine in a perfect‐model framework the skill of the various LIMs to forecast the dynamics of the surface temperatures and provide spatial field reconstructions over the last millennium in a SIR PF. Our results show that the LIMs allow for skillful ensemble forecasts at 1‐year lead‐time based on GCMs dynamical knowledge with best prediction in the tropics and the North Atlantic. The PDA further provides a set of physically consistent spatial fields allowing robust uncertainty quantification related to climate models biases and proxy spatial sampling. Our results indicate that the LIM yields dynamical memory improving climate variability reconstructions and support the use of the LIM as a GCM‐emulator in real reconstruction to propagate large ensembles of particles at low cost in SIR PF. Plain Language Summary: To detect and attribute anthropogenic climate changes, it is necessary to quantify and understand the wide range of the natural climate variability. The length of the instrumental period is relatively short for investigating slow climate features and their regional impacts. The field of paleoclimate data assimilation offers a way to extend the window of observation and provide physically consistent spatial fields by combining information from both CMIP‐class climate models and proxy records documenting the last 2000 years. These approaches are however confronted with limitations due to the use of expensive global climate models over long periods and the very large number of simulations required to correctly describe the space of possible climate states, which increases exponentially with the dimension of the problem. Here we present a new Proxy Data Assimilation strategy that relies on a statistical model called LIM to reproduce the spatiotemporal dynamics of the surface temperatures simulated by costly CMIP‐class climate models. This climate model emulator is then used to simulate large sets of particles that are optimally combined with the information provided by proxy records to deduce skillful spatial climate fields reconstructions and their low‐frequency variability over the Common Era. Key Points: Linear inverse modeling trained as emulators of CMIP‐class models produce realistic climate and skillful ensemble forecasts as compared to the host‐based modelAn online particle filter incorporating LIMs (SIR‐LIM) overcomes the so‐called data assimilation "curse of dimensionality" problemSIR‐LIM reconstructions for the last millennium allow for uncertainty quantifications due to climate model biases and proxy spatial sampling [ABSTRACT FROM AUTHOR]

Convections still need to be parameterized in general circulation models (GCMs) in the coming decades. Performances of GCMs are significantly influenced by the convection schemes used. In contrast to most conventional cloud models that ignore in‐cloud inhomogeneities, a new convective cloud model explicitly considering in‐cloud inhomogeneities is developed. The new model adopts a single bulk plume approach, but divides the plume into a series of interacting sub‐plumes in order to mimic the transition from the plume core to its edges. We implemented the new cloud model in NCAR Community Atmosphere Model version 5.0 (CAM5) and evaluated its performance. Single‐column tests show a stronger mass flux profile and low‐level detrainment with the new model. In global simulations, the long‐lasting wet bias in the tropical free troposphere of CAM5 is alleviated. Spatial pattern and intensity distribution of precipitation is improved, but the global hydrological cycle is over‐enhanced. The diurnal cycle of tropical precipitation is also better captured though the precipitation peak still occurs earlier than the observations and the amplitude of the cycle tends to weaken. MJO simulation is slightly improved with the new scheme. In addition, the new cloud model generated more tropical cyclones, in better agreement with observations. Plain Language Summary: Convection scheme significantly impacts the performance of general circulation models. In nature, physical properties inside convective clouds are highly inhomogeneous. For example, convective cloud top reaches higher at its core but lower at its edges. However, most convection schemes use a bulk homogeneous plume to describe convective clouds. Here, we developed a new cloud model that explicitly considered the in‐cloud inhomogeneity by dividing the plume into a series of sub‐plumes in order to mimic the cloud property transition from the core to its edges. Implementation and test of this new cloud model in NCAR CAM5 are introduced. Compared to the default scheme, a stronger mass flux profile and low‐level detrainment are produced in single column simulations. Tropical precipitation spatial distribution, diurnal cycle, and intensity distribution are improved. More tropical cyclones are generated using the new cloud model, in better agreement with observations. Key Points: A new deep convective cloud model considering in‐cloud inhomogeneity is proposedThe new cloud model is implemented and tested in NCAR CAM5Tropical precipitation distribution and diurnal cycle are improved using the new cloud model [ABSTRACT FROM AUTHOR]

Mesoscale eddies play an important role in transport of heat and biogeochemical tracers in the global ocean circulation. Resolving these energetic eddies, however, is challenging in ocean general circulation models (OGCM) because it requires a horizontal grid spacing of ≲1/8° that is computationally expensive. As a result, we are required to parameterize mesoscale eddy effects on large‐scale ocean flows. In this work, we introduce a new subgrid‐scale (SGS) model that is developed based on a Taylor series expansion of resolved variables to parameterize subgrid mesoscale eddy transports and momentum fluxes in OGCM. We have performed an a priori study to evaluate the performance of our new gradient model using high‐resolution ocean simulations. Our results show that the gradient model well predicts the actual SGS thickness fluxes in the zonal and meridional directions in coarse‐resolution simulations with the grid spacing ≳1/4°. The unresolved kinetic energy at the ocean surface is also skillfully estimated. More importantly, unlike current mesoscale eddy parameterizations, which are mainly developed based on an assumption of flat bottom ocean, our new SGS model can capture the structure of unresolved standing meanders at the ocean surface. We have also developed a dynamic procedure for setting in non‐dimensional parameters in our new parameterization through a non‐ad hoc and tuning‐free method. Overall, this work suggests that implementing the gradient model in OGCM can improve the model accuracy with an affordable computational cost in eddy‐permitting and non‐eddying simulations. Plain Language Summary: Mass, momentum and heat transports in the global ocean circulation can be significantly influenced by mesoscale eddies. In ocean models, however, it is not computationally affordable to resolve these energetic eddies; hence, we need to parameterize their effects on large‐scale resolved motions. In this work, we propose a new generation of subgrid mesoscale eddy parameterizations based on the gradient of resolved velocity and thickness. Our new model is evaluated using coarse‐grained ocean simulations, where both actual and gradient‐based‐parameterized subgrid transports are computed. The new subgrid parameterization can reduce biases in coarse‐resolution ocean models by skillfully representing unresolved volume and momentum fluxes, leading to higher accuracy in the prediction of zonal and meridional transports with a reasonable computational cost. Key Points: A new subgrid mesoscale eddy parameterization based on the gradients of resolved variables are introduced in ocean modelingThe gradient model can skillfully reproduce unresolved volume transport and momentum fluxes in non‐eddying and eddy‐permitting ocean resolutionsUnlike current subgrid mesoscale models, the gradient model can capture the structure of unresolved standing meanders at the ocean surface [ABSTRACT FROM AUTHOR]

In the construction of an ensemble‐based data assimilation system for a complex fully coupled general circulation model (CGCM), the model state errors at initial time of assimilation have an important influence on assimilation quality. In this study, with the Community Earth System Model (CESM) and Data Assimilation Research Testbed (DART), we found that the influence of initial states errors persists throughout a vicious cycle and cannot be automatically remedied via consequent assimilations. As such, two strategies were applied to alleviate the initial state errors, and a reliable assimilation system was developed. Data assimilation experiments using oceanic observations were conducted over the period from 2005 to 2014 to investigate the impact of these different strategies. The evaluation revealed that the assimilation of observation‐derived climatological data is an effective approach to reduce initial state errors and preserve the balance between different variables to the largest extent, which significantly improved the performance of the assimilation system in the investigated time period. It was further found that the developed assimilation system can produce high‐quality oceanic analysis results comparable to the ECDA and GODAS, two widely used reanalysis products. Perspectives toward further improvement of coupled data assimilation are also outlined. Plain Language Summary: Data assimilation is a widely used way to constrain the model with real world observations. Based on Bayesian Theorem, a good initial state is required as the prior distribution estimate for the assimilation. It is still an open question of how to construct the initial states. In this study, we proposed a new approach to focus on this issue. Via comparison of different experiments, the effectiveness of the approach was revealed. The results indicated that the current data assimilation system with the constructed initial state can produce high‐quality initial conditions for seasonal predictions. Key Points: A new approach is proposed to construct the initial state in the development of an ensemble‐based coupled data assimilation systemThe influence of different initial states is investigated and the effectiveness of the new approach is revealedEvaluation of the developed data assimilation system reveals the good performance of the current system [ABSTRACT FROM AUTHOR]