Our understanding of the processes that control the burden and budget of tropospheric ozone has changed dramatically over the last 60 years. Models are the key tools used to understand these changes, and these underscore that there are many processes important in controlling the tropospheric ozone budget. In this critical review, we assess our evolving understanding of these processes, both physical and chemical. We review model simulations from the International Global Atmospheric Chemistry Atmospheric Chemistry and Climate Model Intercomparison Project and Chemistry Climate Modelling Initiative to assess the changes in the tropospheric ozone burden and its budget from 1850 to 2010. Analysis of these data indicates that there has been significant growth in the ozone burden from 1850 to 2000 (approximately 43 ± 9%) but smaller growth between 1960 and 2000 (approximately 16 ± 10%) and that the models simulate burdens of ozone well within recent satellite estimates. The Chemistry Climate Modelling Initiative model ozone budgets indicate that the net chemical production of ozone in the troposphere plateaued in the 1990s and has not changed since then inspite of increases in the burden. There has been a shift in net ozone production in the troposphere being greatest in the northern mid and high latitudes to the northern tropics, driven by the regional evolution of precursor emissions. An analysis of the evolution of tropospheric ozone through the 21st century, as simulated by Climate Model Intercomparison Project Phase 5 models, reveals a large source of uncertainty associated with models themselves (i.e., in the way that they simulate the chemical and physical processes that control tropospheric ozone). This structural uncertainty is greatest in the near term (two to three decades), but emissions scenarios dominate uncertainty in the longer term (2050–2100) evolution of tropospheric ozone. This intrinsic model uncertainty prevents robust predictions of near-term changes in the tropospheric ozone burden, and we review how progress can be made to reduce this limitation.

The first Tropospheric Ozone Assessment Report (TOAR) provides information on present-day distributions and long-term trends of ozone metrics relevant for climate change, human health, and vegetation. However, only few results are available in TOAR for China due to limited long-term ozone observations. Here, we present an integrated analysis of long-term measurements of surface ozone from eight sites distributed in the North China Plain (NCP) and Yangtze River Delta (YRD), the relatively underdeveloped region Northeast China, and the remote regions in Northwest and Southwest China. Trends and present-day values for seven annual and five seasonal ozone metrics were calculated following the TOAR methodologies. We compare the seasonal and diurnal cycles of ozone concentrations as well as the present-day values of ozone among sites and discuss the long-term trends in the ozone metrics. Large and significant increases of ozone are detected at the background site in the NCP, moderate increases at the global baseline site in western China, significant decreases at the northwestern edge of China, and nearly no trend at other sites. Extremely high values of ozone occurred in the NCP and YRD, particularly in warmer seasons. The present-day levels of summer ozone metrics in the NCP are much higher than the thresholds set in TOAR for the highest value groups of ozone metrics. The summer ozone metrics at the Shangdianzi background site in the NCP indicate increases at rates of more than 2%/yr during 2004–2016. In contrast, ozone at the Lin’an background site in the YRD was constant over the period 2006–2016. Our results fill some knowledge gaps in spatiotemporal changes of ozone in China and may be of useful in the assessment of ozone impacts on human health and vegetation.

From the earliest observations of ozone in the lower atmosphere in the 19th century, both measurement methods and the portion of the globe observed have evolved and changed. These methods have different uncertainties and biases, and the data records differ with respect to coverage (space and time), information content, and representativeness. In this study, various ozone measurement methods and ozone datasets are reviewed and selected for inclusion in the historical record of background ozone levels, based on relationship of the measurement technique to the modern UV absorption standard, absence of interfering pollutants, representativeness of the well-mixed boundary layer and expert judgement of their credibility. There are significant uncertainties with the 19th and early 20th-century measurements related to interference of other gases. Spectroscopic methods applied before 1960 have likely underestimated ozone by as much as 11% at the surface and by about 24% in the free troposphere, due to the use of differing ozone absorption coefficients. There is no unambiguous evidence in the measurement record back to 1896 that typical mid-latitude background surface ozone values were below about 20 nmol mol –1 , but there is robust evidence for increases in the temperate and polar regions of the northern hemisphere of 30–70%, with large uncertainty, between the period of historic observations, 1896–1975, and the modern period (1990–2014). Independent historical observations from balloons and aircraft indicate similar changes in the free troposphere. Changes in the southern hemisphere are much less. Regional representativeness of the available observations remains a potential source of large errors, which are difficult to quantify. The great majority of validation and intercomparison studies of free tropospheric ozone measurement methods use ECC ozonesondes as reference. Compared to UV-absorption measurements they show a modest (~1–5% ±5%) high bias in the troposphere, but no evidence of a change with time. Umkehr, lidar, and FTIR methods all show modest low biases relative to ECCs, and so, using ECC sondes as a transfer standard, all appear to agree to within one standard deviation with the modern UV-absorption standard. Other sonde types show an increase of 5–20% in sensitivity to tropospheric ozone from 1970–1995. Biases and standard deviations of satellite retrieval comparisons are often 2–3 times larger than those of other free tropospheric measurements. The lack of information on temporal changes of bias for satellite measurements of tropospheric ozone is an area of concern for long-term trend studies.

