Royal Netherlands Meteorological Institute

Abstract 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 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. © Royal Meteorological Society, 2005. The contributions of N. A. Rayner and R. W. Saunders are Crown copyright.

A suite of climate change indices derived from daily temperature and precipitation data, with a primary focus on extreme events, were computed and analyzed. By setting an exact formula for each index and using specially designed software, analyses done in different countries have been combined seamlessly. This has enabled the presentation of the most up‐to‐date and comprehensive global picture of trends in extreme temperature and precipitation indices using results from a number of workshops held in data‐sparse regions and high‐quality station data supplied by numerous scientists world wide. Seasonal and annual indices for the period 1951–2003 were gridded. Trends in the gridded fields were computed and tested for statistical significance. Results showed widespread significant changes in temperature extremes associated with warming, especially for those indices derived from daily minimum temperature. Over 70% of the global land area sampled showed a significant decrease in the annual occurrence of cold nights and a significant increase in the annual occurrence of warm nights. Some regions experienced a more than doubling of these indices. This implies a positive shift in the distribution of daily minimum temperature throughout the globe. Daily maximum temperature indices showed similar changes but with smaller magnitudes. Precipitation changes showed a widespread and significant increase, but the changes are much less spatially coherent compared with temperature change. Probability distributions of indices derived from approximately 200 temperature and 600 precipitation stations, with near‐complete data for 1901–2003 and covering a very large region of the Northern Hemisphere midlatitudes (and parts of Australia for precipitation) were analyzed for the periods 1901–1950, 1951–1978 and 1979–2003. Results indicate a significant warming throughout the 20th century. Differences in temperature indices distributions are particularly pronounced between the most recent two periods and for those indices related to minimum temperature. An analysis of those indices for which seasonal time series are available shows that these changes occur for all seasons although they are generally least pronounced for September to November. Precipitation indices show a tendency toward wetter conditions throughout the 20th century.

Abstract When surrounded by a transparent emission region, black holes are expected to reveal a dark shadow caused by gravitational light bending and photon capture at the event horizon. To image and study this phenomenon, we have assembled the Event Horizon Telescope, a global very long baseline interferometry array observing at a wavelength of 1.3 mm. This allows us to reconstruct event-horizon-scale images of the supermassive black hole candidate in the center of the giant elliptical galaxy M87. We have resolved the central compact radio source as an asymmetric bright emission ring with a diameter of 42 ± 3 μ as, which is circular and encompasses a central depression in brightness with a flux ratio ≳10:1. The emission ring is recovered using different calibration and imaging schemes, with its diameter and width remaining stable over four different observations carried out in different days. Overall, the observed image is consistent with expectations for the shadow of a Kerr black hole as predicted by general relativity. The asymmetry in brightness in the ring can be explained in terms of relativistic beaming of the emission from a plasma rotating close to the speed of light around a black hole. We compare our images to an extensive library of ray-traced general-relativistic magnetohydrodynamic simulations of black holes and derive a central mass of M = (6.5 ± 0.7) × 10 9 M ⊙ . Our radio-wave observations thus provide powerful evidence for the presence of supermassive black holes in centers of galaxies and as the central engines of active galactic nuclei. They also present a new tool to explore gravity in its most extreme limit and on a mass scale that was so far not accessible.

