Search

Abstract During polar spring, periods of elevated tropospheric bromine drive near complete removal of surface ozone. These events impact the tropospheric oxidative capacity and are an area of active research with multiple approaches for representing the underlying processes in global models. We present a method for parameterizing emissions of molecular bromine (Br2) over the Arctic using satellite retrievals of bromine monoxide (BrO) from the Ozone Monitoring Instrument (OMI). OMI retrieves column BrO with daily near global coverage, and we use the GEOS‐Chem chemical mechanism, run online within the Goddard Earth Observing System Earth System Model to identify hotspots of BrO likely associated with polar processes. To account for uncertainties in modeling background BrO, hotspots are only identified where the difference between OMI and modeled columns exceeds a statistical threshold. The resulting hotspot columns are a lower‐limit for the portion of OMI BrO attributable to bromine explosion events. While these hotspots are correlated with BrO measured in the lower troposphere over the Arctic Ocean, a case study of missing detections of near‐surface BrO is identified. Daily flux of Br2 is estimated from hotspot columns of BrO using internal model parameters. When the emissions are applied, BrO hotspots are modeled with a 5% low bias. The sensitivity of the resulting ozone simulations to the treatment of background uncertainties in the BrO column is demonstrated. While periods of isolated, large (>50%) decreases in surface ozone are modeled, this technique does not simulate the low ozone observed at coastal stations and consistently underestimates ozone loss during March.

The sixth phase of the Coupled Model Intercomparison Project (CMIP6) is the latest modeling effort for general circulation models to simulate and project various aspects of climate change. Many of the general circulation models (GCMs) participating in CMIP6 provide archived output that can be used to calculate effective climate sensitivity (ECS) and forecast future temperature change based on emissions scenarios from several Shared Socioeconomic Pathways (SSPs). Here we use our multiple linear regression energy balance model, the Empirical Model of Global Climate (EM-GC), to simulate and project changes in global mean surface temperature (GMST), calculate ECS, and compare to results from the CMIP6 multi-model ensemble. An important aspect of our study is a comprehensive analysis of uncertainties due to radiative forcing of climate from tropospheric aerosols (AER RF) in the EM-GC framework. We quantify the attributable anthropogenic warming rate (AAWR) from the climate record using the EM-GC and use AAWR as a metric to determine how well CMIP6 GCMs replicate human-driven global warming over the last 40 years. The CMIP6 multi-model ensemble indicates a median value of AAWR over 1975–2014 of 0.221 ∘C per decade (range of 0.151 to 0.299 ∘C per decade; all ranges given here are for 5th and 95th confidence intervals), which is notably faster warming than our median estimate for AAWR of 0.157 ∘C per decade (range of 0.120 to 0.195 ∘C per decade) inferred from the analysis of the Hadley Centre Climatic Research Unit version 5 data record for GMST. Estimates of ECS found using the EM-GC assuming that climate feedback does not vary over time (best estimate 2.33 ∘C; range of 1.40 to 3.57 ∘C) are generally consistent with the range of ECS of 1.5 to 4.5 ∘C given by the IPCC's Fifth Assessment Report. The CMIP6 multi-model ensemble exhibits considerably larger values of ECS (median 3.74 ∘C; range of 2.19 to 5.65 ∘C). Our best estimate of ECS increases to 3.08 ∘C (range of 2.23 to 5.53 ∘C) if we allow climate feedback to vary over time. The dominant factor in the uncertainty for our empirical determinations of AAWR and ECS is imprecise knowledge of AER RF for the contemporary atmosphere, though the uncertainty due to time-dependent climate feedback is also important for estimates of ECS. We calculate the likelihood of achieving the Paris Agreement target (1.5 ∘C) and upper limit (2.0 ∘C) of global warming relative to pre-industrial for seven of the SSPs using both the EM-GC and the CMIP6 multi-model ensemble. In our model framework, SSP1-2.6 has a 53 % probability of limiting warming at or below the Paris target by the end of the century, and SSP4-3.4 has a 64 % probability of achieving the Paris upper limit. These estimates are based on the assumptions that climate feedback has been and will remain constant over time since the prior temperature record can be fit so well assuming constant climate feedback. In addition, we quantify the sensitivity of future warming to the curbing of the current rapid growth of atmospheric methane and show that major near-term limits on the future growth of methane are especially important for achievement of the 1.5 ∘C goal of future warming. We also quantify warming scenarios assuming climate feedback will rise over time, a feature common among many CMIP6 GCMs; under this assumption, it becomes more difficult to achieve any specific warming target. Finally, we assess warming projections in terms of future anthropogenic emissions of atmospheric carbon. In our model framework, humans can emit only another 150±79 Gt C after 2019 to have a 66 % likelihood of limiting warming to 1.5 ∘C and another 400±104 Gt C to have the same probability of limiting warming to 2.0 ∘C. Given the estimated emission of 11.7 Gt C per year for 2019 due to combustion of fossil fuels and deforestation, our EM-GC simulations suggest that the 1.5 ∘C warming target of the Paris Agreement will not be achieved unless carbon and methane emissions are severely curtailed in the next 10 years.

