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.

Upper ocean vertical supply: A neglected primary factor controlling the distribution of neodymium concentrations of open ocean surface waters?

Neodymium (Nd) isotopes are an important geochemical tool to trace the present and past water mass mixing as well as continental inputs. The distribution of Nd concentrations in open ocean surface waters (0�100 m) is generally assumed to be controlled by lateral mixing of Nd from coastal surface currents and by removal through reversible particle scavenging. However, using 228Ra activity as an indicator of coastal water mass influence, surface water Nd concentration data available on key oceanic transects as a whole do not support the above scenario. From a global compilation of available data, we find that more stratified regions are generally associated with low surface Nd concentrations. This implies that upper ocean vertical supply may be an as yet neglected primary factor in determining the basin-scale variations of surface water Nd concentrations. Similar to the mechanism of nutrients supply, it is likely that stratification inhibits vertical supply of Nd from the subsurface thermocline waters and thus the magnitude of Nd flux to the surface layer. Consistently, the estimated required input flux of Nd to the surface layer to maintain the observed concentrations could be nearly two orders of magnitudes larger than riverine/dust flux, and also larger than the model-based estimation on shelf-derived coastal flux. In addition, preliminary results from modeling experiments reveal that the input from shallow boundary sources, riverine input, and release from dust are actually not the primary factors controlling Nd concentrations most notably in the Pacific and Southern Ocean surface waters.

Regulation of Phytoplankton Carbon to Chlorophyll Ratio By Light, Nutrients and Temperature in the Equatorial Pacific Ocean: A Basin-Scale Model

The complex effects of light, nutrients and temperature lead to a variable carbon to chlorophyll (C:Chl) ratio in phytoplankton cells. Using field data collected in the Equatorial Pacific, we derived a new dynamic model with a non-steady C:Chl ratio as a function of irradiance, nitrate, iron, and temperature. The dynamic model is implemented into a basin-scale ocean circulation-biogeochemistry model and tested in the Equatorial Pacific Ocean. The model reproduces well the general features of phytoplankton dynamics in this region. For instance, the simulated deep chlorophyll maximum (DCM) is much deeper in the western warm pool (similar to 100 m) than in the Eastern Equatorial Pacific (similar to 50 m). The model also shows the ability to reproduce chlorophyll, including not only the zonal, meridional and vertical variations, but also the interannual variability. This modeling study demonstrates that combination of nitrate and iron regulates the spatial and temporal variations in the phytoplankton C:Chl ratio in the Equatorial Pacific. Sensitivity simulations suggest that nitrate is mainly responsible for the high C:Chl ratio in the western warm pool while iron is responsible for the frontal features in the C:Chl ratio between the warm pool and the upwelling region. In addition, iron plays a dominant role in regulating the spatial and temporal variations of the C:Chl ratio in the Central and Eastern Equatorial Pacific. While temperature has a relatively small effect on the C:Chl ratio, light is primarily responsible for the vertical decrease of phytoplankton C:Chl ratio in the euphotic zone.

Dispersal Modeling of Fish Early Life Stages: Sensitivity with Application to Atlantic Cod in the Western Gulf of Maine

As an initial step in establishing mechanistic relationships between environmental variability and recruitment in Atlantic cod Gadhus morhua along the coast of the western Gulf of Maine, we assessed transport success of larvae from major spawning grounds to nursery areas with particle tracking using the unstructured grid model FVCOM (finite volume coastal ocean model). In coastal areas, dispersal of early planktonic life stages of fish and invertebrate species is highly dependent on the regional dynamics and its variability, which has to be captured by our models. With state-of-the-art forcing for the year 1995, we evaluate the sensitivity of particle dispersal to the timing and location of spawning, the spatial and temporal resolution of the model, and the vertical mixing scheme. A 3 d frequency for the release of particles is necessary to capture the effect of the circulation variability into an averaged dispersal pattern of the spawning season. The analysis of sensitivity to model setup showed that a higher resolution mesh, tidal forcing, and current variability do not change the general pattern of connectivity, but do tend to increase within-site retention. Our results indicate strong downstream connectivity among spawning grounds and higher chances for successful transport from spawning areas closer to the coast. The model run for January egg release indicates 1 to 19 % within-spawning ground retention of initial particles, which may be sufficient to sustain local populations. A systematic sensitivity analysis still needs to be conducted to determine the minimum mesh and forcing resolution that adequately resolves the complex dynamics of the western Gulf of Maine. Other sources of variability, i.e. large-scale upstream forcing and the biological environment, also need to be considered in future studies of the interannual variability in transport and survival of the early life stages of cod.

