Further improvement and evaluation of nudging in the E3SM Atmosphere Model version 1 (EAMv1): simulations of the mean climate, weather events, and anthropogenic aerosol effects
- 1Pacific Northwest National Laboratory, Richland, Washington, USA
- 2National Center for Atmospheric Research, Boulder, Colorado, USA
Correspondence: Shixuan Zhang (firstname.lastname@example.org) and Kai Zhang (email@example.com)
A previous study on the use of nudging in E3SM Atmosphere Model version 1 (EAMv1) had an unresolved issue; i.e., a simulation nudged to EAMv1's own meteorology showed non-negligible deviations from the free-running baseline simulation over some of the subtropical marine stratocumulus and trade cumulus regions. Here, we demonstrate that the deviations can be substantially reduced by (1) changing where the nudging tendency is calculated in the time integration loop of a nudged EAM simulation so as to improve consistency with the free-running baseline and (2) increasing the frequency of the constraining data so as to better capture strong sub-diurnal variations.
The fact that modification (2) improves the climate representativeness of the nudged simulations has motivated us to investigate whether the use of newer reanalysis products with higher data frequency can help improve nudged hindcast simulations by better capturing the observed weather events. To answer this question, we present simulations conducted at EAMv1's standard horizontal resolution (approximately 1∘) with nudging towards 6-hourly ERA-Interim reanalysis or 6-hourly, 3-hourly, or hourly ERA5 reanalysis. These simulations are evaluated against the climatology of free-running EAMv1 simulations as well as reanalyses, satellite retrievals, and in situ measurements from the Atmospheric Radiation Measurement user facility. For the 1∘ EAMv1 simulations, we recommend using the relocated nudging tendency calculation and the ERA5 reanalysis at 3-hourly or higher frequency.
Simulations used for estimating the anthropogenic aerosol effects often use nudging to help discern signal from noise. The sensitivity of such estimates to the configuration of nudging is investigated in EAMv1, again using the standard 1∘ horizontal resolution. We find that, when estimating the global mean effects, the frequency of constraining data has relatively small impacts, while the choice of nudged variables can change the results substantially. The nudging of air temperature (in addition to horizontal winds) has two non-negligible effects. First, when the constraining data come from reanalysis, the nudging-induced mean bias correction can cause significant changes in the simulated clouds and hence substantially different estimates of the aerosol effects. The impact of the mean bias correction on ice cloud formation has been noted in previous studies and is also seen in EAMv1. For applications like ours, where the preferred configurations of nudging are those capable of providing results consistent with the multi-year free-running simulations, the consequence of the mean bias correction is undesirable. The second important impact of temperature nudging is a significant suppression of adjustments to aerosol forcing, which also causes changes in the estimated aerosol effects. This effect can be seen in simulations nudged to either reanalysis or EAM's own meteorology. These results suggest that nudging horizontal winds but not temperature is a better choice for estimating the anthropogenic aerosol effects.
Nudging (or Newtonian relaxation) is widely used for diagnosing sensitivities of climate simulations to modifications in model formulation and parameters (Lohmann and Hoose, 2009; Zhang et al., 2012; Separovic et al., 2012; Lin et al., 2016) as well as changes in computational methods (e.g., Wan et al., 2014) and external forcing (Kooperman et al., 2012; Zhang et al., 2014). It has been shown that, by constraining the large-scale meteorological conditions (e.g., horizontal winds) toward weather reanalysis or a baseline simulation, nudging can help reduce noise caused by natural variability and hence allow for the detection of signals without long simulations or large ensembles (e.g., Kooperman et al., 2012). However, nudging should be used with care. The configuration of the nudged simulations must be carefully evaluated based on the purpose of the sensitivity experiment. Many studies have shown that the forcing terms introduced by nudging can be sufficiently strong to break the internal balance between the resolved dynamics and parameterized physics (e.g., Jeuken et al., 1996) or to cause significant changes in the model's climate (e.g., Zhang et al., 2014), making the results less useful for interpreting the behavior of the original model.
Sun et al. (2019) evaluated two types of nudged simulations conducted with the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1, Rasch et al., 2019; Xie et al., 2018) at the standard horizontal resolution with approximately 1∘ grid spacing. One type of the simulations was constrained by reanalysis products, and the second type was constrained by meteorological fields written out from a free-running baseline simulation conducted with the same model (hereafter referred to as the “baseline nudging” method). They showed that simulations using baseline nudging closely resembled the free-running simulation for the key meteorological variables evaluated therein, as evidenced by the high spatial and temporal correlations between the nudged and free-running simulations. On the other hand, systematic decreases in the annual mean shortwave cloud radiative forcing (SWCF) were seen in subtropical and tropical regions when nudging was used, with local annual averages as large as 8 W m−2. The discrepancies are inconvenient as they result in inaccuracies in the anthropogenic aerosol effects estimated using baseline nudging.
The study presented here starts with an effort to address these discrepancies. The sequence of calculations related to nudging in EAMv1's time integration loop is reviewed (Sects. 2.2 and 3.1), and the time-step-by-time-step temporal evolution of the model state in the subtropics is analyzed (Sect. 3.2). We demonstrate that the discrepancy issue in 1∘ simulations in Sun et al. (2019) can be substantially alleviated by two revisions of the nudging implementation: first, changing the sequence of calculations in a nudged EAM simulation to improve consistency with the free-running baseline; second, increasing the frequency of constraining data from 6-hourly to 3-hourly to better capture strong sub-diurnal variations. The resulting improvements in climate representativeness are presented in Sect. 3.
Motivated by the improvements, additional simulations and analyses are presented in Sect. 4 to explore the potential benefits of using newer reanalysis products with higher data frequencies in nudged simulations that aim at capturing the observed weather events. In many previous studies (e.g., Telford et al., 2008; Zhang et al., 2014), the reanalysis products used for generating the nudging data were available only four times per day. This was the case, for example, for the ERA-Interim reanalysis (Dee et al., 2011) from the European Centre for Medium-Range Weather Forecasts (ECMWF) as well as the reanalysis of Kanamitsu et al. (2002) from the National Centers for Environmental Prediction and the National Center for Atmospheric Research (NCEP/NCAR). In recent years, reanalysis data with higher temporal frequency have been emerging. For example, MERRA-2 (Gelaro et al., 2017) from the National Aeronautics and Space Administration's Global Modeling and Assimilation Office is available every 3 h, while the ERA5 reanalysis from the ECMWF (Hersbach et al., 2020) has hourly data. On the one hand, using high-frequency reanalysis data for nudging may better constrain a simulation. On the other hand, processing more data before and during a simulation will consume more resources for data processing and storage. Therefore, it is useful to evaluate the benefit of using high-frequency nudging data. Furthermore, since ERA5 is a new reanalysis product that has not been widely used for nudged simulations, it is useful to compare simulations nudged towards ERA5 and ERA-Interim, evaluate hindcast skills of these simulations, and provide a recommendation. For those purposes, we present in Sect. 4 simulations constrained using 6-hourly ERA-Interim reanalysis and 6-hourly, 3-hourly, or hourly ERA5 reanalysis. Hindcast skills of the nudged 1∘ EAMv1 simulations are evaluated against global-scale satellite retrievals of outgoing longwave radiation and precipitation as well as in situ measurements of air temperature, humidity, and horizontal winds from the Atmospheric Radiation Measurement (ARM) user facility. Since one of our primary interests in using nudged simulations is to efficiently estimate the climate impacts of anthropogenic aerosols, we present in Sect. 5 some analysis of the sensitivity of the estimate to the configuration of nudged simulations. Our findings and recommendations are summarized in Sect. 6.
2.1 A brief overview of EAMv1
E3SM is a global Earth system model developed by the U.S. Department of Energy (Golaz et al., 2019). The present study focuses on nudging applications in the E3SM Atmosphere Model version 1 (EAMv1; Rasch et al., 2019; Xie et al., 2018). EAMv1 uses the hydrostatic spectral element (SE) dynamical core on a cubed-sphere mesh (Dennis et al., 2012; Taylor and Fournier, 2010) to solve the equations for large-scale dynamics and tracer transport. The key subgrid-scale physical processes considered in EAMv1 include deep convection (hereafter Deep Cu; Zhang and McFarlane, 1995), turbulence and shallow convection (Golaz et al., 2002; Larson et al., 2002), cloud microphysics (Morrison and Gettelman, 2008; Gettelman and Morrison, 2015; Wang et al., 2014), aerosol life cycle (Liu et al., 2016; Wang et al., 2020), and radiation (Iacono et al., 2008; Mlawer et al., 1997). EAMv1 is interactively coupled with a land model (Oleson et al., 2013).
Figure 1a shows the sequence of dynamics and physics calculations (i.e., the time integration loop) in EAMv1. More detailed descriptions of the time stepping and coupling of physics and dynamics can be found in Zhang et al. (2018) and Wan et al. (2021, 2022). One important feature relevant to the discussion below is that most of the atmospheric processes are numerically coupled using sequential splitting. This means that after a model component (e.g., a parameterization) predicts the rate of change (also called the tendency) of the model state caused by the atmospheric process it represents, the model state will be updated using the predicted tendency before being handed to the next model component (e.g., another parameterization).
The simulations presented in this paper use a horizontal resolution of approximately 1∘ (∼ 110 km). There are 72 layers in the vertical, extending from the Earth's surface to ∼ 0.1 hPa (∼ 64 km). The vertical grid spacing is uneven, with the layer thickness ranging typically from 20 to 100 m near the surface and up to 600 m near the model top.
2.2 Nudging in EAMv1
The nudging implementation in EAMv1 was described and evaluated in Sun et al. (2019), so we only provide a brief introduction here. Nudging constrains the model solution toward prescribed atmospheric conditions for a certain variable by adding a relaxation term to the prognostic equation:
where X in Eq. (1) represents a model state variable like horizontal winds (U, V), temperature (T), or specific humidity (Q). Subscript m refers to the model-predicted value. Subscript p indicates the prescribed field that is taken or derived from either a global weather reanalysis or a free-running simulation using the same model. τ denotes the relaxation timescale. All three quantities, Xm, Xp, and τ, can affect the sign and strength of the nudging-induced forcing.
Pink boxes in Fig. 1a illustrate where the nudging-related calculations occur in the default EAMv1. In a nudged simulation, after the resolved dynamics (see the blue box in the figure) has been calculated, a nudging tendency term in the form of Eq. (1) is calculated for each nudged variable, with Xm being the value of X after the dynamical core. After the entire physics parameterization suite has been calculated, the sum of the parameterization-induced tendencies and the nudging tendencies is passed to the physics–dynamics coupling interface.
It is worth noting that, when an EAM simulation is considered to be a baseline simulation, the dynamical and thermodynamical variables (e.g., U, V, T, Q, and the surface pressure PS) that are archived – and subsequently used in a nudged simulation as the prescribed atmospheric state – are the values saved before the radiation calculation (see the pink dashed box in Fig. 1a). In other words, in the default EAMv1, the Xp on the right-hand side of Eq. (1) is archived before radiation, while the Xm in that same equation corresponds to the model state after the dynamical core. As is discussed in Sect. 3.1, the fact that Xp and Xm correspond to different locations in the time integration loop plays an important role in causing the issue in Sun et al. (2019) that motivated this study.
The EAMv1 simulations presented in this paper are summarized in Table 1. All the simulations involved active atmosphere and land but used prescribed sea surface temperature (SST) and sea ice extension, following the protocol from the Atmospheric Model Intercomparison Project (Gates et al., 1999). The SST and sea ice extension used in this study are weekly data from the National Oceanic and Atmospheric Administration (NOAA) Optimum Interpolation (OI) analysis (Reynolds et al., 2002). Other external forcings, including volcanic aerosols, solar variability, concentrations of greenhouse gases, and anthropogenic emissions of aerosols and their precursors, were prescribed following the World Climate Research Programme (WCRP) Coupled Model Intercomparison Project Phase 6 (CMIP6; Eyring et al., 2016; Hoesly et al., 2018; Feng et al., 2020).
All simulations were performed from 1 October 2009 to 31 December 2010. The first 3 months were discarded as model spinup, and the remaining 1 year of model output was used for analysis. The choice of simulation year was based on convenience, as hourly ERA5 data of 2010 were readily available to us. Sun et al. (2019) have shown that the annual mean cloud radiative forcing and its shortwave and longwave components derived from 1-year nudged simulations are representative of the corresponding longer-term (e.g., 5-year) statistics (see, e.g., Fig. 19 therein).
The anthropogenic aerosol effect we are interested in estimating in this study is the effective radiative forcing (ERF) defined in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, namely, the changes in the TOA radiative fluxes when all physical variables in a climate model are allowed to respond to perturbations except for those concerning the ocean and sea ice (Myhre et al., 2013). Our primary focus is the net TOA flux and its shortwave and longwave components. These are denoted by FNET, FSW, and FLW, respectively, in the remainder of the paper, with positive values indicating fluxes downward (i.e., into the atmosphere). For the readers who have worked with EAM's output, we note that the FSW presented here is EAM's output variable FSNT, while the FLW here is −FLNT, as FLNT in EAM is defined as positive upward. The net flux is calculated as
The changes in FNET, FSW, and FLW caused by anthropogenic aerosols are denoted by ΔFNET, ΔFSW, and ΔFLW, respectively.
We are also interested in the impact of anthropogenic aerosols on the cloud radiative effect (CRE). CRE is defined as the change in a TOA radiative flux caused by the presence of clouds; here we denote the CRE on the net, shortwave, and longwave TOA radiative fluxes by CRENET, CRESW, and CRELW, respectively, with positive values indicating more fluxes into the atmosphere. The CRESW and CRELW presented here are EAM's output variables SWCF and LWCF, respectively, both of which are diagnosed during a simulation by performing the radiation calculations twice (with and without clouds) and then computing the difference. The net CRE is calculated by
The changes in CRENET, CRESW, and CRELW caused by anthropogenic aerosols are denoted by ΔCRENET, ΔCRESW, and ΔCRELW, respectively.
To estimate the anthropogenic aerosol effects (i.e., the Δ quantities) mentioned above, pairs of simulations were conducted. Each pair had an identical experimental setup except that the emissions of aerosols and their precursor gases were set to the values of the year 2010 to represent the present-day (PD) condition in one simulation and the values of the year 1850 to represent the preindustrial (PI) condition in the second simulation. The greenhouse gas concentrations, SST, and sea ice extent were unchanged (i.e., fixed at their year 2010 values). The main differences between PI and PD aerosol emissions included anthropogenic emissions of sulfate, black carbon, organic carbon, primary organic carbon, and the precursors of secondary organic aerosols (applied as yields). Biomass burning emissions were also different between the PD and PI conditions. Dust, sea salt, and marine organic aerosol emissions were calculated online using the surface wind speed and surface properties predicted in each simulation.
Three groups of simulations are presented in this paper. The first group consists of five pairs of 15-month simulations. The first pair are two free-running baseline simulations referred to as CLIM PD and CLIM PI in the remainder of the paper. From the CLIM PD simulation, the before-radiation values of U, V, T, Q, and PS were archived at 1, 3, and 6 h frequencies to constrain some of the subsequent simulations. The other four pairs in group 1 were nudged to 6-hourly temperature output from the CLIM PD simulation but using long relaxation timescales of 10, 10.1, 10.2, and 10.3 d, respectively. These relaxation timescales correspond to values of on the order of 10−6, which resulted in physically insignificant constraints on the simulations. Therefore, the four pairs of nudged simulations can effectively be considered to be free-running although with perturbations introduced to the 3D temperature field that can be used to quantify natural variability in the evolution of the atmospheric state. A similar experimentation strategy has been used by Liu et al. (2018) to generate hindcast ensembles to investigate the radiative forcing of fire-emitted aerosols.
The second group of simulations was nudged to the meteorology archived from the CLIM PD simulation in group 1, regardless of whether the PD or PI emissions were used in the nudged simulations. Nudging was applied at every time step and vertical level using a 6 h relaxation timescale. The simulations labeled DNDG_UV6 and DNDG_UVT6 used the sequence of calculations shown in Fig. 1a (i.e., the default EAMv1), while RNDG_UV6 and RNDG_UVT6 used the revised sequence shown in Fig. 1b and explained in Sect. 3.1. The impact of the revised sequence of calculations is evaluated in Sect. 3.1. The difference between the experiments labeled with “_UV” and “_UVT” is whether only the horizontal winds were nudged (“_UV”) or both winds and temperature were nudged (“_UVT”). The ending number 6 in an experiment name indicates the use of 6-hourly output from CLIM PD. Additional simulations were conducted, also using the revised sequence of calculations but constrained by 3-hourly or 1-hourly output from the CLIM PD simulation (RNDG_UV3 and RNDG_UVT3; RNDG_UV1 and RNDG_UVT1). These simulations are compared to RNDG_UV6 and RNDG_UVT6 in Sect. 3.2 to evaluate the impact of the frequency of the constraining data.
The third group of simulations was nudged toward two reanalysis products, ERA-Interim (Dee et al., 2011) and ERA5 (Hersbach et al., 2020), to assess whether using a newer product (ERA5) and its higher data frequency, instead of the older ERA-Interim at 6 h intervals, can provide nudged hindcast simulations that agree better with the observational data. The reanalysis products were spatially remapped to the cubed-sphere grid and 72 model layers used by EAMv1, following the method used in the Community Earth System Model Version 2 (CESM2; https://ncar.github.io/CAM/doc/build/html/users_guide/physics-modifications-via-the-namelist.html#nudging, last access: 17 August 2021). Topographical differences between EAMv1 and the reanalysis model were taken into account during the vertical interpolation. This group of simulations is compared to global or quasi-global observational data of surface precipitation rate and the top-of-the-atmosphere (TOA) outgoing longwave radiation (OLR) from satellite retrievals (Sect. 4.2) as well as in situ measurements from the ARM user facility (Sect. 4.3).
This section focuses on analyzing the PD simulations listed in group 2 of Table 1. The CLIM PD simulation in group 1 is used as the baseline simulation and referred to as CLIM for brevity.
Before this work, the EAMv1 simulations nudged to 6-hourly output from CLIM were known to show non-negligible differences from CLIM. For example, Fig. 15b in Sun et al. (2019) showed the weakening of 1-year mean SWCF (=CRFSW here) on the order of 2–8 W m−2 in large areas of the subtropical marine and coastal regions when horizontal winds and temperature were both nudged. The differences exceeded 8 W m−2 in some regions over the southeastern Pacific Ocean and South America. Figure 15a in that same paper showed that constraining only the horizontal winds (i.e., no temperature nudging) would remove the discrepancies in most of the subtropical regions, although one would find 4–8 W m−2 of strengthening of the annual mean SWCF close to the coast of Peru. The corresponding discrepancies seen in CRENET and cloud cover are shown in Figs. 2b and C1d in this paper. When winds and temperature are both nudged, we see a substantial number of grid cells in the subtropical Pacific and Atlantic oceans, where the relative differences on the order of 10 % to 20 % are seen in CRENET when compared to the annual mean CRENET in the baseline simulation CLIM (Fig. C2d). Discrepancies of such magnitudes are counterintuitive since the constraining data were generated from the same model driven by the same external forcing. On the other hand, since nudging introduces forcing terms in the form of Eq. (1) to the model's governing equations, any differences between Xm and Xp will lead to deviations from a free-running simulation. Below, we show that such deviations can be significantly reduced by revising the sequence of calculations in nudged simulations, thereby achieving better consistency with the free-running baseline (Sect. 3.1), as well as by increasing the data frequency of the constraining meteorology to better capture higher-frequency variations in time (Sect. 3.2).
3.1 Calculation of the nudging tendency
As mentioned in Sect. 2.2, in EAMv1's nudging implementation before this study, the baseline simulation's atmospheric state was archived before the radiation parameterization, while the nudging-induced forcing (i.e., Eq. 1) was calculated after the dynamical core. Since EAMv1 uses sequential splitting to couple most of the atmospheric processes (see Sect. 2.1), if we use a subscript “DYN” to label the atmospheric state after the dynamical core and a subscript “ARC” to label the atmospheric state being archived, then the old nudging implementation was, effectively,
In our understanding, the first term on the right-hand side of Eq. (3) is the intended nudging tendency, while the second term is inadvertent. Furthermore, the second term can be understood as the total tendency caused by deep convection, turbulence, and stratiform cloud parameterizations scaled by a factor of , where Δt is the physics time step. Since these moist processes are known to strongly affect the atmospheric state, especially temperature and humidity, it is not surprising that nudged simulations using Eq. (2) deviate from their free-running baseline.
When the calculation of the nudging tendency is moved before the radiation parameterization so that Xp from the baseline simulation and Xm in the nudged simulation come from the same location of the time integration cycle (see the schematic in Fig. 1b), we have, as intended,
Sensitivity experiments confirm that using Eq. (4) instead of Eq. (2) significantly reduces discrepancies between the UVT-nudged and free-running simulations, as can be seen by comparing Fig. 2e with Fig. 2d. The annual mean CRENET differences are reduced to within 1 W m−2 for the majority of the grid cells and within 2 W m−2 in the subtropics and tropics, with only a small number of grid cells showing differences between 2 and 5 W m−2. The discrepancies between UV-nudged and free-running simulations are also reduced, although not as significantly (Fig. 2c versus Fig. 2b). The remaining discrepancies are investigated in the next subsection.
3.2 Frequency of the constraining data
Nudged simulations in the literature (e.g., Kooperman et al., 2012; Subramanian and Zhang, 2014; Ma et al., 2014, 2015; Lin et al., 2016; Fast et al., 2016), including our own work (e.g., Zhang et al., 2014; Sun et al., 2019), often used 6-hourly constraining data. The historical reason was that reanalysis data used to be available only four times per day. Such a frequency, on the other hand, can be insufficient for capturing fast variations because of the problem of aliasing.
Figure 3 shows the evolution of lower-troposphere (700 hPa) zonal wind and temperature averaged over the Peruvian stratocumulus region marked by the magenta box in Fig. 2e for a 2 d period starting from 00:00 Z 2 January 2010. In Fig. 3, the black solid lines are time-step-by-time-step output from CLIM where Δt=30 min. The dashed lines are the linearly interpolated time series used in the calculation of nudging tendencies; green, blue, and red correspond to cases in which the constraining data were provided at 1, 3, and 6 h frequencies, respectively. The EAMv1-simulated wind field in the Peruvian stratocumulus region shows prominent 12 h cycles. Linear interpolation of 6-hourly data misses all the local maxima and minima (red line in Fig. 3a), while the interpolation from 3-hourly data provides substantial improvements (blue line in Fig. 3a). The temperature time series in Fig. 3b also shows 12 h variations, although the amplitude is much smaller compared to the diurnal cycle.
Considering the multiscale nature of the atmospheric motions, one can speculate that there are modes of variability that need higher than 3-hourly sampling frequency to avoid aliasing. The sensitivity experiments conducted using 6, 3, and 1 h constraining data (see group 2 of Table 1 and Fig. 4), however, suggest that nudged simulations using 3-hourly data can provide annual mean cloud forcing estimates that agree with CLIM within 1 W m−2 for most grid cells, at least for the 1∘ simulations considered here. In the future, before nudged simulations are conducted at substantially higher resolutions (e.g., 0.25∘ or convection-permitting), it will be useful to find out whether the better-resolved fine-scale motions will require higher frequencies of constraining data.
3.3 Climate representativeness beyond cloud radiative forcing
The investigations discussed in Sect. 3.1 and 3.2 focused on cloud radiative forcing. In Fig. 5, we further evaluate the climate representativeness of the nudged simulations by assessing the annual averages of 20 2D fields that are often examined during model development and tuning. These fields are labeled along the x axis in Fig. 5d and explained in Appendix B. For each of the nudged PD simulations listed in group 2 of Table 1 and each of the 20 fields, we calculated two error metrics with respect to the CLIM PD simulation: one measuring the difference in the global annual mean (Fig. 5a–b) and one measuring the root-mean-square difference in the annually averaged global geographical pattern (Fig. 5c–d).
Consistent with the cloud forcing results shown in Figs. 2 and 4, the revised sequence of calculations and 3 h data frequency have larger impacts on the UVT-nudged simulations than on UV-nudged simulations. Nevertheless, we see a systematic reduction of global mean and pattern errors across all 20 quantities evaluated in Fig. 5 (i.e., yellow bars are substantially shorter than orange bars; green bars are significantly shorter than yellow bars). In simulations RNDG_UVT3 and RNDG_UV3, the errors in global averages are reduced to less than 1 % (green bars in Fig. 5a–b). The errors in geographical patterns are reduced to 2 % or less for the UVT-nudged simulation and 3 % or less for the UV-nudged simulation (green bars in Fig. 5c–d). Comparing Fig. 5c and Fig. 5d, we see lower errors associated with UVT nudging; this is possibly an indication of better consistency between winds and temperature when both are nudged. Further increase in data frequency to 1 h only leads to limited improvements in the simulated geographical patterns. We consistently see the fact that increasing data frequency from 6-hourly to 3-hourly leads to a better agreement of global averages with the free-running simulation, but a further increase to hourly data no longer leads to substantial differences. This can be seen not only in Fig. 5a–b, but also in the additional cloud- and precipitation-related quantities shown in Table S2.
Therefore, for future applications that use 1∘ simulations nudged to the model's own meteorology, we recommend using the revised sequence of calculations depicted in Fig. 1b and 3-hourly constraining data. Future investigations are needed to find out whether nudged simulations at higher spatial resolutions will require more frequent constraining data.
As mentioned in the introduction, a common application of nudging is to force the simulated large-scale meteorological conditions to follow the trajectory of the observed evolution so as to facilitate process-level model evaluation or composite analyses focused on specific types of weather events. In this case, nudged simulations are typically performed using gridded reanalysis products from an operational weather prediction center as the constraining data. The findings from the previous section, especially the conclusion that higher frequency of the constraining data might help better capture important modes of variability, motivated us to evaluate the potential benefits of using more recent reanalysis products such as ERA5 (Hersbach et al., 2020) and MERRA2 (Gelaro et al., 2017). Since ERA5 has the highest data frequency (i.e., hourly) and ERA5 is also known to show better agreement with observations when compared to its predecessor ERA-Interim (Hersbach et al., 2020), we focus on ERA5-constrained simulations in this section and use the sensitivity experiments listed in group 3 of Table 1 to answer the following questions.
What is the impact of nudging on the simulated mean climate? (Sect. 4.1)
Do ERA5-nudged hindcast simulations agree better with observations than the ERA-Interim-nudged simulations? (Sect. 4.2)
How frequently should the nudging data be provided to obtain sufficiently good hindcast skill? (Sect. 4.3)
The discussion in this section focuses on simulations performed under PD forcing conditions.
4.1 Global and regional mean climate
Since the long-term climate simulated by the free-running EAMv1 is known to have non-negligible biases with respect to observational data (Rasch et al., 2019; Xie et al., 2018), nudging towards reanalysis is expected to result in significant changes in the statistical features of the simulated climate. When U and V are nudged to 6-hourly meteorology from ERA-Interim, the annual mean net CRENET can deviate from CLIM by more than −20 W m−2 in the Californian, Peruvian, and Namibian stratocumulus regions (Fig. 6a). When T is also nudged, we see deviations on the order of −10 to −20 W m−2 over the storm tracks and 10 to 40 W m−2 over the trade cumulus regions (Fig. 6b). In terms of global averages, nudging only U and V to 6-hourly ERA-Interim data gives a CRENET very close to the value in CLIM; the shortwave and longwave components deviate from the corresponding values in CLIM by about 0.3 W m−2 (see simulation RNDG_ERAI_UV6 in Table S3). If T is nudged in addition to U and V, the global mean CRENET deviates from the value in CLIM by about −1.7 W m−2, attributable mainly to the longwave component (see simulation RNDG_ERAI_UVT6 in Table S3). These results are consistent with the conclusion from Sun et al. (2019) that ERA-nudged runs differ substantially from CLIM.
The second row of Fig. 6 shows the impact of using ERA5 instead of ERA-Interim while keeping a 6-hourly data frequency. The resulting changes are substantially smaller than the differences between ERA-nudged simulations and CLIM, although we still see some CRENET differences in the subtropics and tropics as large as 10 to 20 W m−2. The relatively small impact of replacing ERA-Interim with ERA5 is expected, as the differences between ERA5 and ERA-Interim are substantially smaller than the differences between either reanalysis and the free-running EAMv1 simulations. (As an example, the annual mean zonal mean pressure–latitude cross sections of air temperature differences are shown in Fig. C3). Increasing the data frequency from 6-hourly to 3-hourly can lead to local changes of 1 to 4 W m−2 in CRENET. These magnitudes are similar to what we have seen in Fig. 4a–d for the simulations nudged to CLIM. Further increasing the data frequency to hourly only introduces negligible changes, again similar to what we have seen in simulations nudged to CLIM (Fig. 4e–f).
A large number of model output variables have been examined in addition to CRE, where we consistently see the differences between ERA-Interim-nudged and ERA5-nudged simulations being substantially smaller than the differences between nudged runs and CLIM, although the magnitudes are non-negligible in some regions. We also consistently see the fact that increasing data frequency from 6-hourly to 3-hourly can lead to discernible changes locally, while a further increase to hourly data no longer leads to substantial differences. The impacts of data frequency on the simulated global averages are generally very small (see Table S3).
4.2 Global and regional weather events
To evaluate the simulation of large-scale weather events, we follow the procedure used for Fig. 5 in Sun et al. (2019) and examine the anomaly correlation between nudged simulations and the observations. Here, an anomaly is defined as the deviation of a simulated or observed quantity from the corresponding (simulated or observed) monthly average at the same geographical location. We first examined the anomaly correlation between the nudged simulations and the corresponding reanalysis (ERA-Interim or ERA5) for temperature, specific humidity, as well as horizontal and vertical winds at various pressure levels. The results were found to be very similar to those presented in Fig. 5 in Sun et al. (2019). ERA-Interim and ERA5-nudged simulations show similar correlations to the corresponding reanalyses (see Fig. S1).
Since the discussion in this section focuses on comparing the hindcast skill of the ERA-Interim-nudged and ERA5-nudged simulations, we present in Figs. 7 and 8 an evaluation against global- and regional-scale satellite retrievals of OLR and surface precipitation rate. Panel (a) in each figure shows the annual average of spatial correlations in different latitude bands; panel (b) in each figure shows the spatially averaged temporal correlations of the anomalies. In Fig. 7, the two upper rows in each panel compare ERA-Interim-nudged and ERA5-nudged simulations that used wind-only nudging, while the lower rows compare simulations that also used temperature nudging. Figure 8 compares ERA5-nudged simulations that used different data frequencies. The EAM-simulated OLR is compared to the NOAA's daily retrievals from the High Resolution Infrared Radiation Sounder (HIRS, Lee et al., 2007) and the Advanced Very High Resolution Radiometer (AVHRR, Stowe et al., 2002). The simulated total precipitation rate is compared to 3-hourly data from the Tropical Rainfall Measuring Mission (TRMM) 3B42V7 product (Huffman et al., 2007; Huffman and Bolvin, 2015) and daily data from the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Climate Data Record (labeled “P-CDR” in the figures here, Ashouri et al., 2015). Further details on the datasets and the comparison procedure can be found in Sun et al. (2019).
The anomaly correlations shown in Fig. 7 indicate that the correlations in the high- and mid-latitude regions are very similar between the ERA5 and ERA-Interim nudged simulations, regardless of whether temperature is constrained. At the low latitudes (20∘ S to 20∘ N), the correlations are higher when ERA5 is used as the constraining data, in terms of both OLR and precipitation and both with or without temperature nudging. Figure 8 indicates that the changes associated with higher data frequency are small for the annual or regional averages shown here.
Figure 9 evaluates the simulated zonal and temporal propagation of the meridionally averaged precipitation rate in boreal spring (March to May) of 2010 over the tropical Pacific Ocean (10∘ S–10∘ N, 60∘ E–90∘ W, upper row) and Northern Hemisphere mid-latitude region (25∘ N–50∘ N, 150∘ E-60∘ W, lower row). Panels (a) and (d) are Hovmöller diagrams plotted from the TRMM data. The bar charts show the correlation between the Hovmöller diagram of TRMM data and the corresponding Hovmöller diagrams plotted from various nudged simulations. Consistent with the anomaly correlations shown in Figs. 7 and 8, in the tropics we see a clear improvement in the simulated propagation of precipitation when ERA5 is used as the constraining data (Fig. 9b), while at the mid-latitudes there are no substantial differences between ERA-Interim-nudged and ERA5-nudged results (Fig. 9e). The impact of frequency of the constraining data is negligible (Fig. 9c, f). The same conclusions can be drawn if we use the root-mean-square error (RMSE) as the evaluation metric (cf. Fig. C4) and if we change the evaluation to a different season (see Figs. S2 and S3).
As an aside, we note that the better precipitation hindcast skills at the mid-latitudes than in the tropics (Fig. 9c versus d) are consistent with the findings in Sun et al. (2019). The impact of constraining temperature appears to be negligible for the 2010 results shown here (Fig. 9c–d, solid fill versus hatching), while Sun et al. (2019) showed better precipitation hindcast skill with additional temperature nudging for spring 2011, especially in the tropics (see Figs. 6 and 7 therein). This suggests that the role of temperature nudging can be case-dependent. Future evaluation in this respect will be useful.
4.3 Comparison to the ARM data
To further assess the hindcast skill of the nudged simulations, we use the radiosonde observations collected by the U.S. Department of Energy's ARM user facility. Radiosonde data are often considered to be reliable high-accuracy measurements (Milrad, 2017) and therefore can provide an objective evaluation of model simulations. Data from three ARM atmospheric observatories are selected to cover different climate regimes, including the Southern Great Plains (SGP) site over the mid-latitude land (https://www.arm.gov/capabilities/observatories/sgp, last access: 13 September 2021), the North Slope of Alaska (NSA) site in the Northern Hemisphere polar region (https://www.arm.gov/capabilities/observatories/nsa, last access: 13 September 2021), and the three central facilities at the Tropical Western Pacific site (Manus, TWPC1; Nauru, TWPC2; Darwin, TWPC3; https://www.arm.gov/capabilities/observatories/twp, last access: 13 September 2021). To our knowledge, radiosonde measurements from these sites were not used in the data assimilation system producing the ERA reanalysis products and hence can be considered to be independent data for the evaluation of the simulations nudged to ERA-Interim or ERA5.
The simulated temperature, relative humidity, and horizontal winds in January 2010 are evaluated against measurements collected in the same time period at SGP (Fig. 10, upper row), NSA (Fig. 10, lower row), and TWP (Fig. 11). The ERA-Interim reanalysis (black dashed lines in the figures) and ERA5 (black solid lines) are also included for comparison. The 6-hourly model output and reanalysis products were horizontally remapped to the locations of ARM sites using bilinear interpolation. For each of the meteorological quantities shown here, the RMSEs between the ERA-nudged simulations (or ERA analyses) and the ARM measurements were calculated with all available vertical profiles at the sites in January 2010. For each of the variables shown in Figs. 10 and 11 (i.e., T, RH, U, or V), the numbers of available vertical profiles at SGP, NSA, TWPC1, TWPC2, and TWPC3 were 121, 63, 49, 57, and 127, respectively. (ARM sites provide data four times a day at SGP and TWPC3 and twice a day at NSA, TWPC1, and TWPC2.) The temporal correlations between EAM simulations or ERA analyses and the ARM measurements are shown in Figs. C5 and C6 in the Appendix.
As expected, reanalyses (black lines in Figs. 10, 11, C5, and C6) show better agreement with the ARM radiosonde data compared to the nudged EAM simulations (colored lines). ERA5 (solid black in Figs. 10 and C5) is in general better than ERA-Interim (dashed black) at the mid-latitude SGP site and the high-latitude NSA site. For the three ARM TWP sites in the tropics, ERA5 is not always better than ERA-Interim. For example, ERA5's zonal wind field (U, solid black in Fig. 10d) shows larger RMSEs below 500 hPa compared to the ERA-Interim (dashed black).
The ERA-nudged simulations show good agreement with the ARM radiosonde measurements at the mid-latitude (SGP) and high-latitude (NSA) sites (Fig. 10). Compared to the ERA-Interim-nudged simulations (colored dashed lines), the ERA5-nudged simulations (colored solid lines) have slightly better hindcast skills. In addition, EAM simulations with temperature nudging (red, orange, and green lines) show overall better hindcast skills, regardless of which ERA product was used as the constraining data. Slightly better hindcast skills for horizontal winds can be obtained by using the 3-hourly ERA5 data (green lines) for nudging instead of using 6-hourly data (red lines). Using 3-hourly (green lines) or hourly (orange lines) constraining data gives very similar results.
Consistent with the experiences reported in the literature (e.g., Jeuken et al., 1996; Sun et al., 2019), weather events in the tropics are less well-constrained by nudging. Compared to the ARM SGP and NSA sites (Figs. 10, C5), the magnitudes of the RMSEs in ERA-nudged simulation at the three ARM TWP sites are in similar ranges (Fig. 11), while the temporal correlations are smaller at the tropical sites, especially for temperature and relative humidity (Fig. C6). At the tropical sites, we do not see systematic improvements when switching from EAM-Interim to EAM5 for the constraining data or when increasing the data frequency, although the UVT nudging (red, orange, and green lines in Figs. 11 and C6) still provides better hindcast skills than the UV nudging (blue lines).
Nudging has been recognized as a useful and computationally efficient technique for estimating the anthropogenic aerosol effects in global climate models (Kooperman et al., 2012; Zhang et al., 2014, 2016; Ghan et al., 2016; Liu et al., 2018). In this section, we evaluate the impact of nudging configuration on the estimated anthropogenic aerosol effects in EAMv1. It is of practical value to identify nudged simulations that are capable of providing estimates consistent with the results from free-running simulations, as the EAM developers have identified the anthropogenic aerosol effects as one of the key aspects that need more attention in the future development and evaluation of the model (Golaz et al., 2019; Zhang et al., 2022). Similarly to previous studies, we calculate the anthropogenic aerosol effects by contrasting a pair of EAMv1 simulations conducted with PD (year 2010) and PI (year 1850) emissions of the anthropogenic aerosols and precursors following the CMIP6 protocol (Eyring et al., 2016; Hoesly et al., 2018; Feng et al., 2020).
5.1 Results from the free-running EAMv1
As explained in Sect. 2.3 and summarized in group 1 of Table 1, we carried out five pairs of 1-year simulations without nudging (the “CLIM” runs) or with very weak nudging (simulations “CLIMp1”–“CLIMp4”). The five-member mean, 1-year mean, globally averaged PD − PI difference in the TOA net radiative flux, ΔFNET, is about −1.7 W m−2. The shortwave component is W m−2, and the longwave component is ΔFLW=0.7 W m−2 (see Table S4). These numbers are consistent with the effective aerosol forcing estimates reported in Sect. 6.1 of Golaz et al. (2019). The PD − PI differences in the shortwave and longwave CRE in our CLIM ensemble are W m−2 and ΔCRELW=0.6 W m−2, respectively (Table S4), which gives a net effect of W m−2.
5.2 Impacts of temperature nudging
In Fig. 12, we compare various configurations of the nudged simulations against the CLIM ensemble in terms of the annual mean PD − PI differences averaged over the globe (left column) or the tropics (20∘ S–20∘ N, right column). All results are normalized by the ensemble mean of CLIM. The thick black whiskers attached to the gray bars indicate the 2 standard deviation ranges of the CLIM ensemble. The non-normalized data can be found in Tables S4 and S5 in the Supplement. All of the nudged simulations shown in the figure used the revised sequence of calculations and 3-hourly constraining data. Two of the nudged simulations were constrained by ERA5 (the colored bars with hatch filling) and the other two by CLIM PD (the colored bars with solid filling). Our first focus is to compare simulations conducted using UV nudging to those using UVT nudging.
Generally speaking, we can expect the nudging of temperature to have three potential impacts.
First, if a model has significant and systematic temperature biases relative to reanalysis, then nudging temperature towards reanalysis will introduce a mean bias correction, which can trigger responses of the model atmosphere in the nudged simulations and consequently cause differences (compared to the free-running simulations) in the aerosol effects.
Second, it is worth noting that the anthropogenic aerosol effect we are trying to estimate here is the ERF, which, by definition, is the change in the TOA radiative fluxes when all physical variables in a climate model are allowed to respond to perturbations except for those concerning the ocean and sea ice (Myhre et al., 2013). Because anthropogenic aerosols have significant cooling effects, when a pair of simulations with PD and PI emissions are both nudged to the same temperature fields, the aerosol-induced temperature differences between PD and PI conditions can be substantially reduced regardless of which constraining data are used. This can then lead to a significant suppression of the simulated aerosol effects and hence inconsistency between the temperature-nudged and free-running simulations.
Third, when the constraining temperature data are provided at a relatively low frequency, we might encounter situations similar to those illustrated in Fig. 3; i.e., the linearly interpolated constraining data can fail to capture high-frequency variations in temperature and hence cause responses in the simulated atmospheric state. The impact of such aliasing is expected to be discernible only if the high-frequency modes of variability have substantial roles in determining the anthropogenic aerosol effects.
We show in Sect. 5.3 that the third effect (aliasing) is indeed small in EAMv1. In the remainder of this subsection, we focus on quantifying the first two effects listed above.
Figure 12a shows the global mean PD − PI differences in the TOA fluxes and CRE, and panel (b) of the figure shows the averages over the tropics. The second row of Fig. 12 shows the PD − PI differences in the global and tropical mean total cloud fraction (ΔCLDTOT), cloud liquid and ice water path (ΔLWP and ΔIWP), and total precipitable water (ΔTMQ). All results shown in this figure have been normalized by the ensemble mean of CLIM. The non-normalized data can be found in Tables S4 and S5. Keeping in mind the ensemble spread of the CLIM simulations, we see that the estimates obtained with UV nudging (pink bars) are consistent with the estimates from CLIM, while the estimates obtained with UVT nudging (green bars) show statistically significant deviations from CLIM. This is true for both the ERA5-nudged and CLIM-nudged simulations.
5.2.1 Mean bias correction
The impact of mean bias correction caused by nudging EAM's temperature to reanalysis can be evaluated by comparing simulations RNDG_ERA5_UVT (hatched teal bars in Fig. 12) and RNDG_UVT (solid teal bars) against CLIM (gray bars and whiskers). Overall, the hatched teal bars show larger deviations than the solid teal bars, suggesting that the mean bias correction has significant impacts.
Among the quantities shown in Fig. 12, the tropical mean ΔFLW, ΔCRELW, ΔIWP, and ΔCLDTOT show the largest reduction when temperature is nudged to ERA5 instead of CLIM (panels b and d). This can be explained by the zonal and annual mean temperature differences (ΔT) and in-cloud ice number concentration differences (ΔICINC) shown in Fig. 13. Figure 13a suggests that EAMv1's climatology, when compared to ERA5, features cold biases on the order of 1–2 K in the upper troposphere over the tropical and mid-latitude regions where small ice crystals are often formed through homogeneous ice nucleation. These small ice crystals are known to have a large impact on the simulated CRE. Nudging EAM's temperature towards ERA5 leads to a warmer base state and weakened homogeneous ice nucleation compared to CLIM (Fig. C7b). Consequently, the PD − PI changes in aerosol and precursor emissions cause substantially smaller ΔICNIC compared to CLIM (Fig. 13d versus Fig. 13c), which explains the significant reduction in ΔFLW and ΔCRELW shown as hatched green bars in Fig. 12a–b. This reasoning is consistent with the finding in Zhang et al. (2014) that temperature nudging in EAMv1's predecessor model CAM5 led to a substantial decrease in the ice cloud amount and a weaker impact of anthropogenic aerosols on longwave radiation. Although nudging temperature towards reanalysis can improve the agreement between the observed and EAM-simulated long-term climate, it significantly changes the simulated cloud amounts and properties as well as the CRE in EAMv1 and hence does not achieve the goal of obtaining estimates of the anthropogenic aerosol effects that are consistent with those in the free-running EAM simulations.
5.2.2 Suppression of temperature changes
The impact of nudging PD and PI simulations with the same temperature data and hence suppressing the ERF can be evaluated by contrasting the two simulations RNDG_UV3 and RNDG_UVT3 (solid bars in pink and teal in Fig. 12). For reference, we show in Fig. 14 the zonal mean PD − PI temperature difference in CLIM and in the simulations nudged to CLIM PD. In the free-running simulations, the increased emissions of anthropogenic aerosols and precursors lead to considerable cooling in the lower troposphere in the Northern Hemisphere mid- and high-latitude regions as well as warming in the upper troposphere (Fig. 14a). These temperature changes are captured by the simulations nudged to U and V from CLIM PD (Fig. 14b), while the simulations with temperature nudging show negligible ΔT (as expected).
Consistent with the ΔT shown in Fig. 14, we see that temperature nudging has considerable impacts on the simulated PD − PI differences in LWP, IWP, and CLDTOT (Fig. 12, lower row, teal versus pink solid bars). The impact on global mean ERF is more clearly seen in the shortwave component (Fig. 12a), possibly due to the cloud changes caused by the lack of cooling (compared to CLIM) in the mid- and high-latitude lower troposphere when temperature is nudged (Fig. 14c versus Fig. 14b).
5.2.3 Combined effects
The NDG_ERA5_UVT3 configuration is affected by both the mean state correction and the suppression of PD − PI temperature changes, and hence it is not surprising that the results (hatched teal bars in Fig. 12) deviate most substantially from CLIM compared to all other nudged simulations shown in the figure. For global averages (Fig. 12a), the ΔFSW estimated by NDG_ERA5_UVT3 is about 25 % lower than CLIM, and the ΔFLW is about 50 % lower than CLIM (see also Table S4). In the tropics (Fig. 12b), there is a 33 % decrease in the estimated ΔFSW and a 63 % decrease in ΔFLW (see also Table S5). Large discrepancies between NDG_ERA5_UVT3 and CLIM are also seen in ΔCRESW and ΔCRELW (Fig. 12, upper row) and in the cloud-related fields (Fig. 12, lower row).
It is worth noting that the discrepancies in ΔFNET might appear to be not as large. For example, we see a 15 % difference in the global mean in Fig. 12a and Table S4 and a 17 % difference in the tropical average (Fig. 12b and Table S5), but the smaller differences in the net fluxes are results of the cancellation of large changes in the shortwave and longwave components caused by temperature nudging.
Overall, our results discussed above suggest that nudging the horizontal winds but not temperature is the preferred simulation configuration for estimating the anthropogenic aerosol effects, which is consistent with the results reported in Zhang et al. (2014). When the PD − PI simulations are constrained by the same temperature data, nudging can suppress the adjustments to aerosol forcing. When the constraining data come from a different model (e.g., reanalysis), one can get an additional effect; i.e., the mean bias corrections can significantly modify the simulated clouds and their responses to aerosol forcing. In both cases, temperature nudging can lead to estimates of the anthropogenic aerosol effects that are significantly different from the results obtained from the free-running simulations, which is undesirable for many model development and evaluation studies.
5.3 Impacts of the frequency of constraining data
In Fig. 15, we evaluate the impacts of the frequency of the constraining data used in nudged simulations. At least for the annual mean effects, the results obtained from simulations using 6-hourly constraining data (orange bars in the figure) are very similar to those obtained using 3-hourly constraining data (blue bars), regardless of whether UV nudging (Fig. 15, upper row) or UVT nudging (Fig. 15, lower row) is used. The small impact of constraining data frequency on global and tropical mean Faer estimates is expected. As shown in Sect. 3.2, the impact of constraining data frequency on present-day simulations is sizable only in limited regions where strong diurnal variations exist. Therefore, using 6-hourly constraining data in nudged simulations is sufficient for estimating the global and annual mean Faer.
Nudging has been widely used in the development and evaluation of global and regional atmospheric models. In this work, we further improved the nudging implementation in EAMv1 compared to the work of Sun et al. (2019) and evaluated the impact on the climate representativeness, the hindcast skill of nudged simulations, and the estimation of anthropogenic aerosol effects.
The study was motivated by an unresolved issue in Sun et al. (2019); i.e., a nudged EAMv1 simulation constrained by EAMv1's own meteorology showed non-negligible local deviations from the baseline, with annually averaged CRESW (SWCF) discrepancies as large as 4–8 W m−2 over some of the subtropical marine stratocumulus and trade cumulus regions. Two reasons were identified. First, EAMv1 writes out meteorological fields (from a baseline simulation) for nudging before the radiation parameterization, but the nudging tendency is calculated at a different location in the time integration loop, i.e., after the dynamical core. This inconsistency introduced an unintended contribution to the nudging tendency that was proportional to the effects of deep convection, shallow convection, and cloud microphysics on the simulated atmosphere (Sect. 3.1). Second, the EAM-simulated winds and temperature in the lower troposphere were found to have high-frequency modes with non-negligible magnitudes. For example, the zonal wind in the Peruvian stratocumulus region was found to have a prominent 12 h cycle. Such variations cannot be properly captured by a 6-hourly sampling frequency, hence resulting in significant aliasing issues with the constraining data used for nudging (Sect. 3.2). We showed that, by moving the calculation of the nudging tendency to the same location as the data output (Fig. 1b) and by increasing the frequency of constraining data to 3-hourly, one could largely remove the discrepancies between a 1∘ free-running EAMv1 simulation and a 1∘ nudged simulation constrained by EAM's own meteorology. Further increasing the data frequency to hourly only provided marginal improvements. For future studies that nudge EAM towards its own meteorology, we recommend using the revised implementation and the 3-hourly constraining data for 1∘ simulations. Whether simulations performed at higher horizontal resolutions can benefit from higher data frequency remains to be investigated. For users of EAM, we have provided in Table A1 the nudging-related namelist settings for two of the simulations discussed in this paper to demonstrate how to turn on the revised sequence of calculations and change the constraining data frequency.
The abovementioned improvements further motivated us to investigate the potential benefits of using the ERA5 reanalysis, which is available at a higher frequency compared to ERA-Interim, for nudged hindcast simulations. In terms of the annual mean fields, there were discernible but small regional changes when switching from ERA-Interim to ERA5 or changing the constraining data frequency when using ERA5. The impacts on global mean climate were found to be small (Sect. 4.1). Satellite retrievals of OLR and precipitation were used to evaluate the model's skill in capturing real weather events. When ERA5 was used instead of ERA-Interim, the simulated OLR and precipitation were significantly improved, especially in the tropics. We also evaluated the nudged simulations using radiosonde measurements from several ARM sites in different climate regimes. At the SGP and NSA sites, the simulated horizontal winds, temperature, and relative humidity were systematically improved when replacing ERA-Interim with ERA5 and when using higher-frequency nudging data. Significant improvements were also seen in the mid- and high-latitude ARM sites. At the tropical sites (TWPC1, TWPC2, and TWPC3), the improvements were not as significant. At SGP and NSA, nudging winds and temperature together was found to further improve the hindcast skill of the simulations. Overall, the good agreement in the simulated and observed meteorological conditions provides a good basis for possible future studies that use ARM measurements to help identify parameterization deficiencies and improve the representation of cloud and aerosol-related atmospheric processes in EAM.
Last but not least, we evaluated the impact of nudging on the estimated anthropogenic aerosol effects. Results show that the frequency of the constraining data has negligible impacts on the estimated global and tropical averages of annual mean TOA fluxes and CRE, while the impact of temperature nudging is large. Similar to conclusions from earlier studies, we recommend nudging the horizontal winds but not the air temperature when attempting to obtain estimates of the aerosol effects that are consistent with the estimates from free-running simulations. The reason is twofold: when the simulations forced by PD and PI emissions are constrained by the same meteorological fields, nudging temperature can strongly suppress the temperature responses to the aerosol perturbations and consequently affect other atmospheric adjustments; when the constraining meteorology comes from a different model (e.g., reanalysis), one can get an additional effect; i.e., the effective mean bias corrections introduced by temperature nudging can significantly modify the simulated clouds and their responses to aerosol perturbations. In both cases, temperature nudging can lead to estimates of the anthropogenic aerosol effects that are significantly different from the results obtained from the free-running simulations, which is undesirable for many model development and evaluation studies. Our results from EAMv1 showed that when temperatures in the PD and PI simulations were both nudged towards a free-running PD simulation, the shortwave component of the aerosol effects on the TOA radiative flux and CRE was significantly underestimated. When the PD and PI simulations were nudged towards the PD temperature from the ERA5 reanalysis, both the shortwave and longwave components of the aerosol effects on TOA radiative flux and CRE were significantly underestimated. While the percentage discrepancies in the net TOA flux and CRE appeared to be considerably smaller, this was the result of the cancellation of large discrepancies in the shortwave and longwave components. In contrast, nudging horizontal winds but not temperature provided estimates that were reasonably consistent with the results from the free-running simulations, regardless of whether the nudged simulations were constrained by ERA5 or EAM's own meteorology.
We note that the 1∘ configuration of EAMv1 was used in this study. Due to the relatively coarse grid spacing in our simulations, the benefits of the high temporal and spatial resolutions of the ERA5 data might not have been fully revealed. As pointed out by Jeuken et al. (1996), the linear temporal interpolation in nudging can become more questionable for higher-resolution simulations as more short timescale processes are resolved. Also, compared to ERA-Interim, ERA5 can provide more accurate meteorological variables at finer spatial scales, so the ERA5-nudged simulation might perform even better at high resolutions than seen here in the 1∘ simulations. The high-resolution configuration of EAMv1 was substantially more expensive and hence not used in this study. The newer version of EAM released in September 2021, EAMv2, has become significantly less expensive thanks to the use of a different physics grid (Hannah et al., 2021) and a more efficient large-scale advection scheme (Bradley et al., 2019). Hence, it will be useful to further explore nudged EAM simulations at higher resolutions.
The physical quantities labeled with short names along the x axis in Fig. 5 are explained below in Table B1. The error metrics shown in Fig. 5a–b are the relative differences in the simulated global mean annual averages. The error metrics shown in Fig. 5c–d are the relative differences in the annual mean global patterns (i.e., geographical distributions). Following Wan et al. (2021), a relative difference in the global mean annual average is normalized by the global mean annual average from CLIM. A relative difference in the annual mean global pattern is defined as the centered root-mean-square (rms) difference between the pattern in the test simulation and the pattern in the reference simulation, normalized by the rms of the pattern in the reference simulation. A “pattern” here represents the annual mean, global, and geographical distribution of a physical quantity.
The data used by figures in the paper can be found on Zenodo at https://doi.org/10.5281/zenodo.6988262 (Zhang and Zhang, 2022). This includes (1) post-processed EAMv1 output, (2) reanalysis and satellite data interpolated to the same grid as in (1), and (3) radiosonde measurements at selected sites. The radiosonde measurements were obtained from the ARM (https://doi.org/10.5439/1021460, ARM user facility, 2002), a U.S. Department of Energy Office of Science user facility managed by the Biological and Environmental Research program. Satellite data used in this paper were obtained from the following sources:
https://psl.noaa.gov/data/gridded/data.olrcdr.interp.html (last access: 29 August 2022) (NOAA PSL, 2011) for the AVHRR OLR data.
The supplement related to this article is available online at: https://doi.org/10.5194/gmd-15-6787-2022-supplement.
KZ and HW initiated this study. HW identified the two issues in the implementation of nudging in the default EAMv1. SZ conducted all the simulations, processed the model output, and carried out the analyses with input from KZ and HW. SZ, KZ and HW wrote the paper and JS provided feedback. All the co-authors contributed to the revisions.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors thank William Collins and the anonymous referee for their careful review and helpful comments. The EAMv1 code (https://doi.org/10.11578/E3SM/dc.20180418.36, E3SM Project, DOE, 2018) was obtained from the E3SM project sponsored by DOE BER. Computing resources were provided by the DOE BER Earth System Modeling program's Compy computer cluster located at Pacific Northwest National Laboratory (PNNL).
This research has been primarily supported by the Laboratory Directed Research and Development Program at Pacific Northwest National Laboratory (PNNL), a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy (DOE) under contract DE-AC05-76RL01830. Kai Zhang has been partially supported by the DOE Office of Science, Biological and Environmental Research (BER), via the Energy Exascale Earth System Model (E3SM) project. Hui Wan has been supported by the DOE’s Scientific Discovery through Advanced Computing (SciDAC) program via a partnership in Earth system model development between BER and Advanced Scientific Computing Research (ASCR).
This paper was edited by Fiona O'Connor and reviewed by William Collins and one anonymous referee.
Ashouri, H., Hsu, K.-L., Sorooshian, S., Braithwaite, D. K., Knapp, K. R., Cecil, L. D., Nelson, B. R., and Prat, O. P.: PERSIANN-CDR: Daily precipitation climate data record from multisatellite observations for hydrological and climate studies, B. Am. Meteorol. Soc., 96, 69–83, 2015. a
Atmospheric Radiation Measurement (ARM) user facility: Balloon-Borne Sounding System (SONDEWNPN). December 2009–January 2011: Southern Great Plains (SGP) Central Facility (C1), North Slope Alaska (NSA) Central Facility (C1), Tropical Western Pacific (TWP) Central Facility (C1, C2, C3), compiled by: Keeler, E. and Kyrouac, J., ARM Data Center [data set], https://doi.org/10.5439/1021460, 2002. a
Bradley, A. M., Bosler, P. A., Guba, O., Taylor, M. A., and Barnett, G. A.: Communication-Efficient Property Preservation in Tracer Transport, SIAM J. Sci. Comput., 41, C161–C193, https://doi.org/10.1137/18M1165414, 2019. a
Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011. a, b
Dennis, J. M., Edwards, J., Evans, K. J., Guba, O., Lauritzen, P. H., Mirin, A. A., St-Cyr, A., Taylor, M. A., and Worley, P. H.: CAM-SE: A scalable spectral element dynamical core for the Community Atmosphere Model, Int. J. High Perform., 26, 74–89, https://doi.org/10.1177/1094342011428142, 2012. a
Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958, https://doi.org/10.5194/gmd-9-1937-2016, 2016. a, b
Fast, J. D., Berg, L. K., Zhang, K., Easter, R. C., Ferrare, R. A., Hair, J. W., Hostetler, C. A., Liu, Y., Ortega, I., Sedlacek III, A., Shilling, J. E., Shrivastava, M., Springston, S. R., Tomlinson, J. M., Volkamer, R., Wilson, J., Zaveri, R. A., and Zelenyuk, A.: Model representations of aerosol layers transported from North America over the Atlantic Ocean during the Two-Column Aerosol Project, J. Geophys. Res.-Atmos., 121, 9814–9848, https://doi.org/10.1002/2016JD025248, 2016. a
Feng, L., Smith, S. J., Braun, C., Crippa, M., Gidden, M. J., Hoesly, R., Klimont, Z., van Marle, M., van den Berg, M., and van der Werf, G. R.: The generation of gridded emissions data for CMIP6, Geosci. Model Dev., 13, 461–482, https://doi.org/10.5194/gmd-13-461-2020, 2020. a, b
Gates, W. L., Boyle, J. S., Covey, C., Dease, C. G., Doutriaux, C. M., Drach, R. S., Fiorino, M., Gleckler, P. J., Hnilo, J. J., Marlais, S. M., Phillips, T. J., Potter, G. L., Santer, B. D., Sperber, K. R., Taylor, K. E., and Williams, D. N.: An Overview of the Results of the Atmospheric Model Intercomparison Project (AMIP I), B. Am. Meteorol. Soc., 80, 29–56, https://doi.org/10.1175/1520-0477(1999)080<0029:AOOTRO>2.0.CO;2, 1999. a
Gelaro, R., McCarty, W., Suárez, M. J., Todling, R., Molod, A., Takacs, L., Randles, C. A., Darmenov, A., Bosilovich, M. G., Reichle, R., Wargan, K., Coy, L., Cullather, R., Draper, C., Akella, S., Buchard, V., Conaty, A., da Silva, A. M., Gu, W., Kim, G.-K., Koster, R., Lucchesi, R., Merkova, D., Nielsen, J. E., Partyka, G., Pawson, S., Putman, W., Rienecker, M., Schubert, S. D., Sienkiewicz, M., and Zhao, B.: The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), J. Climate, 30, 5419–5454, https://doi.org/10.1175/JCLI-D-16-0758.1, 2017. a, b
Gettelman, A. and Morrison, H.: Advanced two-moment bulk microphysics for global models, Part I: Off-line tests and comparison with other schemes, J. Climate, 28, 1268–1287, https://doi.org/10.1175/JCLI-D-14-00102.1, 2015. a
Ghan, S., Wang, M., Zhang, S., Ferrachat, S., Gettelman, A., Griesfeller, J., Kipling, Z., Lohmann, U., Morrison, H., Neubauer, D., Partridge, D. G., Stier, P., Takemura, T., Wang, H., and Zhang, K.: Challenges in constraining anthropogenic aerosol effects on cloud radiative forcing using present-day spatiotemporal variability, P. Natl. Acad. Sci. USA, 113, 5804–5811, https://doi.org/10.1073/pnas.1514036113, 2016. a
Golaz, J.-C., Larson, V., and Cotton, W.: A PDF-Based Model for Boundary Layer Clouds. Part I: Method and Model Description, J. Atmos. Sci., 59, 3540–3551, https://doi.org/10.1175/1520-0469(2002)059<3540:APBMFB>2.0.CO;2, 2002. a
Golaz, J.-C., Caldwell, P. M., Van Roekel, L. P., Petersen, M. R., Tang, Q., Wolfe, J. D., Abeshu, G., Anantharaj, V., Asay-Davis, X. S., Bader, D. C., Baldwin, S. A., Bisht, G., Bogenschutz, P. A., Branstetter, M., Brunke, M. A., Brus, S. R., Burrows, S. M., Cameron-Smith, P. J., Donahue, A. S., Deakin, M., Easter, R. C., Evans, K. J., Feng, Y., Flanner, M., Foucar, J. G., Fyke, J. G., Griffin, B. M., Hannay, C., Harrop, B. E., Hoffman, M. J., Hunke, E. C., Jacob, R. L., Jacobsen, D. W., Jeffery, N., Jones, P. W., Keen, N. D., Klein, S. A., Larson, V. E., Leung, L. R., Li, H.-Y., Lin, W., Lipscomb, W. H., Ma, P.-L., Mahajan, S., Maltrud, M. E., Mametjanov, A., McClean, J. L., McCoy, R. B., Neale, R. B., Price, S. F., Qian, Y., Rasch, P. J., Reeves Eyre, J. E. J., Riley, W. J., Ringler, T. D., Roberts, A. F., Roesler, E. L., Salinger, A. G., Shaheen, Z., Shi, X., Singh, B., Tang, J., Taylor, M. A., Thornton, P. E., Turner, A. K., Veneziani, M., Wan, H., Wang, H., Wang, S., Williams, D. N., Wolfram, P. J., Worley, P. H., Xie, S., Yang, Y., Yoon, J.-H., Zelinka, M. D., Zender, C. S., Zeng, X., Zhang, C., Zhang, K., Zhang, Y., Zheng, X., Zhou, T., and Zhu, Q.: The DOE E3SM Coupled Model Version 1: Overview and Evaluation at Standard Resolution, J. Adv. Model. Earth Sy., 11, 2089–2129, https://doi.org/10.1029/2018MS001603, 2019. a, b, c
Hannah, W. M., Bradley, A. M., Guba, O., Tang, Q., Golaz, J.-C., and Wolfe, J.: Separating Physics and Dynamics Grids for Improved Computational Efficiency in Spectral Element Earth System Models, J. Adv. Model. Earth Sy., 13, e2020MS002419, https://doi.org/10.1029/2020MS002419, 2021. a
Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñz Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., De Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049, https://doi.org/10.1002/qj.3803, 2020. a, b, c, d
Hoesly, R. M., Smith, S. J., Feng, L., Klimont, Z., Janssens-Maenhout, G., Pitkanen, T., Seibert, J. J., Vu, L., Andres, R. J., Bolt, R. M., Bond, T. C., Dawidowski, L., Kholod, N., Kurokawa, J.-I., Li, M., Liu, L., Lu, Z., Moura, M. C. P., O'Rourke, P. R., and Zhang, Q.: Historical (1750–2014) anthropogenic emissions of reactive gases and aerosols from the Community Emissions Data System (CEDS), Geosci. Model Dev., 11, 369–408, https://doi.org/10.5194/gmd-11-369-2018, 2018. a, b
Huffman, G., Bolvin, D., Braithwaite, D., Hsu, K., Joyce, R., and Xie, P.: Integrated Multi-satellitE Retrievals for GPM (IMERG), version 4.4, NASA’s Precipitation Processing Center [data set], https://doi.org/10.5067/TRMM/TMPA/3H/7, 2014. a
Huffman, G. J. and Bolvin, D. T.: TRMM and other data precipitation data set documentation, NASA Global Change Master Directory Doc., NASA, Greenbelt, USA, 44 pp., ftp://precip.gsfc.nasa.gov/pub/trmmdocs/3B42_3B43_doc.pdf (last access: 18 March 2021), 2015. a
Huffman, G. J., Bolvin, D. T., Nelkin, E. J., Wolff, D. B., Adler, R. F., Gu, G., Hong, Y., Bowman, K. P., and Stocker, E. F.: The TRMM multisatellite precipitation analysis (TMPA): Quasi-global, multiyear, combined-sensor precipitation estimates at fine scales, J. Hydrometeorol., 8, 38–55, 2007. a
Iacono, M. J., Delamere, J. S., Mlawer, E. J., Shephard, M. W., Clough, S. A., and Collins, W. D.: Radiative forcing by long–lived greenhouse gases: Calculations with the AER radiative transfer models, J. Geophys. Res.-Atmos., 113, D13103, https://doi.org/10.1029/2008JD009944, 2008. a
Jeuken, A. B. M., Siegmund, P. C., Heijboer, L. C., Feichter, J., and Bengtsson, L.: On the potential of assimilating meteorological analyses in a global climate model for the purpose of model validation, J. Geophys. Res.-Atmos., 101, 16939–16950, https://doi.org/10.1029/96JD01218, 1996. a, b, c
Kanamitsu, M., Ebisuzaki, W., Woollen, J., Yang, S.-K., Hnilo, J. J., Fiorino, M., and Potter, G. L.: NCEP–DOE AMIP-II Reanalysis (R-2), B. Am. Meteorol. Soc., 83, 1631–1644, https://doi.org/10.1175/BAMS-83-11-1631, 2002. a
Kooperman, G. J., Pritchard, M. S., Ghan, S. J., Wang, M., Somerville, R. C. J., and Russell, L. M.: Constraining the influence of natural variability to improve estimates of global aerosol indirect effects in a nudged version of the Community Atmosphere Model 5, J. Geophys. Res.-Atmos., 117, D23204, https://doi.org/10.1029/2012JD018588, 2012. a, b, c, d
Larson, V. E., Golaz, J.-C., and Cotton, W. R.: A PDF-Based Model for Boundary Layer Clouds. Part I: Method and Model Description, J. Atmos. Sci., 59, 3519–3539, https://doi.org/10.1175/1520-0469(2002)059<3519:SSAMVI>2.0.CO;2, 2002. a
Lee, H.-T. and NOAA CDR Program: NOAA Climate Data Record (CDR) of Daily Outgoing Longwave Radiation (OLR), Version 1.2, NOAA National Climatic Data Center [data set], https://doi.org/10.7289/V5SJ1HH2, 2011. a
Lee, H.-T., Gruber, A., Ellingson, R. G., and Laszlo, I.: Development of the HIRS Outgoing Longwave Radiation Climate Dataset, J. Atmos. Ocean. Tech., 24, 2029–2047, https://doi.org/10.1175/2007JTECHA989.1, 2007. a
Lin, G., Wan, H., Zhang, K., Qian, Y., and Ghan, S. J.: Can nudging be used to quantify model sensitivities in precipitation and cloud forcing?, J. Adv. Model. Earth Sy., 8, 1073–1091, https://doi.org/10.1002/2016MS000659, 2016. a, b
Liu, X., Ma, P.-L., Wang, H., Tilmes, S., Singh, B., Easter, R. C., Ghan, S. J., and Rasch, P. J.: Description and evaluation of a new four-mode version of the Modal Aerosol Module (MAM4) within version 5.3 of the Community Atmosphere Model, Geosci. Model Dev., 9, 505–522, https://doi.org/10.5194/gmd-9-505-2016, 2016. a
Liu, Y., Zhang, K., Qian, Y., Wang, Y., Zou, Y., Song, Y., Wan, H., Liu, X., and Yang, X.-Q.: Investigation of short-term effective radiative forcing of fire aerosols over North America using nudged hindcast ensembles, Atmos. Chem. Phys., 18, 31–47, https://doi.org/10.5194/acp-18-31-2018, 2018. a, b
Ma, H.-Y., Chuang, C. C., Klein, S. A., Lo, M.-H., Zhang, Y., Xie, S., Zheng, X., Ma, P.-L., Zhang, Y., and Phillips, T. J.: An improved hindcast approach for evaluation and diagnosis of physical processes in global climate models, J. Adv. Model. Earth Sy., 7, 1810–1827, https://doi.org/10.1002/2015MS000490, 2015. a
Ma, P.-L., Rasch, P. J., Fast, J. D., Easter, R. C., Gustafson Jr., W. I., Liu, X., Ghan, S. J., and Singh, B.: Assessing the CAM5 physics suite in the WRF-Chem model: implementation, resolution sensitivity, and a first evaluation for a regional case study, Geosci. Model Dev., 7, 755–778, https://doi.org/10.5194/gmd-7-755-2014, 2014. a
Milrad, S.: 4 – Upper-Air Observations, in: Synoptic analysis and forecasting: An introductory toolkit, edited by: Milrad, S., Elsevier, Amsterdam, Netherlands, 41–47, https://doi.org/10.1016/B978-0-12-809247-7.00004-1, 2017. a
Mlawer, E. J., Taubman, S. J., Brown, P. D., Iacono, M. J., and Clough, S. A.: Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave, J. Geophys. Res.-Atmos., 102, 16663–16682, https://doi.org/10.1029/97JD00237, 1997. a
Morrison, H. and Gettelman, A.: A New Two-Moment Bulk Stratiform Cloud Microphysics Scheme in the Community Atmosphere Model, Version 3 (CAM3). Part I: Description and Numerical Tests, J. Climate, 21, 3642–3659, https://doi.org/10.1175/2008JCLI2105.1, 2008. a
Myhre, G., Shindell, D., Bréon, F.-M., Collins, W., Fuglestvedt, J., Huang, J., Koch, D., Lamarque, J.-F., Lee, D., Mendoza, B., Nakajima, T., Robock, A., Stephens, G., Takemura, T., and Zhang, H.: Anthropogenic and natural radiative forcing, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Doschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, UK, 659–740, https://doi.org/10.1017/CBO9781107415324.018, 2013. a, b
NOAA PSL: NOAA Interpolated Outgoing Longwave Radiation (OLR), NOAA National Climatic Data Center [data set], https://psl.noaa.gov/data/gridded/data.olrcdr.interp.html (last access: 29 August 2022), 2001. a
Oleson, K. W., Lawrence, M., Bonan, B., Drewniak, B. A., Huang, M., Koven, D., Levis, S., Li, F., Riley, J. P., Subin, M. C., Swenson, S. C., Thornton, E., Bozbiyik, A., Fisher, R. A., Heald, L., Kluzek, E., Lamarque, J., Lawrence, J., Leung, R., Lipscomb, W. H., Muszala, P., Ricciuto, M., Sacks, J., Sun, Y., Tang, J., and Yang, Z.: Technical description of version 4.5 of the Community Land Model (CLM), NCAR Technical Note NCAR/TN-503+STR, NCAR, USA, 420 pp., https://doi.org/10.5065/D6RR1W7M, 2013. a
Rasch, P. J., Xie, S., Ma, P.-L., Lin, W., Wang, H., Tang, Q., Burrows, S. M., Caldwell, P., Zhang, K., Easter, R. C., Cameron-Smith, P., Singh, B., Wan, H., Golaz, J.-C., Harrop, B. E., Roesler, E., Bacmeister, J., Larson, V. E., Evans, K. J., Qian, Y., Taylor, M., Leung, L. R., Zhang, Y., Brent, L., Branstetter, M., Hannay, C., Mahajan, S., Mametjanov, A., Neale, R., Richter, J. H., Yoon, J.-H., Zender, C. S., Bader, D., Flanner, M., Foucar, J. G., Jacob, R., Keen, N., Klein, S. A., Liu, X., Salinger, A., Shrivastava, M., and Yang, Y.: An Overview of the Atmospheric Component of the Energy Exascale Earth System Model, J. Adv. Model. Earth Sy., 11, 2377–2411, https://doi.org/10.1029/2019MS001629, 2019. a, b, c
Reynolds, R. W., Rayner, N. A., Smith, T. M., Stokes, D. C., and Wang, W.: An Improved In Situ and Satellite SST Analysis for Climate, J. Climate, 15, 1609–1625, https://doi.org/10.1175/1520-0442(2002)015<1609:AIISAS>2.0.CO;2, 2002. a
Separovic, L., de Elía, R., and Laprise, R.: Impact of spectral nudging and domain size in studies of RCM response to parameter modification, Clim. Dynam., 38, 1325–1343, https://doi.org/10.1007/s00382-011-1072-7, 2012. a
Sorooshian, S., Hsu, K., Braithwaite, D., Ashouri, H., and NOAA CDR Program: NOAA Climate Data Record (CDR) of Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN-CDR), Version 1 Revision 1, NOAA National Centers for Environmental Information [data set], https://doi.org/10.7289/V51V5BWQ, 2014. a
Stowe, L. L., Jacobowitz, H., Ohring, G., Knapp, K. R., and Nalli, N. R.: The Advanced Very High Resolution Radiometer (AVHRR) Pathfinder Atmosphere (PATMOS) Climate Dataset: Initial Analyses and Evaluations, J. Climate, 15, 1243–1260, https://doi.org/10.1175/1520-0442(2002)015<1243:TAVHRR>2.0.CO;2, 2002. a
Subramanian, A. C. and Zhang, G. J.: Diagnosing MJO hindcast biases in NCAR CAM3 using nudging during the DYNAMO field campaign, J. Geophys. Res.-Atmos., 119, 7231–7253, https://doi.org/10.1002/2013JD021370, 2014. a
Sun, J., Zhang, K., Wan, H., Ma, P.-L., Tang, Q., and Zhang, S.: Impact of Nudging Strategy on the Climate Representativeness and Hindcast Skill of Constrained EAMv1 Simulations, J. Adv. Model. Earth Sy., 11, 3911–3933, https://doi.org/10.1029/2019MS001831, 2019. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q
Telford, P. J., Braesicke, P., Morgenstern, O., and Pyle, J. A.: Technical Note: Description and assessment of a nudged version of the new dynamics Unified Model, Atmos. Chem. Phys., 8, 1701–1712, https://doi.org/10.5194/acp-8-1701-2008, 2008. a
Wan, H., Rasch, P. J., Zhang, K., Qian, Y., Yan, H., and Zhao, C.: Short ensembles: an efficient method for discerning climate-relevant sensitivities in atmospheric general circulation models, Geosci. Model Dev., 7, 1961–1977, https://doi.org/10.5194/gmd-7-1961-2014, 2014. a
Wan, H., Zhang, S., Rasch, P. J., Larson, V. E., Zeng, X., and Yan, H.: Quantifying and attributing time step sensitivities in present-day climate simulations conducted with EAMv1, Geosci. Model Dev., 14, 1921–1948, https://doi.org/10.5194/gmd-14-1921-2021, 2021. a, b
Wan, H., Zhang, K., Rasch, P. J., Larson, V. E., Zeng, X., Zhang, S., and Dixon, R.: CondiDiag1.0: a flexible online diagnostic tool for conditional sampling and budget analysis in the E3SM atmosphere model (EAM), Geosci. Model Dev., 15, 3205–3231, https://doi.org/10.5194/gmd-15-3205-2022, 2022. a
Wang, H., Easter, R. C., Zhang, R., Ma, P.-L., Singh, B., Zhang, K., Ganguly, D., Rasch, P. J., Burrows, S. M., Ghan, S. J., Lou, S., Qian, Y., Yang, Y., Feng, Y., Flanner, M., Leung, L. R., Liu, X., Shrivastava, M., Sun, J., Tang, Q., Xie, S., and Yoon, J.-H.: Aerosols in the E3SM Version 1: New Developments and Their Impacts on Radiative Forcing, J. Adv. Model. Earth Sy., 12, e2019MS001851, https://doi.org/10.1029/2019MS001851, 2020. a
Wang, Y., Liu, X., Hoose, C., and Wang, B.: Different contact angle distributions for heterogeneous ice nucleation in the Community Atmospheric Model version 5, Atmos. Chem. Phys., 14, 10411–10430, https://doi.org/10.5194/acp-14-10411-2014, 2014. a
Xie, S., Lin, W., Rasch, P. J., Ma, P.-L., Neale, R., Larson, V. E., Qian, Y., Bogenschutz, P. A., Caldwell, P., Cameron-Smith, P., Golaz, J.-C., Mahajan, S., Singh, B., Tang, Q., Wang, H., Yoon, J.-H., Zhang, K., and Zhang, Y.: Understanding Cloud and Convective Characteristics in Version 1 of the E3SM Atmosphere Model, J. Adv. Model. Earth Sy., 10, 2618–2644, https://doi.org/10.1029/2018MS001350, 2018. a, b, c
Zhang, G. J. and McFarlane, N. A.: Sensitivity of climate simulations to the parameterization of cumulus convection in the Canadian Climate Centre general circulation model, Atmos. Ocean, 33, 407–446, https://doi.org/10.1080/07055900.1995.9649539, 1995. a
Zhang, K., O'Donnell, D., Kazil, J., Stier, P., Kinne, S., Lohmann, U., Ferrachat, S., Croft, B., Quaas, J., Wan, H., Rast, S., and Feichter, J.: The global aerosol-climate model ECHAM-HAM, version 2: sensitivity to improvements in process representations, Atmos. Chem. Phys., 12, 8911–8949, https://doi.org/10.5194/acp-12-8911-2012, 2012. a
Zhang, K., Wan, H., Liu, X., Ghan, S. J., Kooperman, G. J., Ma, P.-L., Rasch, P. J., Neubauer, D., and Lohmann, U.: Technical Note: On the use of nudging for aerosol–climate model intercomparison studies, Atmos. Chem. Phys., 14, 8631–8645, https://doi.org/10.5194/acp-14-8631-2014, 2014. a, b, c, d, e, f, g
Zhang, K., Rasch, P. J., Taylor, M. A., Wan, H., Leung, R., Ma, P.-L., Golaz, J.-C., Wolfe, J., Lin, W., Singh, B., Burrows, S., Yoon, J.-H., Wang, H., Qian, Y., Tang, Q., Caldwell, P., and Xie, S.: Impact of numerical choices on water conservation in the E3SM Atmosphere Model version 1 (EAMv1), Geosci. Model Dev., 11, 1971–1988, https://doi.org/10.5194/gmd-11-1971-2018, 2018. a
Zhang, K., Zhang, W., Wan, H., Rasch, P. J., Ghan, S. J., Easter, R. C., Shi, X., Wang, Y., Wang, H., Ma, P.-L., Zhang, S., Sun, J., Burrows, S. M., Shrivastava, M., Singh, B., Qian, Y., Liu, X., Golaz, J.-C., Tang, Q., Zheng, X., Xie, S., Lin, W., Feng, Y., Wang, M., Yoon, J.-H., and Leung, L. R.: Effective radiative forcing of anthropogenic aerosols in E3SM version 1: historical changes, causality, decomposition, and parameterization sensitivities, Atmos. Chem. Phys., 22, 9129–9160, https://doi.org/10.5194/acp-22-9129-2022, 2022. a
Zhang, S., Wang, M., Ghan, S. J., Ding, A., Wang, H., Zhang, K., Neubauer, D., Lohmann, U., Ferrachat, S., Takeamura, T., Gettelman, A., Morrison, H., Lee, Y., Shindell, D. T., Partridge, D. G., Stier, P., Kipling, Z., and Fu, C.: On the characteristics of aerosol indirect effect based on dynamic regimes in global climate models, Atmos. Chem. Phys., 16, 2765–2783, https://doi.org/10.5194/acp-16-2765-2016, 2016. a