Community Atmosphere Model Version Research Articles

A Method for Interpreting the Role of Parameterized Turbulence on Global Metrics in the Community Earth System Model

AbstractThe parameterization of subgrid‐scale processes such as boundary layer (PBL) turbulence introduces uncertainty in Earth System Model (ESM) results. This uncertainty can contribute to or exacerbate existing biases in representing key physical processes. This study analyzes the influence of tunable parameters in an experimental version of the Cloud Layers Unified by Binormals (CLUBBX) scheme. CLUBB is the operational PBL parameterization in the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community ESM version 2 (CESM2). We perform the Morris one‐at‐a‐time (MOAT) parameter sensitivity analysis using short‐term (3‐day), initialized hindcasts of CAM6‐CLUBBX with 24 unique initial conditions. Several input parameters modulating vertical momentum flux appear most influential for various regionally‐averaged quantities, namely surface stress and shortwave cloud forcing (SWCF). These parameter sensitivities have a spatial dependence, with parameters governing momentum flux most influential in regions of high vertical wind shear (e.g., the mid‐latitude storm tracks). We next evaluate several experimental 20‐year simulations of CAM6‐CLUBBX with targeted parameter perturbations. We find that parameter perturbations produce similar physical mechanisms in both short‐term and long‐term simulations, but these physical responses can be muted due to nonlinear feedbacks manifesting over time scales longer than 3 days, thus causing differences in how output metrics respond in the long‐term simulations. Analysis of turbulent fluxes in CLUBBX indicates that the influential parameters affect vertical fluxes of heat, moisture, and momentum, providing physical pathways for the sensitivities identified in this study.

Response of Convectively Coupled Kelvin Waves to Surface Temperature Forcing in Aquaplanet Simulations

AbstractThis study investigates changes in the propagation and maintenance of convectively coupled Kelvin waves (KWs) in response to surface warming. We use a set of three aquaplanet simulations made with the Community Atmospheric Model version 6 by varying the sea surface temperature boundary conditions to represent the current climate as well as warmer (+4 K) and cooler (−4 K) climates. Results show that KWs accelerate at the rate of about 7.1%/K and their amplitudes decrease by 4.7%/K. The dampening of KWs with warming is found to be associated with a weakening of the internal thermodynamic feedback between diabatic heating and temperature anomalies that generates KW eddy available potential energy (EAPE). The phase speed of KWs closely matches that of the second baroclinic mode KW in −4 K, while the phase speed of KWs is approximately that of the first baroclinic mode KW in +4 K. Meanwhile, the coupling between the two baroclinic modes weakens with warming. We hypothesize that in −4 K, as the first and second baroclinic modes are strongly coupled, KWs destabilize by positive EAPE generation within the second baroclinic mode and propagate more slowly, following the second baroclinic mode KW phase speed. In +4 K, as the first and second baroclinic modes decouple, KWs are damped by negative EAPE generation within the first baroclinic mode and propagate faster, following the first baroclinic mode KW phase speed.

The Interaction Between Climate Forcing and Feedbacks

AbstractA Perturbed Parameter Ensemble (PPE) with the Community Atmosphere Model version 6 (CAM6) is used to better understand the sensitivity of aerosol forcing and cloud feedbacks to changes in model processes. Aerosol forcing through aerosol‐cloud interactions is mostly negative (a cooling) due to shortwave radiation, while feedbacks are positive or negative in different regions due to contrasting longwave and shortwave effects. Both forcing and feedbacks are related to the mean climate state. Higher magnitude cloud radiative effects generally mean larger magnitude net negative forcing and larger magnitude net positive feedback. Aerosol forcing is broadly related to the susceptibility of clouds to drop number. Feedbacks also related to susceptibility, but to a lesser extent and in different regions to aerosol forcing. Aerosol forcing and cloud feedbacks are anti‐correlated in the CAM6 PPE such that stronger negative forcing is associated with stronger positive feedbacks. Even the processes governing forcing and feedback sensitivity in the PPE are similar. These include the warm rain formation process, ice loss processes and deep convective intensity.

