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.

Optimizing Models of the North Atlantic Spring Bloom Using Physical, Chemical and Bio-Optical Observations from a Lagrangian Float

The North Atlantic spring bloom is one of the main events that lead to carbon export to the deep ocean and drive oceanic uptake of CO(2) from the atmosphere. Here we use a suite of physical, bio-optical and chemical measurements made during the 2008 spring bloom to optimize and compare three different models of biological carbon export. The observations are from a Lagrangian float that operated south of Iceland from early April to late June, and were calibrated with ship-based measurements. The simplest model is representative of typical NPZD models used for the North Atlantic, while the most complex model explicitly includes diatoms and the formation of fast sinking diatom aggregates and cysts under silicate limitation. We carried out a variational optimization and error analysis for the biological parameters of all three models, and compared their ability to replicate the observations. The observations were sufficient to constrain most phytoplankton-related model parameters to accuracies of better than 15 %. However, the lack of zooplankton observations leads to large uncertainties in model parameters for grazing. The simulated vertical carbon flux at 100 m depth is similar between models and agrees well with available observations, but at 600 m the simulated flux is larger by a factor of 2.5 to 4.5 for the model with diatom aggregation. While none of the models can be formally rejected based on their misfit with the available observations, the model that includes export by diatom aggregation has a statistically significant better fit to the observations and more accurately represents the mechanisms and timing of carbon export based on observations not included in the optimization. Thus models that accurately simulate the upper 100 m do not necessarily accurately simulate export to deeper depths.

A Lagrangian Analysis of the Northern Hemisphere Subtropical Jet

This study presents a 5-yr climatology of 7-day back trajectories started from the Northern Hemisphere subtropical jet. These trajectories provide insight into the seasonally and regionally varying angular momentum and potential vorticity characteristics of the air parcels that end up in the subtropical jet. The trajectories reveal preferred pathways of the air parcels that reach the subtropical jet from the tropics and the extratropics and allow estimation of the tropical and extratropical forcing of the subtropical jet. The back trajectories were calculated 7 days back in time and started every 6 h from December 2005 to November 2010 using the Interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-Interim) dataset as a basis. The trajectories were started from the 345-K isentrope in areas where the wind speed exceeded a seasonally varying threshold and where the wind shear was confined to upper levels. During winter, the South American continent, the Indian Ocean, and the Maritime Continent are preferred areas of ascent into the upper troposphere. From these areas, air parcels follow an anticyclonic pathway into the subtropical jet. During summer, the majority of air parcels ascend over the Himalayas and Southeast Asia. Angular momentum is overall well conserved for trajectories that reach the subtropical jet from the deep tropics. In winter and spring, the hemispheric-mean angular momentum loss amounts to approximately 6%; in summer, it amounts to approximately 18%; and in fall, it amounts to approximately 13%. This seasonal variability is confirmed using an independent potential vorticity–based method to estimate tropical and extratropical forcing of the subtropical jet.

Ingestion rate and egg production of copepods collected in the upper 20m in the eastern Mediterranean Sea in April 2008 during SES_GR2

The ingestion on ciliates and phytoplankton dataset is based on samples taken during October 2008 in Northern Aegean Sea, the area influenced by the Black Sea water outflow. A Lagrangian experiment was established and copepod ingestion was estimated from experiments performed at stations according to the different positions of drifters during the cruise. Copepods for the experiments were obtained with slow non-quantitative tows from the upper 20 m layer of the water column using 200 µm mesh size nets fitted with a large non-filtering cod end. For the grazing experiments we used the following copepod species: Clausocalanus furcatus, and Temoraa stylifera according to the relevant reference (Bamstedt et al. 2000). Copepod clearance rates on ciliates were calculated according to Frost equations (Frost 1972). Ingestion rates were calculated by multiplying clearance rates by the initial standing stocks (Bamstedt et al. 2000). The egg production dataset is based on samples taken during October 2008 in Northern Aegean Sea, the area influenced by the Black Sea water outflow. A Lagrangian experiment was established and copepod egg production was estimated from experiments performed at stations according to the different positions of drifters during the cruise. Egg production rates of the dominant calanoid copepods were determined by incubation of fertilised females (eggs female/day) collected in the 0-20m layer. Copepod egg production was measured for the copepods Clausocalanus furcatus, Temora stylifera. On board experiments for the estimation of copepod egg production were taken place. For the estimation of copepod production (mgC/m**2/day), lengths (copepods and eggs) were converted to body carbon (Hopcroft et al., 1998) and production was estimated from biomass and weight-specific egg production rates, by assuming that those rates are representative for juvenile specific growth rates (Berggreen et al., 1988).

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