This Tropospheric Ozone Assessment Report (TOAR) on the current state of knowledge of ozone metrics of relevance to vegetation ( TOAR-Vegetation ) reports on present-day global distribution of ozone at over 3300 vegetated sites and the long-term trends at nearly 1200 sites. TOAR-Vegetation focusses on three metrics over vegetation-relevant time-periods across major world climatic zones: M12, the mean ozone during 08:00–19:59; AOT40, the accumulation of hourly mean ozone values over 40 ppb during daylight hours, and W126 with stronger weighting to higher hourly mean values, accumulated during 08:00–19:59. Although the density of measurement stations is highly variable across regions, in general, the highest ozone values (mean, 2010–14) are in mid-latitudes of the northern hemisphere, including southern USA, the Mediterranean basin, northern India, north, north-west and east China, the Republic of Korea and Japan. The lowest metric values reported are in Australia, New Zealand, southern parts of South America and some northern parts of Europe, Canada and the USA. Regional-scale assessments showed, for example, significantly higher AOT40 and W126 values in East Asia (EAS) than Europe (EUR) in wheat growing areas ( p < 0.05), but not in rice growing areas. In NAM, the dominant trend during 1995–2014 was a significant decrease in ozone, whilst in EUR it was no change and in EAS it was a significant increase. TOAR-Vegetation provides recommendations to facilitate a more complete global assessment of ozone impacts on vegetation in the future, including: an increase in monitoring of ozone and collation of field evidence of the damaging effects on vegetation; an investigation of the effects on peri-urban agriculture and in mountain/upland areas; inclusion of additional pollutant, meteorological and inlet height data in the TOAR dataset; where not already in existence, establishing new region-specific thresholds for vegetation damage and an innovative integration of observations and modelling including stomatal uptake of the pollutant.

The Tropospheric Ozone Assessment Report (TOAR) is an activity of the International Global Atmospheric Chemistry Project. This paper is a component of the report, focusing on the present-day distribution and trends of tropospheric ozone relevant to climate and global atmospheric chemistry model evaluation. Utilizing the TOAR surface ozone database, several figures present the global distribution and trends of daytime average ozone at 2702 non-urban monitoring sites, highlighting the regions and seasons of the world with the greatest ozone levels. Similarly, ozonesonde and commercial aircraft observations reveal ozone’s distribution throughout the depth of the free troposphere. Long-term surface observations are limited in their global spatial coverage, but data from remote locations indicate that ozone in the 21 st century is greater than during the 1970s and 1980s. While some remote sites and many sites in the heavily polluted regions of East Asia show ozone increases since 2000, many others show decreases and there is no clear global pattern for surface ozone changes since 2000. Two new satellite products provide detailed views of ozone in the lower troposphere across East Asia and Europe, revealing the full spatial extent of the spring and summer ozone enhancements across eastern China that cannot be assessed from limited surface observations. Sufficient data are now available (ozonesondes, satellite, aircraft) across the tropics from South America eastwards to the western Pacific Ocean, to indicate a likely tropospheric column ozone increase since the 1990s. The 2014–2016 mean tropospheric ozone burden (TOB) between 60°N–60°S from five satellite products is 300 Tg ± 4%. While this agreement is excellent, the products differ in their quantification of TOB trends and further work is required to reconcile the differences. Satellites can now estimate ozone’s global long-wave radiative effect, but evaluation is difficult due to limited in situ observations where the radiative effect is greatest.

