Influence of the Central American Seaway and Drake Passage on ocean circulation and neodymium isotopes: A model study

The sensitivity of the neodymium isotopic composition (ϵNd) to tectonic rearrangements of seaways is investigated using an Earth System Model of Intermediate Complexity. The shoaling and closure of the Central American Seaway (CAS) is simulated, as well as the opening and deepening of Drake Passage (DP). Multiple series of equilibrium simulations with various intermediate depths are performed for both seaways, providing insight into ϵNd and circulation responses to progressive throughflow evolutions. Furthermore, the sensitivity of these responses to the Atlantic Meridional Overturning Circulation (AMOC) and the neodymium boundary source is examined. Modeled ϵNd changes are compared to sediment core and ferromanganese (Fe-Mn) crust data. The model results indicate that the North Atlantic ϵNd response to the CAS shoaling is highly dependent on the AMOC state, i.e., on the AMOC strength before the shoaling to shallow depths (preclosure). Three scenarios based on different AMOC forcings are discussed, of which the model-data agreement favors a shallow preclosure (Miocene) AMOC (∼6 Sv). The DP opening causes a rather complex circulation response, resulting in an initial South Atlantic ϵNd decrease preceding a larger increase. This feature may be specific to our model setup, which induces a vigorous CAS throughflow that is strongly anticorrelated to the DP throughflow. In freshwater experiments following the DP deepening, ODP Site 1090 is mainly influenced by AMOC and DP throughflow changes, while ODP Site 689 is more strongly influenced by Southern Ocean Meridional Overturning Circulation and CAS throughflow changes. The boundary source uncertainty is largest for shallow seaways and at shallow sites.

Modelling Nd-isotopes with a coarse resolution ocean circulation model: Sensitivities to model parameters and source/sink distributions

The neodymium (Nd) isotopic composition (Nd) of seawater is a quasi-conservative tracer of water mass mixing and is assumed to hold great potential for paleoceanographic studies. Here we present a comprehensive approach for the simulation of the two neodymium isotopes 143Nd, and 144Nd using the Bern3D model, a low resolution ocean model. The high computational efficiency of the Bern3D model in conjunction with our comprehensive approach allows us to systematically and extensively explore the sensitivity of Nd concentrations and Nd to the parametrisation of sources and sinks. Previous studies have been restricted in doing so either by the chosen approach or by computational costs. Our study thus presents the most comprehensive survey of the marine Nd cycle to date. Our model simulates both Nd concentrations as well as Nd in good agreement with observations. Nd covaries with salinity, thus underlining its potential as a water mass proxy. Results confirm that the continental margins are required as a Nd source to simulate Nd concentrations and Nd consistent with observations. We estimate this source to be slightly smaller than reported in previous studies and find that above a certain magnitude its magnitude affects Nd only to a small extent. On the other hand, the parametrisation of the reversible scavenging considerably affects the ability of the model to simulate both, Nd concentrations and Nd. Furthermore, despite their small contribution, we find dust and rivers to be important components of the Nd cycle. In additional experiments, we systematically varied the diapycnal diffusivity as well as the Atlantic-to-Pacific freshwater flux to explore the sensitivity of Nd concentrations and its isotopic signature to the strength and geometry of the overturning circulation. These experiments reveal that Nd concentrations and Nd are comparatively little affected by variations in diapycnal diffusivity and the Atlantic-to-Pacific freshwater flux. In contrast, an adequate representation of Nd sources and sinks is crucial to simulate Nd concentrations and Nd consistent with observations. The good agreement of our results with observations paves the way for the evaluation of the paleoceanographic potential of Nd in further model studies.

Simultaneous Determination of Stable Carbon, Oxygen, and Hydrogen Isotopes in Cellulose