The MEDESS-GIB drifters database

On September 9th 2014, an intensive drifter deployment was carried out in the Strait of Gibraltar. In the frame of the EU MED Program MEDESS-4MS, the MEDESS-GIB experiment consisted of the deployment of 35 satellite tracked drifters, mostly of CODE-type, equipped with temperature sensor sampling at a rate of 30 minutes. Drifters were distributed along and on both sides of the Strait of Gibraltar. The MEDESS-GIB deployment plan was designed as to ensure quasi-synoptic spatial coverage. To this end, 4 boats covering an area of about 680 NM2 in 6 hours were coordinated. As far as authors know, this experiment is the most important exercise in the area in terms of number of drifters released. Collected satellite-tracked data along drifter trajectories have been quality controlled and processed to build the here presented MEDESS-GIB data set.

Ingestion rate and egg production of copepods collected in the upper 20m in the eastern Mediterranean Sea in March and April 2008 during SES_GR1

The ingestion on ciliates and phytoplankton dataset is based on samples taken during April 2008 in Northern Aegean Sea, the area influenced by the Black Sea water outflow. A Lagrangian experiment was established and copepod ingestion was estimated from experiments performed at stations according to the different positions of drifters during the cruise. Copepods for the experiments were obtained with slow non-quantitative tows from the upper 20 m layer of the water column using 200 µm mesh size nets fitted with a large non-filtering cod end. For the grazing experiments we used the following copepod species: Centropages typicus and Calanus helgolandicus according to the relevant reference (Bamstedt et al. 2000). Copepod clearance rates on ciliates were calculated according to Frost equations (Frost 1972). Ingestion rates were calculated by multiplying clearance rates by the initial standing stocks (Bamstedt et al. 2000). The egg production dataset is based on samples taken during April 2008 in Northern Aegean Sea, the area influenced by the Black Sea water outflow. A Lagrangian experiment was established and copepod egg production was estimated from experiments performed at stations according to the different positions of drifters during the cruise. Egg production rates of the dominant calanoid copepods were determined by incubation of fertilised females (eggs female/day) collected in the 0-20m layer. Copepod egg production was measured for the copepods Centropages typicus, Calanus helgolandicus. On board experiments for the estimation of copepod egg production were taken place. For the estimation of copepod production (mgC/ m**2 /day), lengths (copepods and eggs) were converted to body carbon (Hopcroft et al., 1998) and production was estimated from biomass and weight-specific egg production rates, by assuming that those rates are representative for juvenile specific growth rates (Berggreen et al., 1988).

Airport Taxi Planning: Lagrangian Decomposition

The airport taxi planning (TP) module is a decision tool intended to guide airport surface management operations. TP is defined by a flow network optimization model that represents flight ground movements and improves aircraft taxiing routes and schedules during periods of aircraft congestion. TP is not intended to operate as a stand‐alone tool for airport operations management: on the contrary, it must be used in conjunction with existing departing and arriving traffic tools and overseen by the taxi planner of the airport, also known as the aircraft ground controller. TP must be flexible in order to accommodate changing inputs while maintaining consistent routes and schedules already delivered from past executions. Within this dynamic environment, the execution time of TP may not exceed a few minutes. Classic methods for solving binary multi‐commodity flow networks with side constraints are not efficient enough; therefore, a Lagrangian decomposition methodology has been adapted to solve it. We demonstrate TP Lagrangian decomposition using actual data from the Madrid‐Barajas Airport

Lagrangian skeleton on the Gulf of Mexico oil spill

The understanding of the circulation of ocean currents, the exchange of CO2 between atmosphere and oceans, and the influence of the oceans on the distribution of heat on a global scale is key to our ability to predict and assess the future evolution of climate [1, 2]. Global climate change is affecting sea breathing through mechanisms not yet understood.

A Lagrangian descriptor applied to the Kuroshio current

The understanding of the circulation of ocean currents, the exchange of CO2 between atmosphere and oceans, and the in uence of the oceans on the distribution of heat on a global scale is key to our ability to predict and assess the future evolution of climate.

Involutivity of the Hamilton-Cartan equations of a second-order Lagrangian admitting a first-order Hamiltonian formalism

Involutivity of the Hamilton-Cartan equations of a second-order Lagrangian admitting a first-order Hamiltonian formalism