Air-Sea Exchange in the Global Mercury Cycle

We present results from a new global atmospheric mercury model coupled with a mixed layer slab ocean. The ocean model describes the interactions of the mixed layer with the atmosphere and deep ocean, as well as conversion between elemental, divalent, and nonreactive mercury species. Our global mean aqueous concentrations of 0.07 pM elemental, 0.80 pM reactive, and 1.51 pM total mercury agree with observations. The ocean provides a 14.1 Mmol yr−1 source of mercury to the atmosphere, at the upper end of previous estimates. Re-emission of previously deposited mercury constitutes 89% of this flux. Ocean emissions are largest in the tropics and downwind of industrial regions. Midlatitude ocean emissions display a large seasonal cycle induced by biological productivity. Oceans contribute 54% (36%) of surface atmospheric mercury in the Southern (Northern) Hemisphere. We find a large net loss of mercury to the deep ocean (8.7 Mmol yr−1), implying a ∼0.7%/year increase in deep ocean concentrations.Earth and Planetary SciencesEngineering and Applied Science

Remote sensed and in situ constraints on processes affecting tropical tropospheric ozone

We use a global chemical transport model (GEOS-Chem) to evaluate the consistency of satellite measurements of lightning flashes and ozone precursors with in situ measurements of tropical tropospheric ozone. The measurements are tropospheric O₃, NO₂, and HCHO columns from the GOME satellite instrument, lightning flashes from the OTD and LIS satellite instruments, profiles of O₃, CO, and relative humidity from the MOZAIC aircraft program, and profiles of O₃ from the SHADOZ ozonesonde network. We interpret these multiple data sources with our model to better understand what controls tropical tropospheric ozone. Tropical tropospheric ozone is mainly affected by lightning NO_x and convection in the upper troposphere and by surface emissions in the lower troposphere. Scaling the spatial distribution of lightning in the model to the observed flashes improves the simulation of O₃ in the upper troposphere by 5–20 ppbv versus in situ observations and by 1–4 Dobson Units versus GOME retrievals of tropospheric O₃ columns. A lightning source strength of 6±2 Tg N/yr best represents in situ observations from aircraft and ozonesonde. Tropospheric NO₂ and HCHO columns from GOME are applied to provide top-down constraints on emission inventories of NO_x (biomass burning and soils) and VOCs (biomass burning). The top-down biomass burning inventory is larger than the bottom-up inventory by a factor of 2 for HCHO and alkenes, and by a factor of 2.6 for NO_x over northern equatorial Africa. These emissions increase lower tropospheric O₃ by 5–20 ppbv, improving the simulation versus aircraft observations, and by 4 Dobson Units versus GOME observations of tropospheric O₃ columns. Emission factors in the a posteriori inventory are more consistent with a recent compilation from in situ measurements. The ozone simulation using two different dynamical schemes (GEOS-3 and GEOS-4) is evaluated versus observations; GEOS-4 better represents O₃ observations by 5–15 ppbv, reflecting enhanced convective detrainment in the upper troposphere. Heterogeneous uptake of HNO₃ on aerosols reduces simulated O₃ by 5–7 ppbv, reducing a model bias versus in situ observations over and downwind of deserts. Exclusion of HO₂ uptake on aerosols increases O₃ by 5 ppbv in biomass burning regions, reducing a model bias versus MOZAIC aircraft measurements

Chemistry of hydrogen oxide radicals (HO_x) in the Arctic troposphere in spring

We use observations from the April 2008 NASA ARCTAS aircraft campaign to the North American Arctic, interpreted with a global 3-D chemical transport model (GEOS-Chem), to better understand the sources and cycling of hydrogen oxide radicals (HO_x≡H+OH+peroxy radicals) and their reservoirs (HO_y≡HO_x+peroxides) in the springtime Arctic atmosphere. We find that a standard gas-phase chemical mechanism overestimates the observed HO_2 and H_2O_2 concentrations. Computation of HO_x and HO_y gas-phase chemical budgets on the basis of the aircraft observations also indicates a large missing sink for both. We hypothesize that this could reflect HO_2 uptake by aerosols, favored by low temperatures and relatively high aerosol loadings, through a mechanism that does not produce H_2O_2. We implemented such an uptake of HO_2 by aerosol in the model using a standard reactive uptake coefficient parameterization with γ(HO_2) values ranging from 0.02 at 275 K to 0.5 at 220 K. This successfully reproduces the concentrations and vertical distributions of the different HO_x species and HO_y reservoirs. HO_2 uptake by aerosol is then a major HO_x and HO_y sink, decreasing mean OH and HO_2 concentrations in the Arctic troposphere by 32% and 31% respectively. Better rate and product data for HO_2 uptake by aerosol are needed to understand this role of aerosols in limiting the oxidizing power of the Arctic atmosphere

