ALBERT

Description: A four-dimensional variational (4D-var) data assimilation system for inverse modelling of atmospheric methane emissions is presented. The system is based on the TM5 atmospheric transport model. It can be used for assimilating large volumes of measurements, in particular satellite observations and quasi-continuous in-situ observations, and at the same time it enables the optimization of a large number of model parameters, specifically grid-scale emission rates. Furthermore, the variational method allows to estimate uncertainties in posterior emissions. Here, the system is applied to optimize monthly methane emissions over a 1-year time window on the basis of surface observations from the NOAA-ESRL network. The results are rigorously compared with an analogous inversion by Bergamaschi et al. (2007), which was based on the traditional synthesis approach. The posterior emissions as well as their uncertainties obtained in both inversions show a high degree of consistency. At the same time we illustrate the advantage of 4D-Var in reducing aggregation errors by optimizing emissions at the grid scale of the transport model. The full potential of the assimilation system is exploited in Meirink et al. (2008), who use satellite observations of column-averaged methane mixing ratios to optimize emissions at high spatial resolution, taking advantage of the zooming capability of the TM5 model.

Description: Inverse modeling is widely employed to provide "top-down" emission estimates using atmospheric measurements. Here, we analyze the dependence of derived CH4 emissions on the sampling frequency and density of the observational surface network, using the TM5-4DVAR inverse modeling system and synthetic observations. This sensitivity study focuses on Europe. The synthetic observations are created by TM5 forward model simulations. The inversions of these synthetic observations are performed using virtually no knowledge on the a priori spatial and temporal distribution of emissions, i.e. the emissions are derived mainly from the atmospheric signal detected by the measurement network. Using the European network of stations for which continuous or weekly flask measurements are available for 2001, the synthetic experiments can retrieve the "true" annual total emissions for single countries such as France within 20%, and for all North West European countries together within ~5%. However, larger deviations are obtained for South and East European countries due to the scarcity of stations in the measurement network. Upgrading flask sites to stations with continuous measurements leads to an improvement for central Europe in emission estimates. For realistic emission estimates over the whole European domain, however, a major extension of the number of stations in the existing network is required. We demonstrate the potential of an extended network of a total of ~60 European stations to provide realistic emission estimates over the whole European domain.

Description: The remote sensing of the atmospheric greenhouse gases methane (CH4) and carbon dioxide (CO2) in the troposphere from instrumentation aboard satellites is a new area of research. In this manuscript, results obtained from observations of the up-welling radiation in the near-infrared by SCIAMACHY on board ENVISAT are presented. Vertical columns of CH4, CO2 and oxygen (O2) have been retrieved and the (air or) O2-normalised CH4 and CO2 column amounts, the dry air column averaged mixing ratios XCH4 and XCO2 derived. In this manuscript the first results, obtained by using the version 0.4 of the Weighting Function Modified (WFM) DOAS retrieval algorithm applied to SCIAMACHY data, are described and compared with global models. For the set of individual cloud free measurements over land the standard deviation of the difference with respect to the models is in the range ~100–200 ppbv (5–10%) for XCH4 and ~14–32 ppmv (4–9%) for XCO2. The inter-hemispheric difference of the methane mixing ratio, as determined from single day data, is in the range 30–110 ppbv and in reasonable agreement with the corresponding model data (48–71 ppbv). The weak inter-hemispheric difference of the CO2 mixing ratio can also be detected with single day data. The spatiotemporal pattern of the measured and the modelled XCO2 are in reasonable agreement. However, the amplitude of the difference between the maximum and the minimum for SCIAMACHY XCO2 is about ±20 ppmv which is about a factor of four larger than the variability of the model data which is about ±5 ppmv. More studies are needed to explain the observed differences. The XCO2 model field shows low CO2 concentrations beginning of January 2003 over a spatially extended CO2 sink region located in southern tropical/sub-tropical Africa. The SCIAMACHY data also show low CO2 mixing ratios over this area. According to the model the sink region becomes a source region about six months later and exhibits higher mixing ratios. The SCIAMACHY and the model data over this region show a similar time dependence over the period from January to October 2003. These results indicate that for the first time a regional CO2 surface source/sink region has been detected by measurements from space. The interpretation of the SCIAMACHY CO2 and CH4 measurements is difficult, e.g., because the error analysis of the currently implemented retrieval algorithm indicates that the retrieval errors are on the same order as the small greenhouse gas mixing ratio changes that are to be detected.

Description: Inverse modeling is widely employed to provide "top-down" emission estimates using atmospheric measurements. Here, we analyze the dependence of derived CH4 emissions on the sampling frequency and density of the observational surface network, using the TM5-4DVAR inverse modeling system and synthetic observations. This sensitivity study focuses on Europe. The synthetic observations are created by TM5 forward model simulations. The inversions of these synthetic observations are performed using virtually no knowledge on the a priori spatial and temporal distribution of emissions, i.e. the emissions are derived mainly from the atmospheric signal detected by the measurement network. Using the European network of stations for which continuous or weekly flask measurements are available for 2001, the synthetic experiments can retrieve the "true" annual total emissions for single countries such as France within 20%, and for all North West European countries together within ~5%. However, larger deviations are obtained for South and East European countries due to the scarcity of stations in the measurement network. Upgrading flask sites to stations with continuous measurements leads to an improvement for central Europe in emission estimates. For realistic emission estimates over the whole European domain, however, a major extension of the number of stations in the existing network is required. We demonstrate the potential of an extended network of a total of ~60 European stations to provide realistic emission estimates over the whole European domain.

