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]

MATS (Mesospheric Airglow/Aerosol Tomography and Spectroscopy) is a Swedish satellite designed to investigate atmospheric dynamics in the mesosphere and lower thermosphere (MLT). By observing structures in noctilucent clouds over polar regions and oxygen atmospheric-band (A-band) emissions globally, MATS will provide the research community with properties of the MLT atmospheric wave field. Individual A-band images taken by MATS's main instrument, a six-channel limb imager, are transformed through tomography and spectroscopy into three-dimensional temperature fields, within which the wave structures are embedded. To identify wave properties, particularly the gravity wave momentum flux, from the temperature field, smaller-scale perturbations (associated with the targeted waves) must be separated from large-scale background variations using a method of scale separation. This paper investigates the possibilities of employing a simple method based on smoothing polynomials to separate the smaller and larger scales. Using using synthetic tomography data based on the HIAMCM (HIgh Altitude Mechanistic general Circulation Model), we demonstrate that smoothing polynomials can be applied to MLT temperatures to obtain fields corresponding to global-scale separation at zonal wavenumber 18. The simplicity of the method makes it a promising candidate for studying wave dynamics in MATS temperature fields. [ABSTRACT FROM AUTHOR]

All climate models use parameterizations and tuning in order to be accurate. The different parameterizations and tuning processes are the primary source of difference between models. Because models are tuned with present observations of Earth, they may not accurately simulate climates of other planets or palaeoclimate. A model with no adjustable parameter that happens to fit today's observations is probably more universal and should be more appropriate to model palaeoclimate. However, to our knowledge, such a model does not exist or is yet to be developed. This paper aims to improve a parameter-free radiative–convective model that computes a realistic temperature vertical profile to compute the water cycle, giving a value on average tropical precipitation. Although it is known that the radiative transfer constrains the order of magnitude of precipitation, no parameter-free model has yet been able to compute precipitation. Our model finds a precipitation value closer to observations than similar radiative–convective models or some general circulation models (GCMs). [ABSTRACT FROM AUTHOR]

The impacts of climate change increasingly show themselves in many forms in our everyday lives such as heatwaves and droughts. Drought is one of the critical events today for increasing drought frequency. This study focuses on meteorological drought because it directly affects other drought types. Hence, this study focuses on how the future drought conditions will vary under climate change effects in the Mediterranean basin (Turkey). In doing so, this study utilizes precipitation data from three General Circulation Models (GCMs) and three Regional Circulation Models (RCMs). The GCMs are CNRM-CM6, GFDL-CM4, and MPI-ESM1, while the RCMs are (RCA4)-CNRM-CM5, (Reg CM4)-GFDL-ESM2M, and (RCA4)-MPI-ESM-MR. Mitigating biases of the climate models, this study utilizes four statistical downscaling methods (SD), linear scaling (LS), local intensity scaling (LOCI), power transformation (PT), and distribution mapping (DM). Here, the study has two purposes. The main aim of the paper here is to compare the performance of SD methods in improving the representation of observed climate variables in climate models. In addition, the study shows how different methods will affect the spatial drought distribution in the area under the SSP2 4.5 and SSP5 8.5 scenarios. Consequently, the study uses the standardized precipitation index (SPI) and Z-score index (ZSI) to quantify future drought conditions and reaches the following results. The study reveals that mild drought conditions are prevalent in the basin for future periods, and drought indices go down to − 0.55. The study also shows that different SD methods affect the results obtained by each climate model diversely. For example, while the LS method causes the most drought conditions on the results based on CNRM-CM5 and CNRM-CM6, the DM method has a similar impact on outcomes based on GFDL-CM4 and GFDL-ESM2M and causes the most drought conditions. [ABSTRACT FROM AUTHOR]

This paper presents a stochastic approach for the aggregation process rate in the ICOsahedral Nonhydrostatic general circulation model (ICON-AES), which takes subgrid-scale variability into account. This method creates a stochastic parameterization of the process rate by choosing a new specific cloud ice mass at random from a uniform distribution function. This distribution, which is consistent with the model's cloud cover scheme, is evaluated in terms of cloud ice mass variance with a combined satellite retrieval product (DARDAR) from the satellite cloud radar CloudSat and the Cloud–Aerosol Lidar and Infrared Pathfinder Observations (CALIPSO). The global patterns of simulated and observed cloud ice mixing ratio variance are in a good agreement, despite an underestimation in the tropical regions, especially at lower altitudes, and an overestimation in higher latitudes from the modeled variance. Due to this stochastic approach the yearly mean of cloud ice shows an overall decrease. As a result of the nonlinear nature of the aggregation process, the yearly mean of the process rates increases when taking subgrid-scale variability into account. An increased process rate leads to a stronger transformation of cloud ice into snow and therefore to a cloud ice loss. The yearly averaged global mean aggregation rate is more than 20 % higher at selected pressure levels due to the stochastic approach. A strong interaction of aggregation and accretion, however, lowers the effect of cloud ice loss due to a higher aggregation rate. The new stochastic method presented lowers the bias of the aggregation rate. [ABSTRACT FROM AUTHOR]

To adapt to Earth's rapidly changing climate, detailed modelling of thermal stress is needed. Dangerous stress levels are becoming more frequent, longer, and more severe. While traditional measurements of thermal stress have focused on air temperature and humidity, modern measures including radiation and wind speed are becoming widespread. However, projecting such indices has presented a challenging problem, due to the need for appropriate bias correction of multiple variables that vary on hourly timescales. In this paper, we aim to provide a detailed understanding of changing thermal stress patterns incorporating modern measurements, bias correction techniques, and hourly projections to assess the impact of climate change on thermal stress at human scales. To achieve these aims, we conduct a case study of projected thermal stress in central Hobart, Australia for 2040–2059, compared to the historical period 1990–2005. We present the first hourly metre-scale projections of thermal stress driven by multivariate bias-corrected data. We bias correct four variables from six dynamically downscaled General Circulation Models. These outputs drive the Solar and LongWave Environmental Irradiance Geometry model at metre scale, calculating mean radiant temperature and the Universal Thermal Climate Index. We demonstrate that multivariate bias correction can correct means on multiple time scales while accurately preserving mean seasonal trends. Changes in mean air temperature and UTCI by hour of the day and month of the year reveal diurnal and annual patterns in both temporal trends and model agreement. We present plots of future median stress values in the context of historical percentiles, revealing trends and patterns not evident in mean data. Our modelling illustrates a future Hobart that experiences higher and more consistent numbers of hours of heat stress arriving earlier in the year and extending further throughout the day. [ABSTRACT FROM AUTHOR]

