Long-Range Pollution Transport: Trans-Atlantic Mechanisms and Lagrangian Modeling Methods

Over the past several decades, it has become apparent that anthropogenic activities have resulted in the large-scale enhancement of the levels of many trace gases throughout the troposphere. More recently, attention has been given to the transport pathway taken by these emissions as they are dispersed throughout the atmosphere. The transport pathway determines the physical characteristics of emissions plumes and therefore plays an important role in the chemical transformations that can occur downwind of source regions. For example, the production of ozone (O3) is strongly dependent upon the transport its precursors undergo. O3 can initially be formed within air masses while still over polluted source regions. These polluted air masses can experience continued O3 production or O3 destruction downwind, depending on the air mass's chemical and transport characteristics. At present, however, there are a number of uncertainties in the relationships between transport and O3 production in the North Atlantic lower free troposphere. The first phase of the study presented here used measurements made at the Pico Mountain observatory and model simulations to determine transport pathways for US emissions to the observatory. The Pico Mountain observatory was established in the summer of 2001 in order to address the need to understand the relationships between transport and O3 production. Measurements from the observatory were analyzed in conjunction with model simulations from the Lagrangian particle dispersion model (LPDM), FLEX-PART, in order to determine the transport pathway for events observed at the Pico Mountain observatory during July 2003. A total of 16 events were observed, 4 of which were analyzed in detail. The transport time for these 16 events varied from 4.5 to 7 days, while the transport altitudes over the ocean ranged from 2-8 km, but were typically less than 3 km. In three of the case studies, eastward advection and transport in a weak warm conveyor belt (WCB) airflow was responsible for the export of North American emissions into the FT, while transport in the FT was governed by easterly winds driven by the Azores/Bermuda High (ABH) and transient northerly lows. In the fourth case study, North American emissions were lofted to 6-8 km in a WCB before being entrained in the same cyclone's dry airstream and transported down to the observatory. The results of this study show that the lower marine FT may provide an important transport environment where O3 production may continue, in contrast to transport in the marine boundary layer, where O3 destruction is believed to dominate. The second phase of the study presented here focused on improving the analysis methods that are available with LPDMs. While LPDMs are popular and useful for the analysis of atmospheric trace gas measurements, identifying the transport pathway of emissions from their source to a receptor (the Pico Mountain observatory in our case) using the standard gridded model output, particularly during complex meteorological scenarios can be difficult can be difficult or impossible. The transport study in phase 1 was limited to only 1 month out of more than 3 years of available data and included only 4 case studies out of the 16 events specifically due to this confounding factor. The second phase of this study addressed this difficulty by presenting a method to clearly and easily identify the pathway taken by only those emissions that arrive at a receptor at a particular time, by combining the standard gridded output from forward (i.e., concentrations) and backward (i.e., residence time) LPDM simulations, greatly simplifying similar analyses. The ability of the method to successfully determine the source-to-receptor pathway, restoring this Lagrangian information that is lost when the data are gridded, is proven by comparing the pathway determined from this method with the particle trajectories from both the forward and backward models. A sample analysis is also presented, demonstrating that this method is more accurate and easier to use than existing methods using standard LPDM products. Finally, we discuss potential future work that would be possible by combining the backward LPDM simulation with gridded data from other sources (e.g., chemical transport models) to obtain a Lagrangian sampling of the air that will eventually arrive at a receptor.

The role of ocean transport in the uptake of anthropogenic CO2