Modeling postglacial vegetation dynamics of temperate forests on the Olympic Peninsula (WA, USA) with special regard to snowpack

Past and future forest composition and distribution in temperate mountain ranges is strongly influenced by temperature and snowpack. We used LANDCLIM, a spatially explicit, dynamic vegetation model, to simulate forest dynamics for the last 16,000 years and compared the simulation results to pollen and macrofossil records at five sites on the Olympic Peninsula (Washington, USA). To address the hydrological effects of climate-driven variations in snowpack on simulated forest dynamics, we added a simple snow accumulation-and-melt module to the vegetation model and compared simulations with and without the module. LANDCLIM produced realistic present-day species composition with respect to elevation and precipitation gradients. Over the last 16,000 years, simulations driven by transient climate data from an atmosphere-ocean general circulation model (AOGCM) and by a chironomid-based temperature reconstruction captured Late-glacial to Late Holocene transitions in forest communities. Overall, the reconstruction-driven vegetation simulations matched observed vegetation changes better than the AOGCM-driven simulations. This study also indicates that forest composition is very sensitive to snowpack-mediated changes in soil moisture. Simulations without the snow module showed a strong effect of snowpack on key bioclimatic variables and species composition at higher elevations. A projected upward shift of the snow line and a decrease in snowpack might lead to drastic changes in mountain forests composition and even a shift to dry meadows due to insufficient moisture availability in shallow alpine soils.

Modeling the Effects of Size on Patch Dynamics of an Inert Tracer

Mesoscale iron enrichment experiments have revealed that additional iron affects the phytoplankton productivity and carbon cycle. However, the role of initial size of fertilized patch in determining the patch evolution is poorly quantified due to the limited observational capability and complex of physical processes. Using a three-dimensional ocean circulation model, we simulated different sizes of inert tracer patches that were only regulated by physical circulation and diffusion. Model results showed that during the first few days since release of inert tracer, the calculated dilution rate was found to be a linear function with time, which was sensitive to the initial patch size with steeper slope for smaller size patch. After the initial phase of rapid decay, the relationship between dilution rate and time became an exponential function, which was also size dependent. Therefore, larger initial size patches can usually last longer and ultimately affect biogeochemical processes much stronger than smaller patches.

A global seasonal surface ocean climatology of phytoplankton types based on CHEMTAX analysis of HPLC pigments

Much advancement has been made in recent years in field data assimilation, remote sensing and ecosystem modeling, yet our global view of phytoplankton biogeography beyond chlorophyll biomass is still a cursory taxonomic picture with vast areas of the open ocean requiring field validations. High performance liquid chromatography (HPLC) pigment data combined with inverse methods offer an advantage over many other phytoplankton quantification measures by way of providing an immediate perspective of the whole phytoplankton community in a sample as a function of chlorophyll biomass. Historically, such chemotaxonomic analysis has been conducted mainly at local spatial and temporal scales in the ocean. Here, we apply a widely tested inverse approach, CHEMTAX, to a global climatology of pigment observations from HPLC. This study marks the first systematic and objective global application of CHEMTAX, yielding a seasonal climatology comprised of ~1500 1°x1° global grid points of the major phytoplankton pigment types in the ocean characterizing cyanobacteria, haptophytes, chlorophytes, cryptophytes, dinoflagellates, and diatoms, with results validated against prior regional studies where possible. Key findings from this new global view of specific phytoplankton abundances from pigments are a) the large global proportion of marine haptophytes (comprising 32 ± 5% of total chlorophyll), whose biogeochemical functional roles are relatively unknown, and b) the contrasting spatial scales of complexity in global community structure that can be explained in part by regional oceanographic conditions. These publicly accessible results will guide future parameterizations of marine ecosystem models exploring the link between phytoplankton community structure and marine biogeochemical cycles.

Biomarker proxies on sediment cores in the central Arctic Ocean

