Description: The first concerted multi-model intercomparison of halogenated very short-lived substances (VSLS) has been performed, within the framework of the ongoing Atmospheric Tracer Transport Model Intercomparison Project (TransCom). Eleven global models or model variants participated, simulating the major natural bromine VSLS, bromoform (CHBr3) and dibromomethane (CH2Br2), over a 20-year period (1993-2012). The overarching goal of TransCom-VSLS was to provide a reconciled model estimate of the stratospheric source gas injection (SGI) of bromine from these gases, to constrain the current measurement-derived range, and to investigate inter-model differences due to emissions and transport processes.Models ran with standardised idealised chemistry, to isolate differences due to transport, and we investigated the sensitivity of results to a range of VSLS emission inventories. Models were tested in their ability to reproduce the observed seasonal and spatial distribution of VSLS at the surface, using measurements from NOAA’s long-term global monitoring network, and in the tropical troposphere, using recent aircraft measurements - including high altitude observations from the NASA Global Hawk platform. The models generally capture the seasonal cycle of surface CHBr3 and CH2Br2 well, with a strong model-measurement correlation (r ≥ 0.7) and a low sensitivity to the choice of emission inventory, at most sites. In a given model, the absolute model-measurement agreement is highly sensitive to the choice of emissions and inter-model differences are also apparent, even when using the same inventory, highlighting the challenges faced in evaluating such inventories at the global scale. Across the ensemble, most consistency is found within the tropics where most of the models (8 out of 11) achieve optimal agreement to surface CHBr3 observations using the lowest of the three CHBr3 emission inventories tested (similarly, 8 out of 11 models for CH2 Br2). In general, the models are able to reproduce well observations of CHBr3 and CH2 Br2 obtained in the tropical tropopause layer (TTL) at various locations throughout the Pacific. Zonal variability in VSLS loading in the TTL is generally consistent among models, with CHBr3 (and to a lesser extent CH2 Br2) most elevated over the tropical West Pacific during boreal winter. The models also indicate the Asian Monsoon during boreal summer to be an important pathway for VSLS reaching the stratosphere, though the strength of this signal varies considerably among models. We derive an ensemble climatological mean estimate of the stratospheric bromine SGI from CHBr3 and CH2 Br2 of 2.0 (1.2-2.5) ppt, ≫ 57% larger than the best estimate from the most recent World Meteorological Organization (WMO) Ozone Assessment Report. We find no evidence for a long-term, transport-driven trend in the stratospheric SGI of bromine over the simulation period. However, transport-driven inter-annual variability in the annual mean bromine SGI is of the order of a ±5%, with SGI exhibiting a strong positive correlation with ENSO in the East Pacific.

Description: We investigate the contribution of oceanic methyl iodide (CH3I) to the stratospheric iodine budget. Based on CH3I measurements from three tropical ship campaigns and the Lagrangian transport model FLEXPART, we provide a detailed analysis of CH3I transport from the ocean surface to the cold point in the upper tropical tropopause layer (TTL). While average oceanic emissions differ by less than 50% from campaign to campaign, the measurements show much stronger variations within each campaign. A positive correlation between the oceanic CH3I emissions and the efficiency of CH3I troposphere–stratosphere transport has been identified for some cruise sections. The mechanism of strong horizontal surface winds triggering large emissions on the one hand and being associated with tropical convective systems, such as developing typhoons, on the other hand, could explain the identified correlations. As a result of the simultaneous occurrence of large CH3I emissions and strong vertical uplift, localized maximum mixing ratios of 0.6 ppt CH3I at the cold point have been determined for observed peak emissions during the SHIVA (Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere)-Sonne research vessel campaign in the coastal western Pacific. The other two campaigns give considerably smaller maxima of 0.1 ppt CH3I in the open western Pacific and 0.03 ppt in the coastal eastern Atlantic. In order to assess the representativeness of the large local mixing ratios, we use climatological emission scenarios to derive global upper air estimates of CH3I abundances. The model results are compared with available upper air measurements, including data from the recent ATTREX and HIPPO2 aircraft campaigns. In the eastern Pacific region, the location of the available measurement campaigns in the upper TTL, the comparisons give a good agreement, indicating that around 0.01 to 0.02 ppt of CH3I enter the stratosphere. However, other tropical regions that are subject to stronger convective activity show larger CH3I entrainment, e.g., 0.08 ppt in the western Pacific. Overall our model results give a tropical contribution of 0.04 ppt CH3I to the stratospheric iodine budget. The strong variations in the geographical distribution of CH3I entrainment suggest that currently available upper air measurements are not representative of global estimates and further campaigns will be necessary in order to better understand the CH3I contribution to stratospheric iodine.

