TM3DAM: Assimilation of GOME data

TM3DAM: Assimilated ozone fields based on GOME data

Henk Eskes

Royal Netherlands Meteorological Institute (KNMI),
P.O.Box 201, 3730 AE De Bilt, The Netherlands

This page gives a brief introduction to the data assimilation system TM3DAM, version 3, which is used to generate global ozone fields and forecasts, based on GOME total ozone observations.

Last update: March 2003

What is data assimilation?

Data assimilation is a general technique to combine measurements with a model that describe the time evolution of the system. Assimilation improves the model state based on the information contained in the observations. It provides the most probable state of the system, given the uncertainties in the observations and estimated errors of the model forecast.

figure 1, top figure 1, bottom

Figure 1: Top: a collection of 24 hours of GOME data, 30 Nov 1999. Bottom: the assimilated ozone field at 12 GMT, based on the GOME data. Apart from the ozone hole at the South pole (blue), very low ozone values are observed above the North Sea. Analysis based on the model shows that the low ozone values are mainly caused by an exceptional atmospheric flow and a transport of ozone-poor air from the subtropics to higher latitudes. The event lasted only a few days, and after that the ozone layer in the Northern Hemisphere returned to normal values. Ozone columns in Dobson units.

Why data assimilation?

There are several motives for using data assimilation in the study of atmospheric chemisty in general and for GOME data in particular:
To extend the use of GOME data. Data assimilation converts a sequential list of satellite observations into global maps of ozone, filling the gaps in between the observations. See figure 1. These global maps are easy to use for intercomparisons with other observations and with modelled ozone distributions.
To generate a complete data base of three dimensional daily ozone distributions.
To obtain information about the accuracy of the observations, by confronting these observations with the model forecasts, and by comparing the assimilated products with other observations.
To obtain information about the performance of the various aspects in the model and possible problems, through the continuous confrontation of the model with new measurements.
To provide ozone forecasts [Eskes, 2002]. These provide a basis for UV radiation forecasts and are useful, for instance, for the planning of atmospheric measurement campaigns.

The TM3DAM model: introduction

The model provides an approximate description of the evolution of the three-dimensional distribution of ozone in the atmosphere. This evolution is determined by transport and chemistry. Ozone behaves as a tracer and is transported by the horizontal and vertical winds in the atmosphere. Ozone is also a reactive chemical, and the concentration of ozone is determined by the concentrations of other chemicals, emissions and deposition at the surface and the UV radiation from the Sun. The paper of Eskes et al. [2003] provides a more detailed discussion of the model, assimilation approach and references to related work.

The model: transport

The model calculates the horizontal and vertical transport of ozone masses. It is driven by the meteorological fields (wind, surface pressure, temperature) from the European Centre for Medium-Range Weather Forecasts (ECMWF). These fields are updated every 6 hours.

The model divides the atmosphere in 120 longitude by 90 latitude gridboxes and 44 vertical layers from the surface to 0.1 hPa (the mesosphere). The vertical layers are identical to the ECMWF layers in the stratosphere (where the ozone layer is situated). In the trosphere the amount of layers is reduced with respect to the ECMWF model. The numerical implementation of the transport is based on the second moments advection scheme developed by M. Prather. The ozone distributions are provided on a 1.5 by 1 degree longitude-latitude grid.

The model: chemistry

Stratospheric ozone chemistry is described by two parametrisations, one for gas-phase chemistry and one for heterogeneous ozone breakdown. These simplified schemes are fast: a more complete treatment of the chemistry would require the explicit modelling of in the oder of 50 chemical species, which is too computationally expensive for our purpose.

The gas-phase production and loss of ozone in the stratosphere is described by the chemistry parametrisation developed at Météo France. The implementation is discussed in Eskes [2003].

Ozone breakdown as observed in the antarctic and arctic ozone holes are initiated by heterogeneous chemical processes that take place on liquid and solid cloud particles in the stratosphere, in the presence of polar stratospheric clouds. The heterogeneous chemistry is modelled by a simple scheme that was developed by the Centre for Atmospheric Science, Cambridge University. This scheme introduces an additional 3D tracer field which describes the degree of chlorine activation as a result of the psc formation. In the presence of activation ozone is depleted.

GOME observations

The GOME instrument measures the sunlight reflected from the Earth's atmosphere and surface in the spectral range 240-790 nm (UV-visible) [Burrows, 1999]. GOME is part of the ESA ERS-2 polar orbiting satellite. TM3DAM assimilates vertical columns of ozone which are retrieved from the GOME spectra by the fast delivery service at the KNMI [Valks, 2003], or the GOME products of the European Space Agency, provided by the German Aerospace Center [Spurr, 2002]. The pixel size of GOME is 40 by 320 km, and the swath width is 960 km, see figure 1. With this swath width GOME obtains global coverage in three days.