A technological development is described through which the stable carbon-, oxygen-, and nonexchangeable hydrogen-isotopic ratios (δ13C,δ18O,δ2H) are determined on a single carbohydrate (cellulose) sample with precision equivalent to conventional techniques (δ13 C 0.15‰,δ18O 0.30‰,δ2H 3.0‰). This triple-isotope approach offers significant new research opportunities, most notably in physiology and medicine, isotope biogeochem- istry, forensic science, and palaeoclimatology, when isotopic analysis of a common sample is desirable or when sample material is limited.

Measurement of flow harmonics with multi-particle cumulants in Pb+Pb collisions at √sNN=2.76 TeV with the ATLAS detector

ATLAS measurements of the azimuthal anisotropy in lead–lead collisions at √sNN = 2.76 TeV are shown using a dataset of approximately 7μb−1 collected at the LHC in 2010. The measurements are performed for charged particles with transversemomenta 0.5 < pT < 20 GeV and in the pseudorapidity range |η| < 2.5. The anisotropy is characterized by the Fourier coefficients, vn, of the charged-particle azimuthal angle distribution for n = 2–4. The Fourier coefficients are evaluated using multi-particle cumulants calculated with the generating function method. Results on the transverse momentum, pseudorapidity and centrality dependence of the vn coefficients are presented. The elliptic flow, v2, is obtained from the two-, four-, six- and eight-particle cumulants while higher-order coefficients, v3 and v4, are determined with two- and four-particle cumulants. Flow harmonics vn measured with four-particle cumulants are significantly reduced compared to the measurement involving two-particle cumulants. A comparison to vn measurements obtained using different analysis methods and previously reported by the LHC experiments is also shown. Results of measurements of flow fluctuations evaluated with multiparticle cumulants are shown as a function of transverse momentum and the collision centrality. Models of the initial spatial geometry and its fluctuations fail to describe the flow fluctuations measurements.

Measurement of inclusive jet charged-particle fragmentation functions in Pb+Pb collisions at √sNN = 2.76 TeV with the ATLAS detector

Measurements of charged-particle fragmentation functions of jets produced in ultra-relativistic nuclear collisions can provide insight into the modification of parton showers in the hot, dense medium created in the collisions. ATLAS has measured jets in √sNN=2.76 TeV Pb+Pb collisions at the LHC using a data set recorded in 2011 with an integrated luminosity of 0.14 nb−1. Jets were reconstructed using the anti-kt algorithm with distance parameter values R = 0.2, 0.3, and 0.4. Distributions of charged-particle transverse momentum and longitudinal momentum fraction are reported for seven bins in collision centrality for R=0.4 jets with pjetT>100 GeV. Commensurate minimum pT values are used for the other radii. Ratios of fragment distributions in each centrality bin to those measured in the most peripheral bin are presented. These ratios show a reduction of fragment yield in central collisions relative to peripheral collisions at intermediate z values, 0.04≲z≲0.2 and an enhancement in fragment yield for z≲0.04. A smaller, less significant enhancement is observed at large z and large pT in central collisions.

Pb pollution from leaded gasoline in South America in the context of a 2000-year metallurgical history

Exploitation of the extensive polymetallic deposits of the Andean Altiplano in South America since precolonial times has caused substantial emissions of neurotoxic lead (Pb) into the atmosphere; however, its historical significance compared to recent Pb pollution from leaded gasoline is not yet resolved. We present a comprehensive Pb emission history for the last two millennia for South America, based on a continuous, high-resolution, ice core record from Illimani glacier. Illimani is the highest mountain of the eastern Bolivian Andes and is located at the northeastern margin of the Andean Altiplano. The ice core Pb deposition history revealed enhanced Pb enrichment factors (EFs) due to metallurgical processing for silver production during periods of the Tiwanaku/Wari culture (AD 450–950), the Inca empires (AD 1450–1532), colonial times (AD 1532–1900), and tin production at the beginning of the 20th century. After the 1960s, Pb EFs increased by a factor of 3 compared to the emission level from metal production, which we attribute to gasoline-related Pb emissions. Our results show that anthropogenic Pb pollution levels from road traffic in South America exceed those of any historical metallurgy in the last two millennia, even in regions with exceptional high local metallurgical activity.

