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.

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 numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.

Molar gas ratios and stable oxygen and nitrogen isotopes of air entrapped in ice, North America

The molar ratios of atmospheric gases change during dissolution in water due to differences in their relative solubilities. We exploited this characteristic to develop a tool to clarify the origin of ice formations in permafrost regions. Extracted from ice, molar gas ratios can distinguish buried glacier ice from intrasedimental ground ice formed by freezing groundwaters. An extraction line was built to isolate gases from ice by melting and trapping with liquid He, followed by analysis of N2, O2, Ar, 18O-O2 and 15N-N2, by continuous flow mass spectrometry. The method was tested using glacier ice, aufeis ice (river icing) and intrasedimental ground ice from sites in the Canadian Arctic. O2/Ar and N2/Ar ratios clearly distinguish between atmospheric gas in glacial ice and gases from intrasedimental ground ice, which are exsolved from freezing water. 615NN2 and 618OO2 in glacier ice, aufeis ice and intrasedimental ground ice do not show clear distinguishing trends as they are affected by various physical processes during formation such as gravitational settling, excess air addition, mixing with snow pack, and respiration.

Heat flow measurements in the Weddell Sea

From January to March 1987, heat flow measurements were tried at four sites (Sites 689, 690, 695, and 696) during ODP Leg 113, in the Weddell Sea, Antarctica. At Site 690 (Maud Rise), a convex upward shaped temperature vs. depth profile was observed. This profile cannot be explained by steady-state conduction through solid materials only. We conclude that the minimum heat flow value at Site 690 is 45 mW/m2. A prominent bottom simulating reflector (BSR) was observed at 600 mbsf at Site 695. However, the observed temperature is too high to explain the BSR as a gas hydrate. The origin of the BSR remains unknown, although it is probably of biogenic origin as observed in the Bering Sea during DSDP Leg 19. After correcting for the effects of sedimentation, heat flow values at Sites 695 and 696 are 69 and 63 mW/m2, respectively. Furthermore, we compiled heat flow data south of 50°S. In the Weddell Sea region, the eastern part shows relatively low heat flow in comparison with the western part, with the boundary between them at about 15°W longitude.

Heat flow measurements on METEOR cruise M86/5

Mud volcanoes (MV) are sources of mass and energy, transported from deeper levels of the sediment pile to the surface. Together with fluid and gas, thermal energy is emitted through these structures. Therefore heat flow determination is a sensible tool to detect and quantify the amount of convective flow. In the Gulf of Cadiz several mud volcanoes can be found along major tectonic lines (SWIM faults). We employ geothermal measurements to observe the activity of mud volcanoes and possible leakage at the faults apart from pronounced structures.

Heat flow values in the Sea of Okhotsk region

Tectonic structure and anomalous distributions of geophysical fields of the Sea of Okhotsk region are considered; the lack of reliable data on age of the lithosphere beneath basins of various origin in the Sea of Okhotsk is noted. Model calculations based on geological and geophysical data yielded 65 Ma (Cretaceous-Paleocene boundary) age for the Central Okhotsk rise underlain by the continental lithosphere. This estimate agrees with the age (the end of Cretaceous) derived from seismostratigraphic data. A comparative analysis of theoretical and measured heat flows in the Akademii Nauk Rise, underlain by the thinned continental crust, is performed. The analysis points to a higher (by 20%) value of the measured thermal background of the rise, which is consistent with high negative gradient of gravity anomalies in this area. Calculations yielded 36 Ma (Early Oligocene) age and lithosphere thickness of 50 km for the South Okhotsk depression, whose seafloor was formed by processes of back-arc spreading. The estimated age of the depression is supported by kinematic data on the region; the calculated thickness of the lithosphere coincides with the value estimated from data of magnetotelluric sounding here. This indicates that formation time (36 Ma) of the South Okhotsk depression was estimated correctly. Numerical modeling performed for determination of the basement age of rifting basins in the Sea of Okhotsk gave the following estimates: 18 Ma (Early Miocene) for the Deryugin Basin, 12 Ma (Middle Miocene) for the TINRO Basin, and 23 Ma (Late Oligocene) for the West Kamchatka Trough. These estimates agree with formation time (Oligocene-Quaternary) of the sedimentary cover in rifting basins of the Sea of Okhotsk derived from geological and geophysical data. Model temperature estimates are obtained for lithologic and stratigraphic boundaries of the sedimentary cover in the Deryugin and TINRO Basins and the West Kamchatka Trough; the temperature analysis indicates that the latter two structures are promising for oil and hydrocarbon gas generation; the West Kamchatka Trough possesses better reservoir properties compared to the TINRO and Deryugin Basins. The latter is promising for generation of hydrocarbon gas. Paleogeodynamic reconstructions of the Sea of Okhotsk region evolution are obtained for times of 90, 66, and 36 Ma on the base of kinematic, geomagnetic, structural, tectonic, geothermal, and other geological and geophysical data.