Assimilation scheme

The assimilation scheme in TM3DAM is based on the Kalman filter assimilation approach, but the scheme avoids the extreme time and memory consuming aspects of the Kalman filter by fixing the correlations between the model forecast errors [Eskes, 2003]. This correlation is described by a time-independent function of the distance between the two points, and the shape and correlation length are determined from the forecast performance statistics. The variance (diagonal of the forecast error covariance matrix), however, depends on space and time and the evolution of the variance is described by the Kalman filter equations. As a result the scheme produces a detailed error estimate of the analysed ozone fields (see figure 2).

Figure 2: Estimated forecast error distribution at 18 GMT, 21 June 2000. Scale in Dobson units. Three aspects of the Kalman filter are demonstrated by the figure: 1) new GOME observations lead to a reduction of the uncertainty at the swaths. 2) The forecast error increases with time, from left to right in the figure. 3) The forecast error is transported by the wind field, leading to irregular shapes of the swaths. The uncertainty is high at the South pole due to the lack of GOME observations in June.

Forecast skill

The figure below shows the difference between the GOME total ozone observation and the corresponding model forecast for the ozone column at the same position and time (o-f). The figure shows both the bias and the root-mean-square (rms) difference of (o-f) as a function of the latitude, averaged over all observations in March 2000. The rms difference is the sum of three contributions, namely the model forecast error, the observation error (noise) and the representativeness error (mismatch due to the difference between the model grid cells and the field of view of the satellite). The bias between the observations and forecast is smaller than 1 % for all latitudes. Note that the rms is largest for northern midlatitudes. This is related to the large variability of ozone in the Northern hemisphere in winter.

Figure 3: GOME observation minus model forecast statistics for March 2000. Solid curve: root-mean-square difference between the GOME ozone column and the model forecast. Dashed curve: average difference or bias between GOME and the forecast. The bias is smaller than 3 DU for all latitudes. Vertical scale in Dobson units.

A remark about biases

In general all observations have biases, or systematic offsets from the truth. Typical differences between the fast-delivery GOME total ozone product and ground-based stations are of the order of 3-5 percent, depending on latitude, season and solar zenith angle. In numerical weather prediction models it is common practice to apply a bias correction to observations before the assimilation. In TM3DAM no such bias correction is applied. Through the assimilation the model adopts the ozone level as given by the observations, and in this sense the assimilated fields are "higher-level" GOME products. This is demonstrated in figure 3 (and figure 1) which shows that the relative bias between the model and GOME is very small. Given the small bias, the assimilated fields can be compared with independent observations to study the differences with either other satellite instruments or ground-based observations. An example of this is shown in figure 4.

figure 4, left figure 4, right

Figure 4: A comparison between a GOME assimilated ozone field (left) and independent TOMS observations for the same day (right, taken from the TOMS home page of NASA). To minimise the mismatch in time, the assimilated field is shown at 12 o'clock local time. Note that the small scale structures in both plots are very similar. The large-scale differences in ozone values reflect the differences in the instruments, instrument calibration and the retrieval algorithms of total ozone. Date: 31 March 2000. Scale in Dobson units.

References:

Burrows, J.P., M. Weber, M. Buchwitz, V. Rozanov, A. Ladstätter-Weibenmayer, A. Richter, R. Debeek, R. Hoogen, K. Bramstedt, K.-U. Eichmann, M. Eisinger, D. Perner, "The Global Ozone Monitoring Experiment (GOME): Mission concept and first results", J. Atmos. Sciences, 56, 151-175, 1999.

R. Spurr, W. Thomas, D. Loyola, "GOME Level 1 to 2 Algorithms Description", Technical Note ER-TN-DLR-GO-0025, Deutsches Zentrum für Luft und Raumfahrt, Oberpfaffenhofen, Germany, July 31, 2002.

Eskes, H.J., P. F. J. van Velthoven, and H. M. Kelder, "Global ozone forecasting based on ERS-2 GOME observations", Atmos. Chem. Phys., 2, 271-278, 2002

Eskes, H. J., van Velthoven, P. F. J., Valks, P. J. M. and Kelder, H. M., "Assimilation of GOME total ozone satellite observations in a three-dimensional tracer transport model", Quarterly Journal of the Royal Meteorological Society, in press, 2003. An electronic version of the article can be found here.

GOME fast delivery service, KNMI.

Valks, P. J. M., Piters, A. J. M., Lambert, J. C., and Zehner, C, and Kelder, H., "A fast delivery system for the retrieval of near-real time ozone columns from GOME data", International Journal of Remote Sensing, 24, 423-436, 2003.