Earth Systems

In this thesis, I study polycyclic aromatic hydrocarbons (PAHs) and perfluorocarboxylic acids (PFCAs). PAHs are by-products of burning and therefore have important anthropogenic sources in the combustion of fuels, biomass, etc. PFCAs and their atmospheric precursors are used in making firefighting foams, non-stick coatings, and other surfactant applications.

I quantitatively examine the relative importance of uncertainty in emissions and physicochemical properties (including reaction rate constants) to Northern Hemisphere (NH) and Arctic PAH concentrations. NH average concentrations are more sensitive to uncertainty in the atmospheric lifetime than to emissions rate. The largest uncertainty reductions would come from precise experimental determination of PHE, PYR and BaP rate constants for the reaction with OH.

I calculate long-chain PFCA formation theoretical maximum yields for the degradation of precursor species at a representative sample of atmospheric conditions from a three dimensional chemical transport model, finding that atmospheric conditions farther from pollution sources have both higher capacities to form long chain PFCAs and higher uncertainties in those capacities.

I present results from newly developed simulations of atmospheric PFCA formation and fate using the chemical transport model GEOS-Chem, simulating the degradation of fluorotelomer precursors, as well as deposition and transport of the precursors, intermediates and end-products of the PFCA formation chemistry. I compare the model results to remote deposition measurements and find that it reproduces Arctic deposition of PFOA effectively. Given the most recent precursor emission inventory, the atmospheric indirect source of PFOA and PFNA is 10-45 t/yr globally and 0.2-0.7 t/yr to the Arctic.

The climate response of precipitation to the effects of anthropogenic aerosols is a critical while not yet fully understood aspect in climate science. Results of selected models that participated the Coupled Model Intercomparison Project Phase 5 and the data from the Twentieth Century Reanalysis Project suggest that, throughout the tropics and also in the extratropical Northern Hemisphere, aerosols have largely dominated the distribution of precipitation changes in reference to the preindustrial era in the second half of the last century. Aerosol-induced cooling has offset some of the warming caused by the greenhouse gases from the tropics to the Arctic and thus formed the gradients of surface temperature anomaly that enable the revealed precipitation change patterns to occur. Improved representation of aerosol-cloud interaction has been demonstrated as the key factor for models to reproduce consistent distributions of past precipitation change with the reanalysis data.

© 2015 the author

We report ground-based atmospheric measurements and emission estimates for the halons H-1211 (CBrClF2), H-1301 (CBrF3), and H-2402 (CBrF2CBrF2) from the AGAGE (Advanced Global Atmospheric Gases Experiment) and the National Oceanic and Atmospheric Administration global networks. We also include results from archived air samples in canisters and from polar firn in both hemispheres, thereby deriving an atmospheric record of nearly nine decades (1930s to present). All three halons were absent from the atmosphere until ∼1970, when their atmospheric burdens started to increase rapidly. In recent years H-1211 and H-2402 mole fractions have been declining, but H-1301 has continued to grow. High-frequency observations show continuing emissions of H-1211 and H-1301 near most AGAGE sites. For H-2402 the only emissions detected were derived from the region surrounding the Sea of Japan/East Sea. Based on our observations, we derive global emissions using two different inversion approaches. Emissions for H-1211 declined from a peak of 11 kt yr−1 (late 1990s) to 3.9 kt yr−1 at the end of our record (mean of 2013–2015), for H-1301 from 5.4 kt yr−1 (late 1980s) to 1.6 kt yr−1, and for H-2402 from 1.8 kt yr−1 (late 1980s) to 0.38 kt yr−1. Yearly summed halon emissions have decreased substantially; nevertheless, since 2000 they have accounted for ∼30% of the emissions of all major anthropogenic ozone depletion substances, when weighted by ozone depletion potentials.