To provide insight into the planetary boundary layer (PBL) production of ozone (O3) over the North China Plain, the Air chemistry Research in Asia (ARIAs) campaign conducted aircraft measurements of air pollutants over Hebei Province, China, between May and June 2016. We evaluate vertical profiles of trace gas species including O3, nitrogen oxides (NOx), carbon monoxide (CO), and volatile organic compounds (VOCs) and relate to rates of O3 production. This analysis shows measured O3 levels ranged from 45 to 146 ppbv, with the peak median concentration (∼ 92 ppbv) occurring between 1000 and 1500 m. The NOx concentrations exhibited strong spatial and altitudinal variations, with a maximum of 53 ppbv. Ratios of CO∕CO2 indicate the prevalence of low-efficiency combustion from biomass burning and residential coal burning but indicate some success of regional pollution controls compared to earlier studies in China. Concentrations of total measured VOCs reveal alkanes dominate the total measured volume mixing ratio of VOCs (68 %), and sources include vehicular emissions, fuel and solvent evaporation, and biomass burning. Alkanes and alkenes/alkynes are responsible for 74 % of the total VOC reactivity assessed by calculating the OH loss rates, while aromatics contribute the most to the total ozone formation potential (OFP) (43 %) with toluene, m/p-xylene, ethylene, propylene, and i-pentane playing significant roles in the aloft production of O3 in this region. In the PBL below 500 m, box model calculations constrained by measured precursors indicate the peak rate of mean O3 production was ∼ 7 ppbv h−1. Pollution frequently extended above the PBL into the lower free troposphere around 3000 m, where NO2 mixing ratios (∼ 400 pptv) led to net production rates of O3 up to ∼ 3 ppbv h−1; this pollution can travel substantial distances downwind. The O3 sensitivity regime is determined to be NOx-limited throughout the PBL, whereas it is more VOC-limited at low altitudes near urban centers, demonstrating that control of both VOCs and NOx is needed to reduce aloft O3 pollution over Hebei.

Reduced-complexity climate models (RCMs) are critical in the policy and decision making space, and are directly used within multiple Intergovernmental Panel on Climate Change (IPCC) reports to complement the results of more comprehensive Earth system models. To date, evaluation of RCMs has been limited to a few independent studies. Here we introduce a systematic evaluation of RCMs in the form of the Reduced Complexity Model Intercomparison Project (RCMIP). We expect RCMIP will extend over multiple phases, with Phase 1 being the first. In Phase 1, we focus on the RCMs' global-mean temperature responses, comparing them to observations, exploring the extent to which they emulate more complex models and considering how the relationship between temperature and cumulative emissions of CO2 varies across the RCMs. Our work uses experiments which mirror those found in the Coupled Model Intercomparison Project (CMIP), which focuses on complex Earth system and atmosphere–ocean general circulation models. Using both scenario-based and idealised experiments, we examine RCMs' global-mean temperature response under a range of forcings. We find that the RCMs can all reproduce the approximately 1 ∘C of warming since pre-industrial times, with varying representations of natural variability, volcanic eruptions and aerosols. We also find that RCMs can emulate the global-mean temperature response of CMIP models to within a root-mean-square error of 0.2 ∘C over a range of experiments. Furthermore, we find that, for the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP)-based scenario pairs that share the same IPCC Fifth Assessment Report (AR5)-consistent stratospheric-adjusted radiative forcing, the RCMs indicate higher effective radiative forcings for the SSP-based scenarios and correspondingly higher temperatures when run with the same climate settings. In our idealised setup of RCMs with a climate sensitivity of 3 ∘C, the difference for the ssp585–rcp85 pair by 2100 is around 0.23∘C(±0.12 ∘C) due to a difference in effective radiative forcings between the two scenarios. Phase 1 demonstrates the utility of RCMIP's open-source infrastructure, paving the way for further phases of RCMIP to build on the research presented here and deepen our understanding of RCMs.

Future trajectories of the stratospheric trace gas background will alter the rates of bromine- and chlorine-mediated catalytic ozone destruction via changes in the partitioning of inorganic halogen reservoirs and the underlying temperature structure of the stratosphere. The current formulation of the bromine alpha factor, the ozone-destroying power of stratospheric bromine atoms relative to stratospheric chlorine atoms, is invariant with the climate state. Here, we refactor the bromine alpha factor, introducing normalization to a benchmark chemistry–climate state, and formulate Equivalent Effective Stratospheric Benchmark-normalized Chlorine (EESBnC) to reflect changes in the rates of both bromine- and chlorine-mediated ozone loss catalysis with time. We show that the ozone-processing power of the extrapolar stratosphere is significantly perturbed by future climate assumptions. Furthermore, we show that our EESBnC-based estimate of the extrapolar ozone recovery date is in closer agreement with extrapolar ozone recovery dates predicted using more sophisticated 3-D chemistry–climate models than predictions made using equivalent effective stratospheric chlorine (EESC).