Observed OH and HO_2 in the upper troposphere suggest a major source from convective injection of peroxides

ER-2 aircraft observations of OH and HO_2 concentrations in the upper troposphere during the NASA/STRAT campaign are interpreted using a photochemical model constrained by local observations of O_3, H_2O, NO, CO, hydrocarbons, albedo and overhead ozone column. We find that the reaction Q(^(1)D) + H_2O is minor compared to acetone photolysis as a primary source of HO_x (= OH + peroxy radicals) in the upper troposphere. Calculations using a diel steady state model agree with observed HO_x concentrations in the lower stratosphere and, for some flights, in the upper troposphere. However, for other flights in the upper troposphere, the steady state model underestimates observations by a factor of 2 or more. These model underestimates are found to be related to a recent (< 1 week) convective origin of the air. By conducting time-dependent model calculations along air trajectories determined for the STRAT flights, we show that convective injection of CH_3OOH and H_2O_2 from the boundary layer to the upper troposphere could resolve the discrepancy. These injections of HO_x reservoirs cause large HO_x increases in the tropical upper troposphere for over a week downwind of the convective activity. We propose that this mechanism provides a major source of HO_x in the upper troposphere. Simultaneous measurements of peroxides, formaldehyde and acetone along with OH and HO_2 are needed to test our hypothesis

Nested-grid simulation of mercury over North America

We have developed a new nested-grid mercury (Hg) simulation over North America with a 1/2° latitude by 2/3° longitude horizontal resolution employing the GEOS-Chem global chemical transport model. Emissions, chemistry, deposition, and meteorology are self-consistent between the global and nested domains. Compared to the global model (4° latitude by 5° longitude), the nested model shows improved skill at capturing the high spatial and temporal variability of Hg wet deposition over North America observed by the Mercury Deposition Network (MDN) in 2008–2009. The nested simulation resolves features such as higher deposition due to orographic precipitation, land/ocean contrast and and predicts more efficient convective rain scavenging of Hg over the southeast United States. However, the nested model overestimates Hg wet deposition over the Ohio River Valley region (ORV) by 27%. We modify anthropogenic emission speciation profiles in the US EPA National Emission Inventory (NEI) to account for the rapid in-plume reduction of reactive to elemental Hg (IPR simulation). This leads to a decrease in the model bias to −2.3% over the ORV region. Over the contiguous US, the correlation coefficient (<i>r</i>) between MDN observations and our IPR simulation increases from 0.60 to 0.78. The IPR nested simulation generally reproduces the seasonal cycle in surface concentrations of speciated Hg from the Atmospheric Mercury Network (AMNet) and Canadian Atmospheric Mercury Network (CAMNet). In the IPR simulation, annual mean gaseous and particulate-bound Hg(II) are within 140% and 11% of observations, respectively. In contrast, the simulation with unmodified anthropogenic Hg speciation profiles overestimates these observations by factors of 4 and 2 for gaseous and particulate-bound Hg(II), respectively. The nested model shows improved skill at capturing the horizontal variability of Hg observed over California during the ARCTAS aircraft campaign. The nested model suggests that North American anthropogenic emissions account for 10–22% of Hg wet deposition flux over the US, depending on the anthropogenic emissions speciation profile assumed. The modeled percent contribution can be as high as 60% near large point sources in ORV. Our results indicate that the North American anthropogenic contribution to dry deposition is 13–20%

Evolution and stoichiometry of heterogeneous processing in the Antarctic stratosphere

Simultaneous in situ measurements of HCl and ClO have been made for the first time in the southern hemisphere, allowing a systematic study of the processes governing chlorine activation between 15 and 20 km in the 1994 Antarctic winter. Data for several other gases (O_3, NO, NO_y, OH, HO_2, N_(2)O, CH_4, CO, H_(2)O, CFCs), particulates, and meteorological parameters were collected from the ER-2 aircraft out of New Zealand as part of the 1994 Airborne Southern Hemisphere Ozone Experiment/Measurements of Atmospheric Effects of Stratospheric Aircraft (ASHOE/MAESA) campaign. Observations from the ER-2 in the fall (April–May), prior to polar night, show that chlorine activation begins with 60–75% of inorganic chlorine as HCl. By midwinter (July–August), near-total removal of HCl is observed. The wintertime loss of HCl in air recently exposed to extreme temperatures is found to be correlated with high levels of reactive chlorine (ClO and its dimer, Cl_(2)O_2) in the linear fashion expected from the stoichiometry of the heterogeneous reaction of hydrochloric acid with chlorine nitrate on polar stratospheric clouds (PSCs): HCl + ClONO_2 → Cl_2 + HNO_3. To constrain the role of different heterogeneous reactions and PSC types, we have used a photochemical trajectory model which includes heterogeneous sulfate and PSC chemistry. Model calculations of the evolution of reactive gases are compared with the in situ observations. In addition, simultaneous measurements of OH and HO_2 are used as a diagnostic for the occurrence of the heterogeneous reaction HOCl + HCl → Cl_2 + H_(2)O, which contributes to suppressed levels of HO_x inside the vortex. It is shown that the amount of chlorine activation is not strongly dependent on the composition of PSCs. However, HO_x levels exhibit different signatures depending on the type of heterogeneous surfaces that affected chlorine activation. Furthermore, this analysis implies that in the edge region of the Antarctic vortex, the observed near-total removal of HCl can result from latitudinal excursions of air parcels in and out of sunlight during the winter, which photochemically resupply HOCl and ClONO_2 as oxidation partners for HCl

