Evaluation of seven European aerosol optical depth retrieval algorithms for climate analysis

Satellite data are increasingly used to provide observation-based estimates of the effects of aerosols on climate. The Aerosol-cci project, part of the European Space Agency's Climate Change Initiative (CCI), was designed to provide essential climate variables for aerosols from satellite data. Eight algorithms, developed for the retrieval of aerosol properties using data from AATSR (4), MERIS (3) and POLDER, were evaluated to determine their suitability for climate studies. The primary result from each of these algorithms is the aerosol optical depth (AOD) at several wavelengths, together with the Ångström exponent (AE) which describes the spectral variation of the AOD for a given wavelength pair. Other aerosol parameters which are possibly retrieved from satellite observations are not considered in this paper. The AOD and AE (AE only for Level 2) were evaluated against independent collocated observations from the ground-based AERONET sun photometer network and against “reference” satellite data provided by MODIS and MISR. Tools used for the evaluation were developed for daily products as produced by the retrieval with a spatial resolution of 10 × 10 km2 (Level 2) and daily or monthly aggregates (Level 3). These tools include statistics for L2 and L3 products compared with AERONET, as well as scoring based on spatial and temporal correlations. In this paper we describe their use in a round robin (RR) evaluation of four months of data, one month for each season in 2008. The amount of data was restricted to only four months because of the large effort made to improve the algorithms, and to evaluate the improvement and current status, before larger data sets will be processed. Evaluation criteria are discussed. Results presented show the current status of the European aerosol algorithms in comparison to both AERONET and MODIS and MISR data. The comparison leads to a preliminary conclusion that the scores are similar, including those for the references, but the coverage of AATSR needs to be enhanced and further improvements are possible for most algorithms. None of the algorithms, including the references, outperforms all others everywhere. AATSR data can be used for the retrieval of AOD and AE over land and ocean. PARASOL and one of the MERIS algorithms have been evaluated over ocean only and both algorithms provide good results

