ERA-Interim is the latest global atmospheric reanalysis produced by the EuropeanCentre for Medium-Range Weather Forecasts (ECMWF). The ERA-Interim project was conducted in part to prepare for a new atmospheric reanalysis to replace ERA-40, which will extend back to the early part of the twentieth century. This article describes the forecast model, data assimilation method, and input datasets used to produce ERA-Interim, and discusses the performance of the system. Special emphasis is placed on various difficulties encountered in the production of ERA-40, including the representation of the hydrological cycle, the quality of the stratospheric circulation, and the consistency in time of the reanalysed fields. We provide evidence for substantial improvements in each of these aspects. We also identify areas where further work is needed and describe opportunities and objectives for future reanalysis projects at ECMWF.

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] The turbulent mixing in thin ocean surface boundary layers (OSBL), which occupy the upper 100 m or so of the ocean, control the exchange of heat and trace gases between the atmosphere and ocean. Here we show that current parameterizations of this turbulent mixing lead to systematic and substantial errors in the depth of the OSBL in global climate models, which then leads to biases in sea surface temperature. One reason, we argue, is that current parameterizations are missing key surface-wave processes that force Langmuir turbulence that deepens the OSBL more rapidly than steady wind forcing. Scaling arguments are presented to identify two dimensionless parameters that measure the importance of wave forcing against wind forcing, and against buoyancy forcing. A global perspective on the occurrence of waveforced turbulence is developed using re-analysis data to compute these parameters globally. The diagnostic study developed here suggests that turbulent energy available for mixing the OSBL is under-estimated without forcing by surface waves. Wave-forcing and hence Langmuir turbulence could be important over wide areas of the ocean and in all seasons in the Southern Ocean. We conclude that surfacewave-forced Langmuir turbulence is an important process in the OSBL that requires parameterization.

CERA-20C is a coupled reanalysis of the twentieth century which aims to reconstruct the past weather and climate of the Earth system including the atmosphere, ocean, land, ocean waves, and sea ice. This reanalysis is based on the CERA coupled atmosphere-ocean assimilation system developed at ECMWF. CERA-20C provides a 10 member ensemble of reanalyses to account for errors in the observational record as well as model error. It benefited from the prior experience of the retrospective atmospheric analysis ERA-20C. The dynamical model and the data assimilation systems initially developed for NWP had been modified to take into account the evolution of the radiative forcing and the observing system. To limit the impact of changes in the observing system throughout the century, only conventional surface observations have been used in the atmosphere. CERA-20C improves the specification of the background and the observation errors, two key elements to ensure a consistent weighting of the uncertainties across geophysical variables, space, and time. The quality of CERA-20C has been evaluated against other centennial reanalyses and independent observations. Although CERA-20C inherits some limitations of ERA-20C to represent correctly the tropical cyclones in the first part of the century, it shows significant improvements in the troposphere, compared to ERA-20C and 20CRv2c (the twentieth century reanalysis produced by NOAA/CIRES). A preliminary study of the climate variability in CERA-20C has been carried out. CERA-20C improves on the representation of atmosphere-ocean heat fluxes and mean sea level pressure compared to previous uncoupled ocean and atmospheric historical reanalyses performed at ECMWF.

[1] A combined high-resolution atmospheric downscaling and wave hindcast based on the ERA-40 reanalysis covering the Norwegian Sea, the North Sea, and the Barents Sea is presented. The period covered is from September 1957 to August 2002. The dynamic atmospheric downscaling is performed as a series of short prognostic runs initialized from a blend of ERA-40 and the previous prognostic run to preserve the fine-scale surface features from the high-resolution model while maintaining the large-scale synoptic field from ERA-40. The nested WAM wave model hindcast consists of a coarse 50 km model covering the North Atlantic forced with ERA-40 winds and a nested 10-11 km resolution model forced with downscaled winds. A comparison against in situ and satellite observations of wind and sea state reveals significant improvement in mean values and upper percentiles of wind vectors and the significant wave height over ERA-40. Improvement is also found in the mean wave period. ERA-40 is biased low in wind speed and significant wave height, a bias which is not reproduced by the downscaling. The atmospheric downscaling also reproduces polar lows, which cannot be resolved by ERA-40, but the lows are too weak and short-lived as the downscaling is not capable of capturing their full life cycle.