This paper describes and analyzes the Reed–Jablonowski (RJ) tropical cyclone (TC) test case used in the 2016 Dynamical Core Model Intercomparison Project (DCMIP2016). This intermediate-complexity test case analyzes the evolution of a weak vortex into a TC in an idealized tropical environment. Reference solutions from nine general circulation models (GCMs) with identical simplified physics parameterization packages that participated in DCMIP2016 are analyzed in this study at 50 km horizontal grid spacing, with five of these models also providing solutions at 25 km grid spacing. Evolution of minimum surface pressure (MSP) and maximum 1 km azimuthally averaged wind speed (MWS), the wind–pressure relationship, radial profiles of wind speed and surface pressure, and wind composites are presented for all participating GCMs at both horizontal grid spacings. While all TCs undergo a similar evolution process, some reach significantly higher intensities than others, ultimately impacting their horizontal and vertical structures. TCs simulated at 25 km grid spacings retain these differences but reach higher intensities and are more compact than their 50 km counterparts. These results indicate that dynamical core choice is an essential factor in GCM development, and future work should be conducted to explore how specific differences within the dynamical core affect TC behavior in GCMs. [ABSTRACT FROM AUTHOR]

In this paper we describe the implementation of the carbon isotopes 13 C and 14 C (radiocarbon) into the marine biogeochemistry model REcoM3. The implementation is tested in long-term equilibrium simulations where REcoM3 is coupled with the ocean general circulation model FESOM2.1, applying a low-resolution configuration and idealized climate forcing. Focusing on the carbon-isotopic composition of dissolved inorganic carbon (δ13 C DIC and Δ14 C DIC), our model results are largely consistent with reconstructions for the pre-anthropogenic period. Our simulations also exhibit discrepancies, e.g. in upwelling regions and the interior of the North Pacific. Some of these differences are due to the limitations of our ocean circulation model setup, which results in a rather shallow meridional overturning circulation. We additionally study the accuracy of two simplified modelling approaches for dissolved inorganic 14 C, which are faster (15 % and about a factor of five, respectively) than the complete consideration of the marine radiocarbon cycle. The accuracy of both simplified approaches is better than 5 %, which should be sufficient for most studies of Δ14 C DIC. [ABSTRACT FROM AUTHOR]

In this paper we introduce a new high-resolution (1/10 °) data-driven dataset of 3D ocean currents developed by the National Research Council of Italy in the framework of the European Space Agency World Ocean Circulation project: the WOC-NATL3D dataset. The product domain extends over a wide portion of the North Atlantic Ocean from the surface down to 1500 m depth, and the dataset covers the period between 2010 and 2019. To generate this product, a diabatic quasi-geostrophic diagnostic model is applied to data-driven 3D temperature and salinity fields obtained through a deep learning technique, along with ERA5 fluxes and empirical estimates of the horizontal Ekman currents based on input provided by the European Copernicus Marine Service. The assessment of WOC-NATL3D currents is performed by direct validation of the total horizontal velocities with independent drifter estimates at various depths (0, 15 and 1000 m) and by comparing them with existing reanalyses that are obtained through the assimilation of observations into ocean general circulation numerical models. Our estimates of the ageostrophic components of the flow improve the total horizontal velocity reconstruction, being more accurate and closer to observations than model reanalyses in the upper layers, also providing an indirect proof of the reliability of the resulting vertical velocities. The reconstructed WOC-NATL3D currents are freely available at 10.12770/0aa7daac-43e6-42f3-9f95-ef7da46bc702. [ABSTRACT FROM AUTHOR]

Xinjiang, located in the arid region of the northwest, is one of the areas most sensitive to global changes. The Kaidu River Basin, situated in the heart of Xinjiang, is one of the sources of China's largest inland river—the Tarim River. The Kaidu River not only bears the responsibility for supplying water for industrial use and agricultural production and people's daily life in the basin, but also plays a crucial role in ecological water supply to the Tarim River. Studying and analyzing the characteristics and trends of meteorological condition in the future under climate change can provide important references and a basis for a deeper understanding of changes in the hydrological process and water resources in the basin. Therefore, this paper selects seven precipitation bias correction methods and four temperature bias correction methods to adjust the precipitation and temperature output data of eight general circulation models of the Sixth Coupled Model Intercomparison Project (CMIP6) within the Kaidu River Basin. The applicability of different bias correction methods in the study area is evaluated, and based on the corrected future meteorological data and calculated extreme meteorological index, the trends of meteorological data (precipitation, temperature) in the future period (2025–2050) under four SSP scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0, SSP5-8.5) in the Kaidu River Basin are analyzed. The results show that: (1) Different types of bias correction methods have different correction focus and effects; their reflections on evaluation indicators are also different. (2) In the future period (2025–2050), the annual precipitation and average temperature in the Kaidu River Basin are higher than those in the historical period (1975–2014). The average annual temperature shows an upward trend in the future, but the annual precipitation shows a downward trend in the future except for the SSP2-4.5 scenario. (3) Compared with the historical period, the extreme precipitation in the future period under the SSP1-2.6, SSP2-4.5, and SSP5-8.5 scenarios is higher than that in the historical period, and the number of rainless days decreases. In the future, under the SSP1-2.6 and SSP5-8.5 scenarios, the probability of meteorological drought events occurring due to high temperatures in the basin may further increase, while under the SSP2-4.5 scenario, the situation of high temperatures and heavy rain in the basin may continue to increase. [ABSTRACT FROM AUTHOR]

The performance of general circulation models (GCMs) in simulating horizontal winds is important because the distribution and variation in horizontal winds are central to investigating atmospheric dynamic characteristics and processes. Also, horizontal wind data can be used to extract some of the required information on gravity waves, tides, and planetary waves. In this context, the present paper evaluates the capability of the Institute of Atmospheric Physics atmospheric general circulation model high-top version (IAP-HAGCM) in simulating the horizontal winds and tides of the troposphere and lower stratosphere by presenting a climatological and statistical comparison against observations of the powerful Beijing mesosphere–stratosphere–troposphere (MST) radar (39.78°N, 116.95°E) during 2012–2014. The results illustrated that the IAP-HAGCM can successfully reproduce the time–altitude distribution of the monthly mean zonal wind and diurnal tide amplitude, albeit with some underestimation. The mean correlation coefficients and root-mean-square error for the zonal (meridional) winds were 0.94 (0.73) and 6.60 m s−1 (2.90 m s–1), respectively. Additionally, the IAP-HAGCM can capture the temporal variation in both the zonal and meridional winds. It is worth noting that, compared with the seven coupled model intercomparison project phase 6 (CMIP6) models, the IAP-HAGCM performs better in meridional wind simulations below 15 km. However, there are discrepancies in altitudinal ranges with large wind velocities, such as the westerly jet, in the transition region of the troposphere and stratosphere, and in February, April, July, and September. It is suggested that model users should take advantage of the model's simulation ability by combining this information regarding when and where it is optimal with their own research purposes. Moreover, the evaluation results in this paper can also serve as a reference for guiding improvements of the IAP-HAGCM. [ABSTRACT FROM AUTHOR]