Dominant modes of interannual variability in spring compound dry and hot events over Northern Asia and the possible mechanisms

In this study, spatial and temporal variations of spring compound dry and hot events (CDHEs) over northern Asia (NA) during 1950–2020 and the related possible mechanisms are investigated on the interannual timescale. The standardized compound event indicator (SCEI) is used to represent compound dry and hot conditions over NA, which is validated by the observed summer case of CDHEs over western Russia in 2010. The first empirical orthogonal function (EOF1) mode of the SCEI over NA presents a monopole pattern, while the EOF2 mode shows an east-west dipole pattern. Possible mechanism analysis indicates that the CDHEs tend to be modulated by local high-pressure anomalies, accompanied by reduced cloud cover and more solar radiation. The high-pressure anomalies can lead to warm temperature and water vapor divergence, which synergistically contributes to the occurrence of CDHEs. Further analysis shows that the EOF1 mode is jointly affected by the Scandinavian (SCAND) pattern and the Arctic Oscillation (AO); while the EOF2 mode is influenced by the Atlantic-Eurasian (AEA) teleconnection. Moreover, the North Atlantic tripolar pattern of sea surface temperature (SST) can influence the EOF1 mode through changing the AO. And the diagonal tripolar SST pattern over the North Atlantic affects the EOF2 mode via altering the AEA pattern. The influence of the SST patterns is further confirmed by numerical experiments using the Community Atmospheric Model version CAM-5.3.

Multidecadal variations in North Atlantic SSTs modulate the relationship between ENSO and the South Atlantic Subtropical Dipole since 1900

Abstract This study investigates the long-term variability of the relationship between the El Niño South/Oscillation (ENSO) and the South Atlantic Subtropical Dipole (SASD), and their connection to multidecadal variations in the North Atlantic sea surface temperature (MDV-SST). Using a century’s worth of SST and atmospheric data and the Community Atmosphere Model version 5, the study found significant interdecadal variability in the correlation between the Niño3.4 index and the SASD index. This variability is closely linked to the North Atlantic MDV-SST, often interpreted as Atlantic Multidecadal Oscillation (AMO). This is demonstrated by significant ENSO-SASD correlations during warm (positive) phases of MDV-SST (1930–1960 and 2001–2020) and insignificant correlations during cold (negative) phases (1900–1929 and 1961–2000). Through a Gill-type response, MDV-SST excites Rossby wave over the tropical Pacific Ocean, influencing surface wind, SST, convective activities and upper-level zonal wind over the Pacific Ocean. During warm MDV-SST phases, the more eastward positioned Rossby wave source (RWS), triggered by SST and precipitation anomalies over the South Pacific Ocean, along with a stronger, more northward subtropical jet stream, propagates the wavetrain more eastwards into the South Atlantic Ocean, thereby strengthening the SST anomalies in the SASD. Conversely, during the cold MDV-SST phases, the more westward-positioned ENSO-related RWS and a stronger mid-latitude jet stream guide the wavetrain southeastwards into the southeastern Pacific Ocean, exerting less influence on the SST anomalies in the SASD. The Ekman pumping caused by anomalous surface pressure and the associated surface wind field as well as surface turbulent heat flux also affect the SST anomalies in the SASD and the ENSO-SASD relation.

Wave‐Convection Interactions Amplify Convective Parameterization Biases in the South Pacific Convergence Zone