AeroCom phase III multi-model evaluation of the aerosol life cycle and optical properties using ground- and space-based remote sensing as well as surface in situ observations

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 (AOD(f), AOD(c)), Angstrom exponent (AE), dry surface scattering (SCdry), and absorption (AC(dry)) 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 spatiotemporal variations in AOD. We find a much larger underestimate in coarse AOD(c) (similar to-45% +/- 25 %) than in fine AOD(f) (similar to-15% +/- 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AOD(c) 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 AOD(c). 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 AC(dry), respectively. The larger underestimate of SCdry than AC(dry) 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 %.Peer reviewe

Development, Production and Evaluation of Aerosol Climate Data Records from European Satellite Observations (Aerosol_cci)

Producing a global and comprehensive description of atmospheric aerosols requires integration of ground-based, airborne, satellite and model datasets. Due to its complexity, aerosol monitoring requires the use of several data records with complementary information content. This paper describes the lessons learned while developing and qualifying algorithms to generate aerosol Climate Data Records (CDR) within the European Space Agency (ESA) Aerosol_cci project. An iterative algorithm development and evaluation cycle involving core users is applied. It begins with the application-specific refinement of user requirements, leading to algorithm development, dataset processing and independent validation followed by user evaluation. This cycle is demonstrated for a CDR of total Aerosol Optical Depth (AOD) from two subsequent dual-view radiometers. Specific aspects of its applicability to other aerosol algorithms are illustrated with four complementary aerosol datasets. An important element in the development of aerosol CDRs is the inclusion of several algorithms evaluating the same data to benefit from various solutions to the ill-determined retrieval problem. The iterative approach has produced a 17-year AOD CDR, a 10-year stratospheric extinction profile CDR and a 35-year Absorbing Aerosol Index record. Further evolution cycles have been initiated for complementary datasets to provide insight into aerosol properties (i.e., dust aerosol, aerosol absorption).Peer reviewe

Concentrations and radiative forcing of anthropogenic aerosols from 1750 to 2014 simulated with the Oslo CTM3 and CEDS emission inventory

We document the ability of the new-generation Oslo chemistry-transport model, Oslo CTM3, to accurately simulate present-day aerosol distributions. The model is then used with the new Community Emission Data System (CEDS) historical emission inventory to provide updated time series of anthropogenic aerosol concentrations and consequent direct radiative forcing (RFari) from 1750 to 2014. Overall, Oslo CTM3 performs well compared with measurements of surface concentrations and remotely sensed aerosol optical depth. Concentrations are underestimated in Asia, but the higher emissions in CEDS than previous inventories result in improvements compared to observations. The treatment of black carbon (BC) scavenging in Oslo CTM3 gives better agreement with observed vertical BC profiles relative to the predecessor Oslo CTM2. However, Arctic wintertime BC concentrations remain underestimated, and a range of sensitivity tests indicate that better physical understanding of processes associated with atmospheric BC processing is required to simultaneously reproduce both the observed features. Uncertainties in model input data, resolution, and scavenging affect the distribution of all aerosols species, especially at high latitudes and altitudes. However, we find no evidence of consistently better model performance across all observables and regions in the sensitivity tests than in the baseline configuration. Using CEDS, we estimate a net RFari in 2014 relative to 1750 of −0.17 W m−2, significantly weaker than the IPCC AR5 2011–1750 estimate. Differences are attributable to several factors, including stronger absorption by organic aerosol, updated parameterization of BC absorption, and reduced sulfate cooling. The trend towards a weaker RFari over recent years is more pronounced than in the IPCC AR5, illustrating the importance of capturing recent regional emission changes

Concentrations and radiative forcing of anthropogenic aerosols from 1750 to 2014 simulated with the Oslo CTM3 and CEDS emission inventory

We document the ability of the new-generation Oslo chemistry-transport model, Oslo CTM3, to accurately simulate present-day aerosol distributions. The model is then used with the new Community Emission Data System (CEDS) historical emission inventory to provide updated time series of anthropogenic aerosol concentrations and consequent direct radiative forcing (RFari) from 1750 to 2014. Overall, Oslo CTM3 performs well compared with measurements of surface concentrations and remotely sensed aerosol optical depth. Concentrations are underestimated in Asia, but the higher emissions in CEDS than previous inventories result in improvements compared to observations. The treatment of black carbon (BC) scavenging in Oslo CTM3 gives better agreement with observed vertical BC profiles relative to the predecessor Oslo CTM2. However, Arctic wintertime BC concentrations remain underestimated, and a range of sensitivity tests indicate that better physical understanding of processes associated with atmospheric BC processing is required to simultaneously reproduce both the observed features. Uncertainties in model input data, resolution, and scavenging affect the distribution of all aerosols species, especially at high latitudes and altitudes. However, we find no evidence of consistently better model performance across all observables and regions in the sensitivity tests than in the baseline configuration. Using CEDS, we estimate a net RFari in 2014 relative to 1750 of −0.17 W m−2, significantly weaker than the IPCC AR5 2011–1750 estimate. Differences are attributable to several factors, including stronger absorption by organic aerosol, updated parameterization of BC absorption, and reduced sulfate cooling. The trend towards a weaker RFari over recent years is more pronounced than in the IPCC AR5, illustrating the importance of capturing recent regional emission changes

Application of the CALIOP layer product to evaluate the vertical distribution of aerosols estimated by global models: AeroCom phase I results

The CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) Layer data are used to evaluate the vertical distribution of aerosols in 12 global models from AeroCom (Aerosol Comparisons between Observations and Models) simulations prepared for the 4th IPCC assessment report. Annual and seasonal extinction profiles are analysed over 13 sub-continental regions representative of industrial, dust, and biomass burning pollution, from CALIOP (2007-2009) observations and AeroCom (2000) simulations. An extinction mean height diagnostic (Zα) is also calculated over the 0-10 km altitude range to quantitatively assess the models performance. The mean extinction profiles derived from CALIOP layer products provide a robust climatology of the aerosol vertical distribution over large regions and for the seasonal time scale. The mean extinction decreases from the surface to about 6 km and shows consistent regional specificities and low inter-annual variability. While most models are significantly correlated with the observation in terms of Zα, some do better than others and 2 of the 12 models perform particularly well in all seasons. In most regions, notably over the maritime regions, downwind of the continents, most models show higher Zα than observed by CALIOP, whereas over the African and Chinese dust source regions, Zα is underestimated during northern hemisphere spring and summer. The positive model bias in Zα is in mainly due to an overestimate of the extinction above 6 km. Potential factors related to models and observations that contribute to the differences are discussed.JRC.H.2-Air and Climat

Overview of the Norwegian Earth System Model (NorESM2) and key climate response of CMIP6 DECK, historical, and scenario simulations

