omsz
>>
     OMSZ       IDŐJÁRÁS      ÉGHAJLAT      LEVEGŐKÖRNYEZET      ISMERET-TÁR   
OMSZ főoldal |  Szolgálatunkról |  Pályázatok, projektek |  Rendezvények |  Irodalom 
Felbontás: KicsiFelbontás: NormálFelbontás: KözepesFelbontás: NagyobbFelbontás: Nagy Copyright © 
Contents

HOMOGENISATION AND GRIDDING OF A CENTRAL-EUROPEAN
TEMPERATURE DATA SET
Wolfgang Schöner1, Ingeborg Auer1, Reinhard Böhm1, Michele Brunetti2, Maurizio Maugeri3, Olivier Mestre4
1 Central Institute for Meteorology and Geodynamics, Hohe Warte 38 A-1190 Vienna,
Austria, +43 1 36026 2290, wolfgang.schoener@zamg.ac.at
2 Institute of Atmosphere and Ocean Science, Bologna, Italy
3 Dipartimento di Science della Terra, Universita degli studi di Milano, Italy
4 Meteo France SCEM/CBD/DEV, Toulouse, France

ABSTRACT

Within the EC-project ALPCLIM (Environmental and Climate Records from High Elevation Alpine Sites) about 120 Central-European temperature station data of up to 240 years length were homogenised on a monthly base. For the procedure of relative homogeneity testing we used a multiple Craddock-curve intercomparison method, the MASH-test and a Bayesian rule based test of Meteo-France. Our experience shows that high-quality homogenisation does not depend on the statistical method alone but need a good combination of the statistical method, a powerful software system and detailed meta-data information. Comparison of homogenised data with raw data shows that the raw data are biased by an increasing trend of a too warm temperature level back in time. This can be explained by problems with screening and non-standardisation at the begin of temperature measurements, an overall trend due to relocation of climate station from city centres to more rural sites (airports) in the 1940ies and 1950ies as well as for the Austrian climate network by a change in observation time from 1970 to 1971. For easier data processing homogenised station data were interpolated to a 1deg grid for the region 43S/4E to 49S/18E. Comparison of ALPCLIM temperature data with the gridded temperature data of CRU (Climate Research Unit, Univ. East Anglia) confirms results of other studies that the CRU data set is not useful for describing regional climate variability.

INTRODUCTION

High quality temperature data sets (station data or gridded data) of different spatial scale (regional to global) are sparse though they are very important in climate research. Moreover, only the minor number of investigators deal with the problem of data homogeneity. However, a lot of statements about climate change and also the outputs of GCM experiments are based on these data sets. On global scale and also on regional scale the data most widely used are the gridded temperature data of CRU (Climate Research Unit). Though the authors of this data set give some attention to the question of data homogeneity (Jones, 1994) the results of comparisons with high-quality homogenised regional data sets (eg. Moberg and Alexandersson, 1997, Auer a.o. 2000) show some remarkable differences in temporal trends. As a result of these comparison studies we conclude that the CRU temperature data set is able rather to describe hemispheric to global temperature trends than regional temperature variability.

A second important and widely used temperature data set is the station data set of GHCN (Global Historical Climatology Network, Peterson and Vose, 1997). This global data set is not homogenised but passed some quality control procedures (Peterson a.o., 1997). Consequently, for climate variability studies data from this data base have to be homogenised first. Other temperature data sets are on regional scale only. The minor number of them are carefully checked for homogeneity.

Central Europe is in the favourable position of a dense and long term climate network. Especially for air temperature instrumental data go back to the 18th century. In addition, for several countries this long-term data set is supported by detailed meta data information. Within EC-project ALPCLIM the high potential of central European instrumental temperature data should be used for calibration of proxy temperature information of high Alpine ice core data (dO18, dH2). Moreover the larger scale of a gridded temperature data set should be used for investigation of area representativity of the high alpine temperature signal. Following from the above mentioned aims of ALPCLIM and statements concerning homogeneity of available temperature data sets it was necessary within project ALPCLIM to homogenise a central European data set and interpolate it to grid points for the region 43S/4E to 49S/18E.

Figure 1:
The ALPCLIM region. Dots are low-level stations (<1500m), triangles are high-level stations (>1500m)

DATA