AbstractClimate models have long‐standing difficulties simulating the South Pacific Convergence Zone (SPCZ) and its variability. For example, the default Zhang‐McFarlane (ZM) convection scheme in the Community Atmosphere Model version 5 (CAM5) produces too much light precipitation and too little heavy precipitation in the SPCZ, with this bias toward light precipitation even more pronounced in the SPCZ than in the tropics as a whole. Here, we show that implementing a recently developed convection scheme in the CAM5 yields significant improvements in the simulated SPCZ during austral summer and discuss the reasons behind these improvements. In addition to intensifying both mean rainfall and its variability in the SPCZ, the new scheme produces a larger heavy rainfall fraction that is more consistent with observations and state‐of‐the‐art reanalyses. This shift toward heavier, more variable rainfall increases both the magnitude and altitude of diabatic heating associated with convective precipitation, intensifying lower tropospheric convergence and increasing the influence of convection on the upper‐level circulation. Increased diabatic production of potential vorticity in the upper troposphere intensifies the distortion effect exerted by convection on transient Rossby waves that pass through the SPCZ. Weaker distortion effects in simulations using the ZM scheme allow waves to propagate continuously through the region rather than dissipating locally, further reducing updrafts and weakening convection in the SPCZ. Our results outline a dynamical framework for evaluating model representations of tropical–extratropical interactions within the SPCZ and clarify why convective parameterizations that produce “top‐heavy” profiles of deep convective heating better represent the SPCZ and its variability.

Enhanced Aging of Black Carbon under Recent Clean Air Actions and Future Carbon Neutrality Scenario in China.

China has implemented strict emission control measures, but it is unclear how they affect black carbon (BC) aging and light absorption. Here, we use the Community Atmosphere Model version 6 (CAM6) with the four-mode version of the Modal Aerosol Module coupled with machine learning (MAM4-ML) to simulate BC aging during 2011-2018 and 2050/2100 following a carbon neutrality scenario (SSP1-2.6), respectively. During 2011-2018, the mass ratio of coatings to BC (RBC) widely increased (5.4% yr-1) over the east of China. The increased secondary organic aerosol (SOA) coatings dominate (88%) the increased RBC, while the sulfate coatings decrease. The drivers of BC coating changes come from the different magnitudes of emission reductions in secondary aerosol precursors (i.e., volatile organic compounds (VOCs) and SO2) and BC. During 2011-2018, the increased RBC enhances the BC mass absorption cross section (MAC, 0.7% yr-1). In 2050/2100 for SSP1-2.6, emission control leads to further increased RBC (95/145%) and BC MAC (12/17%). For both 2011-2018 and 2050/2100, the enhanced BC MAC partly offsets the declining direct radiative effect (DRE) of BC due to direct emission reduction. As a result, the full impact of direct emission reductions of BC on BC DRE is only 75% for 2011-2018 and 90/94% for 2050/2100.

Diagnosing uncertainties in global biomass burning emission inventories and their impact on modeled air pollutants

Abstract. Large uncertainties persist within current biomass burning (BB) inventories, and the choice of these inventories can substantially impact model results when assessing the influence of BB aerosols on weather and climate. We evaluated discrepancies among BB emission inventories by comparing carbon monoxide (CO) and organic carbon (OC) emissions from seven major BB regions globally between 2013 and 2016. Mainstream bottom-up inventories, including the Fire INventory from NCAR 1.5 (FINN1.5) and Global Fire Emissions Database version 4s (GFED4s), along with the top-down inventories Quick Fire Emissions Dataset 2.5 (QFED2.5) and Visible Infrared Imaging Radiometer (VIIRS-)based Fire Emission Inventory version 0 (VFEI0), were selected for this study. Global CO emissions range from 252 to 336 Tg, with regional disparities reaching up to a 6-fold difference. Dry matter is the primary contributor to the regional variation in CO emissions (50 %–80 %), with emission factors accounting for the remaining 20 %–50 %. Uncertainties in dry matter often arise from biases in calculating bottom fuel consumption and burned area, influenced by vegetation classification methods and fire detection products. In the tropics, peatlands contribute more fuel loads and higher emission factors than grasslands. At high latitudes, increased cloud fraction amplifies the discrepancy in estimated burned area (or fire radiative power) by 20 %. The global OC emissions range from 14.9 to 42.9 Tg, exhibiting higher variability than CO emissions due to the corrected emission factors in QFED2.5, with regional disparities reaching a factor of 8.7. Additionally, we applied these BB emission inventories to the Community Atmosphere Model version 6 (CAM6) and assessed the model performance against observations. Our results suggest that the simulations based on the GFED4s agree best with the Measurements of Pollution in the Troposphere (MOPITT-)retrieved CO. While comparing the simulation with the Moderate Resolution Imaging Spectroradiometer (MODIS) and AErosol RObotic NETwork (AERONET) aerosol optical depth (AOD), our results reveal that there is no global optimal choice for BB inventories. In the high latitudes of the Northern Hemisphere, using GFED4s and QFED2.5 can better capture the AOD magnitude and diurnal variation. In equatorial Asia, GFED4s outperforms other models in representing day-to-day changes, particularly during intense burning. In Southeast Asia, we recommend using the OC emission magnitude from FINN1.5 combined with daily variability from QFED2.5. In the Southern Hemisphere, the latest VFEI0 has performed relatively well. This study has implications for reducing the uncertainties in emissions and improving BB emission inventories in further studies.