A new high-resolution regional climate change ensemble has been established for Europe within the World Climate Research Program Coordinated Regional Downscaling Experiment (EURO-CORDEX) initiative. The first set of simulations with a horizontal resolution of 12.5 km was completed for the new emission scenarios RCP4.5 and RCP8.5 with more simulations expected to follow. The aim of this paper is to present this data set to the different communities active in regional climate modelling, impact assessment and adaptation. The EURO-CORDEX ensemble results have been compared to the SRES A1B simulation results achieved within the ENSEMBLES project. The large-scale patterns of changes in mean temperature and precipitation are similar in all three scenarios, but they differ in regional details, which can partly be related to the higher resolution in EURO-CORDEX. The results strengthen those obtained in ENSEMBLES, but need further investigations. The analysis of impact indices shows that for RCP8.5, there is a substantially larger change projected for temperature-based indices than for RCP4.5. The difference is less pronounced for precipitation-based indices. Two effects of the increased resolution can be regarded as an added value of regional climate simulations. Regional climate model simulations provide higher daily precipitation intensities, which are completely missing in the global climate model simulations, and they provide a significantly different climate change of daily precipitation intensities resulting in a smoother shift from weak to moderate and high intensities.

Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH4 continue to increase, making CH4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO2). The relative importance of CH4 compared to CO2 depends on its shorter atmospheric\nlifetime, stronger warming potential, and variations in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in reducing uncertainties in the atmospheric growth rate arise from the variety of geographically overlapping CH4 sources and from the destruction of CH4 by short-lived hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations).\nFor the 2008–2017 decade, global methane emissions are estimated by atmospheric inversions (a top-down approach) to be 576 TgCH4 yr-1 (range 550–594, corresponding to the minimum and maximum estimates of the model ensemble). Of this total, 359 TgCH4 yr-1 or 60% is attributed to anthropogenic sources, that is emissions caused by direct human activity (i.e. anthropogenic emissions; range 336–376 TgCH4 yr-1 or 50 %–65 %). The mean annual total emission for the new decade (2008–2017) is 29 TgCH4 yr-1 larger than our estimate for the previous decade (2000–2009), and 24 TgCH4 yr-1 larger than the one reported in the previous budget for 2003–2012 (Saunois et al., 2016). Since 2012, global CH4 emissions have been tracking the warmest scenarios assessed by the Intergovernmental Panel on Climate Change. Bottom-up methods suggest almost 30% larger global emissions (737 TgCH4 yr-1, range 594–881) than top-down inversion methods. Indeed, bottom-up estimates for natural sources such as natural wetlands, other inland water systems, and geological sources are higher than top-down estimates. The atmospheric constraints on the top-down budget suggest that at least some of these bottom-up emissions are overestimated. The latitudinal distribution of atmospheric observation-based emissions indicates a predominance of tropical emissions (∼65% of the global budget, <30◦N) compared to mid-latitudes (∼30 %, 30–60◦ N) and high northern latitudes (∼4 %, 60–90◦N). The most important source of uncertainty in the methane budget is attributable to natural emissions, especially those from wetlands and other inland waters.\nSome of our global source estimates are smaller than those in previously published budgets (Saunois et al., 2016; Kirschke et al., 2013). In particular wetland emissions are about 35 TgCH4 yr-1 lower due to improved partition wetlands and other inland waters. Emissions from geological sources and wild animals are also found to be smaller by 7 TgCH4 yr-1 by 8 TgCH4 yr-1, respectively. However, the overall discrepancy between bottom-up and top-down estimates has been reduced by only 5% compared to Saunois et al. (2016), due to a higher estimate of emissions from inland waters, highlighting the need for more detailed research on emissions factors. Priorities for improving the methane budget include (i) a global, high-resolution map of water-saturated soils and inundated areas emitting methane based on a robust classification of different types of emitting habitats; (ii) further development of process-based models for inland-water emissions; (iii) intensification of methane observations at local scales (e.g., FLUXNET-CH4 measurements) and urban-scale monitoring to constrain bottom-up land surface models, and at regional scales (surface networks and satellites) to constrain atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) development of a 3D variational inversion system using isotopic and/or co-emitted species such as ethane to improve source partitioning.\nThe data presented here can be downloaded from https://doi.org/10.18160/GCP-CH4-2019 (Saunois et al.,\n2020) and from the Global Carbon Project