The hydroxyl radical (OH) plays critical roles within the troposphere, such as determining the lifetime of methane (CH4), yet is challenging to model due to its fast cycling and dependence on a multitude of sources and sinks. As a result, the reasons for variations in OH and the resulting methane lifetime (τCH4), both between models and in time, are difficult to diagnose. We apply a neural network (NN) approach to address this issue within a group of models that participated in the Chemistry-Climate Model Initiative (CCMI). Analysis of the historical specified dynamics simulations performed for CCMI indicates that the primary drivers of τCH4 differences among 10 models are the flux of UV light to the troposphere (indicated by the photolysis frequency JO1D), the mixing ratio of tropospheric ozone (O3), the abundance of nitrogen oxides (NOx≡NO+NO2), and details of the various chemical mechanisms that drive OH. Water vapour, carbon monoxide (CO), the ratio of NO:NOx, and formaldehyde (HCHO) explain moderate differences in τCH4, while isoprene, methane, the photolysis frequency of NO2 by visible light (JNO2), overhead ozone column, and temperature account for little to no model variation in τCH4. We also apply the NNs to analysis of temporal trends in OH from 1980 to 2015. All models that participated in the specified dynamics historical simulation for CCMI demonstrate a decline in τCH4 during the analysed timeframe. The significant contributors to this trend, in order of importance, are tropospheric O3, JO1D, NOx, and H2O, with CO also causing substantial interannual variability in OH burden. Finally, the identified trends in τCH4 are compared to calculated trends in the tropospheric mean OH concentration from previous work, based on analysis of observations. The comparison reveals a robust result for the effect of rising water vapour on OH and τCH4, imparting an increasing and decreasing trend of about 0.5 % decade−1, respectively. The responses due to NOx, ozone column, and temperature are also in reasonably good agreement between the two studies.

Abstract Over the last decades, climate science has evolved rapidly across multiple expert domains. Our best tools to capture state‐of‐the‐art knowledge in an internally self‐consistent modeling framework are the increasingly complex fully coupled Earth System Models (ESMs). However, computational limitations and the structural rigidity of ESMs mean that the full range of uncertainties across multiple domains are difficult to capture with ESMs alone. The tools of choice are instead more computationally efficient reduced complexity models (RCMs), which are structurally flexible and can span the response dynamics across a range of domain‐specific models and ESM experiments. Here we present Phase 2 of the Reduced Complexity Model Intercomparison Project (RCMIP Phase 2), the first comprehensive intercomparison of RCMs that are probabilistically calibrated with key benchmark ranges from specialized research communities. Unsurprisingly, but crucially, we find that models which have been constrained to reflect the key benchmarks better reflect the key benchmarks. Under the low‐emissions SSP1‐1.9 scenario, across the RCMs, median peak warming projections range from 1.3 to 1.7°C (relative to 1850–1900, using an observationally based historical warming estimate of 0.8°C between 1850–1900 and 1995–2014). Further developing methodologies to constrain these projection uncertainties seems paramount given the international community's goal to contain warming to below 1.5°C above preindustrial in the long‐term. Our findings suggest that users of RCMs should carefully evaluate their RCM, specifically its skill against key benchmarks and consider the need to include projections benchmarks either from ESM results or other assessments to reduce divergence in future projections.

We infer surface fluxes of bromoform (CHBr3) and dibromoform (CH2Br2) from aircraft observations over the western Pacific using a tagged version of the GEOS-Chem global 3-D atmospheric chemistry model and a maximum a posteriori inverse model. Using GEOS-Chem (GC) as an intermediary, we find that the distribution of a priori ocean emissions of these gases are reasonably consistent with observed atmospheric mole fractions of CHBr3 (r = 0.62) and CH2Br2 (r = 0.38). These a priori emissions result in a positive model bias in CHBr3 peaking in the marine boundary layer, but reproduce observed values of CH2Br2 with no significant bias by virtue of its longer atmospheric lifetime. Using GEOS-Chem, we find that observed variations in atmospheric CHBr3 are determined equally by sources over the western Pacific and those outside the study region, but observed variations in CH2Br2 are determined mainly by sources outside the western Pacific. Numerical closed-loop experiments show that the spatial and temporal distribution of boundary layer aircraft data have the potential to substantially improve current knowledge of these fluxes, with improvements related to data density. Using the aircraft data, we estimate aggregated regional fluxes of 3.6±0.3×108 and 0.7±0.1×108 g month−1 for CHBr3 and CH2Br2 over 130–155°E and 0–12° N, respectively, which represent reductions of 20 %–40 % of the prior inventories by Ordóñez et al. (2012) and substantial spatial deviations from different a priori inventories. We find no evidence to support a robust linear relationship between CHBr3 and CH2Br2 oceanic emissions, as used by previous studies. We find that over regions with dense observation coverage, our choice of a priori inventory does not significantly impact our reported a posteriori flux estimates.

We use the GEOS-Chem global 3-D atmospheric chemistry transport model to interpret atmospheric observations of bromoform (CHBr3) and dibromomethane (CH2Br2) collected during the CAST and CONTRAST aircraft measurement campaigns over the western Pacific, January–February 2014. We use a new linearized, tagged version of CHBr3 and CH2Br2, allowing us to study the influence of emissions from specific geographical regions on observed atmospheric variations. The model describes 32 %–37 % of CHBr3 and 15 %–45 % of CH2Br2 observed variability during CAST and CONTRAST, reflecting model errors in vertical transport. The model has a mean positive bias of 30 % that is larger near the surface, reflecting errors in the poorly constrained prior emission estimates. We find using the model that observed variability of CHBr3 and CH2Br2 is driven by open ocean emissions where there is deep convection. Atmospheric variability above 6 km includes a significant contribution from coastal oceans, but it is still dominated by emissions from the open ocean and by older air masses that originate upwind. In the absence of reliable ocean emission estimates, we use a new physical age-of-air simulation to determine the relative abundance of halogens delivered by CHBr3 and CH2Br2 to the tropical transition layer (TTL). We find that 76 % (92 %) of air masses that originate from the ocean reach the TTL within two (three) atmospheric e-folding lifetimes of CHBr3 and almost all of them reach the TTL within one e-folding lifetime of CH2Br2. Over the duration of CAST and CONTRAST, and over our study region, oceans delivered a mean (range) CHBr3 and CH2Br2 mole fraction of 0.46 (0.13–0.72) and 0.88 (0.71–1.01) pptv, respectively, to the TTL, and a mean (range) Bry mole fraction of 3.14 (1.81–4.18) pptv from source gases to the upper troposphere.