Assessment of spatial and temporal variation in the impacts of ozone on human health, vegetation, and climate requires appropriate metrics. A key component of the Tropospheric Ozone Assessment Report (TOAR) is the consistent calculation of these metrics at thousands of monitoring sites globally. Investigating temporal trends in these metrics required that the same statistical methods be applied across these ozone monitoring sites. The nonparametric Mann-Kendall test (for significant trends) and the Theil-Sen estimator (for estimating the magnitude of trend) were selected to provide robust methods across all sites. This paper provides the scientific underpinnings necessary to better understand the implications of and rationale for selecting a specific TOAR metric for assessing spatial and temporal variation in ozone for a particular impact. The rationale and underlying research evidence that influence the derivation of specific metrics are given. The form of 25 metrics (4 for model-measurement comparison, 5 for characterization of ozone in the free troposphere, 11 for human health impacts, and 5 for vegetation impacts) are described. Finally, this study categorizes health and vegetation exposure metrics based on the extent to which they are determined only by the highest hourly ozone levels, or by a wider range of values. The magnitude of the metrics is influenced by both the distribution of hourly average ozone concentrations at a site location, and the extent to which a particular metric is determined by relatively low, moderate, and high hourly ozone levels. Hence, for the same ozone time series, changes in the distribution of ozone concentrations can result in different changes in the magnitude and direction of trends for different metrics. Thus, dissimilar conclusions about the effect of changes in the drivers of ozone variability (e.g., precursor emissions) on health and vegetation exposure can result from the selection of different metrics.

This study quantifies the present-day global and regional distributions (2010–2014) and trends (2000–2014) for five ozone metrics relevant for short-term and long-term human exposure. These metrics, calculated by the Tropospheric Ozone Assessment Report, are: 4 th highest daily maximum 8-hour ozone (4MDA8); number of days with MDA8 > 70 ppb (NDGT70), SOMO35 (annual Sum of Ozone Means Over 35 ppb) and two seasonally averaged metrics (3MMDA1; AVGMDA8). These metrics were explored at ozone monitoring sites worldwide, which were classified as urban or non-urban based on population and nighttime lights data. Present-day distributions of 4MDA8 and NDGT70, determined predominantly by peak values, are similar with highest levels in western North America, southern Europe and East Asia. For the other three metrics, distributions are similar with North–South gradients more prominent across Europe and Japan. Between 2000 and 2014, significant negative trends in 4MDA8 and NDGT70 occur at most US and some European sites. In contrast, significant positive trends are found at many sites in South Korea and Hong Kong, with mixed trends across Japan. The other three metrics have similar, negative trends for many non-urban North American and some European and Japanese sites, and positive trends across much of East Asia. Globally, metrics at many sites exhibit non-significant trends. At 59% of all sites there is a common direction and significance in the trend across all five metrics, whilst 4MDA8 and NDGT70 have a common trend at ~80% of all sites. Sensitivity analysis shows AVGMDA8 trends differ with averaging period (warm season or annual). Trends are unchanged at many sites when a 1995–2014 period is used; although fewer sites exhibit non-significant trends. Over the longer period 1970–2014, most Japanese sites exhibit positive 4MDA8/SOMO35 trends. Insufficient data exist to characterize ozone trends for the rest of Asia and other world regions.