Using climate-optimized flight trajectories is one essential measure to reduce aviation's climate impact. Detailed knowledge of temporal and spatial climate sensitivity for aviation emissions in the atmosphere is required to realize such a climate mitigation measure. The algorithmic Climate Change Functions (aCCFs) represent the basis for such purposes. This paper presents the first version of the Algorithmic Climate Change Function submodel (ACCF 1.0) within the European Centre HAMburg general circulation model (ECHAM) and Modular Earth Submodel System (MESSy) Atmospheric Chemistry (EMAC) model framework. In the ACCF 1.0, we implement a set of aCCFs (version 1.0) to estimate the average temperature response over 20 years (ATR20) resulting from aviation CO 2 emissions and non-CO 2 impacts, such as NO x emissions (via ozone production and methane destruction), water vapour emissions, and contrail cirrus. While the aCCF concept has been introduced in previous research, here, we publish a consistent set of aCCF formulas in terms of fuel scenario, metric, and efficacy for the first time. In particular, this paper elaborates on contrail aCCF development, which has not been published before. ACCF 1.0 uses the simulated atmospheric conditions at the emission location as input to calculate the ATR20 per unit of fuel burned, per NO x emitted, or per flown kilometre. In this research, we perform quality checks of the ACCF 1.0 outputs in two aspects. Firstly, we compare climatological values calculated by ACCF 1.0 to previous studies. The comparison confirms that in the Northern Hemisphere between 150–300 hPa altitude (flight corridor), the vertical and latitudinal structure of NO x -induced ozone and H 2 O effects are well represented by the ACCF model output. The NO x -induced methane effects increase towards lower altitudes and higher latitudes, which behaves differently from the existing literature. For contrail cirrus, the climatological pattern of the ACCF model output corresponds with the literature, except that contrail-cirrus aCCF generates values at low altitudes near polar regions, which is caused by the conditions set up for contrail formation. Secondly, we evaluate the reduction of NO x -induced ozone effects through trajectory optimization, employing the tagging chemistry approach (contribution approach to tag species according to their emission categories and to inherit these tags to other species during the subsequent chemical reactions). The simulation results show that climate-optimized trajectories reduce the radiative forcing contribution from aviation NO x -induced ozone compared to cost-optimized trajectories. Finally, we couple the ACCF 1.0 to the air traffic simulation submodel AirTraf version 2.0 and demonstrate the variability of the flight trajectories when the efficacy of individual effects is considered. Based on the 1 d simulation results of a subset of European flights, the total ATR20 of the climate-optimized flights is significantly lower (roughly 50 % less) than that of the cost-optimized flights, with the most considerable contribution from contrail cirrus. The CO 2 contribution observed in this study is low compared with the non-CO 2 effects, which requires further diagnosis. [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]

Changes in patterns of meteorological parameters, like precipitations, temperature, wind, etc., are causing significant increases in various extreme events. And these extreme events, i.e., floods, heatwaves, hurricanes, droughts, etc., lead to a shortage of water resources, crop failures, wildfires, and economic losses. However, Global Circulation Models (GCMs) are considered the most important tools for quantifying climate change. Therefore, we selected 20 different GCMs of precipitation in our research, as the frequency of extreme events, like drought and flood, is highly related to changes in precipitation patterns. However, this research introduced a new weighting scheme — MCFSAWS-Ensemble: Monte Carlo Feature Selection Adaptive Weighting Scheme to Ensemble multiple GCMs, whereas, Monte Carlo Feature Selection (MCFS) is one of the most popular algorithms for discovering important variables. However, the proposed weighting scheme (MCFSAWS-Ensemble) is mainly based on two sources. Initially, it evaluates the prior performance of each GCM model to define their relative importance using MCFS. Then, it computes value by value difference between the observed and simulated model. In addition, the application of this paper is based on the monthly time series data of precipitation in the Tibet Plateau region of China. In addition, we used twenty GCMs from the Coupled Model Intercomparison Project Phase 6 (CMIP6) to analyze the implications of the MCFSAWS-Ensemble. Further, we compared the performance of the MCFSAWS-Ensemble scheme with Simple Model Averaging (SMA) through Mean Average Error (MAE) and correlation statistics. The results of this research indicate that the proposed weighting scheme (MCFSAWS-Ensemble) is more accurate than the SMA approach. Consequently, we recommend the use of advanced machine learning algorithms such as MCFS for making accurate multi-model ensembles. [ABSTRACT FROM AUTHOR]

It is shown that a recently developed hybrid modeling approach that combines machine learning (ML) with an atmospheric global circulation model (AGCM) can serve as a basis for capturing atmospheric processes not captured by the AGCM. This power of the approach is illustrated by three examples from a decades‐long climate simulation experiment. The first example demonstrates that the hybrid model can produce sudden stratospheric warming, a dynamical process of nature not resolved by the low resolution AGCM component of the hybrid model. The second and third example show that introducing 6‐hr cumulative precipitation and sea surface temperature (SST) as ML‐based prognostic variables improves the precipitation climatology and leads to a realistic ENSO signal in the SST and atmospheric surface pressure. Plain Language Summary: This paper introduces and tests schemes for efficiently enabling significant expansion of the utility and scope of a recently introduced hybrid modeling technique that combines machine learning with an atmospheric global circulation model (AGCM). Simulation experiments are carried out with an implementation of the approach on a low resolution simplified AGCM. An examination of the simulated atmospheric circulation suggests that the hybrid model can capture dynamical process not captured by the AGCM. Moreover, the addition of precipitation and sea surface temperature (SST) as machine learning predicted physical quantities to the model improves the precipitation climatology and leads to a realistic El Niño‐La Niña signal in the SST and atmospheric surface pressure. Key Points: A hybrid system combining an atmospheric global circulation model (AGCM) with a machine‐learning component can capture processes not captured by the AGCMMachine learning provides a flexible framework to introduce additional prognostic variables into the hybrid modelThe prototype hybrid model tested in the paper is stable and has a realistic climate in decades‐long simulation experiments [ABSTRACT FROM AUTHOR]