Although the permanently to seasonally ice-covered Arctic Ocean is a unique and sensitive component in the Earth's climate system, the knowledge of its long-term climate history remains very limited due to the restricted number of pre-Quaternary sedimentary records. During Polarstern Expedition PS87/2014, we discovered multiple submarine landslides over a distance of >350 km along Lomonosov Ridge. Removal of younger sediments from steep headwalls has led to exhumation of Miocene to early Quaternary sediments close to the seafloor, allowing the retrieval of such old sediments with gravity cores. Multi-proxy biomarker analyses of these gravity cores reveal for the first time that the late Miocene central Arctic Ocean was relatively warm (4-7°C) and ice-free during summer, whereas sea ice occurred during spring and autumn/winter. A comparison of our proxy data with Miocene climate simulations seems to favour relatively high late Miocene atmospheric CO2 concentrations. These new findings from the Arctic region provide new benchmarks for groundtruthing global climate reconstructions and modeling.

Particle flux studies of sediment traps from the northern and equatorial Indian Ocean

The Asian monsoon system governs seasonality and fundamental environmental characteristics in the study area from which two distinct peculiarities are most notable: upwelling and convective mixing in the Arabian Sea and low surface salinity and stratification in the Bay of Bengal due to high riverine input and monsoonal precipitation. The respective oceanography sets the framework for nutrient availability and productivity. Upwelling ensures high nitrate concentration with temporal/spatial Si limitation; freshwater-induced stratification leads to reduced nitrogen input from the subsurface but Si enrichment in surface waters. Ultimately, both environments support high abundance of diatoms, which play a central role in the export of organic matter. It is speculated that, additional to eddy pumping, nitrogen fixation is a source of N in stratified waters and contributes to the low-d15N signal in sinking particles formed under riverine impact. Organic carbon fluxes are best correlated to opal but not to carbonate, which is explained by low foraminiferal carbonate fluxes within the river-impacted systems. This observation points to the necessity of differentiating between carbonate sources for carbon flux modeling. As evident from a compilation of previously published and new data on labile organic matter composition (amino acids and carbohydrates), organic matter fluxes are mainly driven by direct input from marine production, except the site off Pakistan where sedimentary input of (marine) organic matter is dominant during the NE monsoon. The explanation of apparently different organic carbon export efficiency calls for further investigations of, for example, food web structure and water column processes.

Geochemistry of lavas and dikes from the upper oceanic crust of Hess Deep Rift, eastern Pacific Ocean

Strontium isotopes are useful tracers of fluid-rock interaction in marine hydrothermal systems and provide a potential way to quantify the amount of seawater that passes through these systems. We have determined the whole-rock Sr-isotopic compositions of a section of upper oceanic crust that formed at the fast-spreading East Pacific Rise, now exposed at Hess Deep. This dataset provides the first detailed comparison for the much-studied Ocean Drilling Program (ODP) drill core from Site 504B. Whole-rock and mineral Sr concentrations indicate that Sr-exchange between hydrothermal fluids and the oceanic crust is complex, being dependent on the mineralogical reactions occurring; in particular, epidote formation takes up Sr from the fluid increasing the 87Sr/86Sr of the bulk-rock. Calculating the fluid-flux required to shift the Sr-isotopic composition of the Hess Deep sheeted-dike complex, using the approach of Bickle and Teagle (1992, doi:10.1016/0012-821X(92)90221-G) gives a fluid-flux similar to that determined for ODP Hole 504B. This suggests that the level of isotopic exchange observed in these two regions is probably typical for modern oceanic crust. Unfortunately, uncertainties in the modeling approach do not allow us to determine a fluid-flux that is directly comparable to fluxes calculated by other methods.

The CMCC Eddy-permitting Global Ocean Physical Reanalysis (C-GLORS v5, 1980-2014), links to NetCDF files

The CMCC Global Ocean Physical Reanalysis System (C-GLORS) is used to simulate the state of the ocean in the last decades. It consists of a variational data assimilation system (OceanVar), capable of assimilating all in-situ observations along with altimetry data, and a forecast step performed by the ocean model NEMO coupled with the LIM2 sea-ice model. KEY STRENGTHS: - Data are available for a large number of ocean parameters - An extensive validation has been conducted and is freely available - The reanalysis is performed at high resolution (1/4 degree) and spans the last 30 years KEY LIMITATIONS: - Quality may be discontinuos and depend on observation coverage - Uncertainty estimates are simply derived through verification skill scores