The data set used within this study is described by Figure 1 and Figure 2. Data sources are the respective data holders as well as for single series the GHCN data set. The longest series available started in 1750 (Basel, Genf) about 10 series go back to 1800. Most of the temperature series started between 1850 and 1900 (see Figure 2). Finally the ALPCLIM raw data set included 120 series of monthly air temperature data most of them longer than 100 years.

Figure 2:
Temporal evolution of ALPCLIM temperature data set

HOMOGENISATION

The method of homogenisation used at ZAMG (Central Institute of Meteorology, Vienna) was described in detail in several papers (eg. Auer, Böhm and Schöner 1999). The main features of this method are:
- the input of meta data information
- not to assume any series to be homogenous
- to fill gaps after homogenisation procedure
- to use a multiple-series comparison method for break point detection (Multiple Craddock curves intercomparison method and MASH-test; see Auer, Böhm and Schöner, 1999 and Szentimrey, 1999 for details)
- that homogenisation is done independently for regional groups on a monthly base

The western part of French temperature series were homogenised by a method described by Mestre (1999). From initial 120 digitised series 97 were found to be homogenisable. In Figure 3 the difference between the raw data and the homogenised data are shown for the examples of the mean of the Italian and the Austrian stations (mean of all stations of each subgroup respectively). These two subsets were selected for further comparison due to their high quality meta data information of the respective country.

Figure 3:
Homogenised values minus raw values for the two ALPCLIM subsets Italy and Austria for winter (10 to 3, thin line) and summer (4 to 9, bold line)

Several clear features can be seen in Figure 3:
- the period up to about 1860 with highest values for breaks; these high values result on the one hand from the low number of stations for computation of mean value on the other hand from the low level of standardisation of meteorological networks (note that the first congress of WMO with important support for climate network standardisation was held in Vienna in 1873); moreover this first part of temperature measurements is most probably characterised by problems of screening of thermometers
- a conspicuous change in the 1940ies to 1950ies where many stations in Italy and Austria moved from city centre sites to airport or rural sites (for Italy these relocations continued also up to the 1970ies)
- the break in 1971 for the Austrian stations due to the change in observation time (from 7-14-21 to 7-14-19).

This systematic computation of inhomogeneity adjustments clearly indicates that the contribution of inhomogeneity function can not be filtered out by computation of multiple series means. Moreover we can conclude from Figure 3 that temperature increase of non-homogenised climate series of Italy and Austria since 1850 is to weak by about 1 K.

GRIDDING

For easier data handling and for climate modelling purposes ALPCLIM station data were interpolated to a 1 to 1 degree latitude longitude grid for the area of 43S/4E to 49S/18E. For interpolation a Gaussian filter function was selected with filter weight 1 for distance 0km decreasing to filter weight 0.1 for distance 200km. Some additional limiting features should avoid unrealistic interpolation results:
- for each grid point interpolation is truncated at the first year of the longest series within 200km distance
- interpolation is interrupted at the Alpine main divide, at 1500m level of elevation and at the border between coastal and continental regions (see Figure 4).

Figure 4:
Map of ALPCLIM grid points and according region for interpolation from station data

Figure 5 shows the derived ALPCLIM gridded temperature data sets (one for low-level sites (<1500m) and one for high level sites(>1500m)). Data are relative values based on the 20th century mean (period 1901-1998). The ALPCLIM data set is available from the ZAMG free of charge.

Figure 5:
Map of ALPCLIM high level (>1500m) and low level (<1500m) grid with individual starting year of grid point series

As a first application of the ALPCLIM temperature data set we compared its temporal trend with that of the CRU temperature data set. As the CRU data set is based on 5 to 5 degree latitude-longitude grid boxes the ALPCLIM grid points have to be computed to the same area extent first. However, this is possible only for two CRU grid boxes with a coverage of same area of 80% (5-10E/45-50N and 10-15E/45-50N). All other CRU grid boxes have a lower common area coverage with derived ALPCLIM grid boxes. We selected the grid box 10-15E/45-50N for the following intercomparison study.

In Figure 6 the temporal trend of the CRU grid box 10-15E/45-50N is compared to the ALPCLIM 10-15E/45-50N grid box. As an independent measure of air temperature for the same region also a modelled temperature curve is included in the figure (mean East Alpine column temperature derived from air pressure measurements using barometric height equation, see Böhm a.o., 1998 for details).

