Search

Wildfires and agricultural burning generate seemingly increasing smoke aerosol emissions, impacting societal and natural ecosystems. To understand smoke's effects on climate and public health, we analyzed the spatiotemporal distribution of smoke aerosols, focusing on two major light-absorbing components, namely black carbon (BC) and brown carbon (BrC) aerosols. Using NASA's Earth Polychromatic Imaging Camera (EPIC) instrument aboard NOAA's Deep Space Climate Observatory (DSCOVR) spacecraft, we inferred BC and BrC volume fractions and particle mass concentrations based on spectral absorption provided by the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm with 1–2 h temporal resolution and ∼ 10 km spatial resolution over North America and central Africa. Our analyses of regional smoke properties reveal distinct characteristics for aerosol optical depth (AOD) at 443 nm, spectral single-scattering albedo (SSA), aerosol layer height (ALH), and BC and BrC amounts. Smoke aerosols in North America showed extremely high AOD up to 6, with elevated ALH (6–7 km) and significant BrC components up to 250 mg m−2 along the transport paths, whereas the smoke aerosols in central Africa exhibited stronger light absorption (i.e., lower SSA) and lower AOD, resulting in higher-BC mass concentrations and similar BrC mass concentrations than the cases in North America. Seasonal burning source locations in central Africa, following the seasonal shift in the Intertropical Convergence Zone and diurnal variations in smoke amounts, were also captured. A comparison of retrieved AOD443, SSA443, SSA680, and ALH with collocated AERONET and CALIOP measurements shows agreement with RMSE values of 0.2, 0.03–0.04, 0.02–0.04, and 0.8–1.3 km, respectively. An analysis of the spatiotemporal average reveals distinct geographical characteristics in smoke properties closely linked to burning types and meteorological conditions. Forest wildfires over western North America generated smoke with a small-BC volume fraction of 0.011 and a high ALH with large variability (2.2 ± 1.2 km), whereas smoke from wildfires and agricultural burning over Mexico region shows more absorption and low ALH. Smoke from savanna fires over central Africa had the most absorption, with a high-BC volume fraction (0.015) and low ALH with a small variation (1.8 ± 0.6 km) among the analyzed regions. Tropical forest smoke was less absorbing and had a high variance in ALH. We also quantify the estimation uncertainties related to the assumptions of BC and BrC refractive indices. The MAIAC EPIC smoke properties with BC and BrC volume and mass fractions and assessment of the layer height provide observational constraints for radiative forcing modeling and air quality and health studies.

The Goddard Chemistry Aerosol Radiation and Transport (GOCART) model, which controls the sources, sinks, and chemistry of aerosols within the Goddard Earth Observing System (GEOS), recently underwent a major refactoring and update, including a revision of the emissions datasets and the addition of brown carbon. A 4-year benchmark simulation utilizing the new version of the model code, termed GOCART Second Generation (GOCART-2G) and coupled to the Goddard Earth Observing System (GEOS) model, was evaluated using in situ and spaceborne measurements to develop a baseline and prioritize future development. A comparison of simulated aerosol optical depth between GOCART-2G and MODIS retrievals indicates the model captures the overall spatial pattern and seasonal cycle of aerosol optical depth but overestimates aerosol extinction over dusty regions and underestimates aerosol extinction over Northern Hemisphere boreal forests, requiring further investigation and tuning of emissions. This MODIS-based analysis is corroborated by comparisons to MISR and selected AERONET stations; however, discrepancies between the Aqua and Terra satellites indicate there is a diurnal component to biases in aerosol optical depth over southern Asia and northern Africa. Despite the underestimate of aerosol optical depth in biomass burning regions in GEOS, there is an overestimate in the surface mass of organic carbon in the United States, especially during the summer months. Over Europe, GOCART-2G is unable to match the summertime peak in aerosol optical depth, opposing the observed late fall and early spring peaks in surface mass concentration. A comparison of the vertical profile of attenuated backscatter to observations from CALIPSO indicates the GEOS model is capable of capturing the vertical profile of aerosol; however, the mid-troposphere plumes of dust in the North Atlantic and smoke in the southeastern Atlantic are perhaps too low in altitude. The results presented highlight priorities for future development with GOCART-2G, including improvements for dust, biomass burning aerosols, and anthropogenic aerosols.

The atmospheric sulfur cycle plays a key role in air quality, climate, and ecosystems, such as pollution, radiative forcing, new particle formation, and acid rain. In this study, we compare the spatially and temporally resolved measurements from the NASA Atmospheric Tomography (ATom) mission with simulations from five AeroCom III models for four sulfur species (dimethyl sulfide (DMS), sulfur dioxide (SO2), particulate methanesulfonate (MSA), and particulate sulfate (SO4)). We focus on remote regions over the Pacific, Atlantic, and Southern oceans from near the surface to ∼ 12 km altitude range covering all four seasons. In general, the differences among model results can be greater than 1 order of magnitude. Comparing with observations, model-simulated SO2 is generally low, whereas SO4 is generally high. Simulated DMS concentrations near the sea surface exceed observed levels by a factor of 5 in most cases, suggesting potential overestimation of DMS emissions in all models. With GEOS model simulations of tagging emission from anthropogenic, biomass burning, volcanic, and oceanic sources, we find that anthropogenic emissions are the dominant source of sulfate aerosol (40 %–60 % of the total amount) in the ATom measurements at almost all altitudes, followed by volcanic emissions (18 %–32 %) and oceanic sources (16 %–32 %). Similar source contributions can also be derived at broad ocean basins and on monthly scales, indicating the representativeness of ATom measurements for global ocean. Our work presents the first assessment of AeroCom sulfur study using ATom measurements, providing directions for improving sulfate simulations, which remain the largest uncertainty in radiative forcing estimates in aerosol climate models.

Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO2 injection, seasonality of SO2 and black carbon (BC) particulate emissions, and the assumed fraction of SO2 emissions that is injected into the atmosphere as particulate phase sulfate (SO4) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO2, SO4, and BC, with a particularly large relative variation for SO2, which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of −0.35 W m−2), SO2 lifetime over Northern Hemisphere land (maximum difference of 0.8 d), surface SO2 concentration (up to 59 % decrease), and surface sulfate concentration (up to 23 % increase). Emitting SO2 at height consistently increased SO2 and SO4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO4 emission fraction are accurately and consistently represented in global models.

Global models are widely used to simulate biomass burning aerosol (BBA). Exhaustive evaluations on model representation of aerosol distributions and properties are fundamental to assess health and climate impacts of BBA. Here we conducted a comprehensive comparison of Aerosol Comparisons between Observations and Models (AeroCom) project model simulations with satellite observations. A total of 59 runs by 18 models from three AeroCom Phase-III experiments (i.e., biomass burning emissions, CTRL16, and CTRL19) and 14 satellite products of aerosols were used in the study. Aerosol optical depth (AOD) at 550 nm was investigated during the fire season over three key fire regions reflecting different fire dynamics (i.e., deforestation-dominated Amazon, Southern Hemisphere Africa where savannas are the key source of emissions, and boreal forest burning in boreal North America). The 14 satellite products were first evaluated against AErosol RObotic NETwork (AERONET) observations, with large uncertainties found. But these uncertainties had small impacts on the model evaluation that was dominated by modeling bias. Through a comparison with Polarization and Directionality of the Earth’s Reflectances measurements with the Generalized Retrieval of Aerosol and Surface Properties algorithm (POLDER-GRASP), we found that the modeled AOD values were biased by −93 % to 152 %, with most models showing significant underestimations even for the state-of-the-art aerosol modeling techniques (i.e., CTRL19). By scaling up BBA emissions, the negative biases in modeled AOD were significantly mitigated, although it yielded only negligible improvements in the correlation between models and observations, and the spatial and temporal variations in AOD biases did not change much. For models in CTRL16 and CTRL19, the large diversity in modeled AOD was in almost equal measures caused by diversity in emissions, lifetime, and the mass extinction coefficient (MEC). We found that in the AeroCom ensemble, BBA lifetime correlated significantly with particle deposition (as expected) and in turn correlated strongly with precipitation. Additional analysis based on Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) aerosol profiles suggested that the altitude of the aerosol layer in the current models was generally too low, which also contributed to the bias in modeled lifetime. Modeled MECs exhibited significant correlations with the Ångström exponent (AE, an indicator of particle size). Comparisons with the POLDER-GRASP-observed AE suggested that the models tended to overestimate the AE (underestimated particle size), indicating a possible underestimation of MECs in models. The hygroscopic growth in most models generally agreed with observations and might not explain the overall underestimation of modeled AOD. Our results imply that current global models contain biases in important aerosol processes for BBA (e.g., emissions, removal, and optical properties) that remain to be addressed in future research.

The iron-oxide content of dust in the atmosphere and most notably its apportionment between hematite (α-Fe2O3) and goethite (α-FeOOH) are key determinants in quantifying dust's light absorption, its top of atmosphere ultraviolet (UV) radiances used for dust monitoring, and ultimately shortwave dust direct radiative effects (DREs). Hematite and goethite column mass concentrations and iron-oxide mass fractions of total dust mass concentration were retrieved from the Deep Space Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) measurements in the ultraviolet–visible (UV–Vis) channels. The retrievals were performed for dust-identified aerosol plumes over land using aerosol optical depth (AOD) and the spectral imaginary refractive index provided by the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm over six continental regions (North America, North Africa, West Asia, Central Asia, East Asia, and Australia). The dust particles are represented as an internal mixture of non-absorbing host and absorbing hematite and goethite. We use the Maxwell Garnett effective medium approximation with carefully selected complex refractive indices of hematite and goethite that produce mass fractions of iron-oxide species consistent with in situ values found in the literature to derive the hematite and goethite volumetric/mass concentrations from MAIAC EPIC products. We compared the retrieved hematite and goethite concentrations with in situ dust aerosol mineralogical content measurements, as well as with published data. Our data display variations within the published range of hematite, goethite, and iron-oxide mass fractions for pure-mineral-dust cases. A specific analysis is presented for 15 sites over the main dust-source regions. Sites in the central Sahara, Sahel, and Middle East exhibit a greater temporal variability of iron oxides relative to other sites. The Niger site (13.52∘ N, 2.63∘ E) is dominated by goethite over the Harmattan season with a median of ∼ 2 weight percentage (wt %) of iron oxide. The Saudi Arabia site (27.49∘ N, 41.98∘ E) over the Middle East also exhibited a surge of goethite content with the beginning of the shamal season. The Sahel dust is richer in iron oxide than Saharan and northern China dust except in summer. The Bodélé Depression area shows a distinctively lower iron-oxide concentration (∼ 1 wt %) throughout the year. Finally, we show that EPIC data allow the constraining of the hematite refractive index. Specifically, we select 5 out of 13 different hematite refractive indices that are widely variable in published laboratory studies by constraining the iron-oxide mass ratio to the known measured values. The provided climatology of hematite and goethite mass fractions across the main dust regions of Earth will be useful for dust shortwave DRE studies and climate modeling.

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.