Evaluating nine different air-sea flux algorithms coupled with CAM6

The accuracy of air-sea flux algorithms significantly impacts the simulation performance of climate models. In this study, we evaluated nine widely-used air-sea flux algorithms: those used in the BCC model, the CAM model, the ECMWF model, and the WRF model (with three options), as well as three other commonly-used algorithms: COARE, UA, and VM. An offline comparison revealed significant discrepancies in the flux calculations obtained with different algorithms. Subsequently, we conducted nine sensitivity experiments by coupling these algorithms into the Community Atmosphere Model version 6 (CAM6) to assess the extent and manner in which the selection of flux algorithms influenced the simulated climate. The online simulations introduced complexities due to interactions between the atmosphere and surface fluxes, making the results more intricate. The iterative process of the model created negative feedback between wind speed and momentum flux, as well as between air temperature and heat flux, leading to inconsistencies between online and offline test results. Based on the comprehensive evaluation metric Distance between Indices of Simulation and Observation (DISO), a concept of combined or segmented algorithm was introduced, and a segmented algorithm based on wind speed was found resulting in an overall superior simulation performance compared to the previous nine experiments.

Connection between Barents Sea Ice in May and Early Summer Monsoon Rainfall in the South China Sea and Its Possible Mechanism

The impacts of Arctic sea ice on climate in middle and high latitudes have been extensively studied. However, its effects on climate in low latitudes, particularly on summer monsoon rainfall in the South China Sea (SCS), have received limited attention. Thus, this study investigates the connection between the Arctic sea ice concentration (SIC) anomaly and the early summer monsoon rainfall (ESMR) in the SCS and its underlying physical mechanism. The results reveal a significant positive correlation between the Barents Sea (BS) SIC in May and the ESMR in the SCS. When there is more (less) SIC in the Barents Sea (BS) during May, this results in a positive (negative) anomaly of the local turbulent heat flux, which lasts until June. This, in turn, excites an upward (downward) air motion anomaly in the vicinity of the BS, causing a corresponding downward (upward) motion anomaly over the Black Sea. Consequently, this triggers a wave train similar to the Eurasian (SEU) teleconnection, propagating eastward towards East Asia. The SEU further leads to an (a) upward (downward) motion anomaly and weakens (strengthens) the western Pacific subtropical high (WPSH) over the SCS, which is accompanied by a southwest adequate (scarce) water vapor anomaly transporting from the Indian Ocean, resulting in more (less) precipitation in the SCS. Furthermore, the response of ESMR in the SCS to the SIC in the BS is further verified by using the Community Atmosphere Model version 5.3 (CAM5.3). This study introduces novel precursor factors that influence the South China Sea summer monsoon (SCSSM), presenting a new insight for climate prediction in this region, which holds significant implications.