In this paper we introduce a new high-resolution (1/10?) data-driven dataset of 3D ocean currents developed by the National Research Council of Italy in the framework of the European Space Agency World Ocean Circulation project: the WOC-NATL3D dataset. The product domain extends over a wide portion of the North Atlantic Ocean from the surface down to 1500 m depth, and the dataset covers the period between 2010 and 2019. To generate this product, a diabatic quasi-geostrophic diagnostic model is applied to data-driven 3D temperature and salinity fields obtained through a deep learning technique, along with ERA5 fluxes and empirical estimates of the horizontal Ekman currents based on input provided by the European Copernicus Marine Service. The assessment of WOC-NATL3D currents is performed by direct validation of the total horizontal velocities with independent drifter estimates at various depths (0, 15 and 1000 m) and by comparing them with existing reanalyses that are obtained through the assimilation of observations into ocean general circulation numerical models. Our estimates of the ageostrophic components of the flow improve the total horizontal velocity reconstruction, being more accurate and closer-to-observations than model reanalyses in the upper layers, also providing an indirect proof of the reliability of the resulting vertical velocities. The reconstructed WOC-NATL3D currents are freely available at https://www.worldocean circulation.org/Products#/metadata/0aa7daac-43e6-42f3-9f95-ef7da46bc702 (Buongiorno Nardelli, 2022). [ABSTRACT FROM AUTHOR]

The ratio of O to N2 number densities (O/N2) at different altitudes is an important parameter in describing thermospheric neutral composition changes and their effects on the ionosphere during geomagnetic storms. However, storm‐induced vertical variations in O/N2 and its dependence on the O and N2 perturbations are still not fully understood. Here, the Thermosphere/Ionosphere Electrodynamics General Circulation Model simulations were used to investigate the responses of thermospheric composition at different pressure levels to the super geomagnetic storm occurred on November 20 and 21 in 2003. Our analysis shows that the behaviors of O/N2 perturbations on different pressure levels are similar above ∼180 km altitude. In the middle and low thermosphere of below ∼300 km, the storm‐time O/N2 decrease is mainly caused by a large reduction of O number density. However, N2 enhancement plays a vital role in O/N2 decreases in the upper thermosphere. The O/N2 enhancement is mainly attributed to the N2 decreases at all pressure levels. The changes of O and N2 number densities at a constant pressure level can be explained by the perturbations of their mass mixing ratio (mmr) and total mass density (ρ). The regions of the O/N2 decrease are characterized by the O mmr decrease and N2 mmr enhancement, whereas the regions of the O/N2 increase are characterized by the O mmr increase and N2 mmr decrease. The ρ value that shows the decrease globally at most pressure levels during the storm either enhance or reduce the O and N2 perturbations. Plain Language Summary: The column O/column N2 density ratio (∑O/N2) was usually used to describe thermospheric neutral composition responses to geomagnetic storms and the storm effects on ionospheric plasma density. However, thermospheric circulation changed considerably during the storm, resulting in discrepancies in composition at different altitudes. Additionally, the daytime electron density changes during geomagnetic storms are more related to those of local O/N2 at a given altitude, not the ∑O/N2. Therefore, it is important to fully understand the storm‐induced vertical variations in O, N2 and O/N2 perturbations. In this paper, the vertical variations in O/N2 and its dependence on the O and N2 perturbations during the 20–21 November 2003 storm are investigated by the numerical simulations. Our results shows that the behaviors of O and N2 perturbations depend much on the altitude, but those of O/N2 on different pressure levels are similar, especially above ∼180 km. This study helps us better understand the physical process of storm‐time ∑O/N2 variations based on the observations. Key Points: In middle and low thermosphere of below ∼300 km, storm‐time decreases of the ratio of O/N2 volume density are mainly caused by O reductionIn the upper thermosphere, N2 enhancement plays a vital role in the decreases of the ratio of O/N2 volume density during the stormAt all pressure levels, storm‐time increases of the ratio of O/N2 volume density depend more on the N2 decreases [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]