The Amazon experiences fires every year, and the resulting biomass burning aerosols, together with cloud particles, influence the penetration of sunlight through the atmosphere, increasing the ratio of diffuse to direct photosynthetically active radiation (PAR) reaching the vegetation canopy and thereby potentially increasing ecosystem productivity. In this study, we use the NASA Goddard Earth Observing System (GEOS) model with coupled aerosol, cloud, radiation, and ecosystem modules to investigate the impact of Amazon biomass burning aerosols on ecosystem productivity, as well as the role of the Amazon's clouds in tempering this impact. The study focuses on a 7-year period (2010–2016) during which the Amazon experienced a variety of dynamic environments (e.g., La Niña, normal years, and El Niño). The direct radiative impact of biomass burning aerosols on ecosystem productivity – called here the aerosol diffuse radiation fertilization effect – is found to increase Amazonian gross primary production (GPP) by 2.6 % via a 3.8 % increase in diffuse PAR (DFPAR) despite a 5.4 % decrease in direct PAR (DRPAR) on multiyear average during burning seasons. On a monthly basis, this increase in GPP can be as large as 9.9 % (occurring in August 2010). Consequently, the net primary production (NPP) in the Amazon is increased by 1.5 %, or ∼92 Tg C yr−1 – equivalent to ∼37 % of the average carbon lost due to Amazon fires over the 7 years considered. Clouds, however, strongly regulate the effectiveness of the aerosol diffuse radiation fertilization effect. The efficiency of this fertilization effect is the highest in cloud-free conditions and linearly decreases with increasing cloud amount until the cloud fraction reaches ∼0.8, at which point the aerosol-influenced light changes from being a stimulator to an inhibitor of plant growth. Nevertheless, interannual changes in the overall strength of the aerosol diffuse radiation fertilization effect are primarily controlled by the large interannual changes in biomass burning aerosols rather than by changes in cloudiness during the studied period.

This study characterizes a massive African dust intrusion into the Caribbean Basin and southern US in June 2020, which is nicknamed the “Godzilla” dust plume, using a comprehensive set of satellite and ground-based observations (including MODIS, CALIOP, SEVIRI, AERONET, and EPA Air Quality network) and the NASA GEOS global aerosol transport model. The MODIS data record registered this massive dust intrusion event as the most intense episode over the past 2 decades. During this event, the aerosol optical depth (AOD) observed by AERONET and MODIS peaked at 3.5 off the coast of West Africa and 1.8 in the Caribbean Basin. CALIOP observations show that the top of the dust plume reached altitudes of 6–8 km in West Africa and descended to about 4 km altitude over the Caribbean Basin and 2 km over the US Gulf of Mexico coast. The dust intrusion event degraded the air quality in Puerto Rico to a hazardous level, with maximum daily PM10 concentration of 453 µg m−3 recorded on 23 June. The dust intrusion into the US raised the PM2.5 concentration on 27 June to a level exceeding the EPA air quality standard in about 40 % of the stations in the southern US. Satellite observations reveal that dust emissions from convection-generated haboobs and other sources in West Africa were large albeit not extreme on a daily basis. However, the anomalous strength and northern shift of the North Atlantic Subtropical High (NASH) together with the Azores low formed a closed circulation pattern that allowed for accumulation of the dust near the African coast for about 4 d. When the NASH was weakened and wandered back to the south, the dust outflow region was dominated by a strong African easterly jet that rapidly transported the accumulated dust from the coastal region toward the Caribbean Basin, resulting in the record-breaking African dust intrusion. In comparison to satellite observations, the GEOS model reproduced the MODIS observed tracks of the meandering dust plume well as it was carried by the wind systems. However, the model substantially underestimated dust emissions from haboobs and did not lift up enough dust to the middle troposphere for ensuing long-range transport. Consequently, the model largely missed the satellite-observed elevated dust plume along the cross-ocean track and underestimated the dust intrusion into the Caribbean Basin by a factor of more than 4. Modeling improvements need to focus on developing more realistic representations of moist convection, haboobs, and the vertical transport of dust.

Black carbon (BC) emissions play an important role in regional climate change in the Arctic. It is necessary to pay attention to the impact of long-range transport from regions outside the Arctic as BC emissions from local sources in the Arctic were relatively small. The task force Hemispheric Transport of Air Pollution Phase 2 (HTAP2) set up a series of simulation scenarios to investigate the response of BC in a given region to different source regions. This study investigated the responses of Arctic BC concentrations and surface temperature to 20 % anthropogenic emission reductions from six regions in 2010 within the framework of HTAP2 based on ensemble modeling results. Emission reductions from East Asia (EAS) had the most (monthly contributions: 0.2–1.5 ng m−3) significant impact on the Arctic near-surface BC concentrations, while the monthly contributions from Europe (EUR), Middle East (MDE), North America (NAM), Russia–Belarus–Ukraine (RBU), and South Asia (SAS) were 0.2–1.0, 0.001–0.01, 0.1–0.3, 0.1–0.7, and 0.0–0.2 ng m−3, respectively. The responses of the vertical profiles of the Arctic BC to the six regions were found to be different due to multiple transport pathways. Emission reductions from NAM, RBU, EUR, and EAS mainly influenced the BC concentrations in the low troposphere of the Arctic, while most of the BC in the upper troposphere of the Arctic derived from SAS. The response of the Arctic BC to emission reductions in six source regions became less significant with the increase in the latitude. The benefit of BC emission reductions in terms of slowing down surface warming in the Arctic was evaluated by using absolute regional temperature change potential (ARTP). Compared to the response of global temperature to BC emission reductions, the response of Arctic temperature was substantially more sensitive, highlighting the need for curbing global BC emissions.

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM20) is approximately 5000 Tg yr−1, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, the relative contributions of the world's major source regions to the global dust cycle remain poorly constrained. This problem hinders accounting for the potentially large impact of regional differences in dust properties on clouds, the Earth's energy balance, and terrestrial and marine biogeochemical cycles. Here, we constrain the contribution of each of the world's main dust source regions to the global dust cycle. We use an analytical framework that integrates an ensemble of global aerosol model simulations with observationally informed constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth (DAOD). We obtain a dataset that constrains the relative contribution of nine major source regions to size-resolved dust emission, atmospheric loading, DAOD, concentration, and deposition flux. We find that the 22–29 Tg (1 standard error range) global loading of dust with a geometric diameter up to 20 µm is partitioned as follows: North African source regions contribute ∼ 50 % (11–15 Tg), Asian source regions contribute ∼ 40 % (8–13 Tg), and North American and Southern Hemisphere regions contribute ∼ 10 % (1.8–3.2 Tg). These results suggest that current models on average overestimate the contribution of North African sources to atmospheric dust loading at ∼ 65 %, while underestimating the contribution of Asian dust at ∼ 30 %. Our results further show that each source region's dust loading peaks in local spring and summer, which is partially driven by increased dust lifetime in those seasons. We also quantify the dust deposition flux to the Amazon rainforest to be ∼ 10 Tg yr−1, which is a factor of 2–3 less than inferred from satellite data by previous work that likely overestimated dust deposition by underestimating the dust mass extinction efficiency. The data obtained in this paper can be used to obtain improved constraints on dust impacts on clouds, climate, biogeochemical cycles, and other parts of the Earth system.