Vertical distribution of mercury, CO, ozone, and aerosol scattering coefficient in the Pacific Northwest during the spring 2006 INTEX‐B campaign

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/94827/1/jgrd14532.pd

NO(y) partitioning from measurements of nitrogen and hydrogen radicals in the upper troposphere

Recent studies using NO, NO(y), OH and HO2 (HO(X)) observations have postulated acetone and convection of peroxides as significant sources of HO(X) in the upper troposphere (UT). This work focuses on the effect these additional HO(X) sources have on the modeled NO(y) partitioning and comparisons of the modeled NO(x)/NO(y) ratio to observations. The measured NO(x)/NO(y) ratio is usually much higher than predicted regardless of the presence of acetone in the model. The exception occurs for air parcels having low NO(y) and O3 values. For these air parcels the measured NO(x)/NO(y) ratio is much lower than the calculated ratio unless acetone is included in the model. In all cases acetone increases the fraction of NO(y) that is peroxy acetyl nitrate (PAN) from typical values of much less than 0.1 to values as high as 0.35. Including acetone also reduces the scatter in a comparison between modeled and observed NO(x)/NO(y) ratios

Statistical inference of OH concentrations and air mass dilution rates from successive observations of non-methane hydrocarbons in single air masses

Bayesian inference has been used to determine rigorous estimates of hydroxyl radical concentrations () and air mass dilution rates (K) averaged following air masses between linked observations of nonmethane hydrocarbons (NMHCs) spanning the North Atlantic during the Intercontinental Transport and Chemical Transformation (ITCT)-Lagrangian-2K4 experiment. The Bayesian technique obtains a refined (posterior) distribution of a parameter given data related to the parameter through a model and prior beliefs about the parameter distribution. Here, the model describes hydrocarbon loss through OH reaction and mixing with a background concentration at rate K. The Lagrangian experiment provides direct observations of hydrocarbons at two time points, removing assumptions regarding composition or sources upstream of a single observation. The estimates are sharpened by using many hydrocarbons with different reactivities and accounting for their variability and measurement uncertainty. A novel technique is used to construct prior background distributions of many species, described by variation of a single parameter . This exploits the high correlation of species, related by the first principal component of many NMHC samples. The Bayesian method obtains posterior estimates of , K and following each air mass. Median values are typically between 0.5 and 2.0 × 106 molecules cm−3, but are elevated to between 2.5 and 3.5 × 106 molecules cm−3, in low-level pollution. A comparison of estimates from absolute NMHC concentrations and NMHC ratios assuming zero background (the “photochemical clock” method) shows similar distributions but reveals systematic high bias in the estimates from ratios. Estimates of K are ∼0.1 day−1 but show more sensitivity to the prior distribution assumed

Ozone-CO Correlations Determined by the TES Satellite Instrument in Continental Outflow Regions

Collocated measurements of tropospheric ozone (O3) and carbon monoxide (CO) from the Tropospheric Emission Spectrometer (TES) aboard the EOS Aura satellite provide information on O3-CO correlations to test our understanding of global anthropogenic influence on O3. We examine the global distribution of TES O3-CO correlations in the middle troposphere (618 hPa) for July 2005 and compare to correlations generated with the GEOS-Chem chemical transport model and with ICARTT aircraft observations over the eastern United States (July 2004). The TES data show significant O3-CO correlations downwind of polluted continents, with dO3/dCO enhancement ratios in the range 0.4–1.0 mol mol−1 and consistent with ICARTT data. The GEOS-Chem model reproduces the O3-CO enhancement ratios observed in continental outflow, but model correlations are stronger and more extensive. We show that the discrepancy can be explained by spectral measurement errors in the TES data. These errors will decrease in future data releases, which should enable TES to provide better information on O3-CO correlations.Earth and Planetary SciencesEngineering and Applied Science