>We analyse simulations performed for the Chemistry-Climate Model Initiative (CCMI) to estimate the return dates of the stratospheric ozone layer from depletion caused by anthropogenic stratospheric chlorine and bromine. We consider a total of 155 simulations from 20 models, including a range of sensitivity studies which examine the impact of climate change on ozone recovery. For the control simulations (unconstrained by nudging towards analysed meteorology) there is a large spread (±20 DU in the global average) in the predictions of the absolute ozone column. Therefore, the model results need to be adjusted for biases against historical data. Also, the interannual variability in the model results need to be smoothed in order to provide a reasonably narrow estimate of the range of ozone return dates. Consistent with previous studies, but here for a Representative Concentration Pathway (RCP) of 6.0, these new CCMI simulations project that global total column ozone will return to 1980 values in 2049 (with a 1σ uncertainty of 2043–2055). At Southern Hemisphere mid-latitudes column ozone is projected to return to 1980 values in 2045 (2039–2050), and at Northern Hemisphere mid-latitudes in 2032 (2020–2044). In the polar regions, the return dates are 2060 (2055–2066) in the Antarctic in October and 2034 (2025–2043) in the Arctic in March. The earlier return dates in the Northern Hemisphere reflect the larger sensitivity to dynamical changes. Our estimates of return dates are later than those presented in the 2014 Ozone Assessment by approximately 5–17 years, depending on the region, with the previous best estimates often falling outside of our uncertainty range. In the tropics only around half the models predict a return of ozone to 1980 values, around 2040, while the other half do not reach the 1980 value. All models show a negative trend in tropical total column ozone towards the end of the 21st century. The CCMI models generally agree in their simulation of the time evolution of stratospheric chlorine and bromine, which are the main drivers of ozone loss and recovery. However, there are a few outliers which show that the multi-model mean results for ozone recovery are not as tightly constrained as possible. Throughout the stratosphere the spread of ozone return dates to 1980 values between models tends to correlate with the spread of the return of inorganic chlorine to 1980 values. In the upper stratosphere, greenhouse gas-induced cooling speeds up the return by about 10–20 years. In the lower stratosphere, and for the column, there is a more direct link in the timing of the return dates of ozone and chlorine, especially for the large Antarctic depletion. Comparisons of total column ozone between the models is affected by different predictions of the evolution of tropospheric ozone within the same scenario, presumably due to differing treatment of tropospheric chemistry. Therefore, for many scenarios, clear conclusions can only be drawn for stratospheric ozone columns rather than the total column. As noted by previous studies, the timing of ozone recovery is affected by the evolution of N2O and CH4. However, quantifying the effect in the simulations analysed here is limited by the few realisations available for these experiments compared to internal model variability. The large increase in N2O given in RCP 6.0 extends the ozone return globally by ∼ 15 years relative to N2O fixed at 1960 abundances, mainly because it allows tropical column ozone to be depleted. The effect in extratropical latitudes is much smaller. The large increase in CH4 given in the RCP 8.5 scenario compared to RCP 6.0 also lengthens ozone return by ∼ 15 years, again mainly through its impact in the tropics. Overall, our estimates of ozone return dates are uncertain due to both uncertainties in future scenarios, in particular those of greenhouse gases, and uncertainties in models. The scenario uncertainty is small in the short term but increases with time, and becomes large by the end of the century. There are still some model–model differences related to well-known processes which affect ozone recovery. Efforts need to continue to ensure that models used for assessment purposes accurately represent stratospheric chemistry and the prescribed scenarios of ozone-depleting substances, and only those models are used to calculate return dates. For future assessments of single forcing or combined effects of CO2, CH4, and N2O on the stratospheric column ozone return dates, this work suggests that it is more important to have multi-member (at least three) ensembles for each scenario from every established participating model, rather than a large number of individual models.