The second version of the coupled Norwegian Earth System Model (NorESM2) is presented and evaluated. NorESM2 is based on the second version of the Community Earth System Model (CESM2) and shares with CESM2 the computer code infrastructure and many Earth system model components. However, NorESM2 employs entirely different ocean and ocean biogeochemistry models. The atmosphere component of NorESM2 (CAM-Nor) includes a different module for aerosol physics and chemistry, including interactions with cloud and radiation; additionally, CAM-Nor includes improvements in the formulation of local dry and moist energy conservation, in local and global angular momentum conservation, and in the computations for deep convection and air–sea fluxes. The surface components of NorESM2 have minor changes in the albedo calculations and to land and sea-ice models. We present results from simulations with NorESM2 that were carried out for the sixth phase of the Coupled Model Intercomparison Project (CMIP6). Two versions of the model are used: one with lower (∼ 2∘) atmosphere–land resolution and one with medium (∼ 1∘) atmosphere–land resolution. The stability of the pre-industrial climate and the sensitivity of the model to abrupt and gradual quadrupling of CO2 are assessed, along with the ability of the model to simulate the historical climate under the CMIP6 forcings. Compared to observations and reanalyses, NorESM2 represents an improvement over previous versions of NorESM in most aspects. NorESM2 appears less sensitive to greenhouse gas forcing than its predecessors, with an estimated equilibrium climate sensitivity of 2.5 K in both resolutions on a 150-year time frame; however, this estimate increases with the time window and the climate sensitivity at equilibration is much higher. We also consider the model response to future scenarios as defined by selected Shared Socioeconomic Pathways (SSPs) from the Scenario Model Intercomparison Project defined under CMIP6. Under the four scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5), the warming in the period 2090–2099 compared to 1850–1879 reaches 1.3, 2.2, 3.0, and 3.9 K in NorESM2-LM, and 1.3, 2.1, 3.1, and 3.9 K in NorESM-MM, robustly similar in both resolutions. NorESM2-LM shows a rather satisfactory evolution of recent sea-ice area. In NorESM2-LM, an ice-free Arctic Ocean is only avoided in the SSP1-2.6 scenario

Overview of the Norwegian Earth System Model (NorESM2) and key climate response of CMIP6 DECK, historical, and scenario simulations

Abstract. The second version of the coupled Norwegian Earth System Model (NorESM2) is presented and evaluated. NorESM2 is based on the second version of the Community Earth System Model (CESM2) and shares with CESM2 the computer code infrastructure and many Earth system model components. However, NorESM2 employs entirely different ocean and ocean biogeochemistry models. The atmosphere component of NorESM2 (CAM-Nor) includes a different module for aerosol physics and chemistry, including interactions with cloud and radiation; additionally, CAM-Nor includes improvements in the formulation of local dry and moist energy conservation, in local and global angular momentum conservation, and in the computations for deep convection and air–sea fluxes. The surface components of NorESM2 have minor changes in the albedo calculations and to land and sea-ice models. We present results from simulations with NorESM2 that were carried out for the sixth phase of the Coupled Model Intercomparison Project (CMIP6). Two versions of the model are used: one with lower (∼ 2∘) atmosphere–land resolution and one with medium (∼ 1∘) atmosphere–land resolution. The stability of the pre-industrial climate and the sensitivity of the model to abrupt and gradual quadrupling of CO2 are assessed, along with the ability of the model to simulate the historical climate under the CMIP6 forcings. Compared to observations and reanalyses, NorESM2 represents an improvement over previous versions of NorESM in most aspects. NorESM2 appears less sensitive to greenhouse gas forcing than its predecessors, with an estimated equilibrium climate sensitivity of 2.5 K in both resolutions on a 150-year time frame; however, this estimate increases with the time window and the climate sensitivity at equilibration is much higher. We also consider the model response to future scenarios as defined by selected Shared Socioeconomic Pathways (SSPs) from the Scenario Model Intercomparison Project defined under CMIP6. Under the four scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5), the warming in the period 2090–2099 compared to 1850–1879 reaches 1.3, 2.2, 3.0, and 3.9 K in NorESM2-LM, and 1.3, 2.1, 3.1, and 3.9 K in NorESM-MM, robustly similar in both resolutions. NorESM2-LM shows a rather satisfactory evolution of recent sea-ice area. In NorESM2-LM, an ice-free Arctic Ocean is only avoided in the SSP1-2.6 scenario