We present a European land‐only daily high‐resolution gridded data set for precipitation and minimum, maximum, and mean surface temperature for the period 1950–2006. This data set improves on previous products in its spatial resolution and extent, time period, number of contributing stations, and attention to finding the most appropriate method for spatial interpolation of daily climate observations. The gridded data are delivered on four spatial resolutions to match the grids used in previous products as well as many of the rotated pole Regional Climate Models (RCMs) currently in use. Each data set has been designed to provide the best estimate of grid box averages rather than point values to enable direct comparison with RCMs. We employ a three‐step process of interpolation, by first interpolating the monthly precipitation totals and monthly mean temperature using three‐dimensional thin‐plate splines, then interpolating the daily anomalies using indicator and universal kriging for precipitation and kriging with an external drift for temperature, then combining the monthly and daily estimates. Interpolation uncertainty is quantified by the provision of daily standard errors for every grid square. The daily uncertainty averaged across the entire region is shown to be largely dependent on the season and number of contributing observations. We examine the effect that interpolation has on the magnitude of the extremes in the observations by calculating areal reduction factors for daily maximum temperature and precipitation events with return periods up to 10 years.

The Ozone Monitoring Instrument (OMI) flies on the National Aeronautics and Space Administration's Earth Observing System Aura satellite launched in July 2004. OMI is a ultraviolet/visible (UV/VIS) nadir solar backscatter spectrometer, which provides nearly global coverage in one day with a spatial resolution of 13 km/spl times/24 km. Trace gases measured include O/sub 3/, NO/sub 2/, SO/sub 2/, HCHO, BrO, and OClO. In addition, OMI will measure aerosol characteristics, cloud top heights, and UV irradiance at the surface. OMI's unique capabilities for measuring important trace gases with a small footprint and daily global coverage will be a major contribution to our understanding of stratospheric and tropospheric chemistry and climate change. OMI's high spatial resolution is unprecedented and will enable detection of air pollution on urban scale resolution. In this paper, the instrument and its performance will be discussed.

Abstract Indices for climate variability and extremes have been used for a long time, often by assessing days with temperature or precipitation observations above or below specific physically‐based thresholds. While these indices provided insight into local conditions, few physically based thresholds have relevance in all parts of the world. Therefore, indices of extremes evolved over time and now often focus on relative thresholds that describe features in the tails of the distributions of meteorological variables. In order to help understand how extremes are changing globally, a subset of the wide range of possible indices is now being coordinated internationally which allows the results of studies from different parts of the world to fit together seamlessly. This paper reviews these as well as other indices of extremes and documents the obstacles to robustly calculating and analyzing indices and the methods developed to overcome these obstacles. Gridding indices are necessary in order to compare observations with climate model output. However, gridding indices from daily data are not always straightforward because averaging daily information from many stations tends to dampen gridded extremes. The paper describes recent progress in attribution of the changes in gridded indices of extremes that demonstrates human influence on the probability of extremes. The paper also describes model projections of the future and wraps up with a discussion of ongoing efforts to refine indices of extremes as they are being readied to contribute to the IPCC's Fifth Assessment Report. WIREs Clim Change 2011, 2:851–870. doi: 10.1002/wcc.147 This article is categorized under: Paleoclimates and Current Trends > Modern Climate Change

A new global dataset of derived indicators has been compiled to clarify whether frequency and/or severity of climatic extremes changed during the second half of the 20th century. This period provides the best spatial coverage of homogenous daily series, which can be used for calculating the proportion of global land area exhibiting a significant change in extreme or severe weather. The authors chose 10 indicators of extreme climatic events, defined from a larger selection, that could be applied to a large variety of climates. It was assumed that data producers were more inclined to release derived data in the form of annual indicator time series than releasing their original daily observations. The indicators are based on daily maximum and minimum temperature series, as well as daily totals of precipitation, and represent changes in all seasons of the year. Only time series which had 40 yr or more of almost complete records were used. A total of about 3000 indicator time series were extracted from national climate archives and collated into the unique dataset described here. Global maps showing significant changes from one multi-decadal period to another during the interval from 1946 to 1999 were produced. Coherent spatial patterns of statistically significant changes emerge, particularly an increase in warm summer nights, a decrease in the number of frost days and a decrease in intra-annual extreme temperature range. All but one of the temperaturebased indicators show a significant change. Indicators based on daily precipitation data show more mixed patterns of change but significant increases have been seen in the extreme amount derived from wet spells and number of heavy rainfall events. We can conclude that a significant proportion of the global land area was increasingly affected by a significant change in climatic extremes during the second half of the 20th century. These clear signs of change are very robust; however, large areas are still not represented, especially Africa and South America.