We report measurements of bromine monoxide (BrO) and use an observationally constrained chemical box model to infer total gas-phase inorganic bromine (Bry) over the tropical western Pacific Ocean (tWPO) during the CONTRAST field campaign (January–February 2014). The observed BrO and inferred Bry profiles peak in the marine boundary layer (MBL), suggesting the need for a bromine source from sea-salt aerosol (SSA), in addition to organic bromine (CBry). Both profiles are found to be C-shaped with local maxima in the upper free troposphere (FT). The median tropospheric BrO vertical column density (VCD) was measured as 1.6×1013 molec cm−2, compared to model predictions of 0.9×1013 molec cm−2 in GEOS-Chem (CBry but no SSA source), 0.4×1013 molec cm−2 in CAM-Chem (CBry and SSA), and 2.1×1013 molec cm−2 in GEOS-Chem (CBry and SSA). Neither global model fully captures the C-shape of the Bry profile. A local Bry maximum of 3.6 ppt (2.9–4.4 ppt; 95 % confidence interval, CI) is inferred between 9.5 and 13.5 km in air masses influenced by recent convective outflow. Unlike BrO, which increases from the convective tropical tropopause layer (TTL) to the aged TTL, gas-phase Bry decreases from the convective TTL to the aged TTL. Analysis of gas-phase Bry against multiple tracers (CFC-11, H2O ∕ O3 ratio, and potential temperature) reveals a Bry minimum of 2.7 ppt (2.3–3.1 ppt; 95 % CI) in the aged TTL, which agrees closely with a stratospheric injection of 2.6 ± 0.6 ppt of inorganic Bry (estimated from CFC-11 correlations), and is remarkably insensitive to assumptions about heterogeneous chemistry. Bry increases to 6.3 ppt (5.6–7.0 ppt; 95 % CI) in the stratospheric "middleworld" and 6.9 ppt (6.5–7.3 ppt; 95 % CI) in the stratospheric "overworld". The local Bry minimum in the aged TTL is qualitatively (but not quantitatively) captured by CAM-Chem, and suggests a more complex partitioning of gas-phase and aerosol Bry species than previously recognized. Our data provide corroborating evidence that inorganic bromine sources (e.g., SSA-derived gas-phase Bry) are needed to explain the gas-phase Bry budget in the upper free troposphere and TTL. They are also consistent with observations of significant bromide in Upper Troposphere–Lower Stratosphere aerosols. The total Bry budget in the TTL is currently not closed, because of the lack of concurrent quantitative measurements of gas-phase Bry species (i.e., BrO, HOBr, HBr, etc.) and aerosol bromide. Such simultaneous measurements are needed to (1) quantify SSA-derived Bry in the upper FT, (2) test Bry partitioning, and possibly explain the gas-phase Bry minimum in the aged TTL, (3) constrain heterogeneous reaction rates of bromine, and (4) account for all of the sources of Bry to the lower stratosphere.

Regulatory air quality models, such as the Community Multiscale Air Quality model (CMAQ), are used by federal and state agencies to guide policy decisions that determine how to best achieve adherence with National Ambient Air Quality Standards for surface ozone. We use observations of ozone and its important precursor NO2 to test the representation of the photochemistry and emission of ozone precursors within CMAQ. Observations of tropospheric column NO2 from the Ozone Monitoring Instrument (OMI), retrieved by two independent groups, show that the model overestimates urban NO2 and underestimates rural NO2 under all conditions examined for July and August 2011 in the US Northeast. The overestimate of the urban to rural ratio of tropospheric column NO2 for this baseline run of CMAQ (CB05 mechanism, mobile NOx emissions from the National Emissions Inventory; isoprene emissions from MEGAN v2.04) suggests this model may underestimate the importance of interstate transport of NOx. This CMAQ simulation leads to a considerable overestimate of the 2-month average of 8 h daily maximum surface ozone in the US Northeast, as well as an overestimate of 8 h ozone at AQS sites during days when the state of Maryland experienced NAAQS exceedances. We have implemented three changes within CMAQ motivated by OMI NO2 as well as aircraft observations obtained in July 2011 during the NASA DISCOVER-AQ campaign: (a) the modeled lifetime of organic nitrates within CB05 has been reduced by a factor of 10, (b) emissions of NOx from mobile sources has been reduced by a factor of 2, and (c) isoprene emissions have been reduced by using MEGAN v2.10 rather than v2.04. Compared to the baseline simulation, the CMAQ run using all three of these changes leads to considerably better simulation of column NO2 in both urban and rural areas, better agreement with the 2-month average of daily 8 h maximum ozone in the US Northeast, fewer number of false positives of an ozone exceedance throughout the domain, as well as an unbiased simulation of surface ozone at ground-based AQS sites in Maryland that experienced an ozone exceedance during July and August 2007. These modifications to CMAQ may provide a framework for use in studies focused on achieving future adherence to specific air quality standards for surface ozone by reducing emission of NOx from various anthropogenic sectors.

Very short-lived (VSL) bromocarbons are produced at a prodigious rate by ocean biology and these source compounds (SGVSL), together with their inorganic degradation products (PGVSL), are lofted by vigorous convection to the tropical tropopause layer (TTL). Using a state-of-the-art photochemical mechanism within a global model, we calculate annual average stratospheric injection of total bromine due to VSL sources to be 5 pptv (parts per trillion by volume), with ~ 3 pptv entering the stratosphere as PGVSL and ~ 2 pptv as SGVSL. The geographic distribution and partitioning of VSL bromine within the TTL, and its consequent stratospheric injection, is highly dependent on the oceanic flux, the strength of convection and the occurrence of heterogeneous recycling reactions. Our calculations indicate atomic Br should be the dominant inorganic species in large regions of the TTL during daytime, due to the low ozone and cold conditions of this region. We propose the existence of a "tropical ring of atomic bromine" located approximately between 15 and 19 km and between 30° N and 30° S. Daytime Br / BrO ratios of up to ~ 4 are predicted within this inhomogeneous ring in regions of highly convective transport, such as the tropical Western Pacific. Therefore, we suggest that experimental programs designed to quantify the bromine budget of the TTL and the stratospheric injection of VSL biogenic bromocarbons should include a strategy for the measurement of atomic Br during daytime as well as HOBr and BrCl during nighttime.