We compare modeled oceanic carbon uptake in response to pulse CO2 emissions using a suite of global ocean models and Earth system models. In response to a CO2 pulse emission of 590 Pg C (corresponding to an instantaneous doubling of atmospheric CO2 from 278 to 556 ppm), the fraction of CO2 emitted that is absorbed by the ocean is: 37±8%, 56±10%, and 81±4% (model mean ±2σ ) in year 30, 100, and 1000 after the emission pulse, respectively. Modeled oceanic uptake of pulse CO2 on timescales from decades to about a century is strongly correlated with simulated present-day uptake of chlorofluorocarbons (CFCs) and CO2 across all models, while the amount of pulse CO2 absorbed by the ocean from a century to a millennium is strongly correlated with modeled radiocarbon in the deep Southern and Pacific Ocean. However, restricting the analysis to models that are capable of reproducing observations within uncertainty, the correlation is generally much weaker. The rates of surface-to-deep ocean transport are determined for individual models from the instantaneous doubling CO2 simulations, and they are used to calculate oceanic CO2 uptake in response to pulse CO2 emissions of different sizes pulses of 1000 and 5000 Pg C. These results are compared with simulated oceanic uptake of CO2 by a number of models simulations with the coupling of climate-ocean carbon cycle and without it. This comparison demonstrates that the impact of different ocean transport rates across models on oceanic uptake of anthropogenic CO2 is of similar magnitude as that of climate-carbon cycle feedbacks in a single model, emphasizing the important role of ocean transport in the uptake of anthropogenic CO2.

A High-Resolution Record of Atmospheric Dust Composition and Variability Since Ad 1650 from a Mount Everest Ice Core

A Mount Everest ice core analyzed at high resolution for major and trace elements (Sr, Cs, Ba, La, Ce, Pr, Nd, Sm, Eu, Tb, Dy, Ho, Er, Tm, Yb, Lu, Bi, U, Tl, Al, S, Ca, Ti, V, Cr, Mn, Fe, Co) and spanning the period A. D. 1650- 2002 is used to investigate the sources of and variations in atmospheric dust through time. The chemical composition of dust varies seasonally, and peak dust concentrations occur during the winter-spring months. Significant correlations between the Everest dust record and dust observations at stations suggest that the Everest record is representative of regional variations in atmospheric dust loading. Back-trajectory analysis in addition to a significant correlation of Everest dust concentrations and the Total Ozone Mapping Spectrometer (TOMS) aerosol index indicates that the dominant winter sources of dust are the Arabian Peninsula, Thar Desert, and northern Sahara. Factors that contribute to dust generation at the surface include soil moisture and temperature, and the long-range transport of dust aerosols appears to be sensitive to the strength of 500-mb zonal winds. There are periods of high dust concentration throughout the 350-yr Mount Everest dust record; however, there is an increase in these periods since the early 1800s. The record was examined for recent increases in dust emissions associated with anthropogenic activities, but no recent dust variations can be conclusively attributed to anthropogenic inputs of dust.

Atmospheric Soluble Dust Records from a Tibetan Ice Core: Possible Climate Proxies and Teleconnection With the Pacific Decadal Oscillation

In autumn 2005, a joint expedition between the University of Maine and the Institute of Tibetan Plateau Research recovered three ice cores from Guoqu Glacier (33 degrees 34'37.80 '' N, 91 degrees 10'35.3 '' E, 5720 m above sea level) on the northern side of Mt. Geladaindong, central Tibetan Plateau. Isotopes ( delta(18)O), major soluble ions (Na(+), K(+), Mg(2+), Ca(2+), Cl(-), NO(3)(-), SO(4)(2-)), and radionuclide (beta-activity) measurements from one of the cores revealed a 70-year record (1935-2005). Statistical analysis of major ion time series suggests that atmospheric soluble dust species dominate the chemical signature and that background dust levels conceal marine ion species deposition. The soluble dust time series have interspecies relations and common structure (empirical orthogonal function (EOF) 1), suggesting a similar soluble dust source or transport route. Annual and seasonal correlations between the EOF 1 time series and National Centers for Environmental Prediction/National Center for Atmospheric Research reanalysis climate variables (1948-2004) suggest that the Mt. Geladaindong ice core record provides a proxy for local and regional surface pressure. An approximately threefold decrease of soluble dust concentrations in the middle to late 1970s, accompanied by regional increases in pressure and temperature and decreases in wind velocity, coincides with the major 1976-1977 shift of the Pacific Decadal Oscillation (PDO) from a negative to a positive state. This is the first ice core evidence of a potential teleconnection between central Asian atmospheric soluble dust loading and the PDO. Analysis of temporally longer ice cores from Mt. Geladaindong may enhance understanding of the relationship between the PDO and central Asian atmospheric circulation and subsequent atmospheric soluble dust loading.