Exploring covariabilities of the high-summer subtropical upper-level pressure systems in the Northern Hemisphere

The Northern Hemisphere subtropical upper-tropospheric pressure systems, i.e., South Asian high (SAH), mid-Pacific Trough (MPT), and mid-Atlantic Trough (MAT), play key roles in regulating the Eurasian–Pacific–North American climate, but their covariabilities and drivers remain unclear. Here, we reveal two seesaw-like teleconnections between MPT and western SAH (WSAH) and between MPT and MAT during the boreal high summer, respectively, and explore their associations with the tropical Pacific and Atlantic sea surface temperature (SST) anomalies. The MPT-WSAH teleconnection, featured by coherent variations between WSAH and MPT, is coupled with enhanced western Pacific subtropical high and Asian mid-latitude low in the low-level troposphere. In contrast, the MPT-MAT teleconnection features an out-of-phase relationship between MPT and MAT, which is associated with westward expanded western Pacific subtropical high but weakened North Atlantic subtropical high. Both two teleconnections are closely linked to Pacific SST anomalies. However, the former is mainly driven by the hybrid La Niña in the tropical Pacific, where cold SST anomalies spreading eastward since preceding winter. Whereas the latter corresponds to the central Pacific La Niña since late-spring mainly driven by preceding successive SST warming in the tropical North Atlantic. Above conclusions are further confirmed by the simulated results from the numerical experiments, which are performed using the US National Center for Atmospheric Research global Community Atmosphere Model version 4.

Responses of Westerly Jets over Asia versus North America in the Last Glacial Maximum and Their Attributions

Abstract By using the Community Atmosphere Model version 4 (CAM4) coupled to a slab ocean model, the total responses of westerly jets (WJ) over Asia and North America during the Last Glacial Maximum (LGM), as well as the relative influences of different individual LGM forcing on the changes in WJ, are evaluated in this study. The results show that the seasonal variation of Northern Hemisphere WJ in LGM is characterized by distinct regional discrepancies. The WJ is intensified and shifts southward in North America, while it generally weakens and slightly displaces northward in central Asia. Over Japan, the WJ is enhanced and moves southward in summer but is attenuated and shifts northward in winter. The change of WJ has a close relationship with the midtroposphere temperature anomalies which result from the different external forcings. In summer, the ice sheet albedo plays a leading role in wind fields at mid-to-high latitudes of the Northern Hemisphere by causing a zonal wave–like response of temperature. In winter, the ice sheet topography dominates the wind fields across North America to North Atlantic by inducing a northwest–southeast-oriented tripole temperature response. The effects of orbital parameters and greenhouse gases largely contribute to the alteration of upper wind fields over Asia and North Pacific by causing a dipole response of temperature.

Improving BC Mixing State and CCN Activity Representation With Machine Learning in the Community Atmosphere Model Version 6 (CAM6)

AbstractRepresenting mixing state of black carbon (BC) is challenging for global climate models (GCMs). The Community Atmosphere Model version 6 (CAM6) with the four‐mode version of the Modal Aerosol Module (MAM4) represents aerosols as fully internal mixtures with uniform composition within each aerosol mode, resulting in high degree of internal mixing of BC with non‐BC species and large mass ratio of coating to BC (RBC, the mass ratio of non‐BC species to BC in BC‐containing particles). To improve BC mixing state representation, we coupled a machine learning (ML) model of BC mixing state index trained on particle‐resolved simulations to the CAM6 with MAM4 (MAM4‐ML). In MAM4‐ML, we use RBC to partition accumulation mode particles into two new modes, BC‐free particles and BC‐containing particles. We adjust RBC to make the modeled BC mixing state index (χmode) match the one predicted by the ML model (χML). On a global average, the mass fraction of BC‐containing particles in accumulation mode decreases from 100% (MAM4‐default) to 48% (MAM4‐ML). The globally averaged χmode decreases from 78% (MAM4‐default) to 63% (MAM4‐ML, 19% reduction) and agrees well with χML (66%). The RBC decreases by 52% for accumulation mode and better agrees with observations. The hygroscopicity drops by 9% for BC‐containing particles in accumulation mode, leading to a 20% reduction in the BC activation fraction. The surface BC concentration increases most (6.9%) in the Arctic, and the BC burden increases by 4%, globally. Our study highlights the application of the ML model for improving key aerosol processes in GCMs.