We present a new retrieval of stratospheric BrO (bromine monoxide) from channel 2 SCIAMACHY (SCanning Imaging Absorption spectrometer for Atmospheric CHartographY) limb observations. Retrievals are shown to agree with independent balloon observations to within one standard deviation of the retrieval noise. We retrieve BrO profiles for all of April 2008, and apply simulated [BrO]/[Bry] (bromine monoxide : stratospheric inorganic bromine) ratios to estimate the stratospheric Bry loading. We find 23.5 ± 6 ppt Br, suggesting 7 ppt Br from short-lived bromocarbons to be at the high end of the current best estimate (3–8 ppt). The 6 ppt Br uncertainty estimate is dominated by the 21% uncertainty in the simulated [BrO] / [Bry] ratio due to propagation of errors from the underlying chemical kinetics.

Trends in the composition of the lower atmosphere (0–1500 m altitude) and surface air quality over the Baltimore/Washington area and surrounding states were investigated for the period from 1997 to 2011. We examined emissions of ozone precursors from monitors and inventories as well as ambient ground-level and aircraft measurements to characterize trends in air pollution. The US EPA Continuous Emissions Monitoring System (CEMS) program reported substantial decreases in emission of summertime nitrogen oxides (NOx) from power plants, up to ∼80% in the mid-Atlantic States. These large reductions in emission of NOx are reflected in a sharp decrease of ground-level concentrations of NOx starting around 2003. The decreasing trend of tropospheric column CO observed by aircraft is ∼0.8 Dobson unit (DU) per year, corresponding to ∼35 ppbv yr−1 in the lower troposphere (the surface to 1500 m above ground level). Satellite observations of long-term, near-surface CO show a ∼40% decrease over western Maryland between 2000 and 2011; the same magnitude is indicated by aircraft measurements above these regions upwind of the Baltimore/Washington airshed. With decreasing emissions of ozone precursors, the ground-level ozone in the Baltimore/Washington area shows a 0.6 ppbv yr−1 decrease in the past 15 yr. Since photochemical production of ozone is substantially influenced by ambient temperature, we introduce the climate penalty factor (CPF) into the trend analysis of long-term aircraft measurements. After compensating for inter-annual variations in temperature, historical aircraft measurements indicate that the daily net production of tropospheric ozone over the Baltimore/Washington area decreased from ∼20 ppbv day−1 in the late 1990s to ∼7 ppbv day−1 in the early 2010s during ozone season. A decrease in the long-term column ozone is observed as ∼0.2 DU yr−1 in the lowest 1500 m, corresponding to an improvement of ∼1.3 ppbv yr−1. Our aircraft measurements were conducted on days when severe ozone pollution was forecasted, and these results represent the decreasing trend in high ozone events over the past 15 yr. Back trajectory cluster analysis demonstrates that emissions of air pollutants from Ohio and Pennsylvania through Maryland influence the column abundances of downwind ozone in the lower atmosphere. The trends in air pollutants reveal the success of regulations implemented over the past decades and the importance of region-wide emission controls in the eastern United States.

Observed reductions in Earth's surface temperature following explosive volcanic eruptions have been used as a proxy for geoengineering of climate by the artificial enhancement of stratospheric sulfate. Earth cools following major eruptions due to an increase in the reflection of sunlight caused by a dramatic enhancement of the stratospheric sulfate aerosol burden. Significant global cooling has been observed following the four major eruptions since 1900: Santa María, Mount Agung, El Chichón and Mt. Pinatubo, leading IPCC (2007) to state "major volcanic eruptions can, thus, cause a drop in global mean surface temperature of about half a degree Celsius that can last for months and even years". We use a multiple linear regression model applied to the global surface temperature anomaly to suggest that exchange of heat between the atmosphere and ocean, driven by variations in the strength of the Atlantic Meridional Overturning Circulation (AMOC), has been a factor in the decline of global temperature following these eruptions. The veracity of this suggestion depends on whether sea surface temperature (SST) in the North Atlantic, sometimes called the Atlantic Multidecadal Oscillation, but here referred to as Atlantic Multidecadal Variability (AMV), truly represents a proxy for the strength of the AMOC. Also, precise quantification of global cooling due to volcanoes depends on how the AMV index is detrended. If the AMV index is detrended using anthropogenic radiative forcing of climate, we find that surface cooling attributed to Mt. Pinatubo, using the Hadley Centre/University of East Anglia surface temperature record, maximises at 0.14 °C globally and 0.32 °C over land. These values are about a factor of 2 less than found when the AMV index is neglected in the model and quite a bit lower than the canonical 0.5 °C cooling usually attributed to Pinatubo. This result is driven by the high amplitude, low frequency component of the AMV index, demonstrating that reduced impact of volcanic cooling upon consideration of the AMV index is driven by variations in North Atlantic SST that occur over time periods much longer than those commonly associated with major volcanic eruptions. The satellite record of atmospheric temperature from 1978 to present and other century-long surface temperature records are also consistent with the suggestion that volcanic cooling may have been over estimated by about a factor of 2 due to prior neglect of ocean circulation. Our study suggests a recalibration may be needed for the proper use of Mt. Pinatubo as a proxy for geoengineering of climate. Finally, we highlight possible shortcomings in simulations of volcanic cooling by general circulation models, which are also being used to assess the impact of geoengineering of climate via stratospheric sulfate injection.