Addresses both fundamental and applied aspects of ocean waves including the use of wave observations made from satellites. More specifically it describes the WAM model, its scientific basis, its actual implementation, and its many applications. This model has been developed by an international group (the Wave Modelling group), and is based on a detailed physical description of air/sea interactions. It is widely used for wave forecasting for meteorological and oceanographic purposes. The three sections of the volume describe the basic statistical theory and the relevant physical processes; the numerical model and its global and regional applications; and satellite observations, their interpretation and use in data assimilation. Written by leading experts, it is a comprehensive guide and reference for researchers and advanced students in physical oceanography, meteorology, fluid dynamics, coastal engineering and physics.

This paper discusses an important issue related to the implementation and interpretation of the analysis scheme in the ensemble Kalman filter. It is shown that the observations must be treated as random variables at the analysis steps. That is, one should add random perturbations with the correct statistics to the observations and generate an ensemble of observations that then is used in updating the ensemble of model states. Traditionally, this has not been done in previous applications of the ensemble Kalman filter and, as will be shown, this has resulted in an updated ensemble with a variance that is too low.

Abstract We describe the construction of a new version of the Europe‐wide E‐OBS temperature (daily minimum, mean, and maximum values) and precipitation data set. This version provides an improved estimation of interpolation uncertainty through the calculation of a 100‐member ensemble of realizations of each daily field. The data set covers the period back to 1950 and provides gridded fields at a spacing of 0.25 ∘ × 0.25 ∘ in regular latitude/longitude coordinates. As with the original E‐OBS data set, the ensemble version is based on the station series collated as part of the ECA&D initiative. Station density varies significantly over the domain, and over time, and a reliable estimation of interpolation uncertainty in the gridded fields is therefore important for users of the data set. The uncertainty quantified by the ensemble data set is more realistic than the uncertainty estimates in the original version, although uncertainty is still underestimated in data‐sparse regions. The new data set is compared against the earlier version of E‐OBS and against regional gridded data sets produced by a selection of National Meteorological Services. In terms of both climatological averages and extreme values, the new version of E‐OBS is broadly comparable to the earlier version. Nonetheless, users will notice differences between the two E‐OBS versions, especially for precipitation, which arises from the different gridding method used.

The possibility of using shape mixtures of randomly oriented spheroids for modeling desert dust aerosol light scattering is discussed. For reducing calculation time, look‐up tables were simulated for quadrature coefficients employed in the numerical integration of spheroid optical properties over size and shape. The calculations were done for 25 bins of the spheroid axis ratio ranging from ∼0.3 (flattened spheroids) to ∼3.0 (elongated spheroids) and for 41 narrow size bins covering the size parameter range from ∼0.012 to ∼625. The look‐up tables were arranged into a software package, which allows fast, accurate, and flexible modeling of scattering by randomly oriented spheroids with different size and shape distributions. In order to evaluate spheroid model and explore the possibility of aerosol shape identification, the software tool has been integrated into inversion algorithms for retrieving detailed aerosol properties from laboratory or remote sensing polarimetric measurements of light scattering. The application of this retrieval technique to laboratory measurements by Volten et al. (2001) has shown that spheroids can closely reproduce mineral dust light scattering matrices. The spheroid model was utilized for retrievals of aerosol properties from atmospheric radiation measured by AERONET ground‐based Sun/sky‐radiometers. It is shown that mixtures of spheroids allow rather accurate fitting of measured spectral and angular dependencies of observed intensity and polarization. Moreover, it is shown that for aerosol mixtures with a significant fraction of coarse‐mode particles (radii ≥ ∼1 μm), the nonsphericity of aerosol particles can be detected as part of AERONET retrievals. The retrieval results indicate that nonspherical particles with aspect ratios ∼1.5 and higher dominate in desert dust plumes, while in the case of background maritime aerosol spherical particles are dominant. Finally, the potential of using AERONET derived spheroid mixtures for modeling the effects of aerosol particle nonsphericity in other remote sensing techniques is discussed. For example, the variability of lidar measurements (extinction to backscattering ratio and signal depolarization ratio) is illustrated and analyzed. Also, some potentially important differences in the sensitivity of angular light scattering to parameters of nonspherical versus spherical aerosols are revealed and discussed.