Description: A four-dimensional variational (4D-Var) data assimilation system for inverse modelling of atmospheric methane emissions is presented. The system is based on the TM5 atmospheric transport model. It can be used for assimilating large volumes of measurements, in particular satellite observations and quasi-continuous in-situ observations, and at the same time it enables the optimization of a large number of model parameters, specifically grid-scale emission rates. Furthermore, the variational method allows to estimate uncertainties in posterior emissions. Here, the system is applied to optimize monthly methane emissions over a 1-year time window on the basis of surface observations from the NOAA-ESRL network. The results are rigorously compared with an analogous inversion by Bergamaschi et al. (2007), which was based on the traditional synthesis approach. The posterior emissions as well as their uncertainties obtained in both inversions show a high degree of consistency. At the same time we illustrate the advantage of 4D-Var in reducing aggregation errors by optimizing emissions at the grid scale of the transport model. The full potential of the assimilation system is exploited in Meirink et al. (2008), who use satellite observations of column-averaged methane mixing ratios to optimize emissions at high spatial resolution, taking advantage of the zooming capability of the TM5 model.

Description: A new satellite-derived climate dataset – denoted CLARA-A1 ("The CM SAF cLoud, Albedo and RAdiation dataset from AVHRR data") – is described. The dataset covers the 28 yr period from 1982 until 2009 and consists of cloud, surface albedo, and radiation budget products derived from the AVHRR (Advanced Very High Resolution Radiometer) sensor carried by polar-orbiting operational meteorological satellites. Its content, anticipated accuracies, limitations, and potential applications are described. The dataset is produced by the EUMETSAT Climate Monitoring Satellite Application Facility (CM SAF) project. The dataset has its strengths in the long duration, its foundation upon a homogenized AVHRR radiance data record, and in some unique features, e.g. the availability of 28 yr of summer surface albedo and cloudiness parameters over the polar regions. Quality characteristics are also well investigated and particularly useful results can be found over the tropics, mid to high latitudes and over nearly all oceanic areas. Being the first CM SAF dataset of its kind, an intensive evaluation of the quality of the datasets was performed and major findings with regard to merits and shortcomings of the datasets are reported. However, the CM SAF's long-term commitment to perform two additional reprocessing events within the time frame 2013–2018 will allow proper handling of limitations as well as upgrading the dataset with new features (e.g. uncertainty estimates) and extension of the temporal coverage.

Description: The near-infrared spectra measured with the SCIAMACHY instrument on board the ENVISAT satellite suffer from several instrument calibration problems. The effects of three important instrument calibration issues on the retrieved methane (CH4) and carbon monoxide (CO) total columns have been investigated: the effects of the growing ice layer on the near-infrared detectors, the effects of the orbital variation of the instrument dark signal, and the effects of the dead/bad detector pixels. Corrections for each of these instrument calibration issues have been defined. The retrieved CH4 and CO total columns including these corrections show good agreement with CO measurements from the MOPITT satellite instrument and with CH4 model calculations by the chemistry transport model TM3. Using a systematic approach, it is shown that all three instrument calibration issues have a significant effect on the retrieved CH4 and CO total columns. However, the impact on the CH4 total columns is more pronounced than for CO, because of its smaller variability. Results for three different wavelength ranges are compared and show good agreement. The growing ice layer and the orbital variation of the dark signal show a systematic, but time-dependent effect on the retrieved CH4 and CO total columns, whereas the effect of the dead/bad pixels is rather unpredictable: some dead pixels show a random effect, some more systematic, and others no effect at all. The importance of accurate corrections for each of these instrument calibration issues is illustrated using examples where inaccurate corrections lead to a wrong interpretation of the results.

Description: We apply a four-dimensional variational (4D-VAR) data assimilation system to optimize carbon monoxide (CO) emissions for 2003 and 2004 and to reduce the uncertainty of emission estimates from individual sources using the chemistry transport model TM5. The system is designed to assimilate large (satellite) datasets, but in the current study only a limited amount of surface network observations from the National Oceanic and Atmospheric Administration Earth System Research Laboratory (NOAA/ESRL) Global Monitoring Division (GMD) is used to test the 4D-VAR system. By design, the system is capable to adjust the emissions in such a way that the posterior simulation reproduces background CO mixing ratios and large-scale pollution events at background stations. Uncertainty reduction up to 60 % in yearly emissions is observed over well-constrained regions and the inferred emissions compare well with recent studies for 2004. However, with the limited amount of data from the surface network, the system becomes data sparse resulting in a large solution space. Sensitivity studies have shown that model uncertainties (e.g., vertical distribution of biomass burning emissions and the OH field) and the prior inventories used, influence the inferred emission estimates. Also, since the observations only constrain total CO emissions, the 4D-VAR system has difficulties in separating anthropogenic and biogenic sources in particular. The inferred emissions are validated with NOAA aircraft data over North America and the agreement is significantly improved from the prior to posterior simulation. Validation with the Measurements Of Pollution In The Troposphere (MOPITT) instrument version 4 (V4) shows a slight improved agreement over the well-constrained Northern Hemisphere and in the tropics (except for the African continent). However, the model simulation with posterior emissions underestimates MOPITT CO total columns on the remote Southern Hemisphere (SH) by about 10 %. This is caused by a reduction in SH CO sources mainly due to surface stations on the high southern latitudes.