Long-term summer temperature variations in the Pyrenees from detrended stable carbon isotopes

Substantial effort has recently been put into the development of climate reconstructions from tree-ring stable carbon isotopes, though the interpretation of long-term trends retained in such timeseries remains challenging. Here we use detrended δ13C measurements in Pinus uncinata tree-rings, from the Spanish Pyrenees, to reconstruct decadal variations in summer temperature back to the 13th century. The June-August temperature signal of this reconstruction is attributed using decadally as well as annually resolved, 20th century δ13C data. Results indicate that late 20th century warming has not been unique within the context of the past 750 years. Our reconstruction contains greater am-plitude than previous reconstructions derived from traditional tree-ring density data, and describes particularly cool conditions during the late 19th century. Some of these differences, including early warm periods in the 14th and 17th centuries, have been retained via δ13C timeseries detrending - a novel approach in tree-ring stable isotope chronology development. The overall reduced variance in earlier studies points to an underestimation of pre-instrumental summer temperature variability de-rived from traditional tree-ring parameters.

Multiple metamorphic stages within an eclogite-facies terrane (Sesia Zone, Western Alps) revealed by Th-U-Pb petrochronology

Convergent plate margins typically experience a transition from subduction to collision dynamics as massive continental blocks enter the subduction channel. Studies of high-pressure rocks indicate that tectonic fragments are rapidly exhumed from eclogite facies to midcrustal levels, but the details of such dynamics are controversial.To understand the dynamics of a subduction channel we report the results of a petrochronological study from the central Sesia Zone, a key element of the internalWestern Alps.This comprises two polymetamorphic basement complexes (Eclogitic Micaschist Complex and Gneiss Minuti Complex) and a thin, dismembered cover sequence (Scalaro Unit) associated with pre-Alpine metagabbros and metasediments (Bonze Unit). Structurally controlled samples from three of these units (Eclogitic Micaschist Complex and Scalaro-Bonze Units) yield unequivocal petrological and geochronological evidence of two distinct high-pressure stages. Ages (U-Th-Pb) of growth zones in accessory allanite and zircon, combined with inclusion and textural relationships, can be tied to the multi-stage evolution of single samples.Two independent tectono-metamorphic ‘slices’ showing a coherent metamorphic evolution during a given time interval have been recognized: the Fondo slice (which includes Scalaro and Bonze rocks) and the Druer slice (belonging to the Eclogitic Micaschist Complex).The new data indicate separate stages of deformation at eclogite-facies conditions for each recognized independent kilometer-sized tectono-metamorphic slice, between ~85 and 60 Ma, with evidence of intermittent decompression (∆P~0.5 GPa) within only the Fondo slice. The evolution path of the Druer slice indicates a different P-T-time evolution with prolonged eclogite-facies metamorphism between ~85 and 75Ma. Our approach, combining structural, petrological and geochronological techniques, yields field-based constraints on the duration and rates of dynamics within a subduction channel.

Accumulation rates and predominant atmospheric sources of natural and anthropogenic Hg and Pb on the Faroe Islands