This paper evaluated and compared the performance of a statistical downscaling method and a dynamical downscaling method to simulate the spatial–temporal rainfall distribution. Outputs from RegCM4 Regional Climate Model (RCM) and the CanESM2 Atmosphere–Ocean General Circulation Model (AOGCM) were selected for the data scarce Gwadar‐Ormara basin, Pakistan. The evaluation was based on the climatological average and standard deviation for historic (1971–2000) and future (2041–2070) time periods under Representative Concentration Pathways (RCP) 4.5 and 8.5 scenarios. The performance evaluation showed that statistical downscaling is preferred to simulate and project rainfall patterns in the study area. Additionally, the Statistical DownScaling Model (SDSM) showed low R2 values in calibration and validation of the simulations with respect to observed data for the historic period. Overall, SDSM generated satisfactory results in simulating the monthly rainfall cycle of the entire basin. In this study, RegCM4 showed large rainfall errors and missed one rainfall season in the historic period. This study also explored whether the grid‐based rainfall time series of the Asian Precipitation—Highly Resolved Observational Daily Integration Towards Evaluation (APHRODITE) dataset could be used to enlarge and complement the sample of in situ observed rainfall time series. A spatial correlogram was used for observed and APHRODITE rainfall data to assess the consistency between the two data sources, which resulted in rejecting APHRODITE data. For the future time period (2041–2070) under RCPs 4.5 and 8.5 scenarios, rainfall projections did not show significant difference for both downscaling approaches. This may relate to the driving model (CanESM2 AOGCM) and not necessarily suggests poor performance of downscaling; either statistical or dynamical. Hence, the study recommends evaluating a multi‐model ensemble including other GCMs and RCMs for the same area of study. [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 evaluates barotropic tides simulated by a newly developed multi-layer ocean general circulation, ICON-O, and assesses processes and model configurations that can impact the quality of the simulated tides. Such an investigation is crucial for applications addressing internal tides that are much more difficult to evaluate than the barotropic tides. Although not specially tuned for tides and not constrained by any observations, ICON-O is capable of producing the main features of the open-ocean barotropic tides as described by the geographical distributions of amplitude, phase, and amphidromic points. An error analysis shows, however, that the open-ocean tides simulated by ICON-O are less accurate than those simulated by two other ocean general circulation models (OGCMs), especially when not properly adjusting the time step and the parameters used in the time-stepping scheme. Based on a suite of tidal experiments, we show that an increase in horizontal resolution only improves tides in shallow waters. Relevant for using ICON-O with its telescoping grid capacity, we show that spatial inhomogeneity does not deteriorate the quality of the simulated tides. We further show that implementing a parameterization of topographic wave drag improves the quality of the simulated tides in deep ocean independent of the model configuration used, whereas the implementation of a self-attraction and loading (SAL) parameterization in a low-resolution (40 km) version of ICON-O degrades the quality of tides in shallow ocean. Finally, we show that the quality of tides simulated by ICON-O with low resolution (40 km) can be significantly improved by adjusting the time step or the parameters in the time-stepping scheme used for obtaining the model solution. [ABSTRACT FROM AUTHOR]

Among several devastating natural hazards, flooding is a common and serious threat to society causing huge loss of lives, properties, and infrastructure throughout the world. The intensity and frequency of this extreme weather event are expected to increase due to significant changes in the present‐day climate and land use and land cover (LULC) pattern. India has a very systematic and organized structural program and policies but lacks proper implementations, and adverse effect of climate change and the extreme event goes on in society. This paper is an analysis of floods in India and hazards due to climate change and LULC change patterns. Three models, namely "Eco‐biogeography‐based optimization (EBO), Random forest (RF), and Support vector machine (SVM)" were used to obtain the final output to prepare a "Flood susceptibility map". The result was validated through the "Receiver operating characteristics (ROC)" with "Area under curve (AUC)" values. The future rainfall scenario has been estimated by considering the "General circulation models" through different "shared socioeconomic pathways (SSPs)". The values of AUC are 0.915 (EBO), 0.887 (RF), and 0.869 (SVM), respectively. After consideration of different SSPs, the result shows that there is an increasing tendency of flood hazards in the projected period. Among all the employed modelling approaches, the EBO model has notable potential in delineating the possible flood‐prone regions for effective flood planning and management. Decision‐makers can benefit from country‐specific information and regional planner to implement sustainable and long‐term measures to overcome this type of hazardous situation. [ABSTRACT FROM AUTHOR]

The Philippines is highly vulnerable to multiple climate-related hazards due to its geographical location and weak adaptation measures. Floods are the most catastrophic hazards that impact lives, livelihoods, and, consequently, the economy at large. Understanding the ability of the general circulation models to simulate the observed rainfall using the latest state-of-the-art model is essential for reliable forecasting. Based on this background, this paper objectively aims at assessing and ranking the capabilities of the recent Coupled Model Intercomparison Project Phase 6 (CMIP6) models in simulating the observed rainfall over the Philippines. The Global Precipitation Climatology Project (GPCP) v2.3 was used as a proxy to gauge the performance of 11 CMIP6 models in simulating the annual and rainy-season rainfall during 1980–2014. Several statistical metrics (mean, standard deviation, normalized root means square error, percentage bias, Pearson correlation coefficient, Mann–Kendall test, Theil–Sen slope estimator, and skill score) and geospatial measures were assessed. The results show that that CMIP6 historical simulations exhibit satisfactory effectiveness in simulating the annual cycle, though some models display wet/dry biases. The CMIP6 models generally underestimate rainfall on the land but overestimate it over the ocean. The trend analysis shows that rainfall over the country is insignificantly increasing both annually and during the rainy seasons. Notably, most of the models could correctly simulate the trend sign but over/underestimate the magnitude. The CMIP6 historical rainfall simulating models significantly agree on simulating the mean annual cycle but diverge in temporal ability simulation. The performance of the models remarkably differs from one metric to another and among different time scales. Nevertheless, the models may be ranked from the best to the least best at simulating the Philippines' rainfall in the order GFDL, NOR, ACCESS, ENS, MRI, CMCC, NESM, FIO, MIROC, CESM, TAI, and CAN. The findings of this study form a good basis for the selection of models to be used in robust future climate projection and impact studies regarding the Philippines. The climate model developers may use the documented shortcoming of these models and improve their physical parametrization for better performance in the future. [ABSTRACT FROM AUTHOR]

We introduce ModE-Sim (Modern Era SIMulations), a medium-sized ensemble of simulations with the atmospheric general circulation model ECHAM6 in its LR (low-resolution) version (T63; approx. 1.8 ∘ horizontal grid width with 47 vertical levels). At the lower boundary we use prescribed sea surface temperatures and sea ice that reflect observed values while accounting for uncertainties in these. Furthermore we use radiative forcings that also reflect observed values while accounting for uncertainties in the timing and strength of volcanic eruptions. The simulations cover the period from 1420 to 2009. With 60 ensemble members between 1420 and 1850 and 36 ensemble members from 1850 to 2009, ModE-Sim consists of 31 620 simulated years in total. ModE-Sim is suitable for many applications as its various subsets can be used as initial-condition and boundary-condition ensembles to study climate variability. The main intention of this paper is to give a comprehensive description of the experimental setup of ModE-Sim and to provide an evaluation, mainly focusing on the two key variables, 2 m temperature and precipitation. We demonstrate ModE-Sim's ability to represent their mean state, to produce a reasonable response to external forcings, and to sample internal variability. Through the example of heat waves, we show that the ensemble is even capable of capturing certain types of extreme events. [ABSTRACT FROM AUTHOR]

Tidal signatures in sporadic E (Es) layer have been confirmed by observations and simulations. However, the effect of gravity waves (GWs) on the Es layer formation process has not yet been fully understood. In this paper, the modulation of Es layers by GWs is examined through numerical simulations, in which a physics‐based model of Es layer is forced by neutral winds from the High Altitude Mechanistic General Circulation Model that can resolve GWs with horizontal wavelengths longer than 156 km (λh > 156 km). Comparison of the simulation results with and without the GWs (1,350 km > λh > 156 km) forcing reveals that the inclusion of GWs leads to short‐period (1.2–3 hr) density perturbations in Es layers, which are also seen in ground‐based ionosonde observations. At a given time, the metallic ion density at altitudes between 120 and 150 km can either increase (by up to ∼+600%) or reduce (by up to −90%) in response to GW forcing. The relative density perturbations are smaller (by up to 60%) between 90 and 120 km altitude. It is also found that the GW effect on the metallic ion density relates to the longitude, which is mostly explained by the geographical distribution of GWs activity in the mesosphere and lower thermosphere region. The longitudinal variation of the background geomagnetic field plays only a secondary role. Key Points: The modulation of Es layers by gravity waves (GWs) (1,350 km > λh > 156 km) is simulatedGWs cause short‐period (1.2–3 hr) Fe+ density perturbations, leading to production of fine structure in mid‐latitude Es layersLongitudinal dependence in Fe+ density perturbations is mainly explained by the longitudinal distribution of GWs activity in the mesosphere and lower thermosphere region [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]

Climate is closely related to human life, food security and ecosystems. Forecasting future climate provides important information for agricultural production, water resources management and so on. In this paper, historical climate data from 1962–2001 was used at three sites in Tianjin Baodi, Tianjin and Tanggu districts as baseline and the model parameters were calibrated by the Long Ashton Research Station Weather Generator (LARS-WG). 2m-temperatures in 2011–2020 were verified under two scenarios, representative concentration pathway (RCP) 4.5 and RCP8.5 in different atmospheric circulation models with optimal minimum 2m-temperatures at the three sites. From 2031–2050, Tianjin will be using more moderate minimum 2m-temperatures in future simulations. Support vector machines (SVM) were used to optimize the simulated data to obtain more accurate future maximum and minimum 2m-temperatures for the three sites. The results showed that the determinant coefficient of LARS-WG simulation was 0.8 and SVM optimized determinant coefficient was 0.9 which greatly improved the prediction accuracy. The minimum and maximum future 2m-temperatures optimized under European Community Earth System Model (EC-EARTH) were relatively low and the same future 2m-temperatures optimized under Hadley Centre Global Environment Model Earth System (Had-GEM2-ES4) were high especially in the RCP8.5 scenario which simulated 2051–2070 climate. The SVM optimization showed that the maximum and minimum 2m-temperatures were in general agreement with the original simulation values. [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 thermospheric neutral density response to the 7–9 September 2017 storms is investigated based on the Swarm satellite observations and the thermosphere‐ionosphere‐electrodynamic general circulation model (TIEGCM) simulation. The Swarm data depicted a prominent interhemispheric asymmetry (IHA) in the afternoon sector during the second storm, a feature that was yet explained. Driven by realistic high‐latitude electric potential and electron precipitation patterns, the TIEGCM is able to reproduce the observed storm‐time neutral density response. The TIEGCM simulation reveals that the differences in the traveling atmospheric disturbances (TADs) is largely responsible for the observed IHA in the neutral mass density response at low and middle latitudes, whereas the difference in mean molecular mass between the two hemispheres may contribute to the IHA in neutral density at higher latitudes. The IHAs in TADs and mean molecular mass are attributed to the IHA in Joule heating dissipation on the night and dawn sides. Plain Language Summary: During geomagnetic storms, thermospheric neutral density can increase drastically due to intense heat inputs. Enhanced neutral density increases the air drag on satellites which could lower their orbits and lead to faster atmosphere reentering of satellites. In addition, enhanced neutral mass density may cause larger uncertainties in orbit prediction which increases the risk of collisions with other space objects. Therefore, it is essential to improve the understanding of storm‐time neutral density variations to improve satellite operations. Recent satellite observations have revealed that thermospheric neutral mass density response can be interhemispherically asymmetric even under equinox conditions when the background density in the two hemispheres is comparable. This paper presents a case study of such an event based on data analysis and numerical simulations to investigate the physical mechanisms behind this intriguing phenomenon. The event took place on 7–9 September 2017 and the most prominent interhemispheric asymmetry of neutral density occurred in the afternoon sector. Our study found that the interhemispheric asymmetry during this event results from the differences in the generation, propagation and interaction of globally propagating atmospheric waves, namely traveling atmospheric disturbances, as well as the difference in mean molecular mass between the northern and southern hemispheres. Key Points: Thermosphere‐ionosphere‐electrodynamic general circulation model captures the observed interhemispheric asymmetry (IHA) in the neutral density response to the 7–9 September 2017 geomagnetic stormsThe IHA in the neutral density response results from the IHA in traveling atmospheric disturbances and mean molecular massThe observed IHA in the neutral density response is not caused by direct Joule heating dissipated at the same local time [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]

Radiation schemes are physically important but computationally expensive components of weather and climate models. This has spurred efforts to replace them with a cheap emulator based on neural networks (NNs), obtaining large speed-ups, but at the expense of accuracy, energy conservation and generalization. An alternative approach, which is slower but more robust than full emulation, is to use NNs to predict optical properties but keep the radiative transfer equations. Recently, NNs were developed to replace the RRTMGP (Rapid Radiative Transfer Model for General circulation model applications–Parallel) gas optics scheme and shown to be accurate while improving speed. However, the evaluations were based solely on offline radiation computations. In this paper, we describe the implementation and prognostic evaluation of RRTMGP-NN in the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). The new gas optics scheme was incorporated into ecRad, the modular ECMWF radiation scheme. Using two new methods to improve accuracy – a hybrid loss function designed to reduce radiative forcing errors and an early stopping method based on monitoring fluxes and heating rates with respect to a line-by-line benchmark – we train NN models on RRTMGP k distributions with reduced spectral resolutions. Offline evaluation of the new NN gas optics, RRTMGP-NN 2.0, shows a very high level of accuracy for clear-sky fluxes and heating rates. For instance, the RMSE in the shortwave surface downwelling flux is 0.78 W m -2 for RRTMGP and 0.80 W m -2 for RRTMGP-NN in a present-day scenario, while upwelling flux errors are actually smaller for the NN. Because our approach does not affect the treatment of clouds, no additional errors will be introduced for cloudy profiles. RRTMGP-NN closely reproduces radiative forcings for five important greenhouse gases across a wide range of concentrations such as 8×CO2. To assess the impact of different gas optics schemes in the IFS, four 1-year coupled ocean–atmosphere simulations were performed for each configuration. The results show that RRTMGP-NN and RRTMGP produce very similar model climates, with the differences being smaller than those between existing schemes and statistically insignificant for zonal means of single-level quantities such as surface temperature. The use of RRTMGP-NN speeds up ecRad by a factor of 1.5 compared to RRTMGP (the gas optics being almost 3 times faster) and is also faster than the older and less accurate RRTMG, which is used in the current operational cycle of the IFS. [ABSTRACT FROM AUTHOR]

The ocean general circulation model (OGCM) is an essential tool for researching oceanography and atmospheric science. The LASG/IAP climate system ocean model version 3 (LICOM3) is a parallel version of the OGCM. Our goal is to implement and optimize a GPU version of LICOM3 based on compute unified device architecture (CUDA) called LICOM3-CUDA. Considering the characteristics of LICOM3 and CUDA, we design and implement some pivotal optimization methods, including redesigning the numerical algorithms of complicated functions, decoupling data dependency, avoiding memory write conflicts, and optimizing communication. In this paper, we selected two experiments, including 1 ∘ (small-scale) and 0.1 ∘ (large-scale) resolutions to evaluate the performance of LICOM3-CUDA. Under the experimental environment of two Intel Xeon Gold 6148 CPUs and four NVIDIA Quadro GV100s, the LICOM3-CUDA (1 ∘ ) achieves a simulation speed of 114.3 simulation-year-per-day (SYPD). Compare with the performance of LICOM3, the LICOM3-CUDA can run much faster with 6.5 times, and the compute-intensive module achieves over 70 × speedup. In addition, the energy consumption for the simulation year is reduced by 41.3%. As for high-resolution and large-scale simulation, the number of GPUs increased from 96 to 1536 as well as the LICOM3-CUDA (0.1 ∘ ) time consumption decreased from 3261 to 720 seconds with approximately 4.5 × of speedup. [ABSTRACT FROM AUTHOR]

Convolutional neural networks (CNNs) have been widely studied and found to obtain favorable results in statistical downscaling to derive high-resolution climate variables from large-scale coarse general circulation models (GCMs). However, there is a lack of research exploring the predictor selection for CNN modeling. This paper presents an effective and efficient greedy elimination algorithm to address this problem. The algorithm has three main steps: predictor importance attribution, predictor removal, and CNN retraining, which are performed sequentially and iteratively. The importance of individual predictors is measured by a gradient-based importance metric computed by a CNN backpropagation technique, which was initially proposed for CNN interpretation. The algorithm is tested on the CNN-based statistical downscaling of monthly precipitation with 20 candidate predictors and compared with a correlation analysis-based approach. Linear models are implemented as benchmarks. The experiments illustrate that the predictor selection solution can reduce the number of input predictors by more than half, improve the accuracy of both linear and CNN models, and outperform the correlation analysis method. Although the RMSE (root-mean-square error) is reduced by only 0.8%, only 9 out of 20 predictors are used to build the CNN, and the FLOPs (Floating Point Operations) decrease by 20.4%. The results imply that the algorithm can find subset predictors that correlate more to the monthly precipitation of the target area and seasons in a nonlinear way. It is worth mentioning that the algorithm is compatible with other CNN models with stacked variables as input and has the potential for nonlinear correlation predictor selection. [ABSTRACT FROM AUTHOR]

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

Damage from climate change is widespread throughout the world. This change has brought about calamities, the most prevalent of which is the emergence of numerous droughts which are increasingly threatening human lives. In this paper, we studied the spatial and temporal variations of drought under the effect of climate change in the Huaibei Plain, which is a very important agricultural zone in China. Drought has attracted increasing attention in research due to its heavy impact on agriculture, the environment, livelihood, and food security. The SPEI (Standardized Precipitation Evapotranspiration Index) has been used in this study to express and identify drought events in the Huaibei Plain due to climate change. A general circulation model (GCM), HadGEM2-AO, which was the most appropriate for the study area's precipitation simulation, and three Representative Concentration Pathways (RCP), RCP 2.6, RCP 4.5, and RCP 8.5, were used to analyze and compare the drought effect for the baseline (1985–2017) and the future climate scenarios (2025–2090). At 3 and 6 months, the SPEI successfully detects agricultural drought in temporal and spatial variation. However, according to the analysis, more severe agricultural drought events are foreseen in the future than in the baseline because of climate change. SPEI performed better than SPI in detecting drought in the baseline and simulated data due to increased evapotranspiration. Between the SPEI-3 and SPEI-6, the Pearson coefficient correlation reveals a positive association. The Mann-Kendall test was used to cover the two studied periods in order to establish the drought trend. Both decreasing and increasing trends, in different timescales, were detected by Sen's Slope in the baseline and future periods with all RCPs. [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]

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]