Within the framework of the AeroCom (Aerosol Comparisons between Observations and Models) initiative, the state-of-the-art modelling of aerosol optical properties is assessed from 14 global models participating in the phase III control experiment (AP3). The models are similar to CMIP6/AerChemMIP Earth System Models (ESMs) and provide a robust multi-model ensemble. Inter-model spread of aerosol species lifetimes and emissions appears to be similar to that of mass extinction coefficients (MECs), suggesting that aerosol optical depth (AOD) uncertainties are associated with a broad spectrum of parameterised aerosol processes. Total AOD is approximately the same as in AeroCom phase I (AP1) simulations. However, we find a 50 % decrease in the optical depth (OD) of black carbon (BC), attributable to a combination of decreased emissions and lifetimes. Relative contributions from sea salt (SS) and dust (DU) have shifted from being approximately equal in AP1 to SS contributing about 2∕3 of the natural AOD in AP3. This shift is linked with a decrease in DU mass burden, a lower DU MEC, and a slight decrease in DU lifetime, suggesting coarser DU particle sizes in AP3 compared to AP1. Relative to observations, the AP3 ensemble median and most of the participating models underestimate all aerosol optical properties investigated, that is, total AOD as well as fine and coarse AOD (AODf, AODc), Ångström exponent (AE), dry surface scattering (SCdry), and absorption (ACdry) coefficients. Compared to AERONET, the models underestimate total AOD by ca. 21 % ± 20 % (as inferred from the ensemble median and interquartile range). Against satellite data, the ensemble AOD biases range from −37 % (MODIS-Terra) to −16 % (MERGED-FMI, a multi-satellite AOD product), which we explain by differences between individual satellites and AERONET measurements themselves. Correlation coefficients (R) between model and observation AOD records are generally high (R>0.75), suggesting that the models are capable of capturing spatio-temporal variations in AOD. We find a much larger underestimate in coarse AODc (∼ −45 % ± 25 %) than in fine AODf (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AODc bias is likely due to missing DU over continental land masses (particularly over the United States, SE Asia, and S. America), while marine AERONET sites and the AATSR SU satellite data suggest more moderate oceanic biases in AODc. Column AEs are underestimated by about 10 % ± 16 %. For situations in which measurements show AE > 2, models underestimate AERONET AE by ca. 35 %. In contrast, all models (but one) exhibit large overestimates in AE when coarse aerosol dominates (bias ca. +140 % if observed AE < 0.5). Simulated AE does not span the observed AE variability. These results indicate that models overestimate particle size (or underestimate the fine-mode fraction) for fine-dominated aerosol and underestimate size (or overestimate the fine-mode fraction) for coarse-dominated aerosol. This must have implications for lifetime, water uptake, scattering enhancement, and the aerosol radiative effect, which we can not quantify at this moment. Comparison against Global Atmosphere Watch (GAW) in situ data results in mean bias and inter-model variations of −35 % ± 25 % and −20 % ± 18 % for SCdry and ACdry, respectively. The larger underestimate of SCdry than ACdry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SCdry than ambient air AOD is consistent with recent findings that models overestimate scattering enhancement due to hygroscopic growth. The broadly consistent negative bias in AOD and surface scattering suggests an underestimate of aerosol radiative effects in current global aerosol models. Considerable inter-model diversity in the simulated optical properties is often found in regions that are, unfortunately, not or only sparsely covered by ground-based observations. This includes, for instance, the Sahara, Amazonia, central Australia, and the South Pacific. This highlights the need for a better site coverage in the observations, which would enable us to better assess the models, but also the performance of satellite products in these regions. Using fine-mode AOD as a proxy for present-day aerosol forcing estimates, our results suggest that models underestimate aerosol forcing by ca. −15 %, however, with a considerably large interquartile range, suggesting a spread between −35 % and +10 %.

This study presents a multiparameter analysis of aerosol trends over the last 2 decades at regional and global scales. Regional time series have been computed for a set of nine optical, chemical-composition and mass aerosol properties by using the observations from several ground-based networks. From these regional time series the aerosol trends have been derived for the different regions of the world. Most of the properties related to aerosol loading exhibit negative trends, both at the surface and in the total atmospheric column. Significant decreases in aerosol optical depth (AOD) are found in Europe, North America, South America, North Africa and Asia, ranging from −1.2 % yr−1 to −3.1 % yr−1. An error and representativity analysis of the spatially and temporally limited observational data has been performed using model data subsets in order to investigate how much the observed trends represent the actual trends happening in the regions over the full study period from 2000 to 2014. This analysis reveals that significant uncertainty is associated with some of the regional trends due to time and space sampling deficiencies. The set of observed regional trends has then been used for the evaluation of 10 models (6 AeroCom phase III models and 4 CMIP6 models) and the CAMS reanalysis dataset and of their skills in reproducing the aerosol trends. Model performance is found to vary depending on the parameters and the regions of the world. The models tend to capture trends in AOD, the column Ångström exponent, sulfate and particulate matter well (except in North Africa), but they show larger discrepancies for coarse-mode AOD. The rather good agreement of the trends, across different aerosol parameters between models and observations, when co-locating them in time and space, implies that global model trends, including those in poorly monitored regions, are likely correct. The models can help to provide a global picture of the aerosol trends by filling the gaps in regions not covered by observations. The calculation of aerosol trends at a global scale reveals a different picture from that depicted by solely relying on ground-based observations. Using a model with complete diagnostics (NorESM2), we find a global increase in AOD of about 0.2 % yr−1 between 2000 and 2014, primarily caused by an increase in the loads of organic aerosols, sulfate and black carbon.

The radiative forcing of the aerosol–radiation interaction can be decomposed into clear-sky and cloudy-sky portions. Two sets of multi-model simulations within Aerosol Comparisons between Observations and Models (AeroCom), combined with observational methods, and the time evolution of aerosol emissions over the industrial era show that the contribution from cloudy-sky regions is likely weak. A mean of the simulations considered is 0.01±0.1 W m−2. Multivariate data analysis of results from AeroCom Phase II shows that many factors influence the strength of the cloudy-sky contribution to the forcing of the aerosol–radiation interaction. Overall, single-scattering albedo of anthropogenic aerosols and the interaction of aerosols with the short-wave cloud radiative effects are found to be important factors. A more dedicated focus on the contribution from the cloud-free and cloud-covered sky fraction, respectively, to the aerosol–radiation interaction will benefit the quantification of the radiative forcing and its uncertainty range.

The spatial distribution and properties of submicron organic aerosol (OA) are among the key sources of uncertainty in our understanding of aerosol effects on climate. Uncertainties are particularly large over remote regions of the free troposphere and Southern Ocean, where very few data have been available and where OA predictions from AeroCom Phase II global models span 2 to 3 orders of magnitude, greatly exceeding the model spread over source regions. The (nearly) pole-to-pole vertical distribution of non-refractory aerosols was measured with an aerosol mass spectrometer onboard the NASA DC-8 aircraft as part of the Atmospheric Tomography (ATom) mission during the Northern Hemisphere summer (August 2016) and winter (February 2017). This study presents the first extensive characterization of OA mass concentrations and their level of oxidation in the remote atmosphere. OA and sulfate are the major contributors by mass to submicron aerosols in the remote troposphere, together with sea salt in the marine boundary layer. Sulfate was dominant in the lower stratosphere. OA concentrations have a strong seasonal and zonal variability, with the highest levels measured in the lower troposphere in the summer and over the regions influenced by biomass burning from Africa (up to 10 µg sm−3). Lower concentrations (∼0.1–0.3 µg sm−3) are observed in the northern middle and high latitudes and very low concentrations ( µg sm−3) in the southern middle and high latitudes. The ATom dataset is used to evaluate predictions of eight current global chemistry models that implement a variety of commonly used representations of OA sources and chemistry, as well as of the AeroCom-II ensemble. The current model ensemble captures the average vertical and spatial distribution of measured OA concentrations, and the spread of the individual models remains within a factor of 5. These results are significantly improved over the AeroCom-II model ensemble, which shows large overestimations over these regions. However, some of the improved agreement with observations occurs for the wrong reasons, as models have the tendency to greatly overestimate the primary OA fraction and underestimate the secondary fraction. Measured OA in the remote free troposphere is highly oxygenated, with organic aerosol to organic carbon (OA ∕ OC) ratios of ∼2.2–2.8, and is 30 %–60 % more oxygenated than in current models, which can lead to significant errors in OA concentrations. The model–measurement comparisons presented here support the concept of a more dynamic OA system as proposed by Hodzic et al. (2016), with enhanced removal of primary OA and a stronger production of secondary OA in global models needed to provide better agreement with observations.

Horizontal grid resolution has a profound effect on model performances on meteorology and air quality simulations. In contribution to MICS-Asia Phase III, one of whose goals was to identify and reduce model uncertainty in air quality prediction, this study examined the impact of grid resolution on meteorology and air quality simulation over East Asia, focusing on the North China Plain (NCP) region. The NASA Unified Weather Research and Forecasting (NU-WRF) model has been applied with the horizontal resolutions at 45, 15, and 5 km. The results revealed that, in comparison with ground observations, no single resolution can yield the best model performance for all variables across all stations. From a regional average perspective (i.e., across all monitoring sites), air temperature modeling was not sensitive to the grid resolution but wind and precipitation simulation showed the opposite. NU-WRF with the 5 km grid simulated the wind speed best, while the 45 km grid yielded the most realistic precipitation as compared to the site observations. For air quality simulations, finer resolution generally led to better comparisons with observations for O3, CO, NOx, and PM2.5. However, the improvement of model performance on air quality was not linear with the resolution increase. The accuracy of modeled surface O3 of the 15 km grid was greatly improved over the one from the 45 km grid. A further increase in grid resolution to 5 km, however, showed diminished impact on model performance improvement on O3 prediction. In addition, a 5 km resolution grid showed large advantage in better capturing the frequency of high-pollution occurrences. This was important for the assessment of noncompliance with ambient air quality standards, which was key to air quality planning and management. Balancing the modeling accuracy and resource limitation, a 15 km grid resolution was suggested for future MICS-Asia air quality modeling activity if the research region remained unchanged. This investigation also found a large overestimate of ground-level O3 and an underestimate of surface NOx and CO, likely due to missing emissions of NOx and CO.

Emissions and long-range transport of mineral dust and combustion-related aerosol from burning fossil fuels and biomass vary from year to year, driven by the evolution of the economy and changes in meteorological conditions and environmental regulations. This study offers both satellite and model perspectives on the interannual variability and possible trends of combustion aerosol and dust in major continental outflow regions over the past 15 years (2003–2017). The decade-long record of aerosol optical depth (AOD, denoted as τ), separately for combustion aerosol (τc) and dust (τd), over global oceans is derived from the Collection 6 aerosol products of the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard both Terra and Aqua. These MODIS Aqua datasets, complemented by aerosol source-tagged simulations using the Community Atmospheric Model version 5 (CAM5), are then analyzed to understand the interannual variability and potential trends of τc and τd in the major continental outflows. Both MODIS and CAM5 consistently yield a similar decreasing trend of −0.017 to −0.020 per decade for τc over the North Atlantic Ocean and the Mediterranean Sea that is attributable to reduced emissions from North America and Europe, respectively. On the contrary, both MODIS and CAM5 display an increasing trend of +0.017 to +0.036 per decade for τc over the tropical Indian Ocean, the Bay of Bengal, and the Arabian Sea, which reflects the influence of increased anthropogenic emissions from South Asia and the Middle East in the last 2 decades. Over the northwestern Pacific Ocean, which is often affected by East Asian emissions of pollution and dust, the MODIS retrievals show a decreasing trend of −0.021 per decade for τc and −0.012 per decade for τd, which is, however, not reproduced by the CAM5 model. In other outflow regions strongly influenced by biomass burning smoke or dust, both MODIS retrievals and CAM5 simulations show no statistically significant trends; the MODIS-observed interannual variability is usually larger than that of the CAM5 simulation.

Aerosols from biomass burning (BB) emissions are poorly constrained in global and regional models, resulting in a high level of uncertainty in understanding their impacts. In this study, we compared six BB aerosol emission datasets for 2008 globally as well as in 14 regions. The six BB emission datasets are (1) GFED3.1 (Global Fire Emissions Database version 3.1), (2) GFED4s (GFED version 4 with small fires), (3) FINN1.5 (FIre INventory from NCAR version 1.5), (4) GFAS1.2 (Global Fire Assimilation System version 1.2), (5) FEER1.0 (Fire Energetics and Emissions Research version 1.0), and (6) QFED2.4 (Quick Fire Emissions Dataset version 2.4). The global total emission amounts from these six BB emission datasets differed by a factor of 3.8, ranging from 13.76 to 51.93 Tg for organic carbon and from 1.65 to 5.54 Tg for black carbon. In most of the regions, QFED2.4 and FEER1.0, which are based on satellite observations of fire radiative power (FRP) and constrained by aerosol optical depth (AOD) data from the Moderate Resolution Imaging Spectroradiometer (MODIS), yielded higher BB aerosol emissions than the rest by a factor of 2–4. By comparison, the BB aerosol emissions estimated from GFED4s and GFED3.1, which are based on satellite burned-area data, without AOD constraints, were at the low end of the range. In order to examine the sensitivity of model-simulated AOD to the different BB emission datasets, we ingested these six BB emission datasets separately into the same global model, the NASA Goddard Earth Observing System (GEOS) model, and compared the simulated AOD with observed AOD from the AErosol RObotic NETwork (AERONET) and the Multiangle Imaging SpectroRadiometer (MISR) in the 14 regions during 2008. In Southern Hemisphere Africa (SHAF) and South America (SHSA), where aerosols tend to be clearly dominated by smoke in September, the simulated AOD values were underestimated in almost all experiments compared to MISR, except for the QFED2.4 run in SHSA. The model-simulated AOD values based on FEER1.0 and QFED2.4 were the closest to the corresponding AERONET data, being, respectively, about 73 % and 100 % of the AERONET observed AOD at Alta Floresta in SHSA and about 49 % and 46 % at Mongu in SHAF. The simulated AOD based on the other four BB emission datasets accounted for only ∼50 % of the AERONET AOD at Alta Floresta and ∼20 % at Mongu. Overall, during the biomass burning peak seasons, at most of the selected AERONET sites in each region, the AOD values simulated with QFED2.4 were the highest and closest to AERONET and MISR observations, followed closely by FEER1.0. However, the QFED2.4 run tends to overestimate AOD in the region of SHSA, and the QFED2.4 BB emission dataset is tuned with the GEOS model. In contrast, the FEER1.0 BB emission dataset is derived in a more model-independent fashion and is more physically based since its emission coefficients are independently derived at each grid box. Therefore, we recommend the FEER1.0 BB emission dataset for aerosol-focused hindcast experiments in the two biomass-burning-dominated regions in the Southern Hemisphere, SHAF, and SHSA (as well as in other regions but with lower confidence). The differences between these six BB emission datasets are attributable to the approaches and input data used to derive BB emissions, such as whether AOD from satellite observations is used as a constraint, whether the approaches to parameterize the fire activities are based on burned area, FRP, or active fire count, and which set of emission factors is chosen.

Topic 3 of the Model Inter-Comparison Study for Asia (MICS-Asia) Phase III examines how online coupled air quality models perform in simulating wintertime haze events in the North China Plain region and evaluates the importance of aerosol radiative feedbacks. This paper discusses the estimates of aerosol radiative forcing, aerosol feedbacks, and possible causes for the differences among the participating models. Over the Beijing–Tianjin–Hebei (BTH) region, the ensemble mean of estimated aerosol direct radiative forcing (ADRF) at the top of atmosphere, inside the atmosphere, and at the surface are −1.1, 7.7, and −8.8 W m−2 during January 2010, respectively. Subdivisions of direct and indirect aerosol radiative forcing confirm the dominant role of direct forcing. During severe haze days (17–19 January 2010), the averaged reduction in near-surface temperature for the BTH region can reach 0.3–1.6 ∘C. The responses of wind speeds at 10 m (WS10) inferred from different models show consistent declines in eastern China. For the BTH region, aerosol–radiation feedback-induced daytime changes in PM2.5 concentrations during severe haze days range from 6.0 to 12.9 µg m−3 ( %). Sensitivity simulations indicate the important effect of aerosol mixing states on the estimates of ADRF and aerosol feedbacks. Besides, black carbon (BC) exhibits a large contribution to atmospheric heating and feedbacks although it accounts for a small share of mass concentration of PM2.5.

We invert global black carbon (BC), organic carbon (OC) and desert dust (DD) aerosol emissions from POLDER/PARASOL spectral aerosol optical depth (AOD) and aerosol absorption optical depth (AAOD) using the GEOS-Chem inverse modeling framework. Our inverse modeling framework uses standard a priori emissions to provide a posteriori emissions that are constrained by POLDER/PARASOL AODs and AAODs. The following global emission values were retrieved for the three aerosol components: 18.4 Tg yr−1 for BC, 109.9 Tg yr−1 for OC and 731.6 Tg yr−1 for DD for the year 2010. These values show a difference of +166.7 %, +184.0 % and −42.4 %, respectively, with respect to the a priori values of emission inventories used in “standard” GEOS-Chem runs. The model simulations using a posteriori emissions (i.e., retrieved emissions) provide values of 0.119 for global mean AOD and 0.0071 for AAOD at 550 nm, which are +13.3 % and +82.1 %, respectively, higher than the AOD and AAOD obtained using the a priori values of emissions. Additionally, the a posteriori model simulation of AOD, AAOD, single scattering albedo, Ångström exponent and absorption Ångström exponent show better agreement with independent AERONET, MODIS and OMI measurements than the a priori simulation. Thus, this study suggests that using satellite-constrained global aerosol emissions in aerosol transport models can improve the accuracy of simulated global aerosol properties.

Atmospheric sea salt plays important roles in marine cloud formation and atmospheric chemistry. We performed an integrated analysis of NASA GEOS model simulations run with the GOCART aerosol module, in situ measurements from the PALMS and SAGA instruments obtained during the NASA ATom campaign, and aerosol optical depth (AOD) measurements from the AERONET Marine Aerosol Network (MAN) and from MODIS satellite observations to better constrain sea salt in the marine atmosphere. ATom measurements and GEOS model simulations both show that sea salt concentrations over the Pacific and Atlantic oceans have a strong vertical gradient, varying up to 4 orders of magnitude from the marine boundary layer to free troposphere. The modeled residence times suggest that the lifetime of sea salt particles with a dry diameter of less than 3 µm is largely controlled by wet removal, followed by turbulent process. During both boreal summer and winter, the GEOS-simulated sea salt mass mixing ratios agree with SAGA measurements in the marine boundary layer (MBL) and with PALMS measurements above the MBL. However, comparison of AOD from GEOS with AERONET/MAN and MODIS aerosol retrievals indicated that the model underestimated AOD over the oceans where sea salt dominates. The apparent discrepancy of slightly overpredicted concentration and large underpredicted AOD could not be explained by biases in the model RH affecting the particle hygroscopic growth, as modeled RH was found to be comparable to or larger than the in situ measurements. This conundrum could at least partially be explained by the difference in sea salt size distribution; the GEOS simulation has much less sea salt percentage-wise in the smaller particle size range and thus less efficient light extinction than what was observed by PALMS.

Haze has been severely affecting the densely populated areas in China recently. While many of the efforts have been devoted to investigating the impact of local anthropogenic emission, limited attention has been paid to the contribution from long-range transport. In this study, we apply simulations from six participating models supplied through the Task Force on Hemispheric Transport of Air Pollution phase 2 (HTAP2) exercise to investigate the long-range transport impact of Europe (EUR) and Russia–Belarus–Ukraine (RBU) on the surface air quality in eastern Asia (EAS), with special focus on their contributions during the haze episodes in China. The impact of 20 % anthropogenic emission perturbation from the source region is extrapolated by a factor of 5 to estimate the full impact. We find that the full impacts from EUR and RBU are 0.99 µg m−3 (3.1 %) and 1.32 µg m−3 (4.1 %) during haze episodes, while the annual averaged full impacts are only 0.35 µg m−3 (1.7 %) and 0.53 µg m−3 (2.6 %). By estimating the aerosol response within and above the planetary boundary layer (PBL), we find that long-range transport from EUR within the PBL contributes to 22–38 % of the total column density of aerosol response in EAS. Comparison with the HTAP phase 1 (HTAP1) assessment reveals that from 2000 to 2010, the long-range transport from Europe to eastern Asia has decreased significantly by a factor of 2–10 for surface aerosol mass concentration due to the simultaneous emission reduction in source regions and emission increase in the receptor region. We also find the long-range transport from the Europe and RBU regions increases the number of haze events in China by 0.15 % and 0.11 %, and the North China Plain and southeastern China has 1–3 extra haze days (

Understanding the role atmospheric aerosols play in the Earth–atmosphere system is limited by uncertainties in the knowledge of their distribution, composition and sources. In this paper, we use the GEOS-Chem based inverse modelling framework for retrieving desert dust (DD), black carbon (BC) and organic carbon (OC) aerosol emissions simultaneously. Aerosol optical depth (AOD) and aerosol absorption optical depth (AAOD) retrieved from the multi-angular and polarimetric POLDER/PARASOL measurements generated by the GRASP algorithm (hereafter PARASOL/GRASP) have been assimilated. First, the inversion framework is validated in a series of numerical tests conducted with synthetic PARASOL-like data. These tests show that the framework allows for retrieval of the distribution and strength of aerosol emissions. The uncertainty of retrieved daily emissions in error free conditions is below 25.8 % for DD, 5.9 % for BC and 26.9 % for OC. In addition, the BC emission retrieval is sensitive to BC refractive index, which could produce an additional factor of 1.8 differences for total BC emissions. The approach is then applied to 1 year (December 2007 to November 2008) of data over the African and Arabian Peninsula region using PARASOL/GRASP spectral AOD and AAOD at six wavelengths (443, 490, 565, 670, 865 and 1020 nm). Analysis of the resulting retrieved emissions indicates 1.8 times overestimation of the prior DD online mobilization and entrainment model. For total BC and OC, the retrieved emissions show a significant increase of 209.9 %–271.8 % in comparison to the prior carbonaceous aerosol emissions. The model posterior simulation with retrieved emissions shows good agreement with both the AOD and AAOD PARASOL/GRASP products used in the inversion. The fidelity of the results is evaluated by comparison of posterior simulations with measurements from AERONET that are completely independent measurements and more temporally frequent than PARASOL observations. To further test the robustness of our posterior emissions constrained using PARASOL/GRASP, the posterior emissions are implemented in the GEOS-5/GOCART model and the consistency of simulated AOD and AAOD with other independent measurements (MODIS and OMI) demonstrates promise in applying this database for modelling studies.

Ambient air pollution from ozone and fine particulate matter is associated with premature mortality. As emissions from one continent influence air quality over others, changes in emissions can also influence human health on other continents. We estimate global air-pollution-related premature mortality from exposure to PM2.5 and ozone and the avoided deaths due to 20 % anthropogenic emission reductions from six source regions, North America (NAM), Europe (EUR), South Asia (SAS), East Asia (EAS), Russia–Belarus–Ukraine (RBU), and the Middle East (MDE), three global emission sectors, power and industry (PIN), ground transportation (TRN), and residential (RES), and one global domain (GLO), using an ensemble of global chemical transport model simulations coordinated by the second phase of the Task Force on Hemispheric Transport of Air Pollutants (TF HTAP2), and epidemiologically derived concentration response functions. We build on results from previous studies of TF HTAP by using improved atmospheric models driven by new estimates of 2010 anthropogenic emissions (excluding methane), with more source and receptor regions, new consideration of source sector impacts, and new epidemiological mortality functions. We estimate 290 000 (95 % confidence interval (CI): 30 000, 600 000) premature O3-related deaths and 2.8 million (0.5 million, 4.6 million) PM2.5-related premature deaths globally for the baseline year 2010. While 20 % emission reductions from one region generally lead to more avoided deaths within the source region than outside, reducing emissions from MDE and RBU can avoid more O3-related deaths outside of these regions than within, and reducing MDE emissions also avoids more PM2.5-related deaths outside of MDE than within. Our findings that most avoided O3-related deaths from emission reductions in NAM and EUR occur outside of those regions contrast with those of previous studies, while estimates of PM2.5-related deaths from NAM, EUR, SAS, and EAS emission reductions agree well. In addition, EUR, MDE, and RBU have more avoided O3-related deaths from reducing foreign emissions than from domestic reductions. For six regional emission reductions, the total avoided extra-regional mortality is estimated as 6000 (−3400, 15 500) deaths per year and 25 100 (8200, 35 800) deaths per year through changes in O3 and PM2.5, respectively. Interregional transport of air pollutants leads to more deaths through changes in PM2.5 than in O3, even though O3 is transported more on interregional scales, since PM2.5 has a stronger influence on mortality. For NAM and EUR, our estimates of avoided mortality from regional and extra-regional emission reductions are comparable to those estimated by regional models for these same experiments. In sectoral emission reductions, TRN emissions account for the greatest fraction (26–53 % of global emission reduction) of O3-related premature deaths in most regions, in agreement with previous studies, except for EAS (58 %) and RBU (38 %) where PIN emissions dominate. In contrast, PIN emission reductions have the greatest fraction (38–78 % of global emission reduction) of PM2.5-related deaths in most regions, except for SAS (45 %) where RES emission dominates, which differs with previous studies in which RES emissions dominate global health impacts. The spread of air pollutant concentration changes across models contributes most to the overall uncertainty in estimated avoided deaths, highlighting the uncertainty in results based on a single model. Despite uncertainties, the health benefits of reduced intercontinental air pollution transport suggest that international cooperation may be desirable to mitigate pollution transported over long distances.

The Stratospheric Sulfur and its Role in Climate (SSiRC) Interactive Stratospheric Aerosol Model Intercomparison Project (ISA-MIP) explores uncertainties in the processes that connect volcanic emission of sulfur gas species and the radiative forcing associated with the resulting enhancement of the stratospheric aerosol layer. The central aim of ISA-MIP is to constrain and improve interactive stratospheric aerosol models and reduce uncertainties in the stratospheric aerosol forcing by comparing results of standardized model experiments with a range of observations. In this paper we present four co-ordinated inter-model experiments designed to investigate key processes which influence the formation and temporal development of stratospheric aerosol in different time periods of the observational record. The Background (BG) experiment will focus on microphysics and transport processes under volcanically quiescent conditions, when the stratospheric aerosol is controlled by the transport of aerosols and their precursors from the troposphere to the stratosphere. The Transient Aerosol Record (TAR) experiment will explore the role of small- to moderate-magnitude volcanic eruptions, anthropogenic sulfur emissions, and transport processes over the period 1998–2012 and their role in the warming hiatus. Two further experiments will investigate the stratospheric sulfate aerosol evolution after major volcanic eruptions. The Historical Eruptions SO2 Emission Assessment (HErSEA) experiment will focus on the uncertainty in the initial emission of recent large-magnitude volcanic eruptions, while the Pinatubo Emulation in Multiple models (PoEMS) experiment will provide a comprehensive uncertainty analysis of the radiative forcing from the 1991 Mt Pinatubo eruption.

Background: Incidence of infection following total elbow replacement (TER) is recognised to be higher compared to hip or knee arthroplasty. Extensive swelling following TER can complicate the wound healing which might lead to infection. Tranexamic Acid (TXA) is proven to reduce blood loss peri-operatively which might contribute to better healing outcomes. Our aim is to assess the effect of TXA in wound healing following TER.Methods: A retrospective review of a single surgeon case series. 10 patients had TER mainly for complicated elbow fractures, four of them were relatively immune-supressed. All patients had 2 grams of TXA and antibiotics intra-operatively. All were reviewed at two weeks following surgery for wound check and removal of surgical clips.Results: Seven females and three males with a mean age of 81.5 had TER and TXA. The mean level of pre-operative haemoglobin was 134.40 g/l and the mean post-operative level was 122.70g/l. No patient in this series required blood transfusion. At two weeks and six weeks follow-up, all wound healed up with no signs of infection.Conclusion: TXA has been proven to be safe an effective way of reducing peri-operative bleeding. TXA maintains haemostasis after releasing the tourniquet and therefore reduces the swelling and wound complications post-operatively

Topic 3 of the Model Inter-Comparison Study for Asia (MICS-Asia) Phase III examines how online coupled air quality models perform in simulating high aerosol pollution in the North China Plain region during wintertime haze events and evaluates the importance of aerosol radiative and microphysical feedbacks. A comprehensive overview of the MICS-Asia III Topic 3 study design, including descriptions of participating models and model inputs, the experimental designs, and results of model evaluation, are presented. Six modeling groups from China, Korea and the United States submitted results from seven applications of online coupled chemistry–meteorology models. Results are compared to meteorology and air quality measurements, including data from the Campaign on Atmospheric Aerosol Research Network of China (CARE-China) and the Acid Deposition Monitoring Network in East Asia (EANET). The correlation coefficients between the multi-model ensemble mean and the CARE-China observed near-surface air pollutants range from 0.51 to 0.94 (0.51 for ozone and 0.94 for PM2.5) for January 2010. However, large discrepancies exist between simulated aerosol chemical compositions from different models. The coefficient of variation (SD divided by the mean) can reach above 1.3 for sulfate in Beijing and above 1.6 for nitrate and organic aerosols in coastal regions, indicating that these compositions are less consistent from different models. During clean periods, simulated aerosol optical depths (AODs) from different models are similar, but peak values differ during severe haze events, which can be explained by the differences in simulated inorganic aerosol concentrations and the hygroscopic growth efficiency (affected by varied relative humidity). These differences in composition and AOD suggest that future models can be improved by including new heterogeneous or aqueous pathways for sulfate and nitrate formation under hazy conditions, a secondary organic aerosol (SOA) formation chemical mechanism with new volatile organic compound (VOCs) precursors, yield data and approaches, and a more detailed evaluation of the dependence of aerosol optical properties on size distribution and mixing state. It was also found that using the ensemble mean of the models produced the best prediction skill. While this has been shown for other conditions (for example, the prediction of high-ozone events in the US (McKeen et al., 2005)), this is to our knowledge the first time it has been shown for heavy haze events.

An assessment of global particulate nitrate and ammonium aerosol based on simulations from nine models participating in the Aerosol Comparisons between Observations and Models (AeroCom) phase III study is presented. A budget analysis was conducted to understand the typical magnitude, distribution, and diversity of the aerosols and their precursors among the models. To gain confidence regarding model performance, the model results were evaluated with various observations globally, including ground station measurements over North America, Europe, and east Asia for tracer concentrations and dry and wet depositions, as well as with aircraft measurements in the Northern Hemisphere mid-to-high latitudes for tracer vertical distributions. Given the unique chemical and physical features of the nitrate occurrence, we further investigated the similarity and differentiation among the models by examining (1) the pH-dependent NH3 wet deposition; (2) the nitrate formation via heterogeneous chemistry on the surface of dust and sea salt particles or thermodynamic equilibrium calculation including dust and sea salt ions; and (3) the nitrate coarse-mode fraction (i.e., coarse/total). It is found that HNO3, which is simulated explicitly based on full O3-HOx-NOx-aerosol chemistry by all models, differs by up to a factor of 9 among the models in its global tropospheric burden. This partially contributes to a large difference in NO3−, whose atmospheric burden differs by up to a factor of 13. The atmospheric burdens of NH3 and NH4+ differ by 17 and 4, respectively. Analyses at the process level show that the large diversity in atmospheric burdens of NO3−, NH3, and NH4+ is also related to deposition processes. Wet deposition seems to be the dominant process in determining the diversity in NH3 and NH4+ lifetimes. It is critical to correctly account for contributions of heterogeneous chemical production of nitrate on dust and sea salt, because this process overwhelmingly controls atmospheric nitrate production (typically > 80 %) and determines the coarse- and fine-mode distribution of nitrate aerosol.

Atmospheric aerosols from anthropogenic and natural sources reach the polar regions through long-range transport and affect the local radiation balance. Such transport is, however, poorly constrained in present-day global climate models, and few multi-model evaluations of polar anthropogenic aerosol radiative forcing exist. Here we compare the aerosol optical depth (AOD) at 550 nm from simulations with 16 global aerosol models from the AeroCom Phase II model intercomparison project with available observations at both poles. We show that the annual mean multi-model median is representative of the observations in Arctic, but that the intermodel spread is large. We also document the geographical distribution and seasonal cycle of the AOD for the individual aerosol species: black carbon (BC) from fossil fuel and biomass burning, sulfate, organic aerosols (OAs), dust, and sea-salt. For a subset of models that represent nitrate and secondary organic aerosols (SOAs), we document the role of these aerosols at high latitudes.The seasonal dependence of natural and anthropogenic aerosols differs with natural aerosols peaking in winter (sea-salt) and spring (dust), whereas AOD from anthropogenic aerosols peaks in late spring and summer. The models produce a median annual mean AOD of 0.07 in the Arctic (defined here as north of 60° N). The models also predict a noteworthy aerosol transport to the Antarctic (south of 70° S) with a resulting AOD varying between 0.01 and 0.02. The models have estimated the shortwave anthropogenic radiative forcing contributions to the direct aerosol effect (DAE) associated with BC and OA from fossil fuel and biofuel (FF), sulfate, SOAs, nitrate, and biomass burning from BC and OA emissions combined. The Arctic modelled annual mean DAE is slightly negative (−0.12 W m−2), dominated by a positive BC FF DAE in spring and a negative sulfate DAE in summer. The Antarctic DAE is governed by BC FF. We perform sensitivity experiments with one of the AeroCom models (GISS modelE) to investigate how regional emissions of BC and sulfate and the lifetime of BC influence the Arctic and Antarctic AOD. A doubling of emissions in eastern Asia results in a 33 % increase in Arctic AOD of BC. A doubling of the BC lifetime results in a 39 % increase in Arctic AOD of BC. However, these radical changes still fall within the AeroCom model range.

This work reports on accurate, high-resolution measurements of the 25Mg(n,γ)26Mg and 25Mg(n,tot) cross sections in the neutron energy range from thermal to about 300 keV, leading to a significantly improved 25Mg(n,γ)26Mg parametrization. The relevant resonances for n+25Mg were characterized from a combined R-matrix analysis of the experimental data. This resulted in an unambiguous spin/parity assignment of the corresponding excited states in 26Mg. With this information experimental upper limits of the reaction rates for 22Ne(α,n)25Mg and 22Ne(α,γ)26Mg were established, potentially leading to a significantly higher (α,n)/(α,γ) ratio than previously evaluated. The impact of these results has been studied for stellar models in the mass range 2 to 25 M⊙.

Abstract Introduction Data on real-world clinical practice and outcomes of patients with pulmonary arterial hypertension (PAH) associated with congenital heart disease (CHD) are scarce. The OPUS/OrPHeUS studies enrolled patients with PAH newly initiating macitentan, including those with PAH associated with CHD (CHD-PAH). Methods OPUS was a prospective, United States, multicenter, long-term, observational drug registry (April 2014–June 2020). OrPHeUS was a retrospective, United States, multicenter medical chart review (October 2013–March 2017). The characteristics, treatment patterns, safety, and outcomes during macitentan treatment of patients with CHD-PAH and the subgroups Eisenmenger syndrome, left-to-right shunts and small/coincidental CHD were descriptively compared. Results The combined OPUS/OrPHeUS population included 272 (6.1%) patients with CHD-PAH (80 patients with Eisenmenger syndrome; 82 patients with left-to-right shunts and 92 patients with small/coincidental defects). Most patients across the CHD-PAH subgroups were in World Health Organization Functional Class II/III (82.9–94.6%). Macitentan was initiated as combination therapy in 65.0% of patients with CHD-PAH. During follow-up, 81.4% of patients experienced ≥ 1 adverse event (AE), the most common being dyspnea (23.5%), nausea (13.7%), dizziness (12.7%), headache (12.7%) and edema (10.8%). The 1- and 2-year Kaplan–Meier (95% confidence limits) estimates of patients with CHD-PAH being free from hospitalization were 64.5% (57.9, 70.4) and 49.3% (41.9, 56.3); for survival, the 1- and 2-year Kaplan–Meier estimates were 93.5% (89.3, 96.1) and 90.2% (84.9, 93.7). Conclusions Macitentan was used in clinical practice in patients with CHD-PAH and its subgroups, including as combination therapy in the majority of patients. Safety in this population was consistent with the known profile of macitentan. Trial Registration OPsumit® Users Registry (OPUS): NCT02126943; Opsumit® Historical Users cohort (OrPHeUS): NCT03197688; URL www.clinicaltrials.gov .

In the Hemispheric Transport of Air Pollution Phase 2 (HTAP2) exercise, a range of global atmospheric general circulation and chemical transport models performed coordinated perturbation experiments with 20 % reductions in emissions of anthropogenic aerosols, or aerosol precursors, in a number of source regions. Here, we compare the resulting changes in the atmospheric load and vertically resolved profiles of black carbon (BC), organic aerosols (OA) and sulfate (SO4) from 10 models that include treatment of aerosols. We use a set of temporally, horizontally and vertically resolved profiles of aerosol forcing efficiency (AFE) to estimate the impact of emission changes in six major source regions on global radiative forcing (RF) pertaining to the direct aerosol effect, finding values between. 51.9 and 210.8 mW m−2 Tg−1 for BC, between −2.4 and −17.9 mW m−2 Tg−1 for OA and between −3.6 and −10.3 W m−2 Tg−1 for SO4. In most cases, the local influence dominates, but results show that mitigations in south and east Asia have substantial impacts on the radiative budget in all investigated receptor regions, especially for BC. In Russia and the Middle East, more than 80 % of the forcing for BC and OA is due to extra-regional emission reductions. Similarly, for North America, BC emissions control in east Asia is found to be more important than domestic mitigations, which is consistent with previous findings. Comparing fully resolved RF calculations to RF estimates based on vertically averaged AFE profiles allows us to quantify the importance of vertical resolution to RF estimates. We find that locally in the source regions, a 20 % emission reduction strengthens the radiative forcing associated with SO4 by 25 % when including the vertical dimension, as the AFE for SO4 is strongest near the surface. Conversely, the local RF from BC weakens by 37 % since BC AFE is low close to the ground. The fraction of BC direct effect forcing attributable to intercontinental transport, on the other hand, is enhanced by one-third when accounting for the vertical aspect, because long-range transport primarily leads to aerosol changes at high altitudes, where the BC AFE is strong. While the surface temperature response may vary with the altitude of aerosol change, the analysis in the present study is not extended to estimates of temperature or precipitation changes.

Accurate modeling of the scattering and absorption of ultraviolet and visible radiation by aerosols is essential for accurate simulations of atmospheric chemistry and climate. Closure studies using in situ measurements of aerosol scattering and absorption can be used to evaluate and improve models of aerosol optical properties without interference from model errors in aerosol emissions, transport, chemistry, or deposition rates. Here we evaluate the ability of four externally mixed, fixed size distribution parameterizations used in global models to simulate submicron aerosol scattering and absorption at three wavelengths using in situ data gathered during the 2008 Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) campaign. The four models are the NASA Global Modeling Initiative (GMI) Combo model, GEOS-Chem v9-02, the baseline configuration of a version of GEOS-Chem with online radiative transfer calculations (called GC-RT), and the Optical Properties of Aerosol and Clouds (OPAC v3.1) package. We also use the ARCTAS data to perform the first evaluation of the ability of the Aerosol Simulation Program (ASP v2.1) to simulate submicron aerosol scattering and absorption when in situ data on the aerosol size distribution are used, and examine the impact of different mixing rules for black carbon (BC) on the results. We find that the GMI model tends to overestimate submicron scattering and absorption at shorter wavelengths by 10–23 %, and that GMI has smaller absolute mean biases for submicron absorption than OPAC v3.1, GEOS-Chem v9-02, or GC-RT. However, the changes to the density and refractive index of BC in GC-RT improve the simulation of submicron aerosol absorption at all wavelengths relative to GEOS-Chem v9-02. Adding a variable size distribution, as in ASP v2.1, improves model performance for scattering but not for absorption, likely due to the assumption in ASP v2.1 that BC is present at a constant mass fraction throughout the aerosol size distribution. Using a core-shell mixing rule in ASP overestimates aerosol absorption, especially for the fresh biomass burning aerosol measured in ARCTAS-B, suggesting the need for modeling the time-varying mixing states of aerosols in future versions of ASP.

The NOAA National Centers for Environmental Prediction (NCEP) implemented the NOAA Environmental Modeling System (NEMS) Global Forecast System (GFS) Aerosol Component (NGAC) for global dust forecasting in collaboration with NASA Goddard Space Flight Center (GSFC). NGAC Version 1.0 has been providing 5-day dust forecasts at 1° × 1° resolution on a global scale, once per day at 00:00 Coordinated Universal Time (UTC), since September 2012. This is the first global system capable of interactive atmosphere aerosol forecasting at NCEP. The implementation of NGAC V1.0 reflects an effective and efficient transitioning of NASA research advances to NCEP operations, paving the way for NCEP to provide global aerosol products serving a wide range of stakeholders, as well as to allow the effects of aerosols on weather forecasts and climate prediction to be considered.

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.

Atmospheric pollution over South Asia attracts special attention due to its effects on regional climate, water cycle and human health. These effects are potentially growing owing to rising trends of anthropogenic aerosol emissions. In this study, the spatio-temporal aerosol distributions over South Asia from seven global aerosol models are evaluated against aerosol retrievals from NASA satellite sensors and ground-based measurements for the period of 2000–2007. Overall, substantial underestimations of aerosol loading over South Asia are found systematically in most model simulations. Averaged over the entire South Asia, the annual mean aerosol optical depth (AOD) is underestimated by a range 15 to 44% across models compared to MISR (Multi-angle Imaging SpectroRadiometer), which is the lowest bound among various satellite AOD retrievals (from MISR, SeaWiFS (Sea-Viewing Wide Field-of-View Sensor), MODIS (Moderate Resolution Imaging Spectroradiometer) Aqua and Terra). In particular during the post-monsoon and wintertime periods (i.e., October–January), when agricultural waste burning and anthropogenic emissions dominate, models fail to capture AOD and aerosol absorption optical depth (AAOD) over the Indo–Gangetic Plain (IGP) compared to ground-based Aerosol Robotic Network (AERONET) sunphotometer measurements. The underestimations of aerosol loading in models generally occur in the lower troposphere (below 2 km) based on the comparisons of aerosol extinction profiles calculated by the models with those from Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) data. Furthermore, surface concentrations of all aerosol components (sulfate, nitrate, organic aerosol (OA) and black carbon (BC)) from the models are found much lower than in situ measurements in winter. Several possible causes for these common problems of underestimating aerosols in models during the post-monsoon and wintertime periods are identified: the aerosol hygroscopic growth and formation of secondary inorganic aerosol are suppressed in the models because relative humidity (RH) is biased far too low in the boundary layer and thus foggy conditions are poorly represented in current models, the nitrate aerosol is either missing or inadequately accounted for, and emissions from agricultural waste burning and biofuel usage are too low in the emission inventories. These common problems and possible causes found in multiple models point out directions for future model improvements in this important region.

