High-resolution glacial and deglacial record of atmospheric methane by continuous-flow and laser spectrometer analysis along the NEEM ice core

The Greenland NEEM (North Greenland Eemian Ice Drilling) operation in 2010 provided the first opportunity to combine trace-gas measurements by laser spectroscopic instruments and continuous-flow analysis along a freshly drilled ice core in a field-based setting. We present the resulting atmospheric methane (CH4) record covering the time period from 107.7 to 9.5 ka b2k (thousand years before 2000 AD). Companion discrete CH4 measurements are required to transfer the laser spectroscopic data from a relative to an absolute scale. However, even on a relative scale, the high-resolution CH4 data set significantly improves our knowledge of past atmospheric methane concentration changes. New significant sub-millennial-scale features appear during interstadials and stadials, generally associated with similar changes in water isotopic ratios of the ice, a proxy for local temperature. In addition to the midpoint of Dansgaard–Oeschger (D/O) CH4 transitions usually used for cross-dating, sharp definition of the start and end of these events brings precise depth markers (with ±20 cm uncertainty) for further cross-dating with other palaeo- or ice core records, e.g. speleothems. The method also provides an estimate of CH4 rates of change. The onsets of D/O events in the methane signal show a more rapid rate of change than their endings. The rate of CH4 increase associated with the onsets of D/O events progressively declines from 1.7 to 0.6 ppbv yr−1 in the course of marine isotope stage 3. The largest observed rate of increase takes place at the onset of D/O event #21 and reaches 2.5 ppbv yr−1.

Continuous Flow Analysis Method for Determination of Dissolved Reactive Phosphorus in Ice Cores

Phosphorus (P) is an essential macronutrient for all living organisms. Phosphorus is often present in nature as the soluble phosphate ion PO43– and has biological, terrestrial, and marine emission sources. Thus PO43– detected in ice cores has the potential to be an important tracer for biological activity in the past. In this study a continuous and highly sensitive absorption method for detection of dissolved reactive phosphorus (DRP) in ice cores has been developed using a molybdate reagent and a 2-m liquid waveguide capillary cell (LWCC). DRP is the soluble form of the nutrient phosphorus, which reacts with molybdate. The method was optimized to meet the low concentrations of DRP in Greenland ice, with a depth resolution of approximately 2 cm and an analytical uncertainty of 1.1 nM (0.1 ppb) PO43–. The method has been applied to segments of a shallow firn core from Northeast Greenland, indicating a mean concentration level of 2.74 nM (0.26 ppb) PO43– for the period 1930–2005 with a standard deviation of 1.37 nM (0.13 ppb) PO43– and values reaching as high as 10.52 nM (1 ppb) PO43–. Similar levels were detected for the period 1771–1823. Based on impurity abundances, dust and biogenic particles were found to be the most likely sources of DRP deposited in Northeast Greenland.

Solar forcing of North Atlantic surface temperature and salinity over the past millennium

There were several centennial-scale fluctuations in the climate and oceanography of the North Atlantic region over the past 1,000 years, including a period of relative cooling from about AD 1450 to 1850 known as the Little Ice Age1. These variations may be linked to changes in solar irradiance, amplified through feedbacks including the Atlantic meridional overturning circulation2. Changes in the return limb of the Atlantic meridional overturning circulation are reflected in water properties at the base of the mixed layer south of Iceland. Here we reconstruct thermocline temperature and salinity in this region from AD 818 to 1780 using paired δ18O and Mg/Ca ratio measurements of foraminifer shells from a subdecadally resolved marine sediment core. The reconstructed centennial-scale variations in hydrography correlate with variability in total solar irradiance. We find a similar correlation in a simulation of climate over the past 1,000 years. We infer that the hydrographic changes probably reflect variability in the strength of the subpolar gyre associated with changes in atmospheric circulation. Specifically, in the simulation, low solar irradiance promotes the development of frequent and persistent atmospheric blocking events, in which a quasi-stationary high-pressure system in the eastern North Atlantic modifies the flow of the westerly winds. We conclude that this process could have contributed to the consecutive cold winters documented in Europe during the Little Ice Age.

Stratigraphic Variation Within Polar Firn Caused by Differential Accumulation and Ice Flow: Interpretation of a 400 Mhz Short-Pulse Radar Profile from West Antarctica