Numerical modeling of marine radiocarbon reservoir ages during the last deglaciation

A critical problem in radiocarbon dating is the spatial and temporal variability of marine reservoir ages (MRAs). We assessed the MRA evolution during the last deglaciation by numerical modeling, applying a self-consistent iteration scheme in which an existing radiocarbon chronology (derived by Hughen et al., Quat. Sci. Rev., 25, pp. 3216-3227, 2006) was readjusted by transient, 3-D simulations of marine and atmospheric Delta14C. To estimate the uncertainties regarding the ocean ventilation during the last deglaciation, we considered various ocean overturning scenarios which are based on different climatic background states (PD: modern climate, GS: LGM climate conditions). Minimum and maximum MRAs are included in file 'MRAminmax_21-14kaBP.nc'. Three further files include MRAs according to equilibrium simulations of the preindustrial ocean (file 'C14age_preindustrial.nc'; this is an update of our results published in 2005) and of the glacial ocean (files 'C14age_spinupLGM_GS.nc' and 'C14age_spinupLGM_PD.nc').

Three species distribution models for the marine diatom Fragilariopsis kerguelensis of the Southern Ocean, links to NetCDF files

Fragilariopsis kerguelensis, a dominant diatom species throughout the Antarctic Circumpolar Current, is coined to be one of the main drivers of the biological silicate pump. Here, we study the distribution of this important species and expected consequences of climate change upon it, using correlative species distribution modeling and publicly available presence-only data. As experience with SDM is scarce for marine phytoplankton, this also serves as a pilot study for this organism group. We used the maximum entropy method to calculate distribution models for the diatom F. kerguelensis based on yearly and monthly environmental data (sea surface temperature, salinity, nitrate and silicate concentrations). Observation data were harvested from GBIF and the Global Diatom Database, and for further analyses also from the Hustedt Diatom Collection (BRM). The models were projected on current yearly and seasonal environmental data to study current distribution and its seasonality. Furthermore, we projected the seasonal model on future environmental data obtained from climate models for the year 2100. Projected on current yearly averaged environmental data, all models showed similar distribution patterns for F. kerguelensis. The monthly model showed seasonality, for example, a shift of the southern distribution boundary toward the north in the winter. Projections on future scenarios resulted in a moderately to negligibly shrinking distribution area and a change in seasonality. We found a substantial bias in the publicly available observation datasets, which could be reduced by additional observation records we obtained from the Hustedt Diatom Collection. Present-day distribution patterns inferred from the models coincided well with background knowledge and previous reports about F. kerguelensis distribution, showing that maximum entropy-based distribution models are suitable to map distribution patterns for oceanic planktonic organisms. Our scenario projections indicate moderate effects of climate change upon the biogeography of F. kerguelensis.

Consumers mediate the effects of experimental ocean acidification and warming on primary producers

It is well known that ocean acidification can have profound impacts on marine organisms. However, we know little about the direct and indirect effects of ocean acidification and also how these effects interact with other features of environmental change such as warming and declining consumer pressure. In this study, we tested whether the presence of consumers (invertebrate mesograzers) influenced the interactive effects of ocean acidification and warming on benthic microalgae in a seagrass community mesocosm experiment. Net effects of acidification and warming on benthic microalgal biomass and production, as assessed by analysis of variance, were relatively weak regardless of grazer presence. However, partitioning these net effects into direct and indirect effects using structural equation modeling revealed several strong relationships. In the absence of grazers, benthic microalgae were negatively and indirectly affected by sediment-associated microalgal grazers and macroalgal shading, but directly and positively affected by acidification and warming. Combining indirect and direct effects yielded no or weak net effects. In the presence of grazers, almost all direct and indirect climate effects were nonsignificant. Our analyses highlight that (i) indirect effects of climate change may be at least as strong as direct effects, (ii) grazers are crucial in mediating these effects, and (iii) effects of ocean acidification may be apparent only through indirect effects and in combination with other variables (e.g., warming). These findings highlight the importance of experimental designs and statistical analyses that allow us to separate and quantify the direct and indirect effects of multiple climate variables on natural communities.