DACH
Deutsch - Österreichisch - Schweizerische
Meteorologen - Tagung
18. - 21. September 2001, Vienna, Austria

official homepage
download


Statistical modelling of total ozone measurements of the northern hemisphere

Jörg A. Mäder and Johannes Stähelin
Institute for Atmospheric and Climate Science, ETH Zürich

September 12. 2001

Introduction and Goals

Multiple regression models are widely used in ozone research because they allow to separate long-term ozone trends from natural variability by using explanatory variables (commonly used: Quasi Biennial Oscillation (QBO) and solar cycle, see e.g. Staehelin et al. 2001). Several recent studies provided evidence, that a significant part of the long-term total ozone winter trends (particularly over Europe) can be attributed to changes in the atmospheric condition rather than by the increase in the stratospheric concentrations of ozone depleting substances (EESC: Effective Equivalent Stratospheric Chlorine) (Steinbrecht et al., 1998, 2001, Appenzeller et al., 2000). In this study we describe the variability of total ozone of 65 ground sites of the northern hemisphere (minimal length: 20 years) by using a statistical model in which we introduce a monthly factor which allows to treat the 12 calendar months by one single statistical model instead of using 12 separate model. We start from an extended data set of explanatory variables which is sequentially reduced to select those which explain most of the variability.

Method

The method is based on a multiple linear regression model using the following equation:
log O3 = log O3 mean + c1 log P1,m + c2 log P2,m + cx Px,m + ... + log EESC + e
where:
O3 measured monthly total ozone value of the respective month (see Table 1)
O3 mean mean value of the respective calendar month of the entire period
P1, P2, ... Px explanatory variables (see used data) with the respective coefficients c1, c2, ... c3
EESC effective equivalent stratospheric chlorine (stratospheric concentrations of the ozone depleting substances weighted over their reactivities)
m monthly factor
e noise term
The monthly factor m is introduced to describe (i) the seasonal variation of ozone and (ii) the influence of the seasonal variation of the relationship of the explanatory variable on total ozone. We used logarithmic variables if possible to give extreme values less weight. In the procedure we started from 24 explanatory variables of Table 1, thereafter the variables which explain less of the variability were stepwise eliminated using F-values, until all variables except one were eliminated. In the procedure the monthly factor of the worst variable is excluded first before the explanatory variable itself is eliminated.

Used Data

Local variables Zonal variables Regional climate pattern, teleconnection indices
total ozone1 O3 aerosol loading6Aero north arctic oscillation9nao pacific north america9pna
pressure at tropopause2 Ptp Hemispheric pattern east atlantic9ea east pacific9ep
Global variables atlantic oscillation7ao east atlantic jetstream9eaj west pacific9wp
light flux at 10.7cm 3 Lf107 Sea surface pressure east atlantic - west russia9eawr north pacific9np
EESC 4 Chlor jakarta sea level pressure8jak polar eurasi9pe pacific transition9pt
quasi biennale oscillation 5 enso sea surface pres.8sst scandinavia9sca subtropical zonal9sz
positive phase QBOpos asian summer9as tropical - north. hemisphere9tnh
negative phase QBOneg southern oscillation8so

Results and Discussion

We obtain for every station a set of ranking of all variables (see below for a example). The set of variables obtained for the mean of all stations (see Figure 1) shows the following characteristics. Beside the seasonal variation the local quantity of the tropopause pressure (tropopause altitude) explains the largest fraction of the variability, followed by EESC, which describes (an almost linear) ozone trend. The next variables include stratospheric aerosols, describing the effects of large volcanic eruptions which have significantly influenced total ozone (e.g. after the eruptions of Mt. Pinatubo and El Chichon), followed by the Arctic Oscillation (ao), which is a hemispheric variable and which is believed to describe the behaviour of the Arctic stratospheric vortex. The teleconnection pattern representing large scale weather conditions explain a rather small part of the variability in this method, possibly because the combination of the hemispheric variable AO and the local variable of the tropopause pressure already partially describes their effects.
rank at arosa, switzerland
1-1011-2021-3031-4041-49
Ptpnaom:pem:QBOnegpt
mnppnam:Aerom:pt
EESCm:QBOposm:pnaeawreaj
aom:naom:soim:eawrm:eaj
Aerotnhm:epLf107as
sstm:tnhjakm:Lf107sz
eam:npscam:eam:jak
QBOposm:aom:scawpm:as
epsoim:EESCm:wpm:sz
m:PtppeQBOnegm:sst

Conclusion and Outlook

In the next step we will extend the investigation by introducing lower stratospheric temperatures as further variables and we will study the pattern of the long-term trends because the statistical models allow to separate the long-term behaviour caused by changes in the atmospheric conditions from those which have to be attributed to the increase in the concentrations of the ozone depleting substances.

References

Data Sources

  1. /www.msc-smc.ec.gc.ca/woudc
  2. www.cdc.noaa.gov/Datasets/ncep.reanalysis/tropopausec
  3. http://ftp.ngdc.noaa.gov/STP/SOLAR_DATAc
  4. WMO/UNEP, 1999
  5. ftp.ncep.noaa.gov/pub/cpc/wd52dg/data/indices/singa50c
  6. www.giss.nasa.gov/data/strataerc
  7. www.atmos.colostate.edu/ao/Datac
  8. tao.atmos.washington.edu/datac
  9. www.cpc.ncep.noaa.gov/data/teledocc

Download

Poster (.ps): a0 - a4

Figures (.ps): 1 - 2 - 3
Figures (.gif): 1 - 2 - 3
Mail home Stand: Tue 15. Jun 2004