We investigate causes of the stratigraphic variation revealed in a 177 km, 400 MHz short-pulse radar profile of firn from West Antarctica. The profile covers 56 m depth, and its direction was close to those of the ice flow and mean wind. The average, near-surface accumulation rates calculated from the time delays of one radar horizon consistently show minima on leeward slopes and maxima on windward slopes, confirming an earlier study based on stake observations. The stratigraphic variation includes up to 30 m depth variation in individual horizons over tens of km, fold limbs that become progressively steeper with depth, and fold-hinge loci that change direction or propagate down-ice with depth over distances far less than predicted by the ice speeds. We use an accumulation rate model to show how local rate anomalies and the effect of ice speed upon a periodic variation in accumulation rate cause these phenomena, and we reproduce two key features seen in the stratigraphic variations. We conclude that the model provides an explanation of changes in spatial stratigraphy and local measures of accumulation history given the constraints of surface topography, ice and wind velocities, and a general accumulation rate for an area.

Modeling of ore-forming and geoenvironmental systems: Roles of fluid flow and chemical reaction processes

Ore-forming and geoenviromental systems commonly involve coupled fluid flowand chemical reaction processes. The advanced numerical methods and computational modeling have become indispensable tools for simulating such processes in recent years. This enables many hitherto unsolvable geoscience problems to be addressed using numerical methods and computational modeling approaches. For example, computational modeling has been successfully used to solve ore-forming and mine site contamination/remediation problems, in which fluid flow and geochemical processes play important roles in the controlling dynamic mechanisms. The main purpose of this paper is to present a generalized overview of: (1) the various classes and models associated with fluid flow/chemically reacting systems in order to highlight possible opportunities and developments for the future; (2) some more general issues that need attention in the development of computational models and codes for simulating ore-forming and geoenviromental systems; (3) the related progresses achieved on the geochemical modeling over the past 50 years or so; (4) the general methodology for modeling of oreforming and geoenvironmental systems; and (5) the future development directions associated with modeling of ore-forming and geoenviromental systems.

On the interference of Kr during carbon isotope analysis of methane using continuous-flow combustion–isotope ratio mass spectrometry

Stable carbon isotope analysis of methane (delta C-13 of CH4) on atmospheric samples is one key method to constrain the current and past atmospheric CH4 budget. A frequently applied measurement technique is gas chromatography (GC) isotope ratio mass spectrometry (IRMS) coupled to a combustion-preconcentration unit. This report shows that the atmospheric trace gas krypton (Kr) can severely interfere during the mass spectrometric measurement, leading to significant biases in delta C-13 of CH4, if krypton is not sufficiently separated during the analysis. According to our experiments, the krypton interference is likely composed of two individual effects, with the lateral tailing of the doubly charged Kr-86 peak affecting the neighbouring m/z 44 and partially the m/z 45 Faraday cups. Additionally, a broad signal affecting m/z 45 and especially m/z 46 is assumed to result from scattered ions of singly charged krypton. The introduced bias in the measured isotope ratios is dependent on the chromatographic separation, the krypton-to-CH4 mixing ratio in the sample, the focusing of the mass spectrometer as well as the detector configuration and can amount to up to several per mil in delta C-13. Apart from technical solutions to avoid this interference, we present correction routines to a posteriori remove the bias.

A chrysophyte-based quantitative reconstruction of winter severity from varved lake sediments in NE Poland during the past millennium and its relationship to natural climate variability