The goal of the Tropospheric Ozone Assessment Report (TOAR) is to provide the research community with an up-to-date scientific assessment of tropospheric ozone, from the surface to the tropopause. While a suite of observations provides significant information on the spatial and temporal distribution of tropospheric ozone, observational gaps make it necessary to use global atmospheric chemistry models to synthesize our understanding of the processes and variables that control tropospheric ozone abundance and its variability. Models facilitate the interpretation of the observations and allow us to make projections of future tropospheric ozone and trace gas distributions for different anthropogenic or natural perturbations. This paper assesses the skill of current-generation global atmospheric chemistry models in simulating the observed present-day tropospheric ozone distribution, variability, and trends. Drawing upon the results of recent international multi-model intercomparisons and using a range of model evaluation techniques, we demonstrate that global chemistry models are broadly skillful in capturing the spatio-temporal variations of tropospheric ozone over the seasonal cycle, for extreme pollution episodes, and changes over interannual to decadal periods. However, models are consistently biased high in the northern hemisphere and biased low in the southern hemisphere, throughout the depth of the troposphere, and are unable to replicate particular metrics that define the longer term trends in tropospheric ozone as derived from some background sites. When the models compare unfavorably against observations, we discuss the potential causes of model biases and propose directions for future developments, including improved evaluations that may be able to better diagnose the root cause of the model-observation disparity. Overall, model results should be approached critically, including determining whether the model performance is acceptable for the problem being addressed, whether biases can be tolerated or corrected, whether the model is appropriately constituted, and whether there is a way to satisfactorily quantify the uncertainty.

In support of the first Tropospheric Ozone Assessment Report (TOAR) a relational database of global surface ozone observations has been developed and populated with hourly measurement data and enhanced metadata. A comprehensive suite of ozone data products including standard statistics, health and vegetation impact metrics, and trend information, are made available through a common data portal and a web interface. These data form the basis of the TOAR analyses focusing on human health, vegetation, and climate relevant ozone issues, which are part of this special feature. Cooperation among many data centers and individual researchers worldwide made it possible to build the world’s largest collection of in-situ hourly surface ozone data covering the period from 1970 to 2015. By combining the data from almost 10,000 measurement sites around the world with global metadata information, new analyses of surface ozone have become possible, such as the first globally consistent characterisations of measurement sites as either urban or rural/remote. Exploitation of these global metadata allows for new insights into the global distribution, and seasonal and long-term changes of tropospheric ozone and they enable TOAR to perform the first, globally consistent analysis of present-day ozone concentrations and recent ozone changes with relevance to health, agriculture, and climate. Considerable effort was made to harmonize and synthesize data formats and metadata information from various networks and individual data submissions. Extensive quality control was applied to identify questionable and erroneous data, including changes in apparent instrument offsets or calibrations. Such data were excluded from TOAR data products. Limitations of a posteriori data quality assurance are discussed. As a result of the work presented here, global coverage of surface ozone data for scientific analysis has been significantly extended. Yet, large gaps remain in the surface observation network both in terms of regions without monitoring, and in terms of regions that have monitoring programs but no public access to the data archive. Therefore future improvements to the database will require not only improved data harmonization, but also expanded data sharing and increased monitoring in data-sparse regions.

Surface ozone is a greenhouse gas and pollutant detrimental to human health and crop and ecosystem productivity. The Tropospheric Ozone Assessment Report (TOAR) is designed to provide the research community with an up-to-date observation-based overview of tropospheric ozone’s global distribution and trends. The TOAR Surface Ozone Database contains ozone metrics at thousands of monitoring sites around the world, densely clustered across mid-latitude North America, western Europe and East Asia. Calculating regional ozone trends across these locations is challenging due to the uneven spacing of the monitoring sites across urban and rural areas. To meet this challenge we conducted a spatial and temporal trend analysis of several TOAR ozone metrics across these three regions for summertime (April–September) 2000–2014, using the generalized additive mixed model (GAMM). Our analysis indicates that East Asia has the greatest human and plant exposure to ozone pollution among investigating regions, with increasing ozone levels through 2014. The results also show that ozone mixing ratios continue to decline significantly over eastern North America and Europe, however, there is less evidence for decreases of daytime average ozone at urban sites. The present-day spatial coverage of ozone monitors in East Asia (South Korea and Japan) and eastern North America is adequate for estimating regional trends by simply taking the average of the individual trends at each site. However the European network is more sparsely populated across its northern and eastern regions and therefore a simple average of the individual trends at each site does not yield an accurate regional trend. This analysis demonstrates that the GAMM technique can be used to assess the regional representativeness of existing monitoring networks, indicating those networks for which a regional trend can be obtained by simply averaging the trends of all individual sites and those networks that require a more sophisticated statistical approach.