The Multi-angle Imaging SpectroRadiometer (MISR) aboard the NASA Earth Observing System's Terra satellite can provide more reliable aerosol optical depth (AOD) and better constraints on particle size (Ångström exponent, or ANG), sphericity, and single-scattering albedo (SSA) than many other satellite instruments. However, many aerosol mixtures pass the algorithm acceptance criteria, yielding a poor constraint, when the particle-type information in the MISR radiances is low, typically at low AOD. We investigate adding value to the MISR aerosol product under these conditions by filtering the list of MISR-retrieved mixtures based on agreement between the mixture ANG and absorbing AOD (AAOD) values, and simulated aerosol properties from the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model. MISR–GOCART ANG difference and AAOD ratio thresholds for applying GOCART constraints were determined based on coincident AOD, ANG, and AAOD measurements from the AErosol RObotic NETwork (AERONET). The results were validated by comparing the adjusted MISR aerosol optical properties over the contiguous USA between 2006 and 2009 with additional AERONET data. The correlation coefficient (r) between the adjusted MISR ANG derived from this study and AERONET improves to 0.45, compared to 0.29 for the MISR Version 22 standard product. The ratio of the adjusted MISR AAOD to AERONET increases to 0.74, compared to 0.5 for the MISR operational retrieval. These improvements occur primarily when AOD < 0.2 for ANG and AOD < 0.5 for AAOD. Spatial and temporal differences among the aerosol optical properties of MISR V22, GOCART, and the adjusted MISR are traced to (1) GOCART underestimation of AOD and ANG in polluted regions; (2) aerosol mixtures lacking in the MISR Version 22 algorithm climatology; (3) low MISR sensitivity to particle type under some conditions; and (4) parameters and thresholds used in our method.