A monolith representing 5420 14C yr of peat accumulation was collected from a blanket bog at Myrarnar, Faroe Islands. The maximum Hg concentration (498 ng/g at a depth of 4.5 cm) coincides with the maximum concentration of anthropogenic Pb (111 μg/g). Age dating of recent peat accumulation using 210Pb (CRS model) shows that the maxima in Hg and Pb concentrations occur at AD 1954 ± 2. These results, combined with the isotopic composition of Pb in that sample (206Pb/207Pb = 1.1720 ± 0.0017), suggest that coal burning was the dominant source of both elements. From the onset of peat accumulation (ca. 4286 BC) until AD 1385, the ratios Hg/Br and Hg/Se were constant (2.2 ± 0.5 × 10-4 and 8.5 ± 1.8 × 10-3, respectively). Since then, Hg/Br and Hg/Se values have increased, also reaching their maxima in AD 1954. The age date of the maximum concentrations of anthropogenic Hg and Pb in the Faroe Islands is consistent with a previous study of peat cores from Greenland and Denmark (dated using the atmospheric bomb pulse curve of 14C), which showed maximum concentrations in AD 1953. The average rate of atmospheric Hg accumulation from 1520 BC to AD 1385 was 1.27 ± 0.38 μg/m2/yr. The Br and Se concentrations and the background Hg/Br and Hg/Se ratios were used to calculate the average rate of natural Hg accumulation for the same period, 1.32 ± 0.36 μg/m2/yr and 1.34 ± 0.29 μg/m2/yr, respectively. These fluxes are similar to the preanthropogenic rates obtained using peat cores from Switzerland, southern Greenland, southern Ontario, Canada, and the northeastern United States. Episodic volcanic emissions and the continual supply of marine aerosols to the Faroe Islands, therefore, have not contributed significantly to the Hg inventory or the Hg accumulation rates, relative to these other areas. The maximum rate of Hg accumulation was 34 μg/m2/yr. The greatest fluxes of anthropogenic Hg accumulation calculated using Br and Se, respectively, were 26 and 31 μg/m2/yr. The rate of atmospheric Hg accumulation in 1998 (16 μg/m2/yr) is comparable to the values recently obtained by atmospheric transport modeling for Denmark, the Faroe Islands, and Greenland.

The capabilities of ROSINA/DFMS to measure argon isotopes at comet 67P/Churyumov–Gerasimenko

Little is known about the noble gas abundances in comets. These highly volatile atoms are possible tracers of the history of cometary matter including the thermal evolution. They can help quantify the contribution of cometary impacts to terrestrial oceans and help elucidate on the formation history of comets and their role in the formation and evolution of planetary atmospheres. This paper focuses on argon and the capabilities to measure this noble gas with in situ mass spectrometry at comet 67P/Churyumov–Gerasimenko, the target of the European Space Agency׳s spacecraft Rosetta. Argon may have been detected by remote sensing in a single Oort cloud comet but to date nothing is known about the isotopic abundances of argon in comets. Furthermore, no detection of argon in a Jupiter-family comet has been reported. Comet 67P/Churyumov–Gerasimenko belongs to the group of Jupiter-family comets and originates most likely in the Kuiper belt. Onboard Rosetta is ROSINA/DFMS (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis/Double Focusing Mass Spectrometer). DFMS has unprecedented mass resolution and high sensitivity and is designed to measure isotopic ratios including argon (Balsiger et al., 2007). Argon measurements using the DFMS lab model (identical to the flight model) demonstrate this capability. At very least, this mass spectrometer has the resolution and sensitivity to reduce the upper limit on the argon outgassing rate relative to water by more than three orders of magnitude (for 38Ar). Most likely, ROSINA/DFMS will provide the first detection of argon in a Jupiter-family comet together with the first determination of the ³⁶Ar/³⁸Ar ratio at a comet.

Bryozoan stable carbon and hydrogen isotopes: relationships between the isotopic composition of zooids, statoblasts and lake water

We explored the extent to which δ13C and δD values of freshwater bryozoan statoblasts can provide information about the isotopic composition of zooids, bryozoan food and surrounding water. Bryozoan samples were collected from 23 sites and encompassed ranges of nearly 30‰ for δ13C and 100‰ for δD values. δ13C offsets between zooids and statoblasts generally ranged from −3 to +4.5‰, with larger offsets observed in four samples. However, a laboratory study with Plumatella emarginata and Lophopus crystallinus demonstrated that, in controlled settings, zooids had only 0–1.2‰ higher δ13C values than statoblasts, and 1.7‰ higher values than their food. At our field sites, we observed a strong positive correlation between median δ13C values of zooids and median δ13C values of corresponding statoblasts. We also observed a positive correlation between median δD values of zooids and statoblasts for Plumatella, and a positive correlation between median δD values of statoblasts and δD values of lake water for Plumatella and when all bryozoan taxa were examined together. Our results suggest that isotope measurements on statoblasts collected from flotsam or sediment samples can provide information on the feeding ecology of bryozoans and the H isotopic composition of lake water.