Rivers discharge 28 ± 13 Mmol yr−1 of mercury (Hg) to ocean margins, an amount comparable to atmospheric deposition to the global oceans. Most of the Hg discharged by rivers is sequestered by burial of benthic sediment in estuaries or the coastal zone, but some is evaded to the atmosphere and some is exported to the open ocean. We investigate the fate of riverine Hg by developing a new global 3-D simulation for Hg in the Massachusetts Institute of Technology ocean general circulation model. The model includes plankton dynamics and carbon respiration (DARWIN project model) coupled to inorganic Hg chemistry. Results are consistent with observed spatial patterns and magnitudes of surface ocean Hg concentrations. We use observational constraints on seawater Hg concentrations and evasion to infer that most Hg from rivers is sorbed to refractory organic carbon and preferentially buried. Only 6% of Hg discharged by rivers (1.8 Mmol yr−1) is transported to the open ocean on a global basis. This fraction varies from a low of 2.6% in East Asia due to the barrier imposed by the Korean Peninsula and Japanese archipelago, up to 25% in eastern North America facilitated by the Gulf Stream. In the Arctic Ocean, low tributary particle loads and efficient degradation of particulate organic carbon by deltaic microbial communities favor a more labile riverine Hg pool. Evasion of Hg to the Arctic atmosphere is indirectly enhanced by heat transport during spring freshet that accelerates sea ice melt and ice rafting. Discharges of 0.23 Mmol Hg yr−1 from Arctic rivers can explain the observed summer maximum in the Arctic atmosphere, and this magnitude of releases is consistent with recent observations. Our work indicates that rivers are major contributors to Hg loads in the Arctic Ocean.

© 2015 American Chemical Society

We present a numerical model of the ocean that couples a three-stream radiative transfer component with a marine biogeochemical–ecosystem component in a dynamic three-dimensional physical framework. The radiative transfer component resolves the penetration of spectral irradiance as it is absorbed and scattered within the water column. We explicitly include the effect of several optically important water constituents (different phytoplankton functional types; detrital particles; and coloured dissolved organic matter, CDOM). The model is evaluated against in situ-observed and satellite-derived products. In particular we compare to concurrently measured biogeochemical, ecosystem, and optical data along a meridional transect of the Atlantic Ocean. The simulation captures the patterns and magnitudes of these data, and estimates surface upwelling irradiance analogous to that observed by ocean colour satellite instruments. We find that incorporating the different optically important constituents explicitly and including spectral irradiance was crucial to capture the variability in the depth of the subsurface chlorophyll a (Chl a) maximum. We conduct a series of sensitivity experiments to demonstrate, globally, the relative importance of each of the water constituents, as well as the crucial feedbacks between the light field, the relative fitness of phytoplankton types, and the biogeochemistry of the ocean. CDOM has proportionally more importance at attenuating light at short wavelengths and in more productive waters, phytoplankton absorption is relatively more important at the subsurface Chl a maximum, and water molecules have the greatest contribution when concentrations of other constituents are low, such as in the oligotrophic gyres. Scattering had less effect on attenuation, but since it is important for the amount and type of upwelling irradiance, it is crucial for setting sea surface reflectance. Strikingly, sensitivity experiments in which absorption by any of the optical constituents was increased led to a decrease in the size of the oligotrophic regions of the subtropical gyres: lateral nutrient supplies were enhanced as a result of decreasing high-latitude productivity. This new model that captures bio-optical feedbacks will be important for improving our understanding of the role of light and optical constituents on ocean biogeochemistry, especially in a changing environment. Further, resolving surface upwelling irradiance will make it easier to connect to satellite-derived products in the future.

© 2015 the authors

Land management for carbon sequestration offers an opportunity to avoid about 0.5°C of warming if landowners have full economic incentives to participate in a global greenhouse gas mitigation policy. In an energy-only policy aimed at about 550 ppm CO2- eq stabilization the additional 0.5°C of avoided warming brings the world close to staying below the 2°C above preindustrial target. While greater incentives for mitigation from energy and land would be needed to actually meet the 2°C target even holding temperature to 2.2°C or so would be a great improvement over the path we are on now. Even with success at the 21st meeting of the Conference of the Parties (COP21) in Paris in 2015, we would still be heading toward a global mean surface temperature increase is in the range of 1.9 to 2.6°C (central estimate 2.2°C) by mid-century relative to the pre-industrial level (1860–1880 mean), and 3.1 to 5.2°C (central estimate 3.7°C) by 2100. This preliminary assessment awaits a full interpretation of the implications of agreements at COP21 with greater clarity on just what countries pledges mean and how they will be implemented. But this first look estimates that the COP21 pledges would shave about 0.2°C more from warming compared with previous international commitments. From that perspective, a half-degree of avoided warming from land carbon sequestration and avoided deforestation, if we could achieve it, is significant. The relatively small contribution from the COP21 agreement says less about the significance of the commitments—many countries stepped up and offered important goals that will mean changes in their energy systems—but rather more about the challenge of weaning the global economy from its dependence on greenhouse gas emitting activities. It also highlights that no source (or potential sink) can be ignored.