Sharp decline of dust events induces regional wetting over arid and semi-arid Northwest China in the NCAR Community atmosphere model

Multiple lines of observational evidence have indicated a significant wetting over the arid and semi-arid Northwest China (NWC) during recent decades, coinciding with a simultaneous sharp decline of dust events. Although recent studies have attributed NWC wetting to different anthropogenic and natural forcings, the mechanisms are not definitive and the regional wetting has been greatly underestimated in the Coupled Model Intercomparison Project historical simulations. Based on sensitivity experiments with different dust emission amounts using the NCAR Community Atmospheric Model version 5 (CAM5), here we find that decreasing dusts exert significant impacts on mixed-phase clouds through reducing the concentration of ice nucleating particles, increase the NWC precipitation and thus induce regional wetting through enhancing convection precipitation. A possible convection invigoration mechanism whereby the atmospheric vertical temperature gradient and convective instability are strengthened by reduced dusts, leading to convection invigoration and increased precipitation. These results are reinforced by simulations over the dust region in North Africa where mixed-phase and ice clouds are rare and reduced dusts do not increase precipitation. This study highlights the possible mechanism of dust-ice cloud interactions in recent NWC wetting and future regional climate change.

Characteristics and mechanisms of the severe compound cold-wet event in southern China during February 2022

A severe compound cold-wet event occurred in southern China (hereafter referred to as CWESC ) in February 2022, leading to enormous socioeconomic losses. In this study, we proposed a new index to denote the severity of the compound cold-wet event. Based on the multivariate survival method, the CWESC in February 2022 is identified as the severest event during the past six decades. Our results indicate that the CWESC in 2022 is jointly regulated by the La Niña-like SST condition in the tropical Pacific and the warm SST anomalies in the North Atlantic, and a teleconnection in the Northern Hemisphere during winter (hereafter referred to as TNHW) plays the key role. The TNHW pattern originates from the tropical Pacific, and it splits into two routes over the North Atlantic. The northern branch of TNHW propagates via the Arctic and Siberia, causing intensified near-surface northerly wind and partially inducing an anomalous anticyclone over the western North Pacific (WNP). The southern branch of TNHW propagates via the Mediterranean and western Asia, inducing a deepened India–Burma trough and partially inducing the anomalous anticyclone over WNP. The intensified near-surface northerly wind causes enhanced cold advection over southern China, while the deepened India–Burma trough and the anomalous anticyclone over WNP cause increased southerly warm and moist air flow towards southern China, resulting in the CWESC in 2022. Moreover, four groups of numerical experiments forced by tropical Pacific, North Pacific, and North Atlantic SST anomalies are conducted based on the Community Atmosphere Model version 5. The results confirm the important roles of the La Niña-like condition and the warm SST anomalies in the North Atlantic in causing the CWESC in 2022.

Extreme heat waves in June 2021 over Europe regulated by shrinking Eurasian snow cover in mid–high latitudes