The multi-wavelength lidar technique was applied to the study of a smoke event near Washington, DC on 26–28 August 2013. Satellite observations combined with transport model predictions imply that the smoke plume originated mainly from Wyoming/Idaho forest fires and its transportation to Washington, DC took approximately 5 days. The NASA Goddard Space Flight Center (GSFC) multi-wavelength Mie–Raman lidar was used to measure the smoke particle intensive parameters such as extinction and backscatter Ångström exponents together with lidar ratios at 355 and 532 nm wavelengths. For interpretation of the observed vertical profiles of the backscatter Ångström exponents γβ at 355–532 and 532–1064 nm, numerical simulation was performed. The results indicate that, for fine-mode dominant aerosols, the Ångström exponents γβ(355–532) and γβ(532–1064) have essentially different dependence on the particle size and refractive index. Inversion of 3 β + 2 α lidar observations on 27–28 August provided vertical variation of the particle volume, effective radius and the real part of the refractive index through the planetary boundary layer (PBL) and the smoke layer. The particle effective radius decreased with height from approximately 0.27 μm inside the PBL to 0.15 μm in the smoke layer, which was situated above the PBL. Simultaneously the real part of the refractive index in the smoke layer increased to mR &approx; 1.5. The retrievals demonstrate also that the fine mode is predominant in the particle size distribution, and that the decrease of the effective radius with height is due to a shift of the fine mode toward smaller radii.