Chrysophyte cysts are recognized as powerful proxies of cold-season temperatures. In this paper we use the relationship between chrysophyte assemblages and the number of days below 4 °C (DB4 °C) in the epilimnion of a lake in northern Poland to develop a transfer function and to reconstruct winter severity in Poland for the last millennium. DB4 °C is a climate variable related to the length of the winter. Multivariate ordination techniques were used to study the distribution of chrysophytes from sediment traps of 37 low-land lakes distributed along a variety of environmental and climatic gradients in northern Poland. Of all the environmental variables measured, stepwise variable selection and individual Redundancy analyses (RDA) identified DB4 °C as the most important variable for chrysophytes, explaining a portion of variance independent of variables related to water chemistry (conductivity, chlorides, K, sulfates), which were also important. A quantitative transfer function was created to estimate DB4 °C from sedimentary assemblages using partial least square regression (PLS). The two-component model (PLS-2) had a coefficient of determination of View the MathML sourceRcross2 = 0.58, with root mean squared error of prediction (RMSEP, based on leave-one-out) of 3.41 days. The resulting transfer function was applied to an annually-varved sediment core from Lake Żabińskie, providing a new sub-decadal quantitative reconstruction of DB4 °C with high chronological accuracy for the period AD 1000–2010. During Medieval Times (AD 1180–1440) winters were generally shorter (warmer) except for a decade with very long and severe winters around AD 1260–1270 (following the AD 1258 volcanic eruption). The 16th and 17th centuries and the beginning of the 19th century experienced very long severe winters. Comparison with other European cold-season reconstructions and atmospheric indices for this region indicates that large parts of the winter variability (reconstructed DB4 °C) is due to the interplay between the oscillations of the zonal flow controlled by the North Atlantic Oscillation (NAO) and the influence of continental anticyclonic systems (Siberian High, East Atlantic/Western Russia pattern). Differences with other European records are attributed to geographic climatological differences between Poland and Western Europe (Low Countries, Alps). Striking correspondence between the combined volcanic and solar forcing and the DB4 °C reconstruction prior to the 20th century suggests that winter climate in Poland responds mostly to natural forced variability (volcanic and solar) and the influence of unforced variability is low.

An ice-rich flow origin for the banded terrain in the Hellas basin, Mars

The interior of Hellas Basin displays a complex landscape and a variety of geomorphological domains. One of these domains, the enigmatic banded terrain covers much of the northwestern part of the basin. We use high-resolution (CTX and HiRISE) Digital Terrain Models to show that most of the complex viscous flowing behavior exhibited by the banded terrain is controlled by topography and flow-like interactions between neighboring banded terrain. Furthermore, the interior of the basin hosts several landforms suggestive of the presence of near-surface ice, which include polygonal patterns with elongated pits, scalloped depressions, isolated mounds and collapse structures. We suggest that thermal contraction cracking and sublimation of near-surface ice are responsible for the formation and the development of most of the ice-related landforms documented in Hellas. The relatively pristine form, lack of superposed craters, and strong association with the banded terrain, suggest an Amazonian (<3 Ga) age of formation for these landforms. Finally, relatively high surface pressures (above the triple point of water) expected in Hellas and summer-time temperatures often exceeding the melting point of water ice suggest that the basin may have recorded relatively “temperate” climatic conditions compared to other places on Mars. Therefore, the potentially ice-rich banded terrain may have deformed with lower viscosity and stresses compared to other locations on Mars, which may account for its unique morphology.

Investigation of hemodynamics in an in vitro system simulating left ventricular support through the right subclavian artery using 4-dimensional flow magnetic resonance imaging.

OBJECTIVES Left ventricular assist devices are an important treatment option for patients with heart failure alter the hemodynamics in the heart and great vessels. Because in vivo magnetic resonance studies of patients with ventricular assist devices are not possible, in vitro models represent an important tool to investigate flow alterations caused by these systems. By using an in vitro magnetic resonance-compatible model that mimics physiologic conditions as close as possible, this work investigated the flow characteristics using 4-dimensional flow-sensitive magnetic resonance imaging of a left ventricular assist device with outflow via the right subclavian artery as commonly used in cardiothoracic surgery in the recent past. METHODS An in vitro model was developed consisting of an aorta with its supra-aortic branches connected to a left ventricular assist device simulating the pulsatile flow of the native failing heart. A second left ventricular assist device supplied the aorta with continuous flow via the right subclavian artery. Four-dimensional flow-sensitive magnetic resonance imaging was performed for different flow rates of the left ventricular assist device simulating the native heart and the left ventricular assist device providing the continuous flow. Flow characteristics were qualitatively and quantitatively evaluated in the entire vessel system. RESULTS Flow characteristics inside the aorta and its upper branching vessels revealed that the right subclavian artery and the right carotid artery were solely supported by the continuous-flow left ventricular assist device for all flow rates. The flow rates in the brain-supplying arteries are only marginally affected by different operating conditions. The qualitative analysis revealed only minor effects on the flow characteristics, such as weakly pronounced vortex flow caused by the retrograde flow via the brachiocephalic artery. CONCLUSIONS The results indicate that, despite the massive alterations in natural hemodynamics due to the retrograde flow via the right subclavian and brachiocephalic arteries, there are no drastic consequences on the flow in the brain-feeding arteries and the flow characteristics in the ascending and descending aortas. It may be beneficial to adjust the operating condition of the left ventricular assist device to the residual function of the failing heart.