Sub-Neptunes (planets with radii between 2 and 4 R ⊕) are abundant around M-dwarf stars, yet the atmospheric dynamics of these planets is relatively unexplored. In this paper, we aim to provide a basic underpinning of the dry dynamics of general low-mean-molecular-weight, temperate sub-Neptune atmospheres. We use the ExoFMS general circulation model (GCM) with an idealized gray-gas radiation scheme to simulate planetary atmospheres with different levels of instellation and rotation rates, using the atmosphere of K2-18b as our control. We find that the atmospheres of tidally locked (TL), temperate sub-Neptunes have weak horizontal temperature gradients owing to their slow rotation rates and hydrogen-dominated composition. The zonal wind structure is dominated by high-latitude cyclostrophic jets driven by the conservation of angular momentum. At low pressures we observe superrotating equatorial jets, which we propose are driven by a Rossbyâ€"Kelvin instability similar to the type seen in simulations of idealized atmospheres with axisymmetric forcing. By viewing the flow in TL coordinates, we find the predominant overturning circulation to be between the day side and night side, and we derive scaling relations linking the TL stream function and vertical velocities to instellation. Comparing our results to the only other GCM study of K2-18b, we find significant qualitative differences in dynamics, highlighting the need for further collaboration and investigation into the effects of different dynamical cores and physical parameterizations. This paper provides a baseline for studying the dry dynamics of temperate sub-Neptunes, which will be built on in part II with the introduction of moist effects. [ABSTRACT FROM AUTHOR]