Atmospheric black carbon (BC) absorbs solar radiation, and exacerbates global warming through exerting positive radiative forcing (RF). However, the contribution of BC to ongoing changes in global climate is under debate. Anthropogenic BC emissions, and the resulting distribution of BC concentration, are highly uncertain. In particular, long-range transport and processes affecting BC atmospheric lifetime are poorly understood. Here we discuss whether recent assessments may have overestimated present-day BC radiative forcing in remote regions. We compare vertical profiles of BC concentration from four recent aircraft measurement campaigns to simulations by 13 aerosol models participating in the AeroCom Phase II intercomparison. An atmospheric lifetime of BC of less than 5 days is shown to be essential for reproducing observations in remote ocean regions, in line with other recent studies. Adjusting model results to measurements in remote regions, and at high altitudes, leads to a 25% reduction in AeroCom Phase II median direct BC forcing, from fossil fuel and biofuel burning, over the industrial era. The sensitivity of modelled forcing to BC vertical profile and lifetime highlights an urgent need for further flight campaigns, close to sources and in remote regions, to provide improved quantification of BC effects for use in climate policy.

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 variations and trends over different land and ocean regions from 1980 to 2009 are analyzed with the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model and observations from multiple satellite sensors and available ground-based networks. Excluding time periods with large volcanic influence, aerosol optical depth (AOD) and surface concentration over polluted land regions generally vary with anthropogenic emissions, but the magnitude of this association can be dampened by the presence of natural aerosols, especially dust. Over the 30-year period in this study, the largest reduction in aerosol levels occurs over Europe, where AOD has decreased by 40–60% on average and surface sulfate concentrations have declined by a factor of up to 3–4. In contrast, East Asia and South Asia show AOD increases, but the relatively high level of dust aerosols in Asia reduces the correlation between AOD and pollutant emission trends. Over major dust source regions, model analysis indicates that the change of dust emissions over the Sahara and Sahel has been predominantly driven by the change of near-surface wind speed, but over Central Asia it has been largely influenced by the change of the surface wetness. The decreasing dust trend in the North African dust outflow region of the tropical North Atlantic and the receptor sites of Barbados and Miami is closely associated with an increase of the sea surface temperature in the North Atlantic. This temperature increase may drive the decrease of the wind velocity over North Africa, which reduces the dust emission, and the increase of precipitation over the tropical North Atlantic, which enhances dust removal during transport. Despite significant trends over some major continental source regions, the model-calculated global annual average AOD shows little change over land and ocean in the past three decades, because opposite trends in different land regions cancel each other out in the global average, and changes over large open oceans are negligible. This highlights the necessity for regional-scale assessment of aerosols and their climate impacts, as global-scale average values can obscure important regional changes.

