SUMMARY ERA-40 is a re-analysis of meteorological observations from September 1957 to August 2002 produced by the European Centre for Medium-Range Weather Forecasts (ECMWF) in collaboration with many institutions. The observing system changed considerably over this re-analysis period, with assimilable data provided by a succession of satellite-borne instruments from the 1970s onwards, supplemented by increasing numbers of observations from aircraft, ocean-buoys and other surface platforms, but with a declining number of radiosonde ascents since the late 1980s. The observations used in ERA-40 were accumulated from many sources. The first part of this paper describes the data acquisition and the principal changes in data type and coverage over the period. It also describes the data assimilation system used for ERA-40. This benefited from many of the changes introduced into operational forecasting since the mid-1990s, when the systems used for the 15-year ECMWF re-analysis (ERA-15) and the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) re-analysis were implemented. Several of the improvements are discussed. General aspects of the production of the analyses are also summarized.A number of results indicative of the overall performance of the data assimilation system, and implicitly of the observing system, are presented and discussed. The comparison of background (short-range) forecasts and analyses with observations, the consistency of the global mass budget, the magnitude of differences between analysis and background fields and the accuracy of medium-range forecasts run from the ERA-40 analyses are illustrated. Several results demonstrate the marked improvement that was made to the observing system for the * Corresponding author: European Centre for Medium-Range Weather Forecasts, Shinfield Park, Reading RG2 9AX, UK. e-mail: adrian.simmons@ecmwf. southern hemisphere in the 1970s, particularly towards the end of the decade. In contrast, the synoptic quality of the analysis for the northern hemisphere is sufficient to provide forecasts that remain skilful well into the medium range for all years. Two particular problems are also examined: excessive precipitation over tropical oceans and a too strong Brewer-Dobson circulation, both of which are pronounced in later years. Several other aspects of the quality of the re-analyses revealed by monitoring and validation studies are summarized. Expectations that the 'second-generation' ERA-40 re-analysis would provide products that are better than those from the firstgeneration ERA-15 and NCEP/NCAR re-analyses are found to have been met in most cases.

[1] This paper presents the aerosol modeling now part of the ECMWF Integrated Forecasting System (IFS). It includes new prognostic variables for the mass of sea salt, dust, organic matter and black carbon, and sulphate aerosols, interactive with both the dynamics and the physics of the model. It details the various parameterizations used in the IFS to account for the presence of tropospheric aerosols. Details are given of the various formulations and data sets for the sources of the different aerosols and of the parameterizations describing their sinks. Comparisons of monthly mean and daily aerosol quantities like optical depths against satellite and surface observations are presented. The capability of the forecast model to simulate aerosol events is illustrated through comparisons of dust plume events. The ECMWF IFS provides a good description of the horizontal distribution and temporal variability of the main aerosol types. The forecastonly model described here generally gives the total aerosol optical depth within 0.12 of the relevant observations and can therefore provide the background trajectory information for the aerosol assimilation system described in part 2 of this paper.

Breakup of the polar stratospheric vortex in the Northern Hemisphere is an event that is known to be predictable for up to a week or so ahead. This is illustrated using data from the 45-yr ECMWF Re-Analysis (ERA-40) for the sudden warmings of January 1958 and February 1979 and operational ECMWF data for February 2003. It is then shown that a similar level of skill was achieved in operational forecasts for the split of the southern stratospheric vortex in late September 2002. The highly unusual flow conditions nevertheless exposed a computational instability of the forecast model. Analyses and forecasts from reruns using improved versions of the forecasting system are presented. Isentropic maps of potential vorticity and specific humidity provide striking pictures of the advective processes at work. Forecasts as well as analyses are shown to be in good agreement with radiosonde measurements of the temperature changes associated with vortex movement, distortion, and breakup during August and September. Forecasts from 17 September onward capture the remarkable temperature rise of about 60°C recorded at 20 hPa by the Halley radiosonde station as the vortex split. Objective forecast verification and data denial experiments are used to characterize the performance of the observing and data assimilation systems and to infer overall forecast, analysis, and observation accuracy. The observations and analyses from 1957 onward in the ERA-40 archive confirm the extreme nature of the 2002 event. Secondary vortex development by barotropic instability is also discussed; in analyses for early October 2002, the process is active in the breakup of the weaker of the two vortices formed by the late-September split.

SUMMARYA vertical finite-element (FE) discretization designed for the European Centre for Medium-Range Weather Forecasts (ECMWF) model with semi-Lagrangian advection is described. Only non-local operations are evaluated in FE representation, while products of variables are evaluated in physical space. With semi-Lagrangian advection the only non-local vertical operations to be evaluated are vertical integrals. An integral operator is derived based on the Galerkin method using B-splines as basis functions with compact support. Two versions have been implemented, one using piecewise linear basis functions (hat functions) and the other using cubic B-splines. No staggering of dependent variables is employed in physical space, making the method well suited for use with semi-Lagrangian advection.The two versions of the FE scheme are compared to finite-difference (FD) schemes in both the Lorenz and the Charney-Phillips staggering of the dependent variables for the linearized model. The FE schemes give more accurate results than the two FD schemes for the phase speeds of most of the linear gravity waves. Evidence is shown that the FE schemes suffer less from the computational mode than the FD scheme with Lorenz staggering, although temperature and geopotential are held at the same set of levels in the FE scheme too. As a result, the FE schemes reduce the level of vertical noise in forecasts with the full model. They also reduce by about 50% a persistent cold bias in the lower stratosphere present with the FD scheme in Lorenz staggering (i.e. the operational scheme at ECMWF before its replacement by the cubic version of the FE scheme described here) and improve the transport in the stratosphere.

SUMMARYStratospheric humidity analyses produced operationally by the European Centre for Medium-Range Weather Forecasts (ECMWF) are discussed for the period since late January 1996 when the practice of resetting the upperlevel specific humidity to a fixed value at each analysis time was abandoned. Near-tropopause analyses are in reasonable overall agreement with independent observations. Very low humidities occur in conjunction with deep convection and a particularly cold tropopause over the equatorial western Pacific during the northern winter. Drying occurs also in the cold core of the Antarctic polar-night vortex. The lower stratosphere is moistened in the outer tropics and subtropics in summer and autumn, predominantly in the northern hemisphere. Changes associated with the latest occurrence of El Niiio are illustrated.Analysed temperatures near the tropical tropopause are generally in good agreement with corresponding radiosonde measurements, with standard-level biases of the order of 0.5 degC or less. The past two years are the coldest by about 1 degC in a series of tropical mean 100 hPa analyses extending back to 1979. A cooling trend of about 0.6 degC per decade is seen in the global means of the 100 hPa analyses.Moisture is spread zonally and upward from the tropical tropopause as the data assimilation proceeds, but the rate of upward transfer is much faster than observed. Substantial lateral mixing can occur within the stratosphere over the course of a season. Moistening at middle and high latitudes due to mixing with more humid tropospheric air is confined, realistically, to a shallow layer at the base of the stratosphere.The rate of upward transfer of tropical stratospheric moisture is much more realistic in a multi-year simulation using a version of the model that has finer stratospheric resolution than the version used for the operational data assimilation. Temperatures at the tropical tropopause and in the Antarctic polar night are accurately simulated, apart from excessive persistence of cold south-polar temperatures in late winter and early spring. The latter is conducive to drying the model stratosphere; lack of a parametrization of moistening due to methane oxidation is an obvious deficiency in this regard.