We introduce ModE-Sim, a medium size ensemble of simulations with the atmospheric general circulation model ECHAM6 in its LR version (T63/approx. 1.8°horizontal with 47 vertical levels) that covers the period from 1420 to 2009. With 60 ensemble members between 1420 and 1850 and 36 ensemble members from 1850 to 2009 ModE-Sim consists of 31620 simulated years in total. The dataset forms the input for a data assimilation procedure that combines historical climate informations with additional constraints from a climate model to produce a novel gridded 3-dimensional dataset of the modern era. Additionally, ModE-Sim on its own is also suitable for many other applications as its various subsets can be used as initial condition and boundary condition ensemble to study climate variability. The main intention of this paper is to give a comprehensive description of the experimental setup of ModE-Sim and to provide an evaluation of the two key variables 2m-temperature and precipitation. We demonstrate ModE-Sim's ability to represent their mean state, to produce a reasonable response to external forcings and to sample internal variability. At the example of heat waves we show that the ensemble is even capable of capturing extreme events. [ABSTRACT FROM AUTHOR]

Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC + MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC + MF parameterization shows weak sensitivity to the vertical grid resolution and model time step. [ABSTRACT FROM AUTHOR]

General circulation models (GCM) show projections of climate variables that when downscaled can be applied to analyse future behaviour in different areas or places. Using them is possible not just to obtain expected values of climate variables but also to calculate their distributions and use those values to assess the effects of climate change at a local level. However, these calculations depend on the GCM selected. In this paper, daily maximum near-surface air temperatures from 21 climate models under representative concentration pathway (RCP) scenarios RCP 4.5 and RCP 8.5 and historic daily maximum temperatures (1990–2019) from nine cities in southern Spain are used with two objectives: first, to investigate past behaviour broken down into a deterministic part and a stochastic part; second, to compare historical data (2006–2019) with the information extracted from the 21 GCMs based on calculating goodness of fit in the period for both deterministic and stochastic parts. The methodology proposed may be useful in selecting a model or a range of models for use in a specific study. The results show positive historical and future trends in maximum daily temperature for these cities. The GCMs with the best fit for each city in this specific case are also presented. [ABSTRACT FROM AUTHOR]