A climatological analysis of high-precipitation events in Dronning Maud Land, Antarctica, and associated large-scale atmospheric conditions

The link between high precipitation in Dronning Maud Land (DML), Antarctica, and the large-scale atmospheric circulation is investigated using ERA-Interim data for 1979–2009. High-precipitation events are analyzed at Halvfarryggen situated in the coastal region of DML and at Kohnen Station located in its interior. This study further includes a comprehensive comparison of high precipitation in ERA-Interim with precipitation data from the Antarctic Mesoscale Prediction System (AMPS) and snow accumulation measurements from automatic weather stations (AWSs), with the limitations of such a comparison being discussed. The ERA-Interim and AMPS precipitation data agree very well. However, the correspondence between high precipitation in ERA-Interim and high snow accumulation at the AWSs is relatively weak. High-precipitation events at both Halvfarryggen and Kohnen are typically associated with amplified upper level waves. This large-scale atmospheric flow pattern is preceded by the downstream development of a Rossby wave train from the eastern South Pacific several days before the precipitation event. At the surface, a cyclone located over the Weddell Sea is the main synoptic ingredient for high precipitation both at Halvfarryggen and at Kohnen. A blocking anticyclone downstream is not a requirement for high precipitation per se, but a larger share of blocking occurrences during the highest-precipitation days in DML suggests that these blocks strengthen the vertically integrated water vapor transport (IVT) into DML. A strong link between high precipitation and the IVT perpendicular to the local orography suggests that IVT could be used as a “proxy” for high precipitation, in particular over DML's interior.

Benchmark problems for reactive transport modeling of the generation and attenuation of acid rock drainage

Acid rock drainage (ARD) is a problem of international relevance with substantial environmental and economic implications. Reactive transport modeling has proven a powerful tool for the process-based assessment of metal release and attenuation at ARD sites. Although a variety of models has been used to investigate ARD, a systematic model intercomparison has not been conducted to date. This contribution presents such a model intercomparison involving three synthetic benchmark problems designed to evaluate model results for the most relevant processes at ARD sites. The first benchmark (ARD-B1) focuses on the oxidation of sulfide minerals in an unsaturated tailing impoundment, affected by the ingress of atmospheric oxygen. ARD-B2 extends the first problem to include pH buffering by primary mineral dissolution and secondary mineral precipitation. The third problem (ARD-B3) in addition considers the kinetic and pH-dependent dissolution of silicate minerals under low pH conditions. The set of benchmarks was solved by four reactive transport codes, namely CrunchFlow, Flotran, HP1, and MIN3P. The results comparison focused on spatial profiles of dissolved concentrations, pH and pE, pore gas composition, and mineral assemblages. In addition, results of transient profiles for selected elements and cumulative mass loadings were considered in the intercomparison. Despite substantial differences in model formulations, very good agreement was obtained between the various codes. Residual deviations between the results are analyzed and discussed in terms of their implications for capturing system evolution and long-term mass loading predictions.

Gas transport in firn: multiple-tracer characterisation and model intercomparison for NEEM, Northern Greenland