Europe has experienced many extreme heat waves over the past few decades. In this study, the physical processes underlying these long-lasting and wide-ranging heat wave events are investigated based on a case study in Europe in June 2021. Heat waves are associated with barotropic anticyclonic anomalies accompanied by positive geopotential height anomalies locally. These anomalies persist under the conditions of increased meridional air temperature gradients of the mid–upper troposphere in the high latitudes of Eurasia and the formation of the Arctic front jet. The shrinking high-latitude snow cover in April–May favors higher surface temperatures and larger meridional temperature gradients in June in the mid–upper troposphere due to the soil moisture–evaporation–temperature positive feedback process. The summer Arctic front jet is then strengthened, and the mid-latitude westerly winds are weakened. This atmospheric circulation background favors waveguide formation and wave resonance that produces high-amplitude atmospheric waves and the stagnation of ridges in the mid-latitudes. Numerical experiments using the Community Atmosphere Model version 5 verify the proposed physical mechanisms, with the climatic responses in sensitivity experiments to anomalous snowfall rates closely resembling the observational results. Therefore, in June 2021, under the identified atmospheric circulation background and the perturbation of the upstream positive phase of the North Atlantic Oscillation, the large-scale barotropic high pressure and barotropic anticyclonic circulation in the study region tended to be stable and persistent, which is favorable for the production of long-lasting and wide-ranging heat wave events.

An Aerosol Optical Module With Observation‐Constrained Black Carbon Properties for Global Climate Models

AbstractAtmospheric black carbon (BC) aerosols have been long‐lasting uncertain components in environmental and climate studies. Global climate models (GCMs) potentially overestimate BC absorption efficiency due to a lack of consideration of complex BC microphysical and mixing properties. We extract multiple BC properties from observations and develop an aerosol optical module known as Advanced Black Carbon (ABC) in the framework of the Modal Aerosol Model version 4 (MAM4). The ABC module is implemented in the Community Atmosphere Model version 6 (CAM6) and evaluated by in situ and remote sensing observations. CAM6‐ABC addresses the shortcomings of CAM6‐MAM4 in terms of BC microphysical and mixing properties, particularly their size, mixing state and optical simulations. Sensitivity simulations show that the global BC absorption aerosol optical depth at 550 nm simulated by CAM6‐ABC is reduced by ∼29% compared with that in CAM6‐MAM4. The BC absorption enhancement simulated by CAM6‐ABC is reduced from ∼2.6 of the default MAM4 to ∼1.4, which is closer to the observed values (mostly less than 1.5). With improved BC absorption estimation, the biases of aerosol single‐scattering coalbedo simulations are reduced by 18%–69% compared with global Aerosol Robotic Network observations. Moreover, the globally averaged BC direct radiative effect is reduced from 0.37 to 0.28 W/m2 at the top of the atmosphere. Our new scheme alleviates the overestimation of BC absorption in GCMs by constraining BC microphysical and mixing properties when assessing aerosol radiative and climate effects, and it can be easily implemented in most modal‐based aerosol modules of climate models.

The Role of In‐Cloud Wet Removal in Simulating Aerosol Vertical Profiles and Cloud Radiative Forcing

AbstractAmong the physical processes controlling aerosol vertical profiles, in‐cloud wet removal is of utmost importance while its representation in global climate models (GCMs) is crude. In this study, we implement into the Community Atmosphere Model version 6 (CAM6) a physically‐based aerosol wet removal parameterization scheme that explicitly treats aerosol activation, removal and resuspension. Evaluation against in‐situ observations shows that the default scheme substantially overestimates the upper tropospheric black carbon (BC) and sea salt mass. Our physically‐based scheme reduces BC and sea salt mass by a factor of 10 and 1,000, respectively, in better agreement with observations. Also, the new scheme slightly increases number of aerosol particles between 12 nm and 4.8 μm in diameter, thereby mitigating the aerosol number underestimation in the default scheme. Our new scheme reduces the overestimation of coarse‐mode aerosol (0.5–4.8 μm) number. Overall, the aerosol property changes (mass decrease and number increase) reduce the cloud condensation nuclei (CCN) concentration at low supersaturation (i.e., 0.02% and 0.1%), and increase CCN at high supersaturations (i.e., 0.5% and 1%). Consequently, the global annual mean cloud liquid water path increases by 1.89 g m−2 and the ice water path increases by 0.51 g m−2. The global annual mean shortwave, longwave, and net cloud radiative forcing change by −1.06, 0.57, and −0.48 W m−2, respectively. Further improvement is needed to reflect the real physics that the removal efficiencies for aerosol mass and number are disproportionate and to advect cloud‐borne (activated) aerosols for a complete aerosol lifecycle.