Description: Emissions of halogenated very short-lived substances (VSLS) are poorly constrained. However, their inclusion in global models is required to simulate a realistic inorganic bromine (Bry) loading in both the troposphere, where bromine chemistry perturbs global oxidizing capacity, and in the stratosphere, where it is a major sink for ozone (O3). We have performed simulations using a 3-D chemical transport model (CTM) including three top-down and a single bottom-up derived emission inventory of the major brominated VSLS bromoform (CHBr3) and dibromomethane (CH2Br2). We perform the first concerted evaluation of these inventories, comparing both the magnitude and spatial distribution of emissions. For a quantitative evaluation of each inventory, model output is compared with independent long-term observations at National Oceanic and Atmospheric Administration (NOAA) ground-based stations and with aircraft observations made during the NSF HIAPER Pole-to-Pole Observations (HIPPO) project. For CHBr3, the mean absolute deviation between model and surface observation ranges from 0.22 (38%) to 0.78 (115%) parts per trillion (ppt) in the tropics, depending on emission inventory. For CH2Br2, the range is 0.17 (24%) to 1.25 (167%) ppt. We also use aircraft observations made during the 2011 "Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere" (SHIVA) campaign, in the tropical West Pacific. Here, the performance of the various inventories also varies significantly, but overall the CTM is able to reproduce observed CHBr3 well in the free troposphere using an inventory based on observed sea-to-air fluxes. Finally, we identify the range of uncertainty associated with these VSLS emission inventories on stratospheric bromine loading due to VSLS (BryVSLS). Our simulations show BryVSLS ranges from ~ 4.0 to 8.0 ppt depending on the inventory. We report an optimised estimate at the lower end of this range (~ 4 ppt) based on combining the CHBr3 and CH2Br2 inventories which give best agreement with the compilation of observations in the tropics.

Description: We carry out a case study of the transport and chemistry of bromoform and its product gases (PGs) in a sea breeze driven convective episode on 19 November 2011 along the North West coast of Borneo during the "Stratospheric ozone: Halogen Impacts in a Varying Atmosphere" (SHIVA) campaign. We use ground based, ship, aircraft and balloon sonde observations made during the campaign, and a 3-D regional online transport and chemistry model capable of resolving clouds and convection explicitly that includes detailed bromine chemistry. The model simulates the temperature, wind speed, wind direction fairly well for the most part, and adequately captures the convection location, timing, and intensity. The simulated transport of bromoform from the boundary layer up to 12 km compares well to aircraft observations to support our conclusions. The model makes several predictions regarding bromine transport from the boundary layer to the level of convective detrainment (11 to 12 km). First, the majority of bromine undergoes this transport as bromoform. Second, insoluble organic bromine carbonyl species are transported to between 11 and 12 km, but only form a small proportion of the transported bromine. Third, soluble bromine species, which include bromine organic peroxides, hydrobromic acid (HBr), and hypobromous acid (HOBr), are washed out efficiently within the core of the convective column. Fourth, insoluble inorganic bromine species (principally Br2) are not washed out of the convective column, but are also not transported to the altitude of detrainment in large quantities. We expect that Br2 will make a larger relative contribution to the total vertical transport of bromine atoms in scenarios with higher CHBr3 mixing ratios in the boundary layer, which have been observed in other regions. Finally, given the highly detailed description of the chemistry, transport and washout of bromine compounds within our simulations, we make a series of recommendations about the physical and chemical processes that should be represented in 3-D chemical transport models (CTMs) and chemistry climate models (CCMs), which are the primary theoretical means of estimating the contribution made by CHBr3 and other very short-lived substances (VSLS) to the stratospheric bromine budget.

Description: Marine-produced short-lived trace gases such as dibromomethane (CH2Br2), bromoform (CHBr3), methyliodide (CH3I) and dimethyl sulfide (DMS) significantly impact tropospheric and stratospheric chemistry. Describing their marine emissions in atmospheric chemistry models as accurately as possible is necessary to quantify their impact on ozone depletion and Earth's radiative budget. So far, marine emissions of trace gases have mainly been prescribed from emission climatologies, thus lacking the interaction between the actual state of the atmosphere and the ocean. Here we present simulations with the chemistry climate model EMAC (ECHAM5/MESSy Atmospheric Chemistry) with online calculation of emissions based on surface water concentrations, in contrast to directly prescribed emissions. Considering the actual state of the model atmosphere results in a concentration gradient consistent with model real-time conditions at the ocean surface and in the atmosphere, which determine the direction and magnitude of the computed flux. This method has a number of conceptual and practical benefits, as the modelled emission can respond consistently to changes in sea surface temperature, surface wind speed, sea ice cover and especially atmospheric mixing ratio. This online calculation could enhance, dampen or even invert the fluxes (i.e. deposition instead of emissions) of very short-lived substances (VSLS). We show that differences between prescribing emissions and prescribing concentrations (−28 % for CH2Br2 to +11 % for CHBr3) result mainly from consideration of the actual, time-varying state of the atmosphere. The absolute magnitude of the differences depends mainly on the surface ocean saturation of each particular gas. Comparison to observations from aircraft, ships and ground stations reveals that computing the air–sea flux interactively leads in most of the cases to more accurate atmospheric mixing ratios in the model compared to the computation from prescribed emissions. Calculating emissions online also enables effective testing of different air–sea transfer velocity (k) parameterizations, which was performed here for eight different parameterizations. The testing of these different k values is of special interest for DMS, as recently published parameterizations derived by direct flux measurements using eddy covariance measurements suggest decreasing k values at high wind speeds or a linear relationship with wind speed. Implementing these parameterizations reduces discrepancies in modelled DMS atmospheric mixing ratios and observations by a factor of 1.5 compared to parameterizations with a quadratic or cubic relationship to wind speed.