Abstract We present a dataset of daily resolution climatic time series that has been compiled for the European Climate Assessment (ECA). As of December 2001, this ECA dataset comprises 199 series of minimum, maximum and/or daily mean temperature and 195 series of daily precipitation amount observed at meteorological stations in Europe and the Middle East. Almost all series cover the standard normal period 1961–90, and about 50% extends back to at least 1925. Part of the dataset (90%) is made available for climate research on CDROM and through the Internet (at http://www.knmi.nl/samenw/eca ). A comparison of the ECA dataset with existing gridded datasets, having monthly resolution, shows that correlation coefficients between ECA stations and nearest land grid boxes between 1946 and 1999 are higher than 0.8 for 93% of the temperature series and for 51% of the precipitation series. The overall trends in the ECA dataset are of comparable magnitude to those in the gridded datasets. The potential of the ECA dataset for climate studies is demonstrated in two examples. In the first example, it is shown that the winter (October–March) warming in Europe in the 1976–99 period is accompanied by a positive trend in the number of warm‐spell days at most stations, but not by a negative trend in the number of cold‐spell days. Instead, the number of cold‐spell days increases over Europe. In the second example, it is shown for winter precipitation between 1946 and 1999 that positive trends in the mean amount per wet day prevail in areas that are getting drier and wetter. Because of its daily resolution, the ECA dataset enables a variety of empirical climate studies, including detailed analyses of changes in the occurrence of extremes in relation to changes in mean temperature and total precipitation. Copyright © 2002 Royal Meteorological Society.

Abstract We present the first Event Horizon Telescope (EHT) observations of Sagittarius A* (Sgr A*), the Galactic center source associated with a supermassive black hole. These observations were conducted in 2017 using a global interferometric array of eight telescopes operating at a wavelength of λ = 1.3 mm. The EHT data resolve a compact emission region with intrahour variability. A variety of imaging and modeling analyses all support an image that is dominated by a bright, thick ring with a diameter of 51.8 ± 2.3 μ as (68% credible interval). The ring has modest azimuthal brightness asymmetry and a comparatively dim interior. Using a large suite of numerical simulations, we demonstrate that the EHT images of Sgr A* are consistent with the expected appearance of a Kerr black hole with mass ∼4 × 10 6 M ⊙ , which is inferred to exist at this location based on previous infrared observations of individual stellar orbits, as well as maser proper-motion studies. Our model comparisons disfavor scenarios where the black hole is viewed at high inclination ( i > 50°), as well as nonspinning black holes and those with retrograde accretion disks. Our results provide direct evidence for the presence of a supermassive black hole at the center of the Milky Way, and for the first time we connect the predictions from dynamical measurements of stellar orbits on scales of 10 3 –10 5 gravitational radii to event-horizon-scale images and variability. Furthermore, a comparison with the EHT results for the supermassive black hole M87* shows consistency with the predictions of general relativity spanning over three orders of magnitude in central mass.

Some time ago, the continuous ranked probability score (CRPS) was proposed as a new verification tool for (probabilistic) forecast systems. Its focus is on the entire permissible range of a certain (weather) parameter. The CRPS can be seen as a ranked probability score with an infinite number of classes, each of zero width. Alternatively, it can be interpreted as the integral of the Brier score over all possible threshold values for the parameter under consideration. For a deterministic forecast system the CRPS reduces to the mean absolute error. In this paper it is shown that for an ensemble prediction system the CRPS can be decomposed into a reliability part and a resolution/uncertainty part, in a way that is similar to the decomposition of the Brier score. The reliability part of the CRPS is closely connected to the rank histogram of the ensemble, while the resolution/uncertainty part can be related to the average spread within the ensemble and the behavior of its outliers. The usefulness of such a decomposition is illustrated for the ensemble prediction system running at the European Centre for Medium-Range Weather Forecasts. The evaluation of the CRPS and its decomposition proposed in this paper can be extended to systems issuing continuous probability forecasts, by realizing that these can be interpreted as the limit of ensemble forecasts with an infinite number of members.

Evidence that extreme rainfall intensity is increasing at the global scale has strengthened considerably in recent years. Research now indicates that the greatest increases are likely to occur in short-duration storms lasting less than a day, potentially leading to an increase in the magnitude and frequency of flash floods. This review examines the evidence for subdaily extreme rainfall intensification due to anthropogenic climate change and describes our current physical understanding of the association between subdaily extreme rainfall intensity and atmospheric temperature. We also examine the nature, quality, and quantity of information needed to allow society to adapt successfully to predicted future changes, and discuss the roles of observational and modeling studies in helping us to better understand the physical processes that can influence subdaily extreme rainfall characteristics. We conclude by describing the types of research required to produce a more thorough understanding of the relationships between local-scale thermodynamic effects, large-scale atmospheric circulation, and subdaily extreme rainfall intensity.

Global tropospheric ozone distributions, budgets, and radiative forcings from an ensemble of 26 state‐of‐the‐art atmospheric chemistry models have been intercompared and synthesized as part of a wider study into both the air quality and climate roles of ozone. Results from three 2030 emissions scenarios, broadly representing “optimistic,” “likely,” and “pessimistic” options, are compared to a base year 2000 simulation. This base case realistically represents the current global distribution of tropospheric ozone. A further set of simulations considers the influence of climate change over the same time period by forcing the central emissions scenario with a surface warming of around 0.7K. The use of a large multimodel ensemble allows us to identify key areas of uncertainty and improves the robustness of the results. Ensemble mean changes in tropospheric ozone burden between 2000 and 2030 for the 3 scenarios range from a 5% decrease, through a 6% increase, to a 15% increase. The intermodel uncertainty (±1 standard deviation) associated with these values is about ±25%. Model outliers have no significant influence on the ensemble mean results. Combining ozone and methane changes, the three scenarios produce radiative forcings of −50, 180, and 300 mW m −2 , compared to a CO 2 forcing over the same time period of 800–1100 mW m −2 . These values indicate the importance of air pollution emissions in short‐ to medium‐term climate forcing and the potential for stringent/lax control measures to improve/worsen future climate forcing. The model sensitivity of ozone to imposed climate change varies between models but modulates zonal mean mixing ratios by ±5 ppbv via a variety of feedback mechanisms, in particular those involving water vapor and stratosphere‐troposphere exchange. This level of climate change also reduces the methane lifetime by around 4%. The ensemble mean year 2000 tropospheric ozone budget indicates chemical production, chemical destruction, dry deposition and stratospheric input fluxes of 5100, 4650, 1000, and 550 Tg(O 3 ) yr −1 , respectively. These values are significantly different to the mean budget documented by the Intergovernmental Panel on Climate Change (IPCC) Third Assessment Report (TAR). The mean ozone burden (340 Tg(O 3 )) is 10% larger than the IPCC TAR estimate, while the mean ozone lifetime (22 days) is 10% shorter. Results from individual models show a correlation between ozone burden and lifetime, and each model's ozone burden and lifetime respond in similar ways across the emissions scenarios. The response to climate change is much less consistent. Models show more variability in the tropics compared to midlatitudes. Some of the most uncertain areas of the models include treatments of deep tropical convection, including lightning NO x production; isoprene emissions from vegetation and isoprene's degradation chemistry; stratosphere‐troposphere exchange; biomass burning; and water vapor concentrations.

Abstract. The Copernicus Atmosphere Monitoring Service (CAMS) reanalysis is the latest global reanalysis dataset of atmospheric composition produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), consisting of three-dimensional time-consistent atmospheric composition fields, including aerosols and chemical species. The dataset currently covers the period 2003–2016 and will be extended in the future by adding 1 year each year. A reanalysis for greenhouse gases is being produced separately. The CAMS reanalysis builds on the experience gained during the production of the earlier Monitoring Atmospheric Composition and Climate (MACC) reanalysis and CAMS interim reanalysis. Satellite retrievals of total column CO; tropospheric column NO2; aerosol optical depth (AOD); and total column, partial column and profile ozone retrievals were assimilated for the CAMS reanalysis with ECMWF's Integrated Forecasting System. The new reanalysis has an increased horizontal resolution of about 80 km and provides more chemical species at a better temporal resolution (3-hourly analysis fields, 3-hourly forecast fields and hourly surface forecast fields) than the previously produced CAMS interim reanalysis. The CAMS reanalysis has smaller biases compared with most of the independent ozone, carbon monoxide, nitrogen dioxide and aerosol optical depth observations used for validation in this paper than the previous two reanalyses and is much improved and more consistent in time, especially compared to the MACC reanalysis. The CAMS reanalysis is a dataset that can be used to compute climatologies, study trends, evaluate models, benchmark other reanalyses or serve as boundary conditions for regional models for past periods.

