Quantifying and Reducing Uncertainty in the Processes Controlling Tropospheric Ozone and OH

Lead Research Organisation: Lancaster University
Department Name: Lancaster Environment Centre

Abstract

Understanding the behaviour of hydroxyl (OH) radicals in the troposphere is vital for explaining and predicting atmospheric composition change and its impacts on air quality and climate. The observed atmospheric abundance of ozone and methane has increased substantially over the past century due to human activity, and the fates of these gases are strongly coupled through the short-lived OH radical. However, we do not currently understand the relative importance of the different processes and variables that govern the abundance of these gases. State-of-the-art global chemistry-climate models show differences in methane lifetime of almost a factor of two, preventing them from simulating realistically the observed atmospheric build-up of methane or correctly attributing its causes. These models are also unable to reproduce ozone observations from the late 19th century, or more recent ozone trends observed over the past two decades.

This project addresses these weaknesses by using novel statistical approaches to quantify the sensitivity of OH, O3 and CH4 in global models to the processes and inputs that govern them, and by developing new observational constraints to reduce this uncertainty. We will apply tried and tested emulation methods to reproduce the response of computationally-expensive atmospheric models and to permit a more complete and quantitative assessment of process contributions to uncertainty in trace gas abundance. A unique aspect of this project is that we will apply this approach to five different global models to provide a robust assessment of model responses and to identify the cause of model differences for the first time. Our feasibility studies have successfully demonstrated the effectiveness and value of this approach. Using atmospheric composition measurements we will then develop new multi-variable observational constraints that allow us to reduce the uncertainty in key processes by applying Generalised Likelihood Uncertainty Estimation methods in this field for the first time.

Using these constraints, we will quantify the contribution from changing emissions and climate to changes in O3 and CH4 since the preindustrial era. This will permit the first clear source attribution for changes in radiative forcing from O3 and CH4, informing future IPCC assessments. We will identify the factors required to match observed trends, allowing us to explain why current models fail to reproduce observations. We will apply the same techniques to propagate uncertainties in our understanding of processes and emissions to provide formal uncertainties in projected future O3 and CH4 for given emission pathways. This new analysis approach is timely and benefits greatly from our involvement in the international Chemistry-Climate Model Intiative (CCMI) multi-model assessment of past and future atmospheric composition change, allowing us to explain the diversity of model results and to reduce uncertainty in the resulting projections of atmospheric change.

Planned Impact

The wider beneficiaries of this work include atmospheric scientists, statistical researchers, policy makers and the general public.

Atmospheric scientists will benefit from new quantitative insight into the processes and variables governing atmospheric composition and its changes. While advancing scientific understanding, the project will inform model development and promises more confident prediction of future atmospheric composition change. Development of the Hadley Centre's UKCA model will benefit current and future users of this national model and contribute to UK Earth System Model (UK-ESM) development. Our involvement in the international CCMI activity ensures that our results and methods will reach a global audience in the research community in a timely manner, and provides a pathway to inform IPCC climate, WMO ozone and HTAP pollutant transport assessment reports. The emulators developed in the project will be of value themselves as "toy" models allowing emission and process-based scenarios to be explored without the need for computationally demanding model simulations.

Statistical researchers will benefit from new interest in uncertainty quantification from the atmospheric chemistry and composition community, presenting new challenges and encouraging development of new approaches that will have lasting benefits across the disciplines.

The international policy community will benefit from improved and more confident assessment of atmospheric composition change through IPCC, WMO and HTAP assessments and from the improved capability of modelling tools such as the UKCA and GISS CCMs which are used for these assessments. Reduced uncertainty in the atmospheric response to emission changes provides a more solid foundation for climate policy, and improved understanding of the processes driving background ozone permits more focussed air quality policies addressing this air pollutant.

The wider public has a general interest in global change and air quality, and will benefit from our improved capability to represent atmospheric composition indirectly through more confident policy-making that is based on reduced uncertainty in atmospheric responses.

We will engage with the research community through participation in national and international conferences and CCMI workshops, and through publication of our results in the peer-reviewed literature. We will engage with the general public through a project website, blogs and interaction with the media, and will contribute to schools outreach activities through a focus on atmospheric composition and environmental uncertainty. We are already involved in IPCC, WMO and HTAP assessments, and will ensure that key outcomes of the project reach governmental and policy communities through these fora.

We will organise a 2-day workshop in Lancaster in the final months of the project that focusses on uncertainty in atmospheric composition and brings together the atmospheric composition and statistics communities with representatives from the policy-making arena.

Publications

10 25 50

publication icon
Zhao A (2019) Strong Influence of Aerosol Reductions on Future Heatwaves in Geophysical Research Letters