Revisiting the ISN Flow Parameters, using a variable IBEX pointing strategy

The Interstellar Boundary Explorer (IBEX) has observed the interstellar neutral (ISN) gas flow over the past 6 yr during winter/spring when the Earth's motion opposes the ISN flow. Since IBEX observes the interstellar atom trajectories near their perihelion, we can use an analytical model based upon orbital mechanics to determine the interstellar parameters. Interstellar flow latitude, velocity, and temperature are coupled to the flow longitude and are restricted by the IBEX observations to a narrow tube in this parameter space. In our original analysis we found that pointing the spacecraft spin axis slightly out of the ecliptic plane significantly influences the ISN flow vector determination. Introducing the spacecraft spin axis tilt into the analytical model has shown that IBEX observations with various spin axis tilt orientations can substantially reduce the range of acceptable solutions to the ISN flow parameters as a function of flow longitude. The IBEX operations team pointed the IBEX spin axis almost exactly within the ecliptic plane during the 2012-2014 seasons, and about 5° below the ecliptic for half of the 2014 season. In its current implementation the analytical model describes the ISN flow most precisely for the spin axis orientation exactly in the ecliptic. This analysis refines the derived ISN flow parameters with a possible reconciliation between velocity vectors found with IBEX and Ulysses, resulting in a flow longitude lambda∞ = 74.°5 ± 1.°7 and latitude beta∞ = -5.°2 ± 0.°3, but at a substantially higher ISN temperature than previously reported.

Static GPS measurements and elevation model in the vicinity of the EDML deep-drilling site

Interpretation of ice-core records requires accurate knowledge of the past and present surface topography and stress-strain fields. The European Project for Ice Coring in Antarctica (EPICA) drilling site (0.0684° E and 75.0025° S, 2891.7 m) in Dronning Maud Land, Antarctica, is located in the immediate vicinity of a transient and splitting ice divide. A digital elevation model is determined from the combination of kinematic GPS measurements with the GLAS12 data sets from the ICESat satellite. Based on a network of stakes, surveyed with static GPS, the velocity field around the EDML drilling site is calculated. The annual mean velocity magnitude of 12 survey points amounts to 0.74 m/a. Flow directions mainly vary according to their distance from the ice divide. Surface strain rates are determined from a pentagon-shaped stake network with one center point, close to the drilling site. The strain field is characterised by along flow compression, lateral dilatation, and vertical layer thinning.

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.

Planimetric ice-flow convergence and flow-orthonormal strain rate across the Antarctic Ice Sheet, with links to model result files

Fast-flowing ice streams discharge most of the ice from the interior of the Antarctic Ice Sheet coastward. Understanding how their tributary organisation is governed and evolves is essential for developing reliable models of the ice sheet's response to climate change. Despite much research on ice-stream mechanics, this problem is unsolved, because the complexity of flow within and across the tributary networks has hardly been interrogated. Here I present the first map of planimetric flow convergence across the ice sheet, calculated from satellite measurements of ice surface velocity, and use it to explore this complexity. The convergence map of Antarctica elucidates how ice-stream tributaries draw ice from the interior. It also reveals curvilinear zones of convergence along lateral shear margins of streaming, and abundant convergence ripples associated with nonlinear ice rheology and changes in bed topography and friction. Flow convergence on ice-stream tributaries and their feeding zones is markedly uneven, and interspersed with divergence at distances of the order of kilometres. For individual drainage basins as well as the ice sheet as a whole, the range of convergence and divergence decreases systematically with flow speed, implying that fast flow cannot converge or diverge as much as slow flow. I therefore deduce that flow in ice-stream networks is subject to mechanical regulation that limits flow-orthonormal strain rates. These properties and the gridded data of convergence and flow-orthonormal strain rate in this archive provide targets for ice- sheet simulations and motivate more research into the origin and dynamics of tributarization.