Figure 6:
Comparison of ALPCLIM annual temperature series with CRU temperature series for grid box 10-15E/45-50N as well as with modelled East Alpine temperature curve (see Böhm a.o., 1998 for explanation)

The figure shows very clear the too weak trend of CRU temperature series in the period since about 1950. This result of ALPCLIM temperature series is well confirmed by the modelled temperature series. The ALPCLIM data set confirms other investigations mentioned earlier that the CRU temperature data set is not able to describe regional trends.

CONCLUSIONS

The paper deals with the experiences of homogenisation of a central European ALPCLIM temperature data set. Results of homogenisation underlines the importance of this data pre-processing part for climate variability investigations. It came up quite clear from this investigation, that quality of homogenisation does not depend on a powerful statistical method alone but on a good combination of meta data information with efficient software including mathematical background. Temperature data are systematically biased and computation of several stations mean can not filter out the inhomogeneity part of series. Comparison of ALPCLIM temperature data set with CRU temperature data set clearly shows that the CRU data set is not able to describe climate variability for central Europe. Findings of our study emphasis the postulation of further homogenisation work of regional to global climate data sets - not only for other regions but also for other climate elements.

ACKNOWLEDGEMENTS

Temperature data were provided by DWD (Offenbach), Meteo Swiss (Zürich), MHSC (Zagreb), SHMI (Ljubljana), HMS (Budapest), SHMI (Bratislava), HA-APB (Bolzano/Bozen), Soc. Met. Subalpina (Torino), UEA-CRU (Norwich) and NOAA-GHCN (Asheville). The study was supported by EC-project ALPCLIM (contract ENV4-CT97-06389).

REFERENCES

Auer, I, R Böhm and W. Schöner, 1999: ALOCLIM - Austrian - Central European long-term climate - creation of a multiple homogenized long-term climate data-set. In: Proceedings of the 2nd seminar for homogenisation of surface climatological data. Budapest, 9-13 Nov. 1998, WCDMP-No.41, WMO-TD No.962, 47-71

Auer, I, R Böhm and W. Schöner, 2000: Long climatic series from Austria, In: Jones PD a.o. (Ed): 2nd International Climate and History Conference, 7th -11th Sep. 1998, Norwich, in press

Böhm, R, I Auer, W Schöner and M Hagen, 1998: Long alpine barometric time series in different altitude as a measure for 19th/20th century warming. Reprints of the 8th Conference on Mountain Meteorology, Flagstaff, Arizona, AMS, Boston, 72-76.

Jones, PD. 1994. Hemispheric surface air temperature variations: a reanalysis and an update to 1993. J.Climate, 7, 1794-1802.

Mestre O. 1999. Step by step procedures for choosing a model with change-points In: Proceedings of the 2nd seminar for homogenisation of surface climatological data. Budapest, 9-13 Nov. 1998, WCDMP-No.41, WMO-TD No.962, 15-26

Moberg, A and Alexandersson H, 1997: Homogenization of Swedish temperature data. Part II: Homogenised gridded air temperature compared with a subset of global air temperature since 1861. International Journal of Climatology 17, 35-54.

Peterson, TC and RS. Vose, 1997: An overview of the Global Historical Climatology Network temperature data base. Bulletin of the American Meteorological Society, 78, 2837-2849.

Peterson, TC, RS Vose, R Schmoyer, and V Razuvaev, 1997: Quality control of monthly temperature data: The GHCN experience. International Journal of Climatology, submitted.

COST Action 0601 
Városklíma 2011 
SEECOF-2 
Climate variability and climate change 
RCM Workshop 2008 
HIRLAM / AAA Workshop 2007 
The preparation of climate atlas 
17th EGOWS Meeting 
ALADIN / HIRLAM 2005 
ALATNET 2003 
ALADIN / RC LACE 2003 
ALADIN 2002 
6th Seminar for homogenization 
5th Seminar for homogenization 
METEOROLÓGIAI VILÁGNAP 2006
Természeti katasztrófák megelőzése, hatásainak csökkentése

2005 Időjárás, éghajlat, víz és
fenntartható fejlődés
2004 Időjárás, klíma és víz az
információs társadalom korában