© 2016 SciTechnol

Atmosphere–surface exchange of mercury, although a critical component of its global cycle, is currently poorly constrained. Here we use the GEOS-Chem chemical transport model to interpret atmospheric Hg0 (gaseous elemental mercury) data collected during the 2013 summer Nitrogen, Oxidants, Mercury and Aerosol Distributions, Sources and Sinks (NOMADSS) aircraft campaign as well as ground- and ship-based observations in terms of their constraints on the atmosphere–surface exchange of Hg0 over eastern North America. Model–observation comparison suggests that the Northwest Atlantic may be a net source of Hg0, with high evasion fluxes in summer (our best sensitivity simulation shows an average oceanic Hg0 flux of 3.3 ng m-2 h-1 over the Northwest Atlantic), while the terrestrial ecosystem in the summer of the eastern United States is likely a net sink of Hg0 (our best sensitivity simulation shows an average terrestrial Hg0 flux of -0.6 ng m-2 h-1 over the eastern United States). The inferred high Hg0 fluxes from the Northwest Atlantic may result from high wet deposition fluxes of oxidized Hg, which are in turn related to high precipitation rates in this region. We also find that increasing simulated terrestrial fluxes of Hg0 in spring compared to other seasons can better reproduce observed seasonal variability of Hg0 concentration at ground-based sites in eastern North America.

Artificial fertilisation of the ocean has been proposed as a possible geoengineering method for removing carbon dioxide from the atmosphere. The associated increase in marine primary productivity may lead to an increase in emissions of dimethyl sulphide (DMS), the primary source of sulphate aerosol over remote ocean regions, potentially causing direct and cloud-related indirect aerosol effects on climate. This pathway from ocean fertilisation to aerosol induced cooling of the climate may provide a basis for solar radiation management (SRM) geoengineering. In this study, we investigate the transient climate impacts of two emissions scenarios: an RCP4.5 (Representative Concentration Pathway 4.5) control; and an idealised scenario, based on RCP4.5, in which DMS emissions are substantially enhanced over ocean areas. We use mini-ensembles of a coupled atmosphere-ocean configuration of CESM1(CAM5) (Community Earth System Model version 1, with the Community Atmosphere Model version 5). We find that the cooling effect associated with enhanced DMS emissions beneficially offsets greenhouse gas induced warming across most of the world. However, the rainfall response may adversely affect water resources, potentially impacting human livelihoods. These results demonstrate that changes in marine phytoplankton activity may lead to a mixture of positive and negative impacts on the climate.

© 2015 the authors