Although the surface energy budget is essential to determine Earth's climate, site measurements of various radiative components are still too scarce to properly characterize their spatial and temporal variations. This has led to the development of a growing number of surface radiation products, mainly including remotely sensed data, model reanalysis data, and simulations using General Circulation Models (GCMs). This collection of papers introduces new techniques, including the use of machine learning methods for radiation estimation, and evaluates and compares various radiation products, as well as their spatio-temporal variations. These studies show large discrepancies among various products across nearly all radiative parameters in either accuracy or spatio-temporal variations. However, remotely sensed radiation products perform relatively better than others. Despite this, there is an urgent need for further efforts to address these discrepancies and improve the accuracy of these estimates. Even though the major radiative parameters including downward shortwave radiation, net longwave radiation, and albedo, from most products show insignificant long-term variation trends on a global scale, only specific regions, such as the Yunnan-Kweichow Plateau (YKP) and regions with permafrost (i.e. Qinghai-Tibet Plateau and Arctic) and glaciers (i.e. Altai Mountains) exhibit remarkable trends. [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]

In this paper we develop a novel framework aimed to significantly reduce biases related to marine stratocumulus clouds in general circulation models (GCMs) while circumventing excessive computational cost requirements. Our strategy is to increase the horizontal resolution using a regionally refined mesh (RRM) over our region of interest in addition to using the Framework for Improvement by Vertical Enhancement (FIVE) to increase the vertical resolution only for specific physical processes that are important for stratocumulus. We apply the RRM off the coast of Peru in the southeastern Pacific, a region that climatologically contains the most marine stratocumulus in the subtropics. We find that our new modeling framework is able to replicate the results of our high-resolution benchmark simulation with much fidelity, while reducing the computational cost by several orders of magnitude. In addition, this framework is able to greatly reduce the long-standing biases associated with marine stratocumulus in GCMs when compared to the standard-resolution control simulation. [ABSTRACT FROM AUTHOR]

Wind stress and turbulent heat fluxes are the major driving forces that modify the ocean dynamics and thermodynamics. In the Nucleus for European Modelling of the Ocean (NEMO) ocean general circulation model, these turbulent air–sea fluxes (TASFs) can critically impact the simulated ocean characteristics. This paper investigates how the various bulk parameterizations used to calculate turbulent air-sea fluxes in NEMOv4 can lead to substantial differences in the estimation of sea surface temperatures (SSTs). Specifically, we study the contributions of different aspects and assumptions of the bulk parameterizations in driving the SST differences in the NEMO global model configuration at 1/4 ∘ of horizontal resolution. These aspects include the use of the skin temperature instead of the bulk SST in the computation of turbulent heat flux components and the estimation of wind stress and turbulent heat flux components, which vary in each parameterization due to different bulk transfer coefficients. The analysis of a set of short-term sensitivity experiments where the only change is related to one of the aspects of the bulk parameterizations shows that parameterization-related SST differences are primarily sensitive to wind stress differences and to the implementation of skin temperature in the computation of turbulent heat flux components. In addition, in order to highlight the role of SST–turbulent heat flux negative feedback at play in ocean simulations, we compare the TASF differences obtained using the NEMO ocean model with the estimations by , who compared the different bulk parameterizations using prescribed SSTs. Our estimations of turbulent heat flux differences between bulk parameterizations are weaker than those found by. [ABSTRACT FROM AUTHOR]

The newly developed fast Earth system model CLIMBER-X is presented. The climate component of CLIMBER-X consists of a 2.5-D semi-empirical statistical–dynamical atmosphere model, a 3-D frictional–geostrophic ocean model, a dynamic–thermodynamic sea ice model and a land surface model. All the model components are discretized on a regular lat–long grid with a horizontal resolution of 5∘×5∘. The model has a throughput of ∼ 10 000 simulation years per day on a single node with 16 CPUs on a high-performance computer and is designed to simulate the evolution of the Earth system on temporal scales ranging from decades to >100000 years. A comprehensive evaluation of the model performance for the present day and the historical period shows that CLIMBER-X is capable of realistically reproducing many observed climate characteristics, with results that generally lie within the range of state-of-the-art general circulation models. The analysis of model performance is complemented by a thorough assessment of climate feedbacks and model sensitivities to changes in external forcings and boundary conditions. Limitations and applicability of the model are critically discussed. CLIMBER-X also includes a detailed representation of the global carbon cycle and is coupled to an ice sheet model, which will be described in separate papers. CLIMBER-X is available as open-source code and is expected to be a useful tool for studying past climate changes and for the investigation of the long-term future evolution of the climate. [ABSTRACT FROM AUTHOR]

Synoptic eddies (with a period of two to eight days) are active in the South China Sea-Philippine Sea (SCS-PS) and control weather variations. In addition, the intensity and frequency of synoptic eddies may change along with variations in sea surface temperatures (SST). This paper presented the influence of SST in the northwest Indian Ocean on synoptic eddies in the lower troposphere over the SCS-PS in June. Our statistical analysis showed a significant negative correlation between the SST in the northwest Indian Ocean and the synoptic scale eddy kinetic energy (EKE) in the SCS-PS. By analyzing the EKE budget of synoptic eddies, we found that the variation in the synoptic scale EKE over the SCS-PS is mainly due to the change in the monthly zonal wind gradient, which affects the barotropic energy conversion between the monthly mean flow and the synoptic eddies. Additionally, the northwest Indian Ocean SST modulates the monthly flow over the SCS-PS by alternating the strength of the Walker circulation in the west Pacific and Indian Ocean. Finally, the influence of SST in the northwest Indian Ocean on EKE in the SCS-PS was reproduced using the simplified atmospheric general circulation model, SPEEDY. [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]

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]