TIMS measurements of full range of natural Ca isotopes with internally consistent fractionation correction

This study presents static measurements of the Ca isotopic composition of standard reference materials SRM 915 a/b on a Triton Plus™ thermal ionization mass spectrometer with a specially developed Faraday cup array allowing simultaneous measurement of 40Ca and 48Ca. The total amount of Ca in all analyses was kept < 1 µg. With this setup the measurement uncertainties were 0.06 ‰ for 40Ca/44Ca and 0.12 ‰ for 48Ca/40Ca. Measuring all isotopes simultaneously better allows to test the internal consistency of different Ca isotope abundances reported in the literature. The exponential law was observed to correct incompletely instrumental mass fractionation. An improved fractionation correction based on the exponential law is proposed. It changes the 40Ca/44Ca ratio of SRM 915a (corrected relative to 42Ca/44Ca = 0.31221; 48Ca/44Ca = 0.08871) from 47.1635 ± 0.0028 to 47.1649 ± 0.0047. The measurements of SRM 915b were performed with different analytical conditions (runs were prolonged till complete filament load depletion). Even if the 40Ca/44Ca ratio of SRM 915b, when corrected with the simple exponential law, appears different (47.1532 ± 0.0038) from that of SRM 915a, it becomes coincident (47.1613 ± 0.0028) when corrected with a second-order refinement. This supports the use of the improved exponential law to obtain internally consistent Ca isotope ratio for natural samples.

Oxygen isotopes of lake marl at Gerzensee and Leysin (Switzerland), covering the Younger Dryas and two minor oscillations, and their correlation to the GRIP ice core

The ratio of oxygen isotopes is a temperature proxy both in precipitation and in the calcite of lacustrine sediments. The very similar oxygen-isotope records from Greenland ice cores and European lake sediments during the Last Glacial Termination suggest that the drastic climatic changes occurred quasi-simultaneously on an extra-regional, probably hemispheric scale. In order to study temporal relations of the different parameters recorded in lake sediments, for example biotic response times to rapid climatic changes, a precise chronology is required. In unlaminated lake sediments there is not yet available a method to provide a high-resolution chronology, especially for periods with radiocarbon plateaux. Alternatively, an indirect time scale can be constructed by linking the lake stratigraphy with other well-dated climate records. New oxygen-isotope records from Gerzensee and Leysin, with an estimated sampling resolution of between 15 and 40 years, match the Greenlandic isotope record in many details. Under the assumption that the main variations in temperature and thus in oxygen isotopes occurred about simultaneously in Greenland and Switzerland, we have assigned a time scale to the lake sediments of Gerzensee and Leysin by wiggle-matching their stable-isotope records with those of Greenland ice cores, which are among the best dated climatic archives. We estimate a precision of 20 to 100 years during the Last Glacial Termination.

Stable carbon and oxygen isotopes in tree rings show physiological responses of Pericopsis elata to precipitation in the Congo Basin

In equatorial regions, where tree rings are less distinct or even absent, the response of forests to high-frequency climate variability is poorly understood. We measured stable carbon and oxygen isotopes in anatomically distinct, annual growth rings of four Pericopsis elata trees from a plantation in the Congo Basin, to assess their sensitivity to recorded changes in precipitation over the last 50 y. Our results suggest that oxygen isotopes have high common signal strength (EPS = 0.74), and respond to multi-annual precipitation variability at the regional scale, with low δ18O values (28–29‰) during wetter conditions (1960–1970). Conversely, δ13C are mostly related to growth variation, which in a light-demanding species are driven by competition for light. Differences in δ13C values between fast- and slow-growing trees (c. 2‰), result in low common signal strength (EPS = 0.37) and are driven by micro-site conditions rather than by climate. This study highlights the potential for understanding the causes of growth variation in P. elata as well as past hydroclimatic changes, in a climatically complex region characterized by a bimodal distribution in precipitation.