3He/4He and 4He/20Ne ratios in pore water and gas in Pacific deep-sea sediments

We have measured the 3He/3He and 3He/20 Ne ratios of thirty-nine pore water and gas samples in deep-sea sediments collected at twelve sites on the Pacific Ocean bottom during the cruises of Deep Sea Drilling Project Legs 87, 89, 90 and 92. The 3He/4He and 4He/20Ne ratios vary from 0.000000215 to 0.00000165 and from 0.29 to 20, respectively. He in the sample is composed of four components: (1) atmospheric He dissolved in seawater; (2) atmospheric He with mantle-derived He in Pacific bottom water; (3) in situ radiogenic He in the sediment; and (4) crustal He in the basement rock. Assuming that the 20Ne contents are constant with the value of seawater, the depth variations in the 4He/20Ne ratios at five Sites, 583D, 594, 597A, 598A and 504B, may provide useful information on 4He flux at the ocean bottom. The estimated 4He fluxes vary from 2000 to 40000 atoms cm**-2 s**-1 and are one to three orders of magnitude less than those calculated from the excess He in deep ocean water. An overall similarity between the geographical distribution of the 3He/4He ratios and heat flow data is found in the study area, between the East Pacific Rise across the Pacific Ocean and the Japanese Islands. The tendency is well explained by a conventional sea-floor spreading model.

Gas emissions and ROV survey tracks of the Regab pockmark (Northern Congo Fan)

Pockmarks are seafloor depressions commonly associated with fluid escape from the seabed and are believed to contribute noticeably to the transfer of methane into the ocean and ultimately into the atmosphere. They occur in many different areas and geological contexts, and vary greatly in size and shape. Nevertheless, the mechanisms of pockmark growth are still largely unclear. Still, seabed methane emissions contribute to the global carbon budget, and understanding such processes is critical to constrain future quantifications of seabed methane release at local and global scales. The giant Regab pockmark (9°42.6' E, 5°47.8' S), located at 3160 m water depth near the Congo deep-sea channel (offshore southwestern Africa), was investigated with state-of-the-art mapping devices mounted on IFREMER's (French Research Institute for Exploitation of the Sea) remotely operated vehicle (ROV) Victor 6000. ROV-borne micro-bathymetry and backscatter data of the entire structure, a high-resolution photo-mosaic covering 105,000 m2 of the most active area, sidescan mapping of gas emissions, and maps of faunal distribution as well as of carbonate crust occurrence are combined to provide an unprecedented detailed view of a giant pockmark. All data sets suggest that the pockmark is composed of two very distinctive zones in terms of seepage intensity. We postulate that these zones are the surface expression of two fluid flow regimes in the subsurface: focused flow through a fractured medium and diffuse flow through a porous medium. We conclude that the growth of giant pockmarks is controlled by self-sealing processes and lateral spreading of rising fluids. In particular, partial redirection of fluids through fractures in the sediments can drive the pockmark growth in preferential directions.