We use 23 atmospheric chemistry transport models to calculate current and future (2030) deposition of reactive nitrogen (NO y , NH x ) and sulfate (SO x ) to land and ocean surfaces. The models are driven by three emission scenarios: (1) current air quality legislation (CLE); (2) an optimistic case of the maximum emissions reductions currently technologically feasible (MFR); and (3) the contrasting pessimistic IPCC SRES A2 scenario. An extensive evaluation of the present‐day deposition using nearly all information on wet deposition available worldwide shows a good agreement with observations in Europe and North America, where 60–70% of the model‐calculated wet deposition rates agree to within ±50% with quality‐controlled measurements. Models systematically overestimate NH x deposition in South Asia, and underestimate NO y deposition in East Asia. We show that there are substantial differences among models for the removal mechanisms of NO y , NH x , and SO x , leading to ±1 σ variance in total deposition fluxes of about 30% in the anthropogenic emissions regions, and up to a factor of 2 outside. In all cases the mean model constructed from the ensemble calculations is among the best when comparing to measurements. Currently, 36–51% of all NO y , NH x , and SO x is deposited over the ocean, and 50–80% of the fraction of deposition on land falls on natural (nonagricultural) vegetation. Currently, 11% of the world's natural vegetation receives nitrogen deposition in excess of the “critical load” threshold of 1000 mg(N) m −2 yr −1 . The regions most affected are the United States (20% of vegetation), western Europe (30%), eastern Europe (80%), South Asia (60%), East Asia (40%), southeast Asia (30%), and Japan (50%). Future deposition fluxes are mainly driven by changes in emissions, and less importantly by changes in atmospheric chemistry and climate. The global fraction of vegetation exposed to nitrogen loads in excess of 1000 mg(N) m −2 yr −1 increases globally to 17% for CLE and 25% for A2. In MFR, the reductions in NO y are offset by further increases for NH x deposition. The regions most affected by exceedingly high nitrogen loads for CLE and A2 are Europe and Asia, but also parts of Africa.