A new global reanalysis data set of atmospheric composition (AC) for the period 2003–2015 has been produced by the Copernicus Atmosphere Monitoring Service (CAMS). Satellite observations of total column (TC) carbon monoxide (CO) and aerosol optical depth (AOD), as well as several TC and profile observations of ozone, have been assimilated with the Integrated Forecasting System for Composition (C-IFS) of the European Centre for Medium-Range Weather Forecasting. Compared to the previous Monitoring Atmospheric Composition and Climate (MACC) reanalysis (MACCRA), the new CAMS interim reanalysis (CAMSiRA) is of a coarser horizontal resolution of about 110 km, compared to 80 km, but covers a longer period with the intent to be continued to present day. This paper compares CAMSiRA with MACCRA and a control run experiment (CR) without assimilation of AC retrievals. CAMSiRA has smaller biases than the CR with respect to independent observations of CO, AOD and stratospheric ozone. However, ozone at the surface could not be improved by the assimilation because of the strong impact of surface processes such as dry deposition and titration with nitrogen monoxide (NO), which were both unchanged by the assimilation. The assimilation of AOD led to a global reduction of sea salt and desert dust as well as an exaggerated increase in sulfate. Compared to MACCRA, CAMSiRA had smaller biases for AOD, surface CO and TC ozone as well as for upper stratospheric and tropospheric ozone. Finally, the temporal consistency of CAMSiRA was better than the one of MACCRA. This was achieved by using a revised emission data set as well as by applying careful selection and bias correction to the assimilated retrievals. CAMSiRA is therefore better suited than MACCRA for the study of interannual variability, as demonstrated for trends in surface CO.
Exploiting the multitude of satellite observations of atmospheric composition (AC) is a key objective of the Copernicus Atmosphere Monitoring Service (CAMS). For its global component CAMS uses the four-dimensional variational (4D-VAR) data assimilation technique to combine satellite observations with chemistry-aerosol modelling to obtain a gridded continuous representation (analysis) of the mass mixing ratios of atmospheric trace gases and aerosols.
The global CAMS system is built on the heritage of the EU-funded GEMS project (Hollingsworth et al., 2008) and a series of Monitoring Atmospheric Composition and Climate (MACC) projects at the European Centre for Medium-Range Weather Forecasts (ECMWF). During these projects the Integrated Forecasting System (IFS) of ECMWF was extended by modules for atmospheric chemistry, aerosols and greenhouse gases in such a way that the 4D-VAR data assimilation system, which had been developed for the analysis of the meteorological fields, could be used for the assimilation of AC retrievals.
Assimilating satellite AC retrievals into an AC model has advantages over the sole use of the AC retrievals because of their specific limitations. First, only a small subset of the trace gases or only total aerosol is directly observable with an accuracy sufficient to have an impact during the assimilation. Second, AC satellite retrievals have incomplete horizontal coverage because of the orbital cycle, viewing geometry, the presence of clouds and other factors such as surface albedo. Third, the vertical distribution of the trace species can often not or only coarsely be retrieved from the satellite observations, while the measurement sensitivity towards the surface is generally low.
The AC analyses are used (i) to initialise AC model forecasts and (ii) for the retrospective analysis (reanalysis) of AC for air quality and climate studies. The reanalysis of the meteorological fields has been an important activity at ECMWF (ERA-40, Uppala et al., 2005; ERA interim, Dee et al., 2011) and other meteorological centres such as National Centers for Environmental Protection (NCEP) (CFSR; Saha et al., 2010, JMA (JRA-55, JRA-25; Onogi et al., 2007)) and NASA-DAO (MERRA; Rienecker, et al., 2011). An important application of these reanalysis data sets is the estimation of the interannual variability and the trends of climate variables over the last decades up to the present day. The complete spatial and temporal coverage makes the trend analysis of reanalyses more robust and universal than the trend analysis of individual observing systems. However, constructing a data set, which is suited for this purpose, is a complex task because of the developing and changing observing system, which can introduce spurious trends and sudden shifts in the reanalysis data record. Careful quality control of the assimilated observations and techniques (e.g. Dee et al., 2004) to address inter-instrument biases are applied to mitigate this problem.
Most meteorological reanalyses contain stratospheric ozone, but other trace
gases, apart from water vapour, are not included. In the last decade chemical
and aerosol data assimilation has matured (Bocquet et al., 2015) and
dedicated reanalysis data sets for AC have emerged. The
multi-sensor reanalysis of total ozone (van der A et al., 2015) for
1970–2012 used ground-based Brewer observations to inter-calibrate satellite
retrievals. The MERRAero reanalysis (2002–present,
The MACC reanalysis (MACCRA) of reactive gases (Inness et al., 2013) and aerosols for the period 2003–2012 is an AC reanalysis that covers tropospheric and stratospheric reactive gases and aerosols as well as the meteorological fields in one consistent data set. MACCRA has proved to be a realistic data set as shown in several evaluation studies for reactive gases (Elguindi et al., 2010; Inness et al., 2013; Katragkou et al., 2015; Gaudel et al., 2015) and aerosols (Cesnulyte et al., 2014; Cuevas et al., 2015). MACCRA is widely used, for example, as a boundary condition for regional models (Schere et al., 2012; Im et al., 2015; Giordano et al., 2015), to construct trace gas climatologies for the IFS radiation schemes (Bechtold et al., 2009), to estimate aerosol radiative forcing (Bellouin et al., 2013), as input for solar radiation schemes for solar energy applications, and to report the current state of aerosol and CO as part of the climate system (Benedetti et al., 2014; Flemming and Inness, 2014).
CAMS is committed to producing a comprehensive high-resolution AC reanalysis
in the next years. The CAMS interim Reanalysis (CAMSiRA) presented here has
an interim status between MACCRA and this planned analysis data set. It was
produced at a lower horizontal resolution (110 km) than the resolution of
MACCRA (80 km), and the number of archived AC fields was limited to the
aerosol variables and selected chemical species such as ozone, HNO
The reasons for producing CAMSiRA before the more comprehensive reanalysis are as follows: the MACCRA for reactive gases was produced using a coupled system consisting of the IFS and the MOZART-3 (Kinnison et al., 2007) CTM as described in Flemming et al. (2009). This coupled system was replaced by the much more computationally efficient online coupled model C-IFS (Flemming et al., 2015), which uses the chemical mechanism CB05 of the TM5 CTM (Huijnen et al., 2010). With the discontinuation of the coupled system it was not possible to extend the MACC reanalysis to the present day. For the AC monitoring service of CAMS, it is however important to be able to compare the present conditions with previous years in a consistent way. Another motivation for producing CAMSiRA was that the aerosol module used for the MACCRA had undergone upgrades (Morcrette et al., 2011) in recent years. Finally, MACCRA suffered from small but noticeable shifts because of changes in the assimilated observations, the emission data and the bias correction approach. These spurious shifts undermine the usefulness of the MACCRA for the reliable estimation of trends. The lessons learnt from the evaluation of CAMSiRA will feed into the set-up of the planned CAMS reanalysis.
Reanalyses of AC are generally less well-constrained by observations than meteorological reanalyses because of the aforementioned limitations of the AC observations and because of the strong impact of the emissions, which are in many cases not constrained by observations. It is therefore good scientific practice to investigate the impact of the AC assimilation by comparing the AC reanalysis to a control experiment that did not assimilate AC observations. The control run (CR) for CAMSiRA was carried out using the same emission data as well as the meteorological fields produced by CAMSiRA.
The purpose of this paper is firstly to document the model system, the emissions, and the assimilated observations used to produce CAMSiRA and secondly to highlight its differences to the set-up of the MACCRA. Since the emissions are important drivers for variability of AC, a presentation of the totals and the interannual variability of the emission data used in CAMSiRA and the CR is given in a Supplement to the paper.
In the remainder of the paper, CO, aerosol as well as tropospheric and stratospheric ozone of CAMSiRA, the CR and MACCRA are inter-compared and evaluated with independent observations in a separate section for each species. The comparison of CAMSiRA with MACCRA has the purpose of reporting progress of and issues with CAMSiRA for potential users of the data sets. The comparison of CAMSiRA with the CR shows the impact of the data assimilation and is helpful to better understand deficiencies of the C-IFS model and its input data.
Each section starts with a discussion of the spatial differences of CAMSiRA, the CR and MACCRA for the considered species. Next, the temporal variability is investigated using time series of monthly mean values averaged over selected regions. We present global burdens and discuss changes in the speciation of the aerosol fields introduced by the assimilation. Finally, the three data sets are compared with independent observations, which were not used in the assimilation. A summary and recommendations for future AC reanalysis will be given in the last section.
CAMSiRA is a data set of 6-hourly reanalyses of AC for the period 2003–2015. A 3-hourly data set consistent with the AC analysis is available from forecasts linking the analyses. The horizontal resolution is about 110 km on a reduced Gaussian grid (T159) and the vertical discretization uses 60 levels from the surface to a model top of 0.1 hPa. Total columns of carbon monoxide (TC CO) form the Measurements Of Pollution In The Troposphere (MOPITT) instrument, MODIS AOD, and several ozone TC and stratospheric profile retrievals (see Table 2) were assimilated together with meteorological in situ and satellite observations.
The description of MACCRA for reactive gases can be found in Inness et al. (2013). Important commonalities and differences between the two AC reanalyses are given in Table 1.
Important commonalities and differences between MACCRA and CAMSiRA.
Assimilated satellite observations in CAMSiRA.
The control run is a forward simulation of C-IFS in monthly segments. The meteorological simulation is relaxed using the approach by Jung et al. (2008) to the meteorological reanalysis produced by the CAMSiRA. The emission input fields are the same as those used for CAMSiRA.
The model C-IFS is documented and evaluated in Flemming et al. (2015). C-IFS applies the chemical mechanism CB05, which describes tropospheric chemistry with 55 species and 126 reactions. Stratospheric ozone chemistry in C-IFS is parameterized by the “Cariolle-scheme” (Cariolle and Dèquè, 1986; Cariolle and Teyssèdre, 2007). Chemical tendencies for stratospheric and tropospheric ozone are merged at an empirical interface of the diagnosed tropopause height in C- IFS. C-IFS benefits from the detailed cloud and precipitation physics of the IFS for the calculation of wet deposition and lightning NO emission. Wet deposition modelling for the chemical species is based on Jacob (2000) and accounts for the subgrid scale distribution of clouds and precipitation. Dry deposition is modelled using precalculated monthly mean dry deposition velocities following Wesely (1989) with a superimposed diurnal cycle. Surface emissions and dry deposition fluxes are applied as surface boundary conditions of the diffusion scheme. Lightning emissions of NO were calculated based on convective precipitation (Meijer et al., 2001).
The aerosol module (Morcrette et al., 2009) is a bulk–bin scheme simulating
desert dust, sea salt at 80 % relative humidity (RH), hydrophilic and
hydrophobic organic carbon and black carbon as well as sulfate aerosol based
on the the LMDZ model of Laboratoire de Météorologie Dynamique aerosol model (Reddy et al., 2005). Sea salt and desert dust are
represented in three size bins. The radius ranges of the dust bins are
0.030–0.55, 0.55–0.9 and 0.9–20
The aerosol and chemistry modules used to simulate source and sink terms are not coupled. Wet and dry deposition are also modelled with different parameterisations but with the same meteorological input as precipitation fields. Aerosol and chemistry have in common that they are advected and vertically distributed by diffusion and convection in the same way. A proportional mass fixer as described in Diamantakis and Flemming (2014) is applied for all tracers in C-IFS.
This section only references the origin of the emission data. The emitted totals and the linear trends of the anthropogenic, biomass burning, and natural emissions as well as the modelled desert dust and sea salt emissions used in CAMSiRA and the CR are presented in a Supplement.
The anthropogenic surface emissions for the chemical species were taken from the MACCity inventory (Granier et al., 2011), which covers the period 1960–2010. MACCity emissions are based on the ACCMIP (Lamarque et al., 2013) inventory but have improved seasonal variability. The changes from 2000 to 2005 and for 2010 are obtained in the MACCity data using the representative concentration pathway (RCP) scenario version 8.5. For the production of CAMSiRA, the MACCity data set was extended to 2015 by also applying the RCP 8.5 scenario. The anthropogenic CO emissions were increased following Stein et al. (2014). Time series of the anthropogenic CO emissions for Europe, North America, East Asia (see Table 3) and the globe are shown in Fig. S2 of the Supplement.
Coordinates of regions.
Timeline of assimilated AC satellite retrievals from different instruments assimilated in CAMSiRA (see Table 2).
The anthropogenic emissions of organic matter, black carbon and aerosol
precursor SO
The biogenic emissions for the chemical species were simulated offline by
the MEGAN2.1 model (Guenther et al., 2006) for the 2000–2010 period
(MEGAN-MACC, Sindelarova et al., 2014). For the remaining years 2011–2015 a
climatology data set of the MEGAN-MACC data was put together. Natural emissions from
soils and oceans for NO
Daily biomass burning emissions for reactive gases and aerosols were produced by the Global Fire Assimilation System (GFAS) version 1.2, which is based on satellite retrievals of fire radiative power (Kaiser et al., 2012). This is an important difference with respect to the MACCRA, which used early versions of the Global Fire Emissions Database (GFED 3.1) data from 2003 until the end of 2008 and daily GFAS v1.0 data from 2009 to 2012. The GFED 3.1 is on average 20 % lower than GFAS v1.2 (Inness et al., 2013). Time series of the biomass burning CO emissions for tropical Africa, South America, Maritime South East Asia (see Table 3) and the globe are shown in Fig. S3.
C-IFS uses an incremental 4D-VAR algorithm (Courtier et al., 1994), which
minimizes a cost function for selected control variables to combine the model
and the observations in order to obtain the best possible representations of
the atmospheric fields. The mass mixing ratios of O
The assimilation of AOD retrievals uses an observation operator that translates the aerosol mass mixing ratios and humidity fields of C-IFS to the respective AOD (550 nm) values using precomputed optical properties. Total aerosol mass mixing ratio is included in the 4D-VAR cost function and the analysis increments are repartitioned into the individual aerosol components according to their fractional contribution to the total aerosol mass. This is an approximation that is assumed to be only valid over the 12 h of the assimilation window. In reality, the relative fraction of the aerosol components is not conserved during the whole assimilation procedure because of differences in the efficiency of the removal processes. Aerosol components with a longer atmospheric lifetime will retain the change imposed by the increments relatively longer and may thereby change the relative contributions.
In the ECMWF data assimilation system the background error covariance matrix is given in a wavelet formulation (Fisher, 2004, 2006). This allows both spatial and spectral variations in the horizontal and vertical background error covariances. The background errors for AC are constant in time.
The background errors for ozone are the same as the ones used for MACCRA
(Inness et al., 2013). Only the vertical correlations of the ozone background
errors were modified and restricted to
The background error statistics for the chemical species and for total aerosol are univariate in order to minimise the feedback effects of the chemical fields on the meteorological variables. Correlations between the background errors of different chemical species are also not accounted for (Inness et al., 2015).
A further potential interaction between the assimilated species could be introduced by the adjoint and tangent linear representations of the chemical mechanism and the aerosol module as part of the 4D-VAR approach. The applied tangent linear and adjoint formulation of C-IFS only accounts for transport processes and not the sources and sinks of atmospheric composition in this study. Because of this limitation and the lack of aerosol or chemistry and meteorology feedbacks in C-IFS, interaction among species and with the meteorology as part of the assimilation procedure are not represented in CAMSiRA.
Table 2 shows the AC composition data sets for CO, ozone and AOD that were assimilated in CAMSiRA. The timeline of the assimilation for the different retrievals is shown in Fig. 1. CO is assimilated from MOPITT V5 TIR only, whereas the MACCRA assimilated the V4 TIR product and Infrared Atmospheric Sounding Interferometer (IASI) TC CO retrievals after April 2008. The biases between the retrievals (George et al., 2015) of the two instruments in middle and higher latitudes could not be reconciled with the variational bias correction and led to a discontinuity in the time series of CO in MACCRA, which consequently could not be used for trend analyses (see Fig. 4 below). It was therefore decided to only use the MOPITT V5 CO data set in CAMSiRA because they cover the whole period from 2003 to 2015. The MOPITT V5 product has better long-term stability and a smaller SH bias than V4 (Deeter et al., 2013). V4 suffered from a positive temporal bias drift and a positive bias in the Southern Hemisphere (SH).
Additional ozone data sets in CAMSiRA were the Michelson Interferometer for
Passive Atmospheric Sounding (MIPAS) ozone profiles, which were assimilated
from 2005 until the end of the ENVISAT mission in April 2012. After the end
of 2012 the version of the assimilated Microwave Limb Sounder (MLS) data set
changed from V2 to V3.4. Information about the differences between the two
versions can be found in
Averaging kernels (AKs) were used for the calculation of the model's first-guess fields in the observation operators for the MOPITT data. For the ozone retrieval, AKs were not used because they were not provided or did not improve the analysis. For example, the high vertical resolution of the MLS ozone retrievals in the stratosphere made the use of AKs unnecessary.
The AC satellite retrievals were thinned to a horizontal resolution of
1
Variational bias correction (Dee, 2004; McNally et al., 2006; Auligné et al., 2007; Dee and Uppala, 2009) was applied to the MODIS AOD data, as well as to ozone column data from the Ozone Monitoring Instrument (OMI), the SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY) and the Global Ozone Monitoring Experiment 2 (GOME-2). The partial column of the Solar Backscatter Ultraviolet Radiometer-2 (SBUV/2), MLS and MIPAS were used to anchor the bias correction. Experience from the MACC reanalysis had shown that it was important to have an anchor for the bias correction to avoid drifts in the fields (Inness et al., 2013).
Global CTMs tend to underestimate the observed CO values in the Northern Hemisphere (NH) (Shindell et al., 2006), but data assimilation (Inness et al., 2013, 2015; Miyazaki et al., 2015; Gaubert et al., 2016) of satellite retrieval is able to successfully reduce the biases of the simulated CO fields. The correct representation of vertical CO profiles by the assimilation remains a challenge (Gaudel et al., 2015). An important next step will be the correct representation of the global CO trends by means of CO reanalyses such as CAMSiRA.
Figure 2 shows the seasonal mean of TC CO over the period 2003–2015 of CAMSiRA and the differences between the CR and MACCRA (2003–2012). Overall, the assimilation of TC CO into CAMSiRA led to an increase in TC CO in the NH and a decrease in TC CO in the SH and most of the tropics. CAMSiRA was about 2–5 % higher than the CR in the NH and as much as 20 % lower in the SH. The reduction was especially large in the tropical and subtropical outflow regions of the biomass burning regions in South America, central Africa and Maritime South East Asia. The largest reduction in these regions occurred in December–February (DJF). The largest negative bias of the CR with respect to CAMSiRA occurred over NH in DJF and March–May (MAM). Overall the zonal patterns of the biases throughout all seasons were rather uniform, indicating an underestimation of the hemispheric CO gradient in the CR. This could point to deficiencies in the simulation of the global chemical loss and production of CO as well as problems with the large-scale transport. Biases in the amount of the emissions seem to play a smaller role for the problem with the hemispheric gradient.
Average TC CO (10
However, more CO-emission-related differences occurred in September–November (SON) and to a smaller extent in June–August (JJA), when the CR had (i) higher values in the biomass burning regions and the respective outflow regions in Central Africa, Maritime South East Asia, and South America and (ii) lower values in the outflow regions of the emissions in North America and East Asia in the eastern and western North Pacific. This suggests that GFAS biomass burning emissions were too high, whereas the anthropogenic emissions in North America and East Asia were too low. Conversely, the CR had higher values than CAMSiRA in South Asia, which indicates that the anthropogenic emissions are too high in India.
Compared to MACCRA, CAMSiRA was up to 10 % higher in the northern high latitudes and up to 20 % higher above the tropical biomass burning regions and above parts of East Asia. The differences over the biomass burning regions can be attributed to the different biomass burning emission data sets (see Sect. 2.3). Over the oceans in the NH and the tropics, apart from biomass burning outflow regions, CAMSiRA CO is slightly lower (3 %) than MACCRA. The differences in the NH high latitudes are mainly caused by the reduction in MACCRA CO in this region introduced by the assimilation of IASI CO retrieval after 2008 (see also Fig. 4 below).
Figure 3 shows the average zonal mean cross section of the CO mass mixing ratio of CAMSiRA and the relative difference to the CR and MACCRA. The overestimation of the CR in the tropics and SH extratropics was found throughout the troposphere. It was most pronounced in relative terms at about 500 hPa. Stratospheric CO in CAMSiRA was much lower than in MACCRA. This might be an improvement since Gaudel et al. (2015) report an overestimation in the MACCRA over this region. In the upper troposphere CAMSiRA had higher CO than MACCRA, most notably in the tropics and SH where values are up to 40 % higher. CO was lower in the middle and lower troposphere in the SH and higher in the NH. These differences in the vertical distribution might be caused by (i) a more consistent modelling approach of the stratosphere–troposphere exchange with the online coupled C-IFS, (ii) the fact that C-IFS CB05 has a very different chemistry treatment compared to MOZART and (iii) updated background error statistics for CO (see Table 1).
Zonally averaged CO cross section of CAMSiRA (ppb) (2003–2015, left) and relative difference (%) compared to the CR (2003–2015, middle) and MACCRA (2003–2012, right).
Time series of monthly mean CO burden (Tg) over different regions (see Table 3) for the period 2003–2015 from CAMSiRA (red), the CR (blue) and MACCRA (green, 2003–2012).
Figure 4 shows time series of the monthly mean CO burden from CAMSiRA, MACCRA
and the CR for selected areas (see Table 3). The modelled global CO burden (CR)
was reduced by the assimilation by about 3 % at the start and by about
7 % at the end of the period. CAMSiRA showed a stepwise decrease in the
global CO burden from 2008 and 2009, which corresponds to a significant
(95 % confidence level) negative linear trend of
The higher global CO burdens of the CR with respect to CAMSiRA originated mainly
from the tropics and the SH mid-latitudes, which are strongly influenced by
biomass burning emissions in tropical Africa and South America. CO was
reduced by the assimilation in CAMSiRA, especially after the start of the
biomass burning season. The reduction in the biomass burning emissions of
2015 was an exceptional year because the global CO burden reached the highest values in the whole period for both CAMSiRA and the CR despite the overall decadal negative trend. The increase was caused by exceptionally high biomass burning emissions in Indonesia because of El Niño-related dry conditions. The El Niño-controlled interannual variability of CO over Maritime South East Asia was reproduced in a very similar way in CAMSiRA and the CR, but the assimilation reduced the burden by about 1 Tg (10 %).
In the regions of high anthropogenic emissions the temporal variability on a monthly scale was very similar between the CR and CAMSiRA. Both in North America and Europe, the CR underestimated the CO maximum of CAMSiRA in early spring by less than 5 % up to the year 2010, but the biases almost disappeared in later years. This means that the negative total CO trend in these regions was larger in CAMSiRA, which contains the MOPITT observations, than in the CR. It could indicate that the anthropogenic emissions were biased low at the beginning of the period but less so towards the end. Over East Asia the difference between the CR and CAMSiRA was generally very small, indicating a high degree of realism of the emissions in the area. A further explanation for this agreement is the fact that this area covers both the underestimation of CAMSiRA by the CR in NH mid-latitudes and the overestimation in the tropics. Both CAMSiRA and the CR had a negative but not significant trend over East Asia.
Stroden et al. (2016) also find good agreement between MOPITT-based and modelled negative trends for the 2000–2010 period of TC CO over Europe and North America but disagreement in the sign of the trend over eastern China, where their model, using MACCity emissions, simulates a positive trend but MOPITT has a negative trend. Over eastern China the CR (2003–2015) also had a small positive linear trend, whereas CAMSiRA had a negative trend; neither trend was statistically significant. The positive trend over eastern China in the CR was mainly driven by directly emitted CO at the surface. Because of the hemispheric influence, i.e. the hemispheric reduction in CO, the CO trend in the CR over eastern China became negative in the middle troposphere.
In the Arctic, which is influenced by the long-range transport from North America, Europe and Asia (Emmons et al., 2015), no MOPITT observations were assimilated (see Table 2) because of the higher biases of the MOPITT data in this region. The variability of the CR and CAMSiRA CO burden also matched well in this region, but the bias was much reduced after 2012.
The time series of the global CO burden of CAMSiRA and MACCRA agree better than CAMSiRA and the CR. The global burden of MACCRA is slightly lower than in CAMSiRA (1 %) until 2010 but starts to exceed CAMSiRA in 2011 and 2012. Hence, larger differences occur at the beginning and end of the MACCRA period.
The CO burden of MACCRA above the biomass burning regions of South America and tropical Africa was lower than CAMSiRA for the period 2003–2010. This is most likely because of the use of the GFED biomass burning emissions until 2008, which are on average 20 % lower than GFAS, which was used for CAMSiRA . In the years 2011–2012 MACCRA had higher values, which even led to a reversal in the sign of the trend over the two regions in MACCRA in comparison to CAMSiRA. MACCRA and CAMSiRA agreed well above the anthropogenic source regions. Only from 2008 onwards was MACCRA slightly lower, which led to enhanced negative trends.
Over the Arctic, CAMSiRA is higher from 2008 onwards, whereas MACCRA was higher at the start. This is consistent with the respective trends over Europe and North America. All data sets showed a step-like reduction in the CO burden in mid-2008 but it was most pronounced in MACCRA.
MOZAIC (Measurements of OZone, water vapour, carbon monoxide and nitrogen oxides by
in-service AIrbus aircraft) and IAGOS (In-service Aircraft for a Global
Observing System) are subsequent programmes of AC observations mounted
on commercial aircraft. The MOZAIC CO data have an accuracy of
We use the CO profiles obtained during take-off and landing to evaluate the CO fields averaged over airports in different regions from 2003 to 2012. The number of MOZAIC–IAGOS CO profiles fluctuated considerably over the years. They decreased from 2003 to 2014 by about 50 % and certain airports had many more observations than others. Since the aircraft used in MOZAIC were based in Frankfurt, the majority of the CO profiles were observed at this airport. Therefore the observations from Frankfurt dominate the European mean values. Observations from Tokyo and other Japanese cities were the largest contributions to the mean over East Asia. Atlanta, Toronto and Vancouver had the largest number of observations in the North American domain. Windhoek had by far the largest number of observations in tropical Africa and Caracas had the most in South America. The mean over Maritime South East Asia is mainly calculated from observations over Jakarta and Kuala Lumpur in 2005, 2006 and 2012, with an unbalanced coverage of the difference months.
Profiles of the mean relative bias of CAMSiRA, MACCRA and the CR against MOZAIC–IAGOS CO observations for different regions (see Table 3) averaged over the period 2003–2012 are shown in Fig. 5. We discuss here only the annual biases since the seasonal relative biases did not differ to a large extent from the annual relative biases.
Average relative bias (%) in CO of CAMSiRA, MACCRA and the CR compared with MOZAIC–IGAOS flight profiles averaged over different regions (see Table 3) for the period 2003–2012.
All three data sets underestimated the observed CO values throughout the
troposphere in Europe, North America and East Asia. At the surface and the
lower planetary boundary layer up to 900 hPa, i.e. where the highest CO concentrations are
observed, CAMSiRA and the CR had relative biases of about
Whereas the CR had the largest underestimation in the NH, it was generally higher than CAMSiRA and MACCRA in the tropics. This led to better agreement with the MOZAIC observations in South America and tropical Africa but also to an overestimation of 20–30 % in Maritime South East Asia. The limited number of observations in that region makes this result less robust. MACCRA and CAMSiRA showed little difference over South America and tropical Africa. The 10 % negative bias of MACCRA and CAMSiRA in tropical Africa is consistent with the 10 % underestimation of MOPITT v5 against MOZAIC–IAGOS over Windhoek reported by de Laat et al. (2014, their Fig. 3). Over Maritime South East Asia below 700 hPa, CAMSiRA and MACCRA overestimated CO, whereas MACCRA underestimated the observations. This could be the consequence of the different fire emissions and the different chemistry schemes, but the limited number of available profiles makes this result less representative.
NOAA Global Monitoring Division (GMD) network of flask CO surface observations (Novelli and Masarie, 2010) has good global coverage, which also includes the high latitudes of the SH and NH, to observe the background concentrations. The tropical stations represent the maritime background because they are mainly located on islands in the tropical oceans. The station density is higher in North America and Europe. The uncertainty of the NOAA-GMD CO observations is estimated to be 1–3 ppm (Novelli et al., 2003).
We calculated the mean and, for reasons of simplicity, only the linear trend at each station for the period 2003–2014 or 2003–2012 (MACCRA). The overall bias averaged over all stations of CAMSiRA and the CR was 3.0 ppb for the whole period but CAMSiRA had a slightly lower RMSE (13 ppb) than the CR (15 ppb). For the 2003–2012 period MACCRA had a bias of 6 ppb, whereas CAMSiRA and the CR had a bias of 3.1 and 3.9 ppb respectively.
Figure 6 shows the zonal means of the observed averages and the corresponding
model values at each station location as well as the median of the estimated
linear trend from the observations and the model results. The graphs were
constructed by calculating the mean concentrations and median trends of all
stations in 15
Zonal average of mean surface CO in ppb observed at NOAA-GMD
stations (2003–2014), values from CAMSiRA, the CR, and MACCRA (2003–2012)
(left) and zonal median of linear trend in ppb yr
In the SH high and mid-latitudes, the typical observed annual mean surface concentration was 50 ppbv. The background levels started to rise in the southern extratropics and reached a maximum of 145 ppbv in the NH mid-latitudes. The values then decreased to about 130 ppb in the Arctic. The general structure of the zonal variation was well represented by all data sets. The CR overestimated the SH middle and high values by 15 ppb, whereas CAMSiRA and MACCRA had a bias of 7 ppb. In the tropics CAMSiRA had slightly lower (3 ppb) values than the observations, whereas MACCRA and the CR overestimated by about 5 ppb. CAMSiRA had the highest values of all three data sets in the NH mid-latitudes but still underestimated the mean of the observations by 7 ppb. However, the observed means at the station locations in this latitude band varied in a range of about 100 ppb. The CR had a slightly larger underestimation than CAMSiRA. MACCRA underestimated the observations by more than 20 ppb in the middle and high latitudes. The reduction in the NH high latitudes in the CR and CAMSiRA was similar to the observations.
The observations in the SH showed essentially no linear trend in the
2003–2014 period. Starting in the tropics, a negative linear trend gradually
occurred, which reached values of about
In contrast to the assimilation of individual chemical gases, the assimilation of AOD observations is “underdetermined” because different combinations of the aerosol components can lead to the same extinction, i.e. AOD value. A further complicating factor is that each aerosol component has different optical properties, which depend on RH for the hydrophilic components such as sea salt, sulfate and organic matter. The correction of the speciation of the assimilated aerosol mass mixing ratio fields is therefore a big challenge despite good success in reproducing independent AOD observations with the aerosol analysis (Eskes et al., 2015).
In this section the global averages of burdens and AOD are presented. Spatial patterns of AOD will be discussed in Sect. 4.2. Global area-weighted averages of AOD at 550nm and the total global burden in teragrams for the different aerosol components are shown in Fig. 7. The figure also shows the median of the global AOD average and burdens simulated by the models of the AeroCom intercomparison study (Kinne et al., 2006; Textor et al., 2006). The CR had the highest total global average aerosol burden of 46 Tg compared to MACCRA and CAMSiRA, which both had 33 Tg. This number was very similar to the AeroCom median of 29 Tg.
Global average of total AOD (550 nm) and species AOD (left),
global total and species, and burden in teragrams (right) of sea salt (SS), desert
dust (DD), organic matter (OM), black carbon (BC), and sulfate aerosol
(SO
The global sea salt burden was about twice as high in the CR (15.1 Tg) than in CAMSiRA (8.3 Tg), and it was 16.1 Tg for MACCRA. In comparison, the median of the sea salt burden from the AeroCom models is 6.3 Tg. Another study of different emission schemes by Spada et al. (2013) found sea salt burdens in the range from 5.0 to 7.2 Tg. In light of these studies as well as the applied correction by the assimilation in CAMSiRA, the simulated sea salt burden of the CR appears to be too high. The simulated sea salt emissions of C-IFS were within the reported range in the literature (see Supplement). This suggests that the loss processes of sea salt were underestimated in C-IFS in comparison to other models. Conversely, the high sea salt burden of MACCRA was likely caused by an exaggeration of the sea salt emission by an earlier version of the emissions module.
The desert dust burden in the CR was 27 Tg, which was higher than the AeroCom median of 20 Tg. It was largely reduced to 18 Tg by the assimilation in CAMSiRA. MACCRA had an even lower desert dust burden of 12 Tg because of the underestimation of the desert dust emissions scheme used in MACCRA. As in the case of the sea salt, the underestimation of the desert dust loss by deposition and sedimentation may play an important role in the overestimation of dust burden in the CR.
The strongest relative change in the global burden by the assimilation occurred for sulfate, which was 1.2 Tg in the CR but was 4.7 Tg in CAMSiRA and 3.3 Tg in MACCRA. The respective AeroCom median value is 2 Tg. Because of the larger extinction per unit mass of sulfate, this increase in sulfate had a large impact on total AOD, which will be discussed further below.
The organic matter and black carbon burden of the CR (0.2 and 2.0 Tg) was increased by the assimilation to 0.36 and 2.4 Tg respectively. The values agreed reasonably well with the AeroCom median of 0.21 and 1.76 Tg.
In contrast to the global burden, the CR had the lowest global AOD average of 0.13. CAMSiRA and MACCRA had values of 0.16 and 0.18. The values for the CR were close to the median of the AeroCom models (0.12), but the two reanalyses had a higher value than the highest global average AOD value of the AeroCom models of 0.15.
The largest fraction of the CAMSiRA AOD came from sulfate, which was strongly increased by the assimilation. The contribution of sulfate AOD to total AOD was 13 % in the CR and 43 % in CAMSiRA. Sulphate was also the largest AOD contribution in MACCRA. The global average of sulfate AOD of the CR (0.018) was about half of the AeroCom median (0.034), which could suggest an underestimation of the global sulfate burden and AOD in the CR. Conversely, global sulfate AOD of CAMSiRA was 0.06, which was higher than the highest value of the AeroCom model ensemble (0.051).
As already discussed for the respective burdens, global desert dust AOD and sea salt AOD were strongly reduced in CAMSiRA compared to the CR. In the CR, sea salt and desert dust AOD contributed about 30 % each to the total AOD, whereas in CAMSiRA the contribution was reduced to 15 and 19 % respectively. The reduction of sea salt by the assimilation was reasonable since the sea salt burden was above the reported range by Textor (2006) and Spada et al. (2012). However, the reduction in sea salt was compensated for by the increase in sulfate, which became the most important contribution to total AOD over many parts of the ocean.
The global sea salt burden of MACCRA was higher than in CAMSiRA but similar to the CR. However, a different distribution of the mass within the size classes meant that the resulting sea salt AOD of MACCRA was 20 % higher than the CR. MACCRA had the lowest desert dust burden but differences in the size distribution towards smaller particles meant that the resulting AOD was slightly higher than the CR and 20 % higher than CAMSiRA. Black carbon and organic matter AOD and burden were similar among CAMSiRA, the CR and MACCRA.
Figure 8 shows the annual mean of total AOD and AOD for desert dust, sea salt, sulfate, black carbon and organic matter for the period 2003–2015 from CAMSiRA and the differences compared to the CR and MACCRA (2003–2012). The global maxima of the total AOD (> 0.5) in CAMSiRA were found over areas of desert dust emissions such as the Sahara, the Arabian Peninsula and the deserts of central Asia. High emissions of black carbon and organic matter from biomass burning sources in tropical Africa and anthropogenic sources in eastern China and northern India also contributed to AOD maxima on the global scale.
Total average AOD (row 1, scale max 1.0), AOD of desert dust (row 2, 1.0), sea salt (row 3, 0.5), sulfate (row 4, 0.5), organic matter (row 5, 0.5), and black carbon (row 6, 0.11) of CAMSiRA (average 2003–2015, left) and differences compared with the CR (average 2003–2015, middle) and MACCRA (average 2003–2012, right).
Zonally averaged total aerosol mass mixing ratio
(10
The increase in the global average AOD in CAMSiRA with respect to the CR by the
assimilation (see Sect. 4.1) occurred in most parts of the globe, in
particular over the areas of industrial activity in North America, Europe and
East Asia (20–30 %) as well as in the polar regions
(> 50 %), where AOD is generally low. The differences between
the CR and CAMSiRA, although varying in magnitude, exhibit similar spatial
patterns in all seasons, with the largest differences occurring throughout the NH
in MAM. As discussed in Sect. 4.1, the increase is mostly caused by a
widespread increase in sulfate AOD. Sulphate AOD increased more strongly in relative
terms over the oceans and higher latitudes. In areas of higher
modelled sulfate AOD, such as North America, Europe, northern Asia and
the Arctic, the contribution to total AOD changed from 40 to 90 %, which
made sulfate by far the most abundant aerosol species in these areas as
well as over the Antarctic, which seems unrealistic given that the global
SO
Time series of monthly mean AOD over the whole globe (land or sea points) and for different regions (see Table 3) for the period 2003–2015 from CAMSiRA (red), the CR (blue) and MACCRA (green, 2003–2012).
The identified reduction of global desert dust in CAMSiRA with respect to the CR was mainly confined to the main desert dust region, where AOD was reduced by by 0.2. Since total AOD was dominated by desert dust, total AOD was strongly reduced in these regions, whereas total AOD of CAMSiRA was always higher than the CR in other parts of the globe. The largest relative reduction of desert dust AOD occurred in the remote outflow regions from Australia, tropical Africa and Eurasia. The reduction of desert dust occurred throughout all seasons, with the largest reduction in JJA.
The strongest reduction in sea salt occurred in CAMSiRA compared to the CR and occurred over the oceans proportional to the sea salt AOD. Because of the increase in sulfate, the sea salt reduction led only to a small reduction of total AOD over the area of the highest sea salt emissions in the North Atlantic in DJF and over the Southern Ocean in JJA and MAM. The contribution of sea salt AOD to total AOD over most of the ocean was changed from more than 80 % in the CR to 50 % in CAMSiRA in middle and high latitudes of the SH and to 30 % over the rest of the maritime area by the assimilation.
Black carbon and organic matter AOD were reduced in CAMSiRA over tropical Africa where biomass burning is the largest source on the global scale, and the CO biomass burning emissions were also too high. The black carbon and organic matter AOD values were higher in CAMSiRA away from the sources where values are generally low. The differences in black carbon and organic matter AOD between CAMSiRA and the CR showed a strong reduction directly over the areas of intense fire emission in tropical Africa and the boreal forest of the NH and an increase in the adjacent outflow regions. This could indicate that the GFAS emissions, as in the case of CO (see Sect. 3.1), were too high but the atmospheric residence times of the aerosol species were too short.
Compared to CAMSiRA, MACCRA AOD values were up to 50 % (
Time series of monthly mean bias of AERONET AOD observations averaged over the whole globe (top left), Europe (top right), North America (middle left), Africa (middle right), South East Asia (bottom left), and South America (bottom right) for CAMSiRA (red), the CR (blue), and MACCRA (green).
Conversely, sea salt AOD over all oceans was much higher in MACCRA than CAMSiRA and it even exceeded the high sea salt AOD of the CR. Despite the higher sea salt AOD, the total AOD of MACCRA over the oceans was lower than in CAMSiRA because of the overall smaller sulfate AOD in maritime regions.
In the regions of boreal fire emissions, MACCRA AOD was lower during the JJA fire season as well as in the South American fire season in SON. For the rest of the globe, CAMSiRA was about 0.05 lower than MACCRA, which meant a large relative reduction (> 50 %), in particular over the oceans.
The differences between MACCRA and CAMSiRA can mainly be explained with the changes in the underlying modelling approach and the emissions since the same MODIS AOD retrievals were assimilated in both reanalyses. Differences in the background error statistics may have contributed to the differences between MACCRA and CAMSiRA, particularly in the high latitudes.
Figure 9 shows a zonally averaged cross section of the total aerosol mixing ratio of CAMSiRA and its relative differences compared to the CR and MACCRA. The highest zonal average occurred over the Southern Ocean because of the continuous sea salt production and over the latitudes of the regions with large desert dust and anthropogenic emissions. Despite the mostly higher AOD values, CAMSiRA had lower mass mixing ratios than the CR throughout the troposphere with the largest relative differences occurring over the SH mid-latitudes and in the region of intense convection in the tropics. This is related to a change in the speciation, which was discussed in Sect. 4.1. CAMSiRA had up to 90 % higher values in the stratosphere and Antarctica. The higher aerosol mixing ratios of CAMSiRA in the upper troposphere were dominated by sulfate aerosol. In relative terms, MACCRA mixing ratios were considerably higher than CAMSiRA throughout the troposphere with the exception of the NH extratropical mid-troposphere, which was caused by the lower dust emissions in MACCRA, and the SH and tropical stratosphere, which was caused by high sulfate concentrations in CAMSiRA.
Time series of monthly mean AOD from AERONET observations (light blue dots), MODIS retrievals (brown dots) and from CAMSiRA (red), the CR (blue) and MACCRA (green) at Nauru (left) and Lake Argyle (right).
Figure 10 shows time series of average AOD from CAMSiRA, the CR and MACCRA for different regions. To better distinguish the impact of sea salt, the regional AOD is averaged over land points only. The global average AOD time series are shown separately for land and sea points.
The CR and CAMSiRA did not have any significant (95 % confidence level) trends in AOD over the whole globe or any of the considered regions. There was a good agreement between CAMSiRA and the CR in their interannual variability with respect to specific years with higher maxima over South and North America as well as over Maritime South East Asia and northern Africa. This demonstrates that despite biases the model was able to reproduce the variability related to fire emissions and wind-driven desert dust suspension. A large relative difference between the CR and CAMSiRA occurred in the Arctic. The CAMSiRA and MACCRA AOD values were almost twice as high as the CR and had a much more pronounced seasonality.
In contrast to the lack of significant trends in the CR and CAMSiRA, MACCRA had
a significant positive trend over all sea points, leading to an increase over
10 years that was as large as the seasonal variation over all sea points.
Averaged over all land points, the seasonal variation is much larger than
over sea. The agreement in AOD in the monthly means time series was generally
high, but MACCRA also showed a significant increasing trend, which was not
present in the other two data sets. Most of this trend in MACCRA was caused
by dust AOD, which increased by 3.7 % yr
Even if the CR and CAMSiRA did not show significant trends in total AOD, sulfate
AOD of CAMSiRA increased significantly by 0.55 % yr
The AOD at 550 nm was evaluated with observations of the AERONET network. The AERONET is a network of about 400 stations measuring spectral AOD aerosol with ground-based sun photometers (Holben et al., 1998). The stations are mostly located over land, with a high number of stations situated in North America and Europe. The global number of stations contributing observations for the evaluation increased from about 60 in 2003 to about 250 in 2014 before it was largely reduced to only a couple of stations at the end of 2015.
Figure 11 shows time series of the monthly biases of CAMSiRA, MACCRA and the CR
for the globe and different regions. Over North America, an area with a high
density of AERONET stations, the CR underestimated AOD by 0.05 on
average. Conversely, the two analyses overestimated AOD by about 0.02,
but CAMSiRA has marginally smaller biases than MACCRA. In South America a
similar pattern was found. However, the average underestimation of the CR and
overestimation of CAMSiRA and MACCRA were
Seasonally averaged TC ozone (DU) from CAMSiRA (left), difference between CAMSIRA and the CR (middle), and difference between CAMSiRA and MACCRA (right, 2003–2012, different scale) for the seasons DJF (row 1), MAM (row 2), JJA (row 3), and SON (row 4).
Over Europe, CAMSiRA had the smallest biases and MACCRA slightly
overestimated,
whereas the CR underestimated the observations. The bias of the CR was
MACCRA had the lowest biases over South East Asia because of small biases in northern India and Indochina. In this area, as almost everywhere, MACCRA was higher than CAMSiRA and the CR. CAMSiRA underestimated the observations in this region by about 0.05. The underestimation by the CR was bigger and showed a pronounced seasonal cycle. The largest negative biases occurred at the time of the seasonal minimum in DJF.
The performance for desert dust and sea salt was more difficult to evaluate
in a robust way with AERONET stations because only few stations are available
in these regions. The average bias over Africa showed a strong reduction by the
assimilation of
the CR peak values, which occurred because of desert dust outbreaks. A good example of the successful reduction of dust by the
assimilations was Lake Argyle (16.11
Zonally averaged ozone partial pressure (mPa) of CAMSiRA (2003–2015, left) and relative difference (%) compared with the CR (2003–2015, middle) and MACCRA (2003–2012).
The AERONET AOD observations over the oceans generally show an overestimation
of all runs, in particular for MACCRA. The bias of the MODIS retrievals with
respect to AERONET (Shi et al., 2011) may be a reason for this
overestimation. The comparison with AOD observations at Mauna Loa Station
(19.54
The experience from the assimilation of TC and stratospheric profile retrievals (Inness et al., 2013; van der A et al., 2015; Levefer et al., 2015) shows that these observations are sufficient to constrain stratospheric ozone in the reanalysis. Because almost the same ozone retrievals were assimilated in CAMSiRA as in MACCRA (see Table 2), most of the differences in the ozone analyses can be attributed to differences in the ozone simulation of the assimilating model. For CAMSiRA, the Cariolle parameterization (Cariolle and Teyssèdre, 2007) of stratospheric ozone chemistry and the chemical mechanism CB05 for the troposphere were used. The tropospheric and stratospheric chemical scheme of the MOZART CTM (Kinnison et al., 2007) was used for MACCRA.
Monthly ozone TC (DU) area averaged over different regions (see Table 3) from CAMSiRA (black), the CR (blue) and MACCRA (green) for 2003–2015.
Figure 13 shows the seasonal average TC ozone from CAMSiRA and the difference between this data set and the CR and MACCRA. The differences between CAMSiRA and the CR had a meridional pattern. The assimilation in CAMSiRA increased the total ozone columns in the tropics and subtropics by up to 25 DU (8 %) and it decreased them by 50–70 DU in the NH middle and high latitudes. The largest reduction occurred in DJF and MAM. Over Antarctica the assimilation also led to lower values in austral winter (JJA), when TC ozone was reduced by up to 30 DU.
CAMSiRA was about 3–5 DU (1 %) lower than MACCRA across the globe. Larger differences of up to 10 DU (2 %) were located mainly over tropical land areas. Their shape suggest that they were partially caused by differences in tropospheric ozone (see Sect. 6.1). On the seasonal scale, CAMSiRA was about 10 DU lower over Antarctica and the Arctic in the spring seasons MAM and SON respectively.
Figure 14 shows the average ozone partial pressure cross section of CAMSiRA and the relative differences with the CR and MACCRA. The tropospheric part of the figure will be discussed in Sect. 6.1. The overestimation of the CR in the high latitudes of the NH and SH was located predominately in the middle and upper stratosphere at around 20 hPa. The underestimation in the tropics had the largest values at around 50 hPa.
In the lower and middle stratosphere, i.e. from 70 to 20 hPa, CAMSiRA and MACCRA differed by less than 5 %. Larger differences occurred above 10 hPa where MACCRA was up to 30 % higher than CAMSiRA, which will be discussed in more detail in Sect. 5.5.
Figure 15 shows area-weighted averages of the monthly TCs for the whole globe, the tropics, southern and NH mid-latitudes, Antarctica and the Arctic.
In the tropics, CAMSiRA had a significant (95 % confidence level) trend
of
Time series of monthly mean bias in DU against WOUDC Dobson sun photometers for the globe (top left), the tropics (top right), NH mid-latitudes (middle left), SH mid-latitudes (middle right), the Arctic (bottom left) and Antarctica (bottom right) for CAMSiRA (red), the CR (blue) and MACCRA (green).
The seasonal range, i.e. the difference between annual maximum and minimum, of TC ozone in CAMSiRA increased from 10 DU in the tropics to up 150 DU in the Arctic and 100 DU in Antarctica. As already mentioned in Sect. 5.1, the CR was 20 % higher than CAMSiRA in NH mid-latitudes and Antarctica. However, the interannual variability agreed reasonably well between CAMSiRA and the CR in the southern and middle hemispheric high and mid-latitudes. For example, the reduced Arctic ozone spring in 2011 (Manney et al., 2011) and the year-to-year differences in mid-latitudes found in CAMSiRA were well reproduced by the CR.
The ozone hole in austral spring is the most important feature of seasonal variability over Antarctica. Despite its simplicity, the Cariolle scheme in the CR reproduced the ozone loss during the ozone hole periods with respect to minimum value and interannual variability of TC ozone very well without assimilating any observations. The years with the deepest ozone holes (2015, 2003 and 2006) and with the shallowest ozone holes (2011,2013 and 2004) were the same in CAMSiRA and the CR. Conversely, the CR overestimated the average TC ozone during the Antarctic winter by about 30 DU.
There was generally good agreement between CAMSiRA and MACCRA over all parts of the globe, but MACCRA was on average about 5–10 DU (2 %) higher than CAMSiRA. The strong positive trend of MACCRA in the tropics together with a significant positive trend in the SH mid-latitudes led to increasing differences in the global average at the end of the MACC period. A larger difference between MACCRA and CAMSiRA occurred in winter (JJA) over Antarctica, when MACCRA was up to 25 DU lower than CAMSIRA. The depth of the ozone hole was slightly deeper in CAMSiRA than in MACCRA.
Mean relative bias of CAMSiRA (red), MACCRA (green) and the CR (blue) compared with ozone sondes in the Arctic (top left), NH mid-latitudes (top middle), tropics (top right), SH mid-latitudes (bottom left) and Antarctica (bottom middle) for the period 2003–2012.
Ozone TCs are observed from the ground with Dobson, Brewer, point filter and
Fourier transform infrared (FTIR) spectrometers. The Dobson instruments provide the longest and best
spatial coverage and we use this data set to evaluate the TC of CAMSiRA,
MACCRA and the CR. The Dobson instruments of the World Ozone and Ultraviolet Radiation Data Centre (WOUDC) network are well
calibrated and their precision is 1 % (Basher, 1982). Factors that
influence the accuracy of the Dobson spectrometer are the temperature
dependency of the ozone absorption coefficient and the presence of SO
Figure 16 shows time series of the monthly bias against the Dobson photometer
observations for different regions. Observations of about 50–60 stations
were available until 2013, but the number of stations dropped steadily to
about 10 stations at the end of 2015. CAMSiRA overestimated the observations
in the tropics and the mid-latitudes of both hemispheres by 2 DU
on average, whereas the mean bias of MACCRA was about 5 DU larger. In Antarctica and the
Arctic the biases showed a more pronounced seasonal cycle, mostly between
The biases of MACCRA increased in the tropics and the SH mid-latitudes from 2003 to 2008, whereas CAMSiRA and CR did not show an obvious change in the biases until 2012. The variability of the bias of CAMSiRA was amplified at the start of 2013 in the NH. Since this change in the bias is not seen at individual stations until the end of 2015, we conclude that the change is caused by the reduction in the number of stations available after 2013. It is not caused by the change of the assimilated MLS data set version (from V2 to V3.4) because this already took place at the beginning of 2013 (see Table 2).
The biases of the CR were much larger than the ones of CAMSiRA, and they had a strong seasonal cycle. In the tropics the CR underestimated the TC by 10 DU in DJF and 0 DU in MAM. The NH biases were positive and varied between 20 and 50 DU and in the Arctic between 20 and 70 DU. Over Antarctica the CR overestimated the observation by 40–60 DU in JJA, but the bias was close to zero or even slightly negative during the time of the ozone hole.
Monthly mean ozone profiles (mPa) at Neumayer Station from ozone sondes of CAMSiRA (red), MACCRA (green) and the CR (blue) for August to November (2003–20012).
The global network of ozone sondes is the most comprehensive independent data
set for the evaluation of the three-dimensional ozone fields from the surface to about
10 hPa, which is the level with the highest stratospheric ozone volume
mixing ratios. The observation error of the sondes is about
Figure 17 shows the profiles of the relative biases of CAMSiRA, MACCRA and the CR over the tropics, Antarctica, the Arctic and the NH and SH mid-latitudes for the period 2003–2012. All available observations were included in the average.
In the tropics, CAMSiRA had a relative bias of mostly below 10 % in most levels of the stratosphere. MACCRA strongly underestimated the ozone sondes (up to 30 %) in the lower stratosphere, but the relative bias of MACCRA was similar or slightly smaller than the bias of CAMSiRA in most parts of the stratosphere, i.e. in the pressure range from 70 to 20 hPa. The CR underestimated the ozone sondes by up to 20 % in the stratosphere up to 30 hPa. The largest underestimation of the CR occurred in the lower and mid-stratosphere, where the maximum in ozone partial pressure is located. In the upper stratosphere above 20 hPa, where the maximum of ozone volume mixing ratio is located, the relative biases of all data sets were smaller than in the levels below. The CR had almost no bias, whereas MACCRA overestimated by up to 10 %.
Over the Arctic and NH mid-latitudes, CAMSiRA and MACCRA agreed well with the sondes in the whole stratosphere, with relative biases below 5 %. The absolute biases of CAMSiRA were slightly smaller than the biases of MACCR, in particular in the lower stratosphere and upper troposphere. The CR overestimated the ozone observations by up to 25 % in the stratosphere and upper troposphere over the Arctic and up to 20 % in the NH mid-latitudes. The relative biases of the CR tended to be slightly smaller in the mid-stratosphere (50 hPa) than in the upper and lower stratosphere.
Over SH mid-latitudes and Antarctica, the annual biases in the stratosphere were slightly smaller in CAMSIRA than in MACCRA, but for both reanalyses they were below 10 %. As over the Arctic, the absolute tropospheric biases, with the exception of the surface values, were smaller in MACCRA since CAMSiRA showed an underestimation of about 10 %. The CR had a stronger underestimation in the lower and upper stratosphere.
Since the process of the ozone-hole formation cannot easily be demonstrated with
annual means, Fig. 18 shows the monthly mean profile from August to November
over Neumayer Station (70.7
In contrast, the CR showed a strong overestimation in August in the middle and lower stratosphere. Ozone in the upper stratosphere in September was underestimated in the CR because of an exaggerated depletion, whereas ozone was overestimated in the lower stratosphere. In the following months the CR ozone remained too high in the lower stratosphere and too low in the upper troposphere, but the resulting TCs matched the observations in a reasonable way (see Fig. 16)
Cross sections (50–0.3 hPa) of the relative biases of zonally averaged ozone (%) of CAMSiRA (left), the CR (middle) and MACCRA (right) compared with the GOZCARDS product (GOZ) for the period 2005–2012.
Ozone sondes do not provide accurate measurements above 10 hPa. The ozone bias profiles shown in Fig. 17 indicate higher values of MACCRA in the upper stratosphere and mesosphere, i.e. from above 10 hPa to the model top of 0.1 hPa. Although the ozone mass in this region is relatively small, the high values of the mixing ratios have a large impact on the radiative transfer and the associated heating rates. To investigate the biases in that region we used the GOZCARDS (Global OZone Chemistry And Related trace gas Data records for the Stratosphere) product (Froidevaux et al., 2015). It consists of merged SAGE I, SAGE II, HALOE, UARS and Aura MLS, and ACE-FTS data from late 1979 to 2012. SAGE II is used as the primary reference in the merging procedure for the instruments. For most of the CAMSiRA period, i.e. from 2004 onwards, Aura MLS and ACE-FTS are the dominating instruments in the upper stratosphere. Tegtmeier et al. (2013) showed that ozone retrievals from various instruments show a considerable spread in the upper stratosphere. ACE-FTS is biased high (5–10 %) above 10 hPa and biased low (5–10 %) below 10 hPa compared with the median of various retrievals.
Seasonally averaged ozone at 850 hPa (ppb) from CAMSiRA (left) difference between CAMSIRA and the CR (middle) and CAMSiRA and MACCRA (right, 2003–2012) for the seasons DJF (row 1), MAM (row 2), JJA (row 3) and SON (row 4).
Figure 19 shows cross sections of the GOZCARDS product and relative bias of
CAMSiRA, MACCRA and the CR in the vertical range from 50 to 0.3 hPa. In the region
from 10 to 5 hPa, MACCRA had a positive bias of 10–15 % in the tropics and
mid-latitudes, which has already been reported in Inness et al. (2013). About
half of the 10 DU higher TCs in MACCRA compared to CAMSiRA were caused by
this overestimation in the levels above 10 hPa. The biases of CAMSiRA in
that region were smaller and vary between 2.5 and
The CR had very similar biases to CAMSiRA above 5 hPa in the tropics and mid-latitudes. This means that the assimilation of observations already had little influence in this region even if no increments were added during the CAMSiRA assimilation above 1 hPa. Below 10 hPa the cross section of the bias shows the already-discussed strong overestimation of the CR in the middle and higher latitudes, which was largest in relative terms at around 20–15 hPa, and the underestimation in the tropics, which was largest at around 50 hPa.
Monthly ozone volume mixing ratios at 850, 500 and 200 hPa over different regions (see Table 3) from CAMSiRA (red), the CR (blue) and MACCRA (green) for 2003–2015.
Correcting tropospheric ozone by the assimilation of TC and stratospheric ozone profiles remains a challenge because the observations are dominated by the higher stratospheric mixing ratios (Wagner et al., 2015). The modelled ozone fields as well as the specification of the vertical background error correlation therefore have a large impact on the analysed tropospheric ozone fields (Inness et al., 2015).
Time series of seasonal mean ozone bias in ppb in the pressure ranges 950–700, 700–400 and 400–300 hPa compared with ozone sondes at Ny-Ålesund, De Bilt, Huntsville, Hong Kong Observatory, Nairobi and Neumayer Station for CAMSiRA (red), the CR (blue) and MACCRA (green).
We focus the discussion of the seasonal spatial patterns of monthly mean tropospheric ozone mole fraction to the 850 hPa pressure level values, but we also discuss tropospheric ozone at 500 and 200 hPa in Sect. 6.2 and comparisons with ozone sondes for different tropospheric layers in Sect. 6.3. Figure 20 shows the seasonal means of CAMSiRA and the differences compared with the CR and MACCRA at 850 hPa. Extratropical NH ozone values of CAMSiRA were mostly in the range from 35 to 55 ppb. The maximum was MAM, when values were about 20 ppb higher than the seasonal minimum in DJF. Regional maxima of over 60 ppb were situated over East Asia and the Arabian Peninsula. JJA was the season when the highest values occurred over the areas of the regional maxima. In this season additional regional maxima occurred over tropical Africa. The SH values were generally below 35 ppb. The seasonal maximum was in austral spring (SON) and the minimum was in austral summer and late autumn (SON).
Average diurnal cycle of ozone at EMEP AirBase stations in Europe (black) for the seasons MAM (top left), JJA (top right), SON (bottom left), and DJF (bottom right) for CAMSiRA (red), the CR (blue), and MACCRA (green).
The CR was about 2–4 ppb higher than CAMSiRA in most parts of the globe. Only in the higher latitudes of the SH as well as over the biomass burning regions in Africa, South America and Maritime South East Asia, was CAMSiRA up to 4 ppb lower than the CR. The biggest large-scale reduction by the assimilation in the NH occurred in DJF and the biggest increase occurred in the SH in SON. The largest absolute increases of CAMSiRA of up to 10 ppb occurred over the southern end of the Arabian Peninsula at the time of the seasonal maximum in JJA. This was the only local maximum in CAMSiRA that was increased by the assimilation.
Tropospheric ozone was the only considered species for which the differences between CAMSiRA and MACCRA were larger than the differences between CAMSiRA and the CR. This indicates the importance of the chemistry model parameterization and the limitations of the data assimilation in this respect. In the extratropics of the NH and SH, CAMSiRA was 2–5 ppb lower than MACCRA, with an increasing difference towards the poles. The largest difference occurred in the NH summer in JJA. CAMSiRA was up to 10 ppb lower than MACCRA over the continents in the tropics. Conversely, CAMSiRA had higher values than MACCRA over the tropical oceans; this was true for the Sahara as well as at the location of the strong maximum over the Arabian Peninsula, which was not present in MACCRA. The strong land–sea contrast in the differences could be caused by (i) a different efficiency of deposition over the oceans, (ii) the discussed differences in biomass burning emissions or (iii) the differences in the chemical treatment (e.g. the isoprene degradation scheme).
The vertical distribution (see Fig. 14) of the mean ozone partial pressure in the troposphere shows that CAMSiRA was lower than the CR in the whole troposphere, except for in the tropical upper troposphere, where it was up to 10 % higher, as well as below 500 hPa in the SH troposphere. Compared to MACCRA, CAMSiRA was up to 20 % higher in the middle and upper troposphere in the tropics and subtropics but increasingly lower towards the surface.
Estimating and understanding tropospheric ozone trends has been widely studied in the literature, as reviewed in Cooper et al. (2014) and Monks et al. (2015). Factors that influence the interannual variability and trends of tropospheric ozone are changes in anthropogenic and biomass burning emissions, the stratosphere–troposphere exchange and the variability of the meteorological fields. The observed trends vary strongly because these different factors are not uniform in space and time. Trends are often confined to specific seasons or levels. Positive trends are more common than negative trends and are found over Europe and North America during spring (Cooper et al., 2014).
Figure 21 shows time series of average ozone volume mixing rations over
selected regions and pressure levels at 850, 500 and 200 hPa. It is beyond
the scope of this paper to investigate the robustness of the trends in CAMSiRA
in detail. However, it is worth noting that there were only positive trends in
the considered region at 850, 500 and 200 hPa in CAMSiRA. The trends varied
between 0 and 1.1 % yr
The time series in Fig. 21 show that the higher values in the NH of the CR with
respect to CAMSiRA occurred in the entire troposphere. In the lower and mid-troposphere, CAMSiRA was lower than the CR, especially during the seasonal minimum.
In the tropics, the CR and CAMSiRA agreed well at 850 hpa, the CR was slightly
higher at 500 hPa and about 5 ppb lower than CAMSiRA at 200 hPa. At this
level CAMSiRA had a significant trend of 0.95 % yr
A more detailed inspection of the time series shows that from the start of 2013 the CR and CAMSiRA agree to a higher degree than before in the middle and upper part of the troposphere in the NH. The agreement is most likely caused by a reduced correction by the assimilation in the NH troposphere in this period. In early 2013 the assimilated MLS ozone retrieval switched from version V2 to the NRT V3.4 product (see Table 2), which had different levels and observation errors. The discontinuation of MIPAS in spring 2012 does not seem to be the reason for this behaviour.
The year-to-year variability of tropospheric ozone from MACCRA often did not resemble that of CAMSiRA. In the NH at 850 hPa (most prominently seen in the Arctic), MACCRA had increasing values until 2008, after which they dropped to the values of CAMSiRA. This drift of MACCRA and the associated negative trends are not realistic (as confirmed in Sect. 6.3). They were caused by applying the variational bias correction scheme to MLS data in MACCRA (see Inness at al., 2013 for more details). The agreement between CAMSiRA and MACCRA increases with increasing height in the extratropics, but in the tropics MACCRA showed a much stronger trend at 200 hPa than CAMSiRA.
Average seasonal cycle of surface ozone at EMEP AirBase stations (left) and at European ozone sonde sites in the pressure range (950–700 hPa) for CAMSiRA (red), the CR (blue) and MACCRA (green).
Figure 22 shows time series of seasonal biases in pressure ranges representing the lower, middle and upper troposphere from six different ozone sonde sites. The selected stations had at least one observation for each month of the 2003–2015 period and are examples for Europe (De Bilt), North America (Huntsville), the tropics (Nairobi), the Arctic (Ny-Ålesund) and Antarctica (Neumayer Station). To represent South Asia we chose Hong Kong Observatory, which had complete cover from 2003 to 2012. These individual time series depend on the specific characteristics of the individual stations and are therefore less representative than the averages over the gridded data sets shown in Sect. 6.2.
In the lower troposphere (950–700 hPa) over De Bilt, Huntsville and Nairobi,
the CR and CAMSiRA had seasonal biases mostly in the range of
In the middle troposphere the absolute biases of CAMSiRA and the CR were of the same magnitude but of different signs. In the upper troposphere the CR overestimated the observations by about 10 ppb, whereas the bias of CAMSiRA remained below 5 ppb. The overestimation of the CR is likely caused by the influence of the stratosphere where the CR was too high (see Sect. 5.4). Over Nairobi the biases of the CR and CAMSiRA were very similar in all levels, but CAMSiRA had overall lower biases in the lower troposphere. In the pressure range 400–300 hPa in the tropics, the impact of stratospheric biases on the CR is weaker because of the higher tropopause height in this region.
The biases for all three data sets at Ny-Ålesund, Huntsville and Hong Kong Observatory showed a pronounced seasonality in the middle and upper troposphere. At Huntsville the spring maximum was especially overestimated, i.e. it occurred 2–3 months too early. At Ny-Ålesund the overestimation was caused by too-high values in summer and autumn. Over Hong Kong Observatory the pronounced observed spring maximum was not well reproduced.
As already discussed in Sect. 6.2, the characteristics of the bias of CAMSiRA changed at the start of 2013 mainly in the upper parts of the NH troposphere but also throughout the troposphere over higher latitudes. In this period the CAMSiRA biases resemble the bias of the CR much more, which often means an increase in the average values, which could cause a spurious enhancement of positive trends.
At Neumayer Station CAMSiRA increased in a step-wise manner already at the start of 2012, which changed the bias from an underestimation to a slight overestimation together with an increased seasonality. This behaviour could have been caused by the discontinuation of MIPAS in spring 2012 (see Table 2). Although the MIPAS retrievals were only stratospheric profiles, the assimilation combined with total column retrievals could have triggered a correction in the troposphere (Flemming et al., 2011).
MACCRA had a less stable bias than CAMSiRA. In the lower and mid-troposphere, biases from 2006 to 2008 were much higher than in the rest of the period, when they resembled the biases of CAMSiRA and the CR more. This confirms that the discussed interannual variability of MACCRA seems less realistic than that of the CR and CAMSiRA.
It should be noted that both MACCRA and CAMSiRA suffered from larger-than-typical negative biases in the NH in the first half of 2003, which can probably be explained by biases in the initial conditions and the short spin-up period of only 1 month.
The AirBase and European Monitoring an Evaluation Programme (EMEP) databases host operational air quality observations
from different national European networks. All EMEP stations are located in
rural areas, while AirBase stations are designed to monitor pollution at
different scales. Stations of the rural regime can capture the larger-scale
signal, in particular for O
Figure 23 shows the average diurnal cycle for each season of the observed values and CAMSiRA, the CR and MACCRA. The CR and CAMSiRA were very similar and matched the shape of the observed diurnal cycle well. However, there was a constant bias of about 5 ppb in MAM and DJF. The CR had slightly smaller biases than CAMSiRA in JJA in the afternoon. MACCRA had a larger diurnal range because the daytime values were higher than the ones of CAMSiRA. This meant smaller day-time biases in MAM and DJF and hence a smaller seasonal bias for MACCRA. However, it also led to a considerable (10 ppb) daytime overestimation in JJA and a smaller overestimation in SON, as well as a less-well fit with the shape of the observed diurnal cycle in all seasons.
The winter and spring underestimation of CAMSiRA and the CR has already been reported in Flemming et al. (2015). To investigate the possible causes of this seasonal bias, Fig. 24 shows the average seasonal cycle at the surface at the EMEP AirBase stations and in the lower troposphere (950–750 hPa) over ozone sonde stations. The differences between CAMSiRA, the CR and MACCRA were more pronounced in the lower troposphere than at the surface. This indicates again that the assimilation has little influence on the surface values. The CR matched the observations in the lower troposphere well in all seasons apart from SON, when it overestimated. MACCRA had biases similar to the CR but overestimated in JJA and especially over southern Europe, as shown in Katragkou et al. (2015). CAMSiRA underestimated throughout the year with the exception of SON. Since the patterns of the seasonal biases were different in the lower troposphere and at the surface, we conclude that the winter and springtime bias at the surface is not predominately caused by tropospheric biases. It is more likely that the simulation of surface processes such as dry deposition and titration by freshly emitted NO are the reasons for this bias at the surface.
CAMSiRA is a new reanalysis data set of aerosol, CO and ozone for the period 2003–2015. It was produced by assimilating satellite retrievals of AOD and TC CO as well as TC and stratospheric ozone profile retrievals from various sensors in C-IFS using the ECMWF 4D-VAR approach. A similar set of observations was assimilated in MACCRA, a previous reanalysis data set for the period 2003–2012. A control run with C-IFS (CR) without the assimilation of AC observations was carried out to infer the impact of the assimilated observations.
Compared to its predecessor MACCRA, CAMSiRA had smaller biases of surface and lower-tropospheric CO, as shown by the comparison with MOZAIC–IAGOS CO profiles and NOAA-GMD CO flask observations. However, MACCRA had lower CO biases in the NH middle and upper troposphere with respect to the MOZAIC–IAGOS CO profiles. The biases of TC ozone compared with the WOUDC Dobson sun photometers were reduced from 5–10 DU in MACCRA to 0–5 DU in CAMSiRA. The biases of CAMSiRA compared with AERONET AOD observations were lower in most parts of the globe, with the exception of South East Asia. A larger improvement was the elimination of the positive bias of upper stratospheric ozone in MACCRA, as shown by the comparison with the GOZCARDS ozone product. CAMSiRA also had a better agreement with the shape of the mean observed diurnal cycle of AirBase ground-level ozone observations in Europe in all seasons, but winter and spring seasonal values were still underestimated by 5 ppb. We attribute all the aforementioned differences between CAMSiRA and MACCRA, which were mainly improvements, to the change in the assimilating model, which was the coupled system IFS–MOZART for MACCRA and C-IFS with updated aerosol parameterisations for CAMSiRA.
Progress achieved by changes in the assimilated observations was a noteworthy
improvement of the temporal consistency of the tropospheric CO and ozone
fields in CAMSiRA. The assimilation of IASI CO in MACCRA from 2008 onwards
had led to a decrease in the TC CO values because of the biases against the
MOPITT data set, which was assimilated during the whole period. Consequently,
the MACCRA CO fields in the middle and high latitudes of both hemispheres
showed strong negative trends which were not in agreement with linear trends
estimated from CO flask surface observations. Conversely, the linear
trends of CAMSiRA agreed well with the observed trends, which were close to
zero in the SH and reached values of about 2 ppb yr
A discontinuity in the upper- and mid-tropospheric ozone field was noted for CAMSiRA after January of 2013 and was due to a change in version of the assimilated MLS ozone retrievals. Although this change in CAMSiRA did not mean an increase in the bias, it has to be considered when trends of tropospheric ozone fields are to be calculated from the CAMSiRA data set.
The AOD in CAMSiRA was about 0.01 lower than MACCRA in most parts of the globe, mainly because of a 50 % lower burden of sea salt in CAMSiRA. CAMSiRA had higher AOD values over the desert-dust-emitting regions in northern Africa, and the global desert dust burden was higher in CAMSiRA. CAMSiRA had 25 % higher AOD contribution from sulfate than MACCRA, which is currently under scrutiny.
The comparison with the CR showed that the assimilation led to a clear improvement for CO, AOD and TC ozone as well as stratospheric and upper-tropospheric ozone.
The assimilation of MOPITT CO increased the values in the NH mid-latitudes more in the beginning of the period, which could indicate a stronger underestimation of the anthropogenic emissions in this period as well as an overestimation of the trend in the emissions. The tropical and SH values were reduced by the assimilation, which may indicate an overestimation of the biomass burning emissions in this region. However, the rather zonally homogeneous CO differences between the CR and CAMSiRA suggest that not only biases in the fire emissions but also in the CO lifetime, chemical production and CO transport need to be investigated further.
The Cariolle scheme for stratospheric ozone, which was used in C-IFS,
suffered from a large overestimation of NH mid-and high latitude
stratospheric ozone (40–60 DU) and an underestimation in the tropics
(
The assimilation had little impact on the ozone values at the surface
and in the lower troposphere, where the biases of CAMSiRA were sometimes
slightly more negative than those of the CR. The small influence could be explained by the
fact, that dry deposition velocities and important ozone precursors such as
NO
CAMSiRA had AOD values about 0.05 higher than the CR, apart from the desert dust emission regions, where the assimilation strongly reduced the modelled values. CAMSiRA tended to slightly overestimate the AERONET AOD observations and the CR tended to underestimate them, but the overall biases of CAMSiRA were smaller.
Despite moderate differences in AOD, the CR and CAMSiRA had considerable differences in the aerosol speciation. The global annual sea salt burden by C-IFS in the CR of 15 Tg was considerably higher than the result of other modelling studies (Textor et al., 2006; Spada et al., 2012). Less efficient loss processes may have played a large role in this overestimation. The assimilation strongly reduced the sea salt burden in CAMSiRA to about half of the value in the CR. The global desert dust burden was also reduced by 25 % by the assimilation, leading to lower total AOD values over the desert dust emission regions of the Sahara, Australia and middle Asia. Despite the fact that CAMSiRA had a 30 % smaller global aerosol burden than the CR, its average global AOD was about 10 % higher than that of the CR. This was caused by a strong increase in sulfate in CAMSiRA. The optical properties and assumed size distribution of sulfate make extinction more efficient for the same amount of mass. Sulphate became the dominant contribution to AOD in the regions away from the main aerosol emissions. The strong contribution of sulfate may have partly compensated for the inadequate representation of other secondary aerosols in C-IFS. However, its magnitude and spread over the whole globe seems excessive. It might be caused by the lack of strong loss processes in the free troposphere as well as biases in the assimilated observations over the open oceans. As the CR underestimates the assimilated AOD, the aerosol mass is increased during the assimilation, initially by the same relative amount for all components. However, a longer lifetime of sulfate causes a longer lasting change compared to the other aerosol species, which made sulfate the dominating aerosol. This distortion of the speciation can not be corrected for by the assimilated MODIS AOD retrievals, which do not contain information about the speciation.
CAMSiRA is a considerable improvement over MACCRA, especially with respect to
the temporal consistency. To further improve on this important aspect, one
should make sure that consistent input emission data sets and assimilated
observations are used. Changes in the assimilated observations, such as the
version change of the MLS data after 2012, should be avoided. The use of MEGAN-simulated biogenic emissions for the whole period is advisable even if no
related jumps were detected in this study. To ensure consistency between the
aerosols and chemistry components, the same SO
Since improvements to lower-tropospheric ozone by assimilating current
satellite observations are difficult to achieve, emphasis needs to be put on
the improved simulation of chemistry and dry deposition. The assimilation of
tropospheric ozone column retrievals as well as of tropospheric NO
One prospect is to enable the correction of emissions based on observations of atmospheric composition with the C-IFS data assimilation system. This could also improve the analysis of tropospheric ozone since ozone precursor emissions would be corrected. An intermediate step in this direction is to better account for the emission uncertainty in the model background error statistics.
The high sulfate burden introduced by the assimilation can be avoided by (i) the introduction of more intensive loss processes in the free troposphere, (ii) an increase in the organic matter to better represent non-accounted SOA components and (iii) changes to the vertical structure of the background errors to avoid the accumulation of aerosol mass away from the surface. In general, any modelling improvements for a better speciation will reflect in a more realistic aerosol analysis and a better exploitation of the available observations. If possible the latest reprocessed MODIS AOD data set should be used (collection 6).
In CAMSiRA and MACCRA the aerosol and chemistry schemes were independent. A better coupling between the two and the meteorological simulation is desirable. For example, the use of aerosol to modulate photolysis rates and heterogeneous uptake of aerosol as well as simulating the impact on aerosols and ozone within the radiation transfer calculation of IFS will be important next steps.
The CAMSiRA, CR and MACCRA data are freely available. Please contact copernicus-support@ecmwf.int.
The authors declare that they have no conflict of interest.
CAMS is funded by the European Union's Copernicus Programme. The GOZCARDS data were obtained from the NASA Goddard Earth Science Data and Information Services Center. We are grateful to the World Ozone and Ultraviolet Radiation Data Centre (WOUDC) for providing ozone sonde and Dobson photometer observations. We thank the Global Atmospheric Watch programme for the provision of CO and ozone surface observations. We thank the European Environmental Agency for providing access to European ozone observations in the AirBase database. We also thank the MOZAIC (Measurements of OZone, water vapour, carbon monoxide and nitrogen oxides by in-service AIrbus aircraft) and IAGOS (In-Service Aircraft for a Global Observing System) programmes for providing CO profile observations. Edited by: C. H. Song Reviewed by: two anonymous referees