Long-term response of the climate system to anthropogenic forcing was investigated with the MIT Earth System Model of intermediate complexity version 2.2 (MESM2.2). The MESM2.2 consists of a 2D (zonally averaged) atmospheric model coupled to an anomaly diffusing ocean model. Climate sensitivity of the MESM can be varied using a cloud adjustment technique and rate of oceanic heat uptake can be varied by changing effective diffusion coefficient. An ensemble of four hundred simulations was carried out for the period 1860-2005 using historical forcing. Values of climate sensitivity, rate of ocean heat uptake, and the strength of the aerosol forcing were drawn from the Libardoni and Forest (2013) distribution presented in the IPCC AR5. A 400-member ensemble was carried out for each of four different RCP scenarios from the year 2006 to the year 2500. By the end of the 21st century (2081-2100), the ensemble mean of surface air temperature increases, relative to 1986-2005 period, by 1.2, 1.8, 2.2 and 3.3oC for RCP26, RCP4.5, RCP6.0 and RCP8.5, respectively. Corresponding numbers for the ensemble of the CMPI5 models are 1.0, 1.8, 2.2 and 3.7oC. In spite of the forcing being fixed beyond year 2150 for RCP4.5 and RCP6.0 and beyond 2250 for RCP8.5, surface air temperature keeps rising until the end of 25th century under these scenarios. The upper bound of the 90% probability interval increases significantly more than the mean. For the RCP4.5 scenario, the mean value of possible SAT change increases by 1.6oC from the end of the 21st century to the end of the 25th century, while the value of the 95th percentile increases by 3.2oC. Corresponding numbers for RCP6.0 and RCP8.5 are 3.6 and 10.2oC for the medians and 7.0 and 14.5oC for the 95th percentiles, respectively. Such changes in the shape of probability distributions with time indicate an increase in the probability that surface warming will exceed a given value. For example, the probability of exceeding 3oC warming under the RCP4.5 scenario increases from 2.5% at the end of 21st century to 32% and 50% at the end of 23rd and 25th centuries, respectively. For the RCP2.6 scenario, in which radiative forcing peaks in the year 2070 before decreasing back to the 1990s level by the year 2300, the ensemble mean surface air temperature is still about 0.5oC above present at the end of the simulation. Obtained results show that in spite of large differences in radiative forcing between different RCP scenarios, uncertainties in the climate system characteristics defining climate system response make a significant contribution into overall uncertainty in possible climate change during the next few centuries. Comparison with simulations carried under SRES scenarios also will be presented

High frequency, in situ observations from 11 globally distributed sites for the period 1994–2014 and archived air measurements dating from 1978 onward have been used to determine the global growth rate of 1,1-difluoroethane (HFC-152a, CH3CHF2). These observations have been combined with a range of atmospheric transport models to derive global emission estimates in a top-down approach. HFC-152a is a greenhouse gas with a short atmospheric lifetime of about 1.5 years. Since it does not contain chlorine or bromine, HFC-152a makes no direct contribution to the destruction of stratospheric ozone and is therefore used as a substitute for the ozone depleting chlorofluorocarbons (CFCs) and hydrochlorofluorocarbons (HCFCs). The concentration of HFC-152a has grown substantially since the first direct measurements in 1994, reaching a maximum annual global growth rate of 0.84"¯±"¯0.05"¯ppt"¯yr−1 in 2006, implying a substantial increase in emissions up to 2006. However, since 2007, the annual rate of growth has slowed to 0.38"¯±"¯0.04"¯ppt"¯yr−1 in 2010 with a further decline to an annual average rate of growth in 2013–2014 of −0.06"¯±"¯0.05"¯ppt"¯yr−1. The annual average Northern Hemisphere (NH) mole fraction in 1994 was 1.2"¯ppt rising to an annual average mole fraction of 10.1"¯ppt in 2014. Average annual mole fractions in the Southern Hemisphere (SH) in 1998 and 2014 were 0.84 and 4.5"¯ppt, respectively. We estimate global emissions of HFC-152a have risen from 7.3"¯±"¯5.6"¯Gg"¯yr−1 in 1994 to a maximum of 54.4"¯±"¯17.1"¯Gg"¯yr−1 in 2011, declining to 52.5"¯±"¯20.1"¯Gg"¯yr−1 in 2014 or 7.2"¯±"¯2.8"¯Tg-CO2"¯eq"¯yr−1. Analysis of mole fraction enhancements above regional background atmospheric levels suggests substantial emissions from North America, Asia, and Europe. Global HFC emissions (so called “bottom up” emissions) reported by the United Nations Framework Convention on Climate Change (UNFCCC) are based on cumulative national emission data reported to the UNFCCC, which in turn are based on national consumption data. There appears to be a significant underestimate ( > "¯20"¯Gg) of “bottom-up” reported emissions of HFC-152a, possibly arising from largely underestimated USA emissions and undeclared Asian emissions.

© 2016 the authors

Pages

Subscribe to Earth Systems