Influence of ENSO on Tropospheric Ozone Variability in East Asia

AbstractAbundance of tropospheric ozone is determined by the interplay among emissions, chemistry, and climate dependent transport. It is important to identify the essential factors that modulate ozone in China and Korea because of increasing ozone trends in this area. In this study, we investigated the correlations of tropospheric ozone in East Asia with El Niño Southern Oscillations (ENSO). The model simulated ozone from 2000 to 2020 with Community Atmosphere Model version 6 with chemistry (CAM6‐chem), ECMWF Atmospheric Composition Reanalysis data, and ozone profiles measured in South Korea, Japan, and Hong Kong. These data were utilized to understand the ozone–ENSO relationship in this region. We found that tropospheric ozone in East Asia is significantly correlated with the wintertime Niño 3.4 index in El Niño/La Niña developing summer and decaying spring, with correlation coefficient up to 0.9. Particularly tropospheric ozone originated from stratosphere is strongly correlated with the Niño 3.4 index. Ozone concentrations in China and the Korean Peninsula are ∼5 ppb higher during El Niño developing summer compared to La Niña in a broad area, because of anomalous cyclonic/anticyclonic flows in the extratropics. In El Niño/La Niña decaying spring, the opposite occurs in the same region. In the decaying spring, stronger variability was observed in the subtropical regions of the lower to mid troposphere, indicating higher ozone concentrations during El Niño compared to La Niña. This variability is associated with the fluctuations in the Walker circulation, large‐scale sinking/rising motions, and the accompanying biomass burning events. The model reproduced the ozone variability associated with ENSO in the ozonesonde observations.

Insights of warm-cloud biases in Community Atmospheric Model 5 and 6 from the single-column modeling framework and Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) observations

Abstract. There has been a growing concern that most climate models predict precipitation that is too frequent, likely due to lack of reliable subgrid variability and vertical variations in microphysical processes in low-level warm clouds. In this study, the warm-cloud physics parameterizations in the singe-column configurations of NCAR Community Atmospheric Model version 6 and 5 (SCAM6 and SCAM5, respectively) are evaluated using ground-based and airborne observations from the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) field campaign near the Azores islands during 2017–2018. The 8-month single-column model (SCM) simulations show that both SCAM6 and SCAM5 can generally reproduce marine boundary layer cloud structure, major macrophysical properties, and their transition. The improvement in warm-cloud properties from the Community Atmospheric Model 5 and 6 (CAM5 to CAM6) physics can be found through comparison with the observations. Meanwhile, both physical schemes underestimate cloud liquid water content, cloud droplet size, and rain liquid water content but overestimate surface rainfall. Modeled cloud condensation nuclei (CCN) concentrations are comparable with aircraft-observed ones in the summer but are overestimated by a factor of 2 in winter, largely due to the biases in the long-range transport of anthropogenic aerosols like sulfate. We also test the newly recalibrated autoconversion and accretion parameterizations that account for vertical variations in droplet size. Compared to the observations, more significant improvement is found in SCAM5 than in SCAM6. This result is likely explained by the introduction of subgrid variations in cloud properties in CAM6 cloud microphysics, which further suppresses the scheme's sensitivity to individual warm-rain microphysical parameters. The predicted cloud susceptibilities to CCN perturbations in CAM6 are within a reasonable range, indicating significant progress since CAM5 which produces an aerosol indirect effect that is too strong. The present study emphasizes the importance of understanding biases in cloud physics parameterizations by combining SCM with in situ observations.