A focus of the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) mission was examination of bromine photochemistry in the spring time high latitude troposphere based on aircraft and satellite measurements of bromine oxide (BrO) and related species. The NASA DC-8 aircraft utilized a chemical ionization mass spectrometer (CIMS) to measure BrO and a mist chamber (MC) to measure soluble bromide. We have determined that the MC detection efficiency to molecular bromine (Br2), hypobromous acid (HOBr), bromine oxide (BrO), and hydrogen bromide (HBr) as soluble bromide (Br−) was 0.9±0.1, 1.06+0.30/−0.35, 0.4±0.1, and 0.95±0.1, respectively. These efficiency factors were used to estimate soluble bromide levels along the DC-8 flight track of 17 April 2008 from photochemical calculations constrained to in situ BrO measured by CIMS. During this flight, the highest levels of soluble bromide and BrO were observed and atmospheric conditions were ideal for the space-borne observation of BrO. The good agreement (R2 = 0.76; slope = 0.95; intercept = −3.4 pmol mol−1) between modeled and observed soluble bromide, when BrO was above detection limit (>2 pmol mol−1) under unpolluted conditions (NO−1), indicates that the CIMS BrO measurements were consistent with the MC soluble bromide and that a well characterized MC can be used to derive mixing ratios of some reactive bromine compounds. Tropospheric BrO vertical column densities (BrOVCD) derived from CIMS BrO observations compare well with BrOTROPVCD from OMI on 17 April 2008.

We derive tropospheric column BrO during the ARCTAS and ARCPAC field campaigns in spring 2008 using retrievals of total column BrO from the satellite UV nadir sensors OMI and GOME-2 using a radiative transfer model and stratospheric column BrO from a photochemical simulation. We conduct a comprehensive comparison of satellite-derived tropospheric BrO column to aircraft in-situ observations of BrO and related species. The aircraft profiles reveal that tropospheric BrO, when present during April 2008, was distributed over a broad range of altitudes rather than being confined to the planetary boundary layer (PBL). Perturbations to the total column resulting from tropospheric BrO are the same magnitude as perturbations due to longitudinal variations in the stratospheric component, so proper accounting of the stratospheric signal is essential for accurate determination of satellite-derived tropospheric BrO. We find reasonably good agreement between satellite-derived tropospheric BrO and columns found using aircraft in-situ BrO profiles, particularly when satellite radiances were obtained over bright surfaces (albedo >0.7), for solar zenith angle

We describe a method of evaluating systematic errors in measurements of total column dry-air mole fractions of CO2 (XCO2) from space, and we illustrate the method by applying it to the v2.8 Atmospheric CO2 Observations from Space retrievals of the Greenhouse Gases Observing Satellite (ACOS-GOSAT) measurements over land. The approach exploits the lack of large gradients in XCO2 south of 25° S to identify large-scale offsets and other biases in the ACOS-GOSAT data with several retrieval parameters and errors in instrument calibration. We demonstrate the effectiveness of the method by comparing the ACOS-GOSAT data in the Northern Hemisphere with ground truth provided by the Total Carbon Column Observing Network (TCCON). We use the observed correlation between free-tropospheric potential temperature and XCO2 in the Northern Hemisphere to define a dynamically informed coincidence criterion between the ground-based TCCON measurements and the ACOS-GOSAT measurements. We illustrate that this approach provides larger sample sizes, hence giving a more robust comparison than one that simply uses time, latitude and longitude criteria. Our results show that the agreement with the TCCON data improves after accounting for the systematic errors, but that extrapolation to conditions found outside the region south of 25° S may be problematic (e.g., high airmasses, large surface pressure biases, M-gain, measurements made over ocean). A preliminary evaluation of the improved v2.9 ACOS-GOSAT data is also discussed.

Key kinetic parameters governing the partitioning of chlorine species in the Antarctic polar stratosphere were retrieved from 28 days of chlorine monoxide (ClO) microwave radiometer measurements made during the late winter/early spring of 2005 at Scott Base (77.85° S, 166.75° E). During day-time the loss of the ClO dimer chlorine peroxide (ClOOCl) occurs mainly by photolysis. Some time after sunrise, a photochemical equilibrium is established and the ClO/ClOOCl partitioning is determined by the ratio of the photolysis frequency, J, and the dimer formation rate, kf. The values of J and kf from laboratory studies remain uncertain to a considerable extent, and as a complement to these ongoing studies, the goal of this work is to provide a constraint on that uncertainty based on observations of ClO profiles in the Antarctic. First an optimal estimation technique was used to derive J/kf ratios for a range of Keq values. The optimal estimation forward model was a photochemical box model that takes J, kf, and Keq as inputs, together with a priori profiles of activated chlorine (ClOx = ClO+2×ClOOCl), profiles of ozone, temperature, and pressure. JPL06 kinetics are used as a priori in the optimal estimation and for all other chemistry in the forward model. Using the more recent JPL09 kinetics results in insignificant differences in the retrieved value of J/kf. A complementary approach was used to derive the optimal kinetic parameters; the full parameter space of J, kf, Keq and ClOx was sampled to find the minimum in differences between measured and modelled ClO profiles. Furthermore, values of Keq up to 2.0 times larger than recommended by JPL06 were explored to test the sensitivity of the J/kf ratio to changes in Keq. The results show that the retrieved J/kf ratios bracket the range of 1.23 to 1.97 times the J/kf value recommended by JPL06 over the range of Keq values considered. The retrieved J/kf ratios lie in the lower half of the large uncertainty range of J/kf recommended by JPL06 and towards the upper portion of the smaller uncertainty range recommended by JPL09.

