Atmospheric aerosol particle concentrations are strongly affected by various wet processes, including below and in-cloud wet scavenging and in-cloud aqueous-phase oxidation. We studied how wet ...scavenging and cloud processes affect particle concentrations and composition during transport to a rural boreal forest site in northern Europe. For this investigation, we employed air mass history analysis and observational data. Long-term particle number size distribution (â¼15 years) and composition measurements (â¼8 years) were combined with air mass trajectories with relevant variables from reanalysis data. Some such variables were rainfall rate, relative humidity, and mixing layer height. Additional observational datasets, such as temperature and trace gases, helped further evaluate wet processes along trajectories with mixed effects models.
This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. ...Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA / OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a−1 (range 34–144 Tg a−1) and the median SOA source strength (natural and anthropogenic) is 19 Tg a−1 (range 13–121 Tg a−1). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a−1 (range 16–121 Tg a−1), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a−1; range 13–20 Tg a−1, with one model at 37 Tg a−1). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6–2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8–9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a−1 (range 28–209 Tg a−1), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model–observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model–measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to −0.62 (−0.51) based on the comparison against OC (OA) urban data of all models at the surface, −0.15 (+0.51) when compared with remote measurements, and −0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. However, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately.
Aerosol-induced absorption of shortwave radiation can modify the climate through local atmospheric heating, which affects lapse rates, precipitation, and cloud formation. Presently, the total amount ...of aerosol absorption is poorly constrained, and the main absorbing aerosol species (black carbon (BC), organic aerosols (OA), and mineral dust) are diversely quantified in global climate models. As part of the third phase of the Aerosol Comparisons between Observations and Models (AeroCom) intercomparison initiative (AeroCom phase III), we here document the distribution and magnitude of aerosol absorption in current global aerosol models and quantify the sources of intermodel spread, highlighting the difficulties of attributing absorption to different species. In total, 15 models have provided total present-day absorption at 550 nm (using year 2010 emissions), 11 of which have provided absorption per absorbing species. The multi-model global annual mean total absorption aerosol optical depth (AAOD) is 0.0054 (0.0020 to 0.0098; 550 nm), with the range given as the minimum and maximum model values. This is 28 % higher compared to the 0.0042 (0.0021 to 0.0076) multi-model mean in AeroCom phase II (using year 2000 emissions), but the difference is within 1 standard deviation, which, in this study, is 0.0023 (0.0019 in Phase II). Of the summed component AAOD, 60 % (range 36 %–84 %) is estimated to be due to BC, 31 % (12 %–49 %) is due to dust, and 11 % (0 %–24 %) is due to OA; however, the components are not independent in terms of their absorbing efficiency. In models with internal mixtures of absorbing aerosols, a major challenge is the lack of a common and simple method to attribute absorption to the different absorbing species. Therefore, when possible, the models with internally mixed aerosols in the present study have performed simulations using the same method for estimating absorption due to BC, OA, and dust, namely by removing it and comparing runs with and without the absorbing species. We discuss the challenges of attributing absorption to different species; we compare burden, refractive indices, and density; and we contrast models with internal mixing to models with external mixing. The model mean BC mass absorption coefficient (MAC) value is 10.1 (3.1 to 17.7) m2 g−1 (550 nm), and the model mean BC AAOD is 0.0030 (0.0007 to 0.0077). The difference in lifetime (and burden) in the models explains as much of the BC AAOD spread as the difference in BC MAC values. The difference in the spectral dependency between the models is striking. Several models have an absorption Ångstrøm exponent (AAE) close to 1, which likely is too low given current knowledge of spectral aerosol optical properties. Most models do not account for brown carbon and underestimate the spectral dependency for OA.
Here we present for the first time a proof of concept for an emulation-based method that uses a large-eddy simulations (LESs) to present sub-grid cloud processes in a general circulation model (GCM). ...We focus on two key variables affecting the properties of shallow marine clouds: updraft velocity and precipitation formation. The LES is able to describe these processes with high resolution accounting for the realistic variability in cloud properties. We show that the selected emulation method is able to represent the LES outcome with relatively good accuracy and that the updraft velocity and precipitation emulators can be coupled with the GCM practically without increasing the computational costs. We also show that the emulators influence the climate simulated by the GCM but do not consistently improve or worsen the agreement with observations on cloud-related properties, although especially the updraft velocity at cloud base is better captured. A more quantitative evaluation of the emulator impacts against observations would, however, have required model re-tuning, which is a significant task and thus could not be included in this proof-of-concept study. All in all, the approach introduced here is a promising candidate for representing detailed cloud- and aerosol-related sub-grid processes in GCMs. Further development work together with increasing computing capacity can be expected to improve the accuracy and the applicability of the approach in climate simulations.
The vertical profile of aerosol is important for its radiative effects, but weakly constrained by observations on the global scale, and highly variable among different models. To investigate the ...controlling factors in one particular model, we investigate the effects of individual processes in HadGEM3-UKCA and compare the resulting diversity of aerosol vertical profiles with the inter-model diversity from the AeroCom Phase II control experiment. In this way we show that (in this model at least) the vertical profile is controlled by a relatively small number of processes, although these vary among aerosol components and particle sizes. We also show that sufficiently coarse variations in these processes can produce a similar diversity to that among different models in terms of the global-mean profile and, to a lesser extent, the zonal-mean vertical position. However, there are features of certain models' profiles that cannot be reproduced, suggesting the influence of further structural differences between models. In HadGEM3-UKCA, convective transport is found to be very important in controlling the vertical profile of all aerosol components by mass. In-cloud scavenging is very important for all except mineral dust. Growth by condensation is important for sulfate and carbonaceous aerosol (along with aqueous oxidation for the former and ageing by soluble material for the latter). The vertical extent of biomass-burning emissions into the free troposphere is also important for the profile of carbonaceous aerosol. Boundary-layer mixing plays a dominant role for sea salt and mineral dust, which are emitted only from the surface. Dry deposition and below-cloud scavenging are important for the profile of mineral dust only. In this model, the microphysical processes of nucleation, condensation and coagulation dominate the vertical profile of the smallest particles by number (e.g. total CN >3 nm), while the profiles of larger particles (e.g. CN>100 nm) are controlled by the same processes as the component mass profiles, plus the size distribution of primary emissions. We also show that the processes that affect the AOD-normalised radiative forcing in the model are predominantly those that affect the vertical mass distribution, in particular convective transport, in-cloud scavenging, aqueous oxidation, ageing and the vertical extent of biomass-burning emissions.
We studied the potential of using machine learning to downscale global-scale climate model output towards ground station data. The aim was to simultaneously analyze both city-level air quality and ...regional- and global-scale radiative forcing values for anthropogenic aerosols. As the city-level air pollution values are typically underestimated in global-scale models, we used a machine learning approach to downscale fine particulate (PM.sub.2.5) concentrations towards measured values. We first simulated the global climate with the aerosol-climate model ECHAM-HAMMOZ and corrected the PM.sub.2.5 values for the Indian megacity New Delhi.
Biomass burning plumes are frequently transported over the southeast Atlantic (SEA) stratocumulus deck during the southern African fire season (June–October). The plumes bring large amounts of ...absorbing aerosols and enhanced moisture, which can trigger a rich set of aerosol–cloud–radiation interactions with climatic consequences that are still poorly understood. We use large-eddy simulation (LES) to explore and disentangle the individual impacts of aerosols and moisture on the underlying stratocumulus clouds, the marine boundary layer (MBL) evolution, and the stratocumulus-to-cumulus transition (SCT) for three different meteorological situations over the southeast Atlantic during August 2017. For all three cases, our LES shows that the SCT is driven by increased sea surface temperatures and cloud-top entrainment as the air is advected towards the Equator. In the LES model, aerosol indirect effects, including impacts on drizzle production, have a small influence on the modeled cloud evolution and SCT, even when aerosol concentrations are lowered to background concentrations. In contrast, local semi-direct effects, i.e., aerosol absorption of solar radiation in the MBL, cause a reduction in cloud cover that can lead to a speed-up of the SCT, in particular during the daytime and during broken cloud conditions, especially in highly polluted situations. The largest impact on the radiative budget comes from aerosol impacts on cloud albedo: the plume with absorbing aerosols produces a total average 3 d of simulations. We find that the moisture accompanying the aerosol plume produces an additional cooling effect that is about as large as the total aerosol radiative effect. Overall, there is still a large uncertainty associated with the radiative and cloud evolution effects of biomass burning aerosols. A comparison between different models in a common framework, combined with constraints from in situ observations, could help to reduce the uncertainty.
In this paper we present simulations of the effect of nitric acid (HNO3) on cloud processing of aerosol particles. Sulfuric acid (H2SO4) production and incloud coagulation are both affected by ...condensed nitric acid as nitric acid increases the number of cloud droplets, which will lead to smaller mean size and higher total surface area of droplets. As a result of increased cloud droplet number concentration (CDNC), the incloud coagulation rate is enhanced by a factor of 1–1.3, so that the number of interstitial particles reduces faster. In addition, sulfuric acid production occurs in smaller particles and so the cloud processed aerosol size distribution is dependent on the HNO3 concentration. This affects both radiative properties of aerosol particles and the formation of cloud droplets during a sequence of cloud formation-evaporation events. It is shown that although the condensation of HNO3 increases the number of cloud droplets during the single updraft, it is possible that presence of HNO3 can actually decrease the cloud droplet number concentration after several cloud cycles when also H2SO4 production is taken into account.
We carried out a closure study of aerosol–cloud interactions during stratocumulus formation using a large eddy simulation model UCLALES–SALSA (University of California Los Angeles large eddy ...simulation model–sectional aerosol module for large applications) and observations from the 2020 cloud sampling campaign at Puijo SMEAR IV (Station for Measuring Ecosystem–Atmosphere Relations) in Kuopio, Finland. The unique observational setup combining in situ and cloud remote sensing measurements allowed a closer look into the aerosol size–composition dependence of droplet activation and droplet growth in turbulent boundary layer driven by surface forcing and radiative cooling. UCLALES–SALSA uses spectral bin microphysics for aerosols and hydrometeors, and incorporates a full description of their interactions into the turbulent-convective radiation-dynamical model of stratocumulus. Based on our results, the model successfully described the probability distribution of updraught velocities and consequently the size dependency of aerosol activation into cloud droplets, and further recreated the size distributions for both interstitial aerosol and cloud droplets. This is the first time such a detailed closure is achieved not only accounting for activation of cloud droplets in different updraughts, but also accounting for processes evaporating droplets and drizzle production through coagulation–coalescence. We studied two cases of cloud formation, one diurnal (24 September 2020) and one nocturnal (31 October 2020), with high and low aerosol loadings, respectively. Aerosol number concentrations differ more than 1 order of magnitude between cases and therefore, lead to cloud droplet number concentration (CDNC) values which range from less than 100 cm−3 up to 1000 cm−3. Different aerosol loadings affected supersaturation at the cloud base, and thus the size of aerosol particles activating to cloud droplets. Due to higher CDNC, the mean size of cloud droplets in the diurnal high aerosol case was lower. Thus, droplet evaporation in downdraughts affected more the observed CDNC at Puijo altitude compared to the low aerosol case. In addition, in the low aerosol case, the presence of large aerosol particles in the accumulation mode played a significant role in the droplet spectrum evolution as it promoted the drizzle formation through collision and coalescence processes. Also, during the event, the formation of ice particles was observed due to subzero temperature at the cloud top. Although the modelled number concentration of ice hydrometeors was too low to be directly measured, the retrieval of hydrometeor sedimentation velocities with cloud radar allowed us to assess the realism of modelled ice particles. The studied cases are presented in detail and can be further used by the cloud modellers to test and validate their models in a well-characterized modelling setup. We also provide recommendations on how increasing amount of information on aerosol properties could improve the understanding of processes affecting cloud droplet number and liquid water content in stratiform clouds.
The number of cloud droplets formed at the cloud base depends on both the properties of aerosol particles and the updraft velocity of an air parcel at the cloud base. As the spatial scale of updrafts ...is too small to be resolved in global atmospheric models, the updraft velocity is commonly parameterised based on the available turbulent kinetic energy. Here we present alternative methods through parameterising updraft velocity based on high-resolution large-eddy simulation (LES) runs in the case of marine stratocumulus clouds. First we use our simulations to assess the accuracy of a simple linear parameterisation where the updraft velocity depends only on cloud top radiative cooling. In addition, we present two different machine learning methods (Gaussian process emulation and random forest) that account for different boundary layer conditions and cloud properties. We conclude that both machine learning parameterisations reproduce the LES-based updraft velocities at about the same accuracy, while the simple approach employing radiative cooling only produces on average lower coefficient of determination and higher root mean square error values. Finally, we apply these machine learning methods to find the key parameters affecting cloud base updraft velocities.