Though many global aerosols models prognose surface deposition, only a few models have been used to directly simulate the radiative effect from black carbon (BC) deposition to snow and sea ice. Here, we apply aerosol deposition fields from 25 models contributing to two phases of the Aerosol Comparisons between Observations and Models (AeroCom) project to simulate and evaluate within-snow BC concentrations and radiative effect in the Arctic. We accomplish this by driving the offline land and sea ice components of the Community Earth System Model with different deposition fields and meteorological conditions from 2004 to 2009, during which an extensive field campaign of BC measurements in Arctic snow occurred. We find that models generally underestimate BC concentrations in snow in northern Russia and Norway, while overestimating BC amounts elsewhere in the Arctic. Although simulated BC distributions in snow are poorly correlated with measurements, mean values are reasonable. The multi-model mean (range) bias in BC concentrations, sampled over the same grid cells, snow depths, and months of measurements, are −4.4 (−13.2 to +10.7) ng g−1 for an earlier phase of AeroCom models (phase I), and +4.1 (−13.0 to +21.4) ng g−1 for a more recent phase of AeroCom models (phase II), compared to the observational mean of 19.2 ng g−1. Factors determining model BC concentrations in Arctic snow include Arctic BC emissions, transport of extra-Arctic aerosols, precipitation, deposition efficiency of aerosols within the Arctic, and meltwater removal of particles in snow. Sensitivity studies show that the model–measurement evaluation is only weakly affected by meltwater scavenging efficiency because most measurements were conducted in non-melting snow. The Arctic (60–90° N) atmospheric residence time for BC in phase II models ranges from 3.7 to 23.2 days, implying large inter-model variation in local BC deposition efficiency. Combined with the fact that most Arctic BC deposition originates from extra-Arctic emissions, these results suggest that aerosol removal processes are a leading source of variation in model performance. The multi-model mean (full range) of Arctic radiative effect from BC in snow is 0.15 (0.07–0.25) W m−2 and 0.18 (0.06–0.28) W m−2 in phase I and phase II models, respectively. After correcting for model biases relative to observed BC concentrations in different regions of the Arctic, we obtain a multi-model mean Arctic radiative effect of 0.17 W m−2 for the combined AeroCom ensembles. Finally, there is a high correlation between modeled BC concentrations sampled over the observational sites and the Arctic as a whole, indicating that the field campaign provided a reasonable sample of the Arctic.

The land surface plays a crucial role in regulating water and energy fluxes at the land–atmosphere (L–A) interface and controls many processes and feedbacks in the climate system. Land cover and vegetation type remains one key determinant of soil moisture content that impacts air temperature, planetary boundary layer (PBL) evolution, and precipitation through soil-moisture–evapotranspiration coupling. In turn, it will affect atmospheric chemistry and air quality. This paper presents the results of a modeling study of the effect of land cover on some key L–A processes with a focus on air quality. The newly developed NASA Unified Weather Research and Forecast (NU-WRF) modeling system couples NASA's Land Information System (LIS) with the community WRF model and allows users to explore the L–A processes and feedbacks. Three commonly used satellite-derived land cover datasets – i.e., from the US Geological Survey (USGS) and University of Maryland (UMD), which are based on the Advanced Very High Resolution Radiometer (AVHRR), and from the Moderate Resolution Imaging Spectroradiometer (MODIS) – bear large differences in agriculture, forest, grassland, and urban spatial distributions in the continental United States, and thus provide an excellent case to investigate how land cover change would impact atmospheric processes and air quality. The weeklong simulations demonstrate the noticeable differences in soil moisture/temperature, latent/sensible heat flux, PBL height, wind, NO2/ozone, and PM2.5 air quality. These discrepancies can be traced to associate with the land cover properties, e.g., stomatal resistance, albedo and emissivity, and roughness characteristics. It also implies that the rapid urban growth may have complex air quality implications with reductions in peak ozone but more frequent high ozone events.

We use the NASA GEOS-5 transport model with tagged tracers to investigate the contributions of different regional sources of CO and black carbon (BC) to their concentrations in the Western Arctic (i.e., 50–90° N and 190–320° E) in spring and summer 2008. The model is evaluated by comparing the results with airborne measurements of CO and BC from the NASA Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) field campaigns to demonstrate the strengths and limitations of our simulations. We also examine the reliability of tagged CO tracers in characterizing air mass origins using the measured fossil fuel tracer of dichloromethane and the biomass burning tracer of acetonitrile. Our tagged CO simulations suggest that most of the enhanced CO concentrations (above background level from CH4 production) observed during April originate from Asian anthropogenic emissions. Boreal biomass burning emissions and Asian anthropogenic emissions are of similar importance in July domain wise, although the biomass burning CO fraction is much larger in the area of the ARCTAS field experiments. The fraction of CO from Asian anthropogenic emissions is larger in spring than in summer. European sources make up no more than 10% of CO levels in the campaign domain during either period. Comparisons of CO concentrations along the flight tracks with regional averages from GEOS-5 show that the along-track measurements are representative of the concentrations within the large domain of the Western Arctic in April but not in July.