The column abundance of NO3 was measured over Table Mountain Facility, CA (34.4° N, 117.7° W) from May 2003 through September 2004, using lunar occultation near full moon with a grating spectrometer. The NO3 column retrieval was performed with the differential optical absorption spectroscopy (DOAS) technique using both the 623 and 662 nm NO3 absorption bands. Other spectral features such as Fraunhofer lines and absorption from water vapor and oxygen were removed using solar spectra obtained at different airmass factors. We observed a seasonal variation, with nocturnally averaged NO3 columns between 5−7 × 1013 molec cm−2 during October through March, and 5−22 × 1013 molec cm−2 during April through September. A subset of the data, with diurnal variability vastly different from the temporal profile obtained from one-dimensional stratospheric model calculations, clearly has boundary layer contributions; this was confirmed by simultaneous long-path DOAS measurements. However, even the NO3 columns that did follow the modeled time evolution were often much larger than modeled stratospheric partial columns constrained by realistic temperatures and ozone concentrations. This discrepancy is attributed to substantial tropospheric NO3 in the free troposphere, which may have the same time dependence as stratospheric NO3.

Chemical ozone loss in winter 1991–1992 is recalculated based on observations of the HALOE satellite instrument, Version 19, ER-2 aircraft measurements and balloon data. HALOE satellite observations are shown to be reliable in the lower stratosphere below 400 K, at altitudes where the measurements are most likely disturbed by the enhanced sulfate aerosol loading, as a result of the Mt.~Pinatubo eruption in June 1991. Significant chemical ozone loss (13–17 DU) is observed below 380 K from Kiruna balloon observations and HALOE satellite data between December 1991 and March 1992. For the two winters after the Mt. Pinatubo eruption, HALOE satellite observations show a stronger extent of chemical ozone loss towards lower altitudes compared to other Arctic winters between 1991 and 2003. In spite of already occurring deactivation of chlorine in March 1992, MIPAS-B and LPMA balloon observations indicate that chlorine was still activated at lower altitudes, consistent with observed chemical ozone loss occurring between February and March and April. Large chemical ozone loss of more than 70 DU in the Arctic winter 1991–1992 as calculated in earlier studies is corroborated here.

Among the major factors controlling ozone loss in the polar vortices in winter/spring is the kinetics of the ClO dimer catalytic cycle. Here, we propose a strategy to test and improve our understanding of these kinetics by comparing and combining information on the thermal equilibrium between ClO and Cl2O2, the rate of Cl2O2 formation, and the Cl2O2 photolysis rate from laboratory experiments, theoretical studies and field observations. Concordant with a number of earlier studies, we find considerable inconsistencies of some recent laboratory results with rate theory calculations and stratospheric observations of ClO and Cl2O2. The set of parameters for which we find the best overall consistency – namely the ClO/Cl2O2 equilibrium constant suggested by Plenge et al. (2005), the Cl2O2 recombination rate constant reported by Nickolaisen et al. (1994) and Cl2O2 photolysis rates based on absorption cross sections in the range between the JPL 2006 assessment and the laboratory study by Burkholder et al. (1990) – is not congruent with the latest recommendations given by the JPL and IUPAC panels and does not represent the laboratory studies currently regarded as the most reliable experimental values. We show that the incorporation of new Pope et al. (2007) Cl2O2 absorption cross sections into several models, combined with best estimates for other key parameters (based on either JPL and IUPAC evaluations or on our study), results in severe model underestimates of observed ClO and observed ozone loss rates. This finding suggests either the existence of an unknown process that drives the partitioning of ClO and Cl2O2, or else some unidentified problem with either the laboratory study or numerous measurements of atmospheric ClO. Our mechanistic understanding of the ClO/Cl2O2 system is grossly lacking, with severe implications for our ability to simulate both present and future polar ozone depletion.

Airborne in situ observations of molecules with a wide range of lifetimes (methane, nitrous oxide, reactive nitrogen, ozone, chlorinated halocarbons, and halon-1211), used in a tropical tracer model, show that mid-latitude air is entrained into the tropical lower stratosphere within about 13.5 months; transport is faster in the reverse direction. Because exchange with the tropics is slower than global photochemical models generally assume, ozone at mid-latitudes appears to be more sensitive to elevated levels of industrial chlorine than is currently predicted. Nevertheless, about 45 percent of air in the tropical ascent region at 21 kilometers is of mid-latitude origin, implying that emissions from supersonic aircraft could reach the middle stratosphere.