Air was sampled from the porous firn layer at the NEEM site in Northern Greenland. We use an ensemble of ten reference tracers of known atmospheric history to characterise the transport properties of the site. By analysing uncertainties in both data and the reference gas atmospheric histories, we can objectively assign weights to each of the gases used for the depth-diffusivity reconstruction. We define an objective root mean square criterion that is minimised in the model tuning procedure. Each tracer constrains the firn profile differently through its unique atmospheric history and free air diffusivity, making our multiple-tracer characterisation method a clear improvement over the commonly used single-tracer tuning. Six firn air transport models are tuned to the NEEM site; all models successfully reproduce the data within a 1σ Gaussian distribution. A comparison between two replicate boreholes drilled 64 m apart shows differences in measured mixing ratio profiles that exceed the experimental error. We find evidence that diffusivity does not vanish completely in the lock-in zone, as is commonly assumed. The ice age- gas age difference (1 age) at the firn-ice transition is calculated to be 182+3−9 yr. We further present the first intercomparison study of firn air models, where we introduce diagnostic scenarios designed to probe specific aspects of the model physics. Our results show that there are major differences in the way the models handle advective transport. Furthermore, diffusive fractionation of isotopes in the firn is poorly constrained by the models, which has consequences for attempts to reconstruct the isotopic composition of trace gases back in time using firn air and ice core records.

Evaluating atmospheric methane inversion model results for Pallas, northern Finland

A state-of-the-art inverse model, CarbonTracker Data Assimilation Shell (CTDAS), was used to optimize estimates of methane (CH4) surface fluxes using atmospheric observations of CH4 as a constraint. The model consists of the latest version of the TM5 atmospheric chemistry-transport model and an ensemble Kalman filter based data assimilation system. The model was constrained by atmospheric methane surface concentrations, obtained from the World Data Centre for Greenhouse Gases (WDCGG). Prior methane emissions were specified for five sources: biosphere, anthropogenic, fire, termites and ocean, of which bio-sphere and anthropogenic emissions were optimized. Atmospheric CH 4 mole fractions for 2007 from northern Finland calculated from prior and optimized emissions were compared with observations. It was found that the root mean squared errors of the posterior esti - mates were more than halved. Furthermore, inclusion of NOAA observations of CH 4 from weekly discrete air samples collected at Pallas improved agreement between posterior CH 4 mole fraction estimates and continuous observations, and resulted in reducing optimized biosphere emissions and their uncertainties in northern Finland.

Impact of oceanic circulation changes on atmospheric δ¹³ CO₂

δ¹³ CO₂ measured in Antarctic ice cores provides constraints on oceanic and terrestrial carbon cycle processes linked with millennial-scale changes in atmospheric CO₂. However, the interpretation of δ¹³ CO₂ is not straight-forward. Using carbon isotope-enabled versions of the LOVECLIM and Bern3D models, we perform a set of sensitivity experiments in which the formation rates of North Atlantic Deep Water (NADW), North Pacific Deep Water (NPDW), Antarctic Bottom Water (AABW), and Antarctic Intermediate Water (AAIW) are varied. We study the impact of these circulation changes on atmospheric δ¹³ CO₂ as well as on the oceanic δ¹³ CO₂ distribution. In general, we find that the formation rates of AABW, NADW, NPDW, and AAIW are negatively correlated with changes in δ¹³ CO₂: namely, strong oceanic ventilation decreases atmospheric δ¹³ CO₂. However, since large-scale oceanic circulation reorganizations also impact nutrient utilization and the Earth’s climate, the relationship between atmospheric δ¹³ CO₂ levels and ocean ventilation rate is not unequivocal. In both models atmospheric δ¹³ CO₂ is very sensitive to changes in AABW formation rates: increased AABW formation enhances the transport of low δ¹³ CO₂ waters to the surface and decreases atmospheric δ¹³ CO₂. By contrast, the impact of NADW changes on atmospheric δ¹³ CO₂ is less robust and might be model dependent. This results from complex interplay between global climate, carbon cycle, and the formation rate of NADW, a water body characterized by relatively high δ¹³ CO₂.

Surface in situ measurements of atmospheric CO, CH4, and CO2 combined with selected meteorological measurements at the University of Wollongong, Australia