publication icon
Zhao A (2019) Climate Forcing and Response to Greenhouse Gases, Aerosols, and Ozone in CESM1 in Journal of Geophysical Research: Atmospheres

 
Description We have developed and applied new approaches to exploring uncertainty in atmospheric chemistry models through the use of Gaussian Process Emulation, and published a summary of these approaches and their application. We have used these approaches across a number of models to demonstrate that they can provide new insight into the origin of model differences, and to identify model weaknesses. We have also used them to perform the first calibration of atmospheric chemistry models.
Exploitation Route The techniques developed show substantial promise in allowing us to identify weaknesses in our current understanding as represented in models, and we are keen to explore these further in future work.
Sectors Environment,Other

 
Title Data and R code for "Calibrating a global atmospheric chemistry transport model using Gaussian process emulation and ground-level concentrations of ozone and carbon monoxide" 
Description These are a sample of the csv and R files used to obtain the results for: "Calibrating a global atmospheric chemistry transport model using Gaussian process emulation and ground-level concentrations of ozone and carbon monoxide" V33, V34 and V35 refer to using the synthetic data for O3 only, CO only, and O3 and CO combined. V39 to V41 is the same as above but using the reanalysis data. R scripts are shown for the MCMC algorithm (Figures 7-10) and also for carrying out the global sensitivity analysis (Figures 4-6). Not all R scripts are stored here because there is a lot of replication in the code used. E.g. when running the MCMC algorithm using the synthetic CO data only this is almost identical to the MCMC algorithm code when using the synthetic O3 data only. Similar comments can be made for the csv files. Please e-mail the creator of these R scripts (edryan803@gmail.com) if further explanation is needed. 
Type Of Material Database/Collection of data 
Year Produced 2021 
Provided To Others? Yes  
Impact Data and R code are fully documented and made available; there have been 130 downloads to date, but evidence of substantial impact is not yet clear. 
URL https://zenodo.org/record/4537614
 
Title Fast sensitivity analysis methods for computationally expensive models: Model data 
Description Processed model output from the FRSGC/UCI atmospheric chemistry transport model used in the paper by Ryan et al. entitled: "Fast sensitivity analysis methods for computationally expensive models with multi-dimensional output". 
Type Of Material Database/Collection of data 
Year Produced 2017 
Provided To Others? Yes  
Impact N/A 
URL https://zenodo.org/record/1038670
 
Title Fast sensitivity analysis methods for computationally expensive models: R code 
Description R code to carry out the different global sensitivity analysis methods described in the paper by Ryan et al., 2018, and code to carry out the emulator diagnostic checks. 
Type Of Material Computer model/algorithm 
Year Produced 2017 
Provided To Others? Yes  
Impact The algorithms described here have been adapted for use in new and unrelated projects to explore sensitivity in urban air quality models. 
URL https://zenodo.org/record/1038668
 
Title Global Sensitivity Analysis of Tropospheric Ozone and OH: Budgets from three global chemistry-climate models 
Description This dataset contains monthly mean global atmospheric distributions of ozone mixing ratio and methane chemical loss rate (the largest sink for atmospheric methane is the hydroxyl radical, OH) from 105 model runs with three independent global chemistry climate models. The models include the Frontier Research System for Global Change version of the University of California, Irvine Chemical Transport Model (FRSGC/UCI CTM), the Goddard Institute for Space Studies Global Climate Model (GISS GCM) and the Community Atmosphere Model with Chemistry (CAM-Chem). All three models performed the same simulations for a one-year period (broadly representative of 2001 meteorology) under standardised conditions (40 TgN/yr surface NOx emissions, 5 TgN/yr lightning emissions, 500 TgC/yr biogenic isoprene emissions, 1760 ppb tropospheric methane). An ensemble of 105 simulations was performed that included a control run (run 0), a set of runs that were used to build Gaussian Process emulators (runs 1-80), and additional runs that were used to evaluate the emulators (runs 81-104). A spin-up of six months was performed for each run, and monthly mean model results were archived for the following 12 months. The meteorological conditions used were the same in each simulation. The data supports the exploration of the sensitivity of tropospheric ozone and the chemical lifetime of methane in the troposphere (a proxy for the hydroxyl radical, OH) to eight variables: (1) NOx emissions from all surface sources (range: 30-50 TgN/yr), (2) NO emissions from lightning (range: 2-8 TgN/yr), (3) biogenic isoprene emissions (range: 200-800 TgC/yr), (4) dry deposition rates of all deposited species (range: +/- 80%), (5) wet deposition rates of all soluble species (range: +/- 80%), (6) atmospheric humidity as used in the chemistry scheme only (range: +/- 50%), (7) cloud optical depth (range: factor of 10), and (8) turbulent mixing in the planetary boundary layer (range: factor of 100). The design of the ensemble runs used a Latin Hypercube method to sample this eight-dimensional parameter space to achieve optimal coverage with only 80 simulations. A separate design was used to select an additional 24 simulations to evaluate the emulators built from the standard 80 runs. The specifications for each simulation (total annual emission rate for surface NOx, biogenic isoprene and lightning NO, and the scaling factors applied to native model dry and wet deposition rates, humidity, cloud optical depth, and boundary layer diffusion coefficient) are provided with the dataset. Wild, O.; Voulgarakis, A.; Lamarque, J.-F. (2020): Global Sensitivity Analysis of Tropospheric Ozone and OH: Budgets from three global chemistry-climate models. Centre for Environmental Data Analysis, https://catalogue.ceda.ac.uk/uuid/d5afa10e50b44229b079c7c5a036e660 
Type Of Material Database/Collection of data 
Year Produced 2020 
Provided To Others? Yes  
Impact Database fully documented and made available to others; no records of direct impact yet. 
URL https://catalogue.ceda.ac.uk/uuid/d5afa10e50b44229b079c7c5a036e660
 
Description Workshop on Next Generation of Surrogate Modelling in Environmental Science (Lancaster) 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Other audiences
Results and Impact A 3-day workshop hosted at Lancaster involving a broad range of participants from across the field of applied statistics (academic and company-based) with interests in applications to environmental phenomenon. The aim was to build a community of researchers from the UK and abroad interested in applying and extending the techniques developed in the project, and the workshop involved both presentations and break-out discussion groups. The workshop generated a lot of interest and a range of suggestions for new project ideas and directions involving new collaborations from the participants and beyond.
Year(s) Of Engagement Activity 2018