The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) examined the short-lived drivers of climate change in current climate models. Here we evaluate the 10 ACCMIP models that included aerosols, 8 of which also participated in the Coupled Model Intercomparison Project phase 5 (CMIP5). The models reproduce present-day total aerosol optical depth (AOD) relatively well, though many are biased low. Contributions from individual aerosol components are quite different, however, and most models underestimate east Asian AOD. The models capture most 1980–2000 AOD trends well, but underpredict increases over the Yellow/Eastern Sea. They strongly underestimate absorbing AOD in many regions. We examine both the direct radiative forcing (RF) and the forcing including rapid adjustments (effective radiative forcing; ERF, including direct and indirect effects). The models' all-sky 1850 to 2000 global mean annual average total aerosol RF is (mean; range) −0.26 W m−2; −0.06 to −0.49 W m−2. Screening based on model skill in capturing observed AOD yields a best estimate of −0.42 W m−2; −0.33 to −0.50 W m−2, including adjustment for missing aerosol components in some models. Many ACCMIP and CMIP5 models appear to produce substantially smaller aerosol RF than this best estimate. Climate feedbacks contribute substantially (35 to −58%) to modeled historical aerosol RF. The 1850 to 2000 aerosol ERF is −1.17 W m−2; −0.71 to −1.44 W m−2. Thus adjustments, including clouds, typically cause greater forcing than direct RF. Despite this, the multi-model spread relative to the mean is typically the same for ERF as it is for RF, or even smaller, over areas with substantial forcing. The largest 1850 to 2000 negative aerosol RF and ERF values are over and near Europe, south and east Asia and North America. ERF, however, is positive over the Sahara, the Karakoram, high Southern latitudes and especially the Arctic. Global aerosol RF peaks in most models around 1980, declining thereafter with only weak sensitivity to the Representative Concentration Pathway (RCP). One model, however, projects approximately stable RF levels, while two show increasingly negative RF due to nitrate (not included in most models). Aerosol ERF, in contrast, becomes more negative during 1980 to 2000. During this period, increased Asian emissions appear to have a larger impact on aerosol ERF than European and North American decreases due to their being upwind of the large, relatively pristine Pacific Ocean. There is no clear relationship between historical aerosol ERF and climate sensitivity in the CMIP5 subset of ACCMIP models. In the ACCMIP/CMIP5 models, historical aerosol ERF of about −0.8 to −1.5 W m−2 is most consistent with observed historical warming. Aerosol ERF masks a large portion of greenhouse forcing during the late 20th and early 21st century at the global scale. Regionally, aerosol ERF is so large that net forcing is negative over most industrialized and biomass burning regions through 1980, but remains strongly negative only over east and southeast Asia by 2000. Net forcing is strongly positive by 1980 over most deserts, the Arctic, Australia, and most tropical oceans. Both the magnitude of and area covered by positive forcing expand steadily thereafter.

Simulated multi-model "diversity" in aerosol direct radiative forcing estimates is often perceived as a measure of aerosol uncertainty. However, current models used for aerosol radiative forcing calculations vary considerably in model components relevant for forcing calculations and the associated "host-model uncertainties" are generally convoluted with the actual aerosol uncertainty. In this AeroCom Prescribed intercomparison study we systematically isolate and quantify host model uncertainties on aerosol forcing experiments through prescription of identical aerosol radiative properties in twelve participating models. Even with prescribed aerosol radiative properties, simulated clear-sky and all-sky aerosol radiative forcings show significant diversity. For a purely scattering case with globally constant optical depth of 0.2, the global-mean all-sky top-of-atmosphere radiative forcing is −4.47 Wm−2 and the inter-model standard deviation is 0.55 Wm−2, corresponding to a relative standard deviation of 12%. For a case with partially absorbing aerosol with an aerosol optical depth of 0.2 and single scattering albedo of 0.8, the forcing changes to 1.04 Wm−2, and the standard deviation increases to 1.01 W−2, corresponding to a significant relative standard deviation of 97%. However, the top-of-atmosphere forcing variability owing to absorption (subtracting the scattering case from the case with scattering and absorption) is low, with absolute (relative) standard deviations of 0.45 Wm−2 (8%) clear-sky and 0.62 Wm−2 (11%) all-sky. Scaling the forcing standard deviation for a purely scattering case to match the sulfate radiative forcing in the AeroCom Direct Effect experiment demonstrates that host model uncertainties could explain about 36% of the overall sulfate forcing diversity of 0.11 Wm−2 in the AeroCom Direct Radiative Effect experiment. Host model errors in aerosol radiative forcing are largest in regions of uncertain host model components, such as stratocumulus cloud decks or areas with poorly constrained surface albedos, such as sea ice. Our results demonstrate that host model uncertainties are an important component of aerosol forcing uncertainty that require further attention.

We report on the AeroCom Phase II direct aerosol effect (DAE) experiment where 16 detailed global aerosol models have been used to simulate the changes in the aerosol distribution over the industrial era. All 16 models have estimated the radiative forcing (RF) of the anthropogenic DAE, and have taken into account anthropogenic sulphate, black carbon (BC) and organic aerosols (OA) from fossil fuel, biofuel, and biomass burning emissions. In addition several models have simulated the DAE of anthropogenic nitrate and anthropogenic influenced secondary organic aerosols (SOA). The model simulated all-sky RF of the DAE from total anthropogenic aerosols has a range from −0.58 to −0.02 Wm−2, with a mean of −0.27 Wm−2 for the 16 models. Several models did not include nitrate or SOA and modifying the estimate by accounting for this with information from the other AeroCom models reduces the range and slightly strengthens the mean. Modifying the model estimates for missing aerosol components and for the time period 1750 to 2010 results in a mean RF for the DAE of −0.35 Wm−2. Compared to AeroCom Phase I (Schulz et al., 2006) we find very similar spreads in both total DAE and aerosol component RF. However, the RF of the total DAE is stronger negative and RF from BC from fossil fuel and biofuel emissions are stronger positive in the present study than in the previous AeroCom study. We find a tendency for models having a strong (positive) BC RF to also have strong (negative) sulphate or OA RF. This relationship leads to smaller uncertainty in the total RF of the DAE compared to the RF of the sum of the individual aerosol components. The spread in results for the individual aerosol components is substantial, and can be divided into diversities in burden, mass extinction coefficient (MEC), and normalized RF with respect to AOD. We find that these three factors give similar contributions to the spread in results.

The challenge of using satellite observations to retrieve aerosol properties in a cloudy environment is to prevent contamination of the aerosol signal from clouds, while maintaining sufficient aerosol product yield to satisfy specific applications. We investigate aerosol retrieval availability at different instrument pixel resolutions using the standard MODIS aerosol cloud mask applied to MODIS data and supplemented with a new GOES-R cloud mask applied to GOES data for a domain covering North America and surrounding oceans. Aerosol product availability is not the same as the cloud free fraction and takes into account the techniques used in the MODIS algorithm to avoid clouds, reduce noise and maintain sufficient numbers of aerosol retrievals. The inherent spatial resolution of each instrument, 0.5×0.5 km for MODIS and 1×1 km for GOES, is systematically degraded to 1×1, 2×2, 1×4, 4×4 and 8×8 km resolutions and then analyzed as to how that degradation would affect the availability of an aerosol retrieval, assuming an aerosol product resolution at 8×8 km. The analysis is repeated, separately, for near-nadir pixels and those at larger view angles to investigate the effect of pixel growth at oblique angles on aerosol retrieval availability. The results show that as nominal pixel size increases, availability decreases until at 8×8 km 70% to 85% of the retrievals available at 0.5 km, nadir, have been lost. The effect at oblique angles is to further decrease availability over land but increase availability over ocean, because sun glint is found at near-nadir view angles. Finer resolution sensors (i.e., 1×1, 2×2 or even 1×4 km) will retrieve aerosols in partly cloudy scenes significantly more often than sensors with nadir views of 4×4 km or coarser. Large differences in the results of the two cloud masks designed for MODIS aerosol and GOES cloud products strongly reinforce that cloud masks must be developed with specific purposes in mind and that a generic cloud mask applied to an independent aerosol retrieval will likely fail.

We use aircraft observations obtained during the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) mission to examine the distributions and source attributions of O3 and NOy in the Arctic and sub-Arctic region. Using a number of marker tracers, we distinguish various air masses from the background troposphere and examine their contributions to NOx, O3, and O3 production in the Arctic troposphere. The background Arctic troposphere has a mean O3 of ~60 ppbv and NOx of ~25 pptv throughout spring and summer with CO decreasing from ~145 ppbv in spring to ~100 ppbv in summer. These observed mixing ratios are not notably different from the values measured during the 1988 ABLE-3A and the 2002 TOPSE field campaigns despite the significant changes in emissions and stratospheric ozone layer in the past two decades that influence Arctic tropospheric composition. Air masses associated with stratosphere-troposphere exchange are present throughout the mid and upper troposphere during spring and summer. These air masses, with mean O3 concentrations of 140–160 ppbv, are significant direct sources of O3 in the Arctic troposphere. In addition, air of stratospheric origin displays net O3 formation in the Arctic due to its sustainable, high NOx (75 pptv in spring and 110 pptv in summer) and NOy (~800 pptv in spring and ~1100 pptv in summer). The air masses influenced by the stratosphere sampled during ARCTAS-B also show conversion of HNO3 to PAN. This active production of PAN is the result of increased degradation of ethane in the stratosphere-troposphere mixed air mass to form CH3CHO, followed by subsequent formation of PAN under high NOx conditions. These findings imply that an adequate representation of stratospheric NOy input, in addition to stratospheric O3 influx, is essential to accurately simulate tropospheric Arctic O3, NOx and PAN in chemistry transport models. Plumes influenced by recent anthropogenic and biomass burning emissions observed during ARCTAS show highly elevated levels of hydrocarbons and NOy (mostly in the form of NOx and PAN), but do not contain O3 higher than that in the Arctic tropospheric background except some aged biomass burning plumes sampled during spring. Convection and/or lightning influences are negligible sources of O3 in the Arctic troposphere but can have significant impacts in the upper troposphere in the continental sub-Arctic during summer.