Wollongong, Australia is an urban site at the intersection of anthropogenic, biomass burning, biogenic and marine sources of atmospheric trace gases. The location offers a valuable opportunity to study drivers of atmospheric composition in the Southern Hemisphere. Here, a record of surface carbon monoxide (CO), methane (CH4) and carbon dioxide (CO2) was measured with an in situ Fourier transform infrared trace gas analyser between April 2011 and August 2014. Clean air was found to arrive at Wollongong in approximately 10% of air masses. Biomass burning influence was evident in the average annual cycle of clean air CO during austral spring. A significant negative short-term trend was found in clean air CO (-1.5 nmol/mol/a), driven by a reduction in northern Australian biomass burning. Significant short-term positive trends in clean air CH4 (5.4 nmol/mol/a) and CO2 (1.9 ?mol/mol/a) were consistent with the long-term global average trends. Polluted Wollongong air was investigated using wind-direction/wind-speed clustering, which revealed major influence from local urban and industrial sources from the south. High values of CH4, with anthropogenic DCH4/DCO2 enhancement ratio signatures, originated from the northwest, in the direction of local coal mining. A pollution climatology was developed for the region using back trajectory analysis and DO3/DCO enhancement ratios. Ozone production environments in austral spring and summer were associated with anticyclonic meteorology on the east coast of Australia, while ozone depletion environments in autumn and winter were associated with continental transport, or fast moving trajectories from southern latitudes. This implies the need to consider meteorological conditions when developing policies for controlling air quality.

Modelling Saharan dust transport into the Mediterranean basin with CMAQ

The need for a better quantification of the influence of Saharan dust transport processes on the air quality modelling in the Mediterranean basin led to the formulation of a dust emission module (DEM) integrated into the Air Quality Risk Assessment System for the Iberian Peninsula (SERCA). This paper is focused on the formulation of DEM based on the GOCART aerosol model, along with its integration and execution into the air quality model. It also addresses the testing of the module and its evaluation by contrasting results against satellite products such as MODIS and CALIPSO and ground-level observations of aerosol optical thickness (AOT) and concentration levels of PM10 for different periods in July 2007. DEM was found capable of reproducing the spatial (horizontal and vertical) and temporal profiles of Saharan dust outbreaks into the Mediterranean basin and the Atlantic coast of Africa. Moreover, it was observed that its combination with CMAQ increased the correlation degree between observed and modelled PM10 concentrations at the selected monitoring locations. DEM also enhanced CMAQ capabilities to reproduce observed AOT, although significant underestimations remain. The implementation of CMAQ + DEM succeeded in capturing Saharan dust transport into the Iberian Peninsula, with contributions up to 25 and 14 μg m−3 in 1 h and 24 h average PM10 respectively. The general improvement of total PM10 predictions in Spain are however moderate. The analysis of model performance for the main PM components points out that remaining PM10 underestimation is due to dust local sources missing in the inventories and misrepresentation of organic aerosol processes, which constitutes the main areas for future improvement of CMAQ capabilities to simulate particulate matter within SERCA.

A methodology to compute emission projections from road transport (EmiTRANS)

Atmospheric emissions from road transport have increased all around the world during the last decades more rapidly than from other pollution sources. For instance, they contribute to more than 25% of total CO, CO2, NOx, and fine particle emissions in most of the European countries. This situation shows the importance of road transport when complying with emission ceilings and air quality standards applied to these pollutants. This paper presents a modelling system to perform atmospheric emission projections (simultaneously both air quality pollutants and greenhouse gases) from road transport including the development of a tailored software tool (EmiTRANS) as a planning tool. The methodology has been developed with two purposes: 1) to obtain outputs used as inputs to the COPERT4 software to calculate emission projections and 2) to summarize outputs for policy making evaluating the effect of emission abatement measures for a vehicle fleet. This methodology has been applied to the calculation of emission projections in Spain up to 2020 under several scenarios, including a sensitivity analysis useful for a better interpretation and confidence building on the results. This case study demonstrates the EmiTRANS applicability to a country, and points out the need for combining both technical and non-technical measures (such as behavioural changes or demand management) to reduce emissions, indirectly improving air quality and contributing to mitigate climate change.