Spatial heterogeneity and lake morphology affect diffusive greenhouse gas emission estimates of lakes

Most estimates of diffusive flux (F) of methane (CH4) and carbon dioxide (CO2) from lakes are based on single-point flux chamber measurements or on piston velocity (k) modeled from wind speed and single-point measurements of surface water gas concentrations (C-aq). We analyzed spatial variability of F of CH4 and CO2, as well as C-aq and k in 22 European lakes during late summer. F and k were higher in the lake centers, leading to considerable bias when extrapolating single-point chamber measurements to whole-lake estimates. The ratio of our empirical k estimates to wind speed-modeled k was related to lake size and shape, suggesting a lake morphology effect on the relationship between wind speed and k. This indicates that the error inherent to established wind speed models can be reduced by determining k and C-aq at multiple sites on lakes to calibrate wind speed-modeled k to the local system.

Geophysical and Hydrological Evaluation of Two Bog Complexes in a Northern Peatland: Implications for the Distribution of Biogenic Gases at the Basin Scale

Ground penetrating radar (GPR) was used to determine peat basin geometry and the spatial distribution of free-phase biogenic gasses in two separate units of a northern peatland (Central and Southern Unit of Caribou Bog, Maine). The Central Unit is characterized by a deep basin structure (15 m maximum depth) and a raised (eccentric) bog topographic profile (up to 2 m topographic variation). Here numerous regions of electromagnetic (EM) wave scattering are considered diagnostic of the presence of extensive free-phase biogenic gas. In contrast, the Southern Unit is shallower (8 m maximum depth) and has a slightly convex upwards bog profile (less than 1 m topographic variation), and areas of EM wave scattering are notably absent. The biogenic gas zones interpreted from GPR in the Central Unit are associated with: (1) wooded heath vegetation at the surface, (2) open pools at the surface, (3) high water table elevations near the center of the basin, and (4) a region of overpressure (at approximately 5 m depth) immediately below the zone of free-phase gas accumulation. The latter suggests (1) a transient pressure head associated with low hydraulic conductivity resulting from the biogenic gasses themselves or confining layers in the peat that restrict both gas release and groundwater flow and/or (2) overpressure in the peat column as a result of the gas buildup itself. In contrast, the Southern Unit, where zones of EM scattering are absent, is characterized by: (1) predominantly shrub vegetation, (2) a lack of open pools, (3) only minor variations (less than 1 m) in water table elevation throughout the entire unit; and (4) generally upward groundwater flow throughout the basin. The results illustrate the nonuniformity of free-phase biogenic gas distribution at the peat basin scale and provide insights into the processes and controls associated with CH4 and CO2 accumulation in peatlands.

From planetesimals to planets: volatile molecules

Context. Solar and extrasolar planets are the subject of numerous studies aiming to determine their chemical composition and internal structure. In the case of extrasolar planets, the composition is important as it partly governs their potential habitability. Moreover, observational determination of chemical composition of planetary atmospheres are becoming available, especially for transiting planets. Aims. The present works aims at determining the chemical composition of planets formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data for models of planet formation and evolution, and future interpretation of chemical composition of solar and extrasolar planets. Methods. We have developed a model that computes the composition of ices in planets in different stellar systems with the use of models of ice and planetary formation. Results. We provide the chemical composition, ice/rock mass ratio and C:O molar ratio for planets in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, we produce a wide variety of planetary chemical compositions as a function of the mass of the disk and distance to the star. The volatile species incorporated in planets are mainly composed of H2O, CO, CO2, CH3OH, and NH3. Icy or ocean planets have systematically higher values of molecular abundances compared to giant and rocky planets. Gas giant planets are depleted in highly volatile molecules such as CH4, CO, and N2 compared to icy or ocean planets. The ice/rock mass ratio in icy or ocean and gas giant planets is, respectively, equal at maximum to 1.01 ± 0.33 and 0.8 ± 0.5, and is different from the usual assumptions made in planet formation models, which suggested this ratio to be 2–3. The C:O molar ratio in the atmosphere of gas giant planets is depleted by at least 30% compared to solar value.

An inter-regional assessment of concentrations and δ13C values of methane and dissolved inorganic carbon in small European lakes

Methane (CH4) and carbon dioxide emissions from lakes are relevant for assessing the greenhouse gas output of wetlands. However, only few standardized datasets describe concentrations of these gases in lakes across different geographical regions. We studied concentrations and stable carbon isotopic composition (δ13C) of CH4 and dissolved inorganic carbon (DIC) in 32 small lakes from Finland, Sweden, Germany, the Netherlands, and Switzerland in late summer. Higher concentrations and δ13C values of DIC were observed in calcareous lakes than in lakes on non-calcareous areas. In stratified lakes, δ13C values of DIC were generally lower in the hypolimnion due to the degradation of organic matter (OM). Unexpectedly, increased δ13C values of DIC were registered above the sediment in several lakes. This may reflect carbonate dissolution in calcareous lakes or methanogenesis in deepwater layers or in the sediments. Surface water CH4 concentrations were generally higher in western and central European lakes than in Fennoscandian lakes, possibly due to higher CH4 production in the littoral sediments and lateral transport, whereas CH4 concentrations in the hypolimnion did not differ significantly between the regions. The δ13C values of CH4 in the sediment suggest that δ13C values of biogenic CH4 are not necessarily linked to δ13C values of sedimentary OM but may be strongly influenced by OM quality and methanogenic pathway. Our study suggests that CH4 and DIC cycling in small lakes differ between geographical regions and that this should be taken into account when regional studies on greenhouse gas emissions are upscaled to inter-regional scales.

Microthermometry and Raman analysis of fluid inclusions in different generations of quartz veins of the Athabasca Basin

The Athabasca Basin (Canada) contains the highest grade unconformity-type uranium deposits in the world. Underlying the Athabasca Group sedimentary rocks of the Dufferin Lake zone are variably graphitic pelitic schists (VGPS), altered to chlorite and hematite (Red/Green Zone: RGZ), and locally bleached near the unconformity during paleoweathering and/or later fluid interaction, leading to a loss of graphite near the unconformity. Fluid inclusions were examined in different generations of quartz veins, using microthermometry and Raman analysis, to characterize and compare the different fluids that interacted with the RGZ and the VGPS. In the VGPS, CH4-, N2- and CO2-rich fluids circulated. CH4- and N2-rich fluids could be the result of the breakdown of graphite to CH4/CO2, whereas N2-rich fluid is interpreted to be the result of breakdown of feldspars/micas to NH4+/N2. In the RGZ, highly saline fluids interpreted to be basinally derived have been recorded. The circulation of the two types of fluids (carbonic and brines) occurred at two different distinct events: 1) during the retrograde metamorphism of the basement rocks before the deposition of the Athabasca Basin for the carbonic fluids, and 2) after the deposition of the Athabasca Basin for the brines. Thus, in addition to possibly be related to graphite depletion in the RGZ, the brines can be linked to uranium mineralization.

Isotopic composition of interstitial fluids and origin of methane at DSDP Leg 84 Holes

CH4 and CO2 species in pore fluids from slope sediments off Guatemala show extreme 13C-enrichment (d13C of -41 and +38 per mil, respectively) compared with the typical degree of 13C-enrichment in pore fluids of DSDP sediments (d13C of - 60 and + 10 per mil). These unusual isotopic compositions are believed to result from microbial decomposition of organic matter, and possibly from additional isotopic fractionation associated with the formation of gas hydrates. In addition to the isotopic fractionation displayed by CH4 and CO2, the pore water exhibits a systematic increase in d18O with decrease in chlorinity. As against seawater d18O values of 0 and chlorinity of 19 per mil, the water collected from decomposed gas hydrate from Hole 570 had a d18O of + 3.0 per mil and chlorinity of 9.5 per mil. The isotopic compositions of pore-fluid constituents change gradually with depth in Hole 568 and discontinuously with depth in Hole 570.

Simulated change in atmospheric CO2 and D14C around 14.6 kyr BP, at the onset of the northern hemispheric warming into the Bølling/Allerød

One of the most abrupt and yet unexplained past rises in atmospheric CO2 (10 p.p.m.v. in two centuries) occurred in quasi-synchrony with abrupt northern hemispheric warming into the Bølling/Allerød, 14,600 years ago. Here we use a U/Th-dated record of atmospheric D14C from Tahiti corals to provide an independent and precise age control for this CO2 rise. We also use model simulations to show that the release of old (nearly 14C-free) carbon can explain these changes in CO2 and D14C. The D14C record provides an independent constraint on the amount of carbon released (125 Pg C). We suggest, in line with observations of atmospheric CH4 and terrigenous biomarkers, that thawing permafrost in high northern latitudes could have been the source of carbon, possibly with contribution from flooding of the Siberian continental shelf during meltwater pulse 1A. Our findings highlight the potential of the permafrost carbon reservoir to modulate abrupt climate changes via greenhouse-gas feedbacks.

Carbon isotopic composition and dissolved CO2 concentration in pore waters of ODP Leg 174A sites

Microbially mediated redox diagenetic processes in marine sediments are reflected in the amount and carbon isotopic composition of dissolved CO2 and CH4 (Claypool and Kaplan, 1974). Oxidation of organic matter gives rise to dissolved CO2 with about the same 13C/12C ratio as the starting organic matter. Subsequent reduction of CO2 to form CH4 involves a large (~70) kinetic isotopic effect, resulting in significant 13C depletion in the CH4, and 13C enrichment in the residual CO2. Ocean Drilling Program Leg 174A (offshore New Jersey) presented an opportunity to study these processes in shelf and upper slope sediments. Holes 1071A-1071D, 1071F, and 1072A were drilled on the shelf in water depths of 88.0-98.1 m. Hole 1073A was drilled on the slope in 639.4 m of water. Pore-water samples were collected for analysis at all three sites, whereas gas samples could only be obtained from Hole 1073A on the slope.

Isotopic ratios of methane, ethane, carbon dioxide, inorganic carbon and organic matter from ODP Leg 164 sites

The isotopic characteristics of CH4 (d13C values range from -101.3 per mil to -61.1 per mil PDB, and dD values range from -256 per mil to -136 per mil SMOW) collected during Ocean Drilling Program (ODP) Leg 164 indicate that the CH4 was produced by microbial CO2 reduction and that there is not a significant contribution of thermogenic CH4 to the sampled sediment gas from the Blake Ridge. The isotopic values of CO2 (d13C range -20.6 per mil to +1.24 per mil PDB) and dissolved inorganic carbon (DIC; d13C range -37.7 per mil to +10.8 per mil PDB) have parallel profiles with depth, but with an offset of 12.5 per mil. Distinct downhole variations in the carbon isotopic composition of CH4 and CO2 cannot be explained by closed-system fractionation where the CO2 is solely derived from the locally available sedimentary organic matter (d13C -2.0 per mil ± 1.4 per mil PDB) and the CH4 is derived from CO2 reduction. The observed isotopic profiles reflect the combined effects of upwards gas migration and decreased microbial activity with depth.

Analysis and compound of carbon sequestration and serpentinization from the Iberian Margin and the Northern Apennine

Fluid circulation in peridotite-hosted hydrothermal systems influences the incorporation of carbon into the oceanic crust and its long-term storage. At low to moderate temperatures, serpentinization of peridotite produces alkaline fluids that are rich in CH4 and H2. Upon mixing with seawater, these fluids precipitate carbonate, forming an extensive network of calcite veins in the basement rocks, while H2 and CH4 serve as an energy source for microorganisms. Here, we analyzed the carbon geochemistry of two ancient peridotite-hosted hydrothermal systems: 1) ophiolites cropping out in the Northern Apennines, and 2) calcite-veined serpentinites from the Iberian Margin (Ocean Drilling Program (ODP) Legs 149 and 173), and compare them to active peridotite-hosted hydrothermal systems such as the Lost City hydrothermal field (LCHF) on the Atlantis Massif near the Mid-Atlantic Ridge (MAR). Our results show that large amounts of carbonate are formed during serpentinization of mantle rocks exposed on the seafloor (up to 9.6 wt.% C in ophicalcites) and that carbon incorporation decreases with depth. In the Northern Apennine serpentinites, serpentinization temperatures decrease from 240 °C to < 150 °C, while carbonates are formed at temperatures decreasing from ~ 150 °C to < 50 °C. At the Iberian Margin both carbonate formation and serpentinization temperatures are lower than in the Northern Apennines with serpentinization starting at ~ 150 °C, followed by clay alteration at < 100 °C and carbonate formation at < 19-44 °C. Comparison with various active peridotite-hosted hydrothermal systems on the MAR shows that the serpentinites from the Northern Apennines record a thermal evolution similar to that of the basement of the LCHF and that tectonic activity on the Jurassic seafloor, comparable to the present-day processes leading to oceanic core complexes, probably led to formation of fractures and faults, which promoted fluid circulation to greater depth and cooling of the mantle rocks. Thus, our study provides further evidence that the Northern Apennine serpentinites host a paleo-stockwork of a hydrothermal system similar to the basement of the LCHF. Furthermore, we argue that the extent of carbonate uptake is mainly controlled by the presence of fluid pathways. Low serpentinization temperatures promote microbial activity, which leads to enhanced biomass formation and the storage of organic carbon. Organic carbon becomes dominant with increasing depth and is the principal carbon phase at more than 50-100 m depth of the serpentinite basement at the Iberian Margin. We estimate that annually 1.1 to 2.7 × 1012 g C is stored within peridotites exposed to seawater, of which 30-40% is fixed within the uppermost 20-50 m mainly as carbonate. Additionally, we conclude that alteration of oceanic lithosphere is an important factor in the long-term global carbon cycle, having the potential to store carbon for millions of years.

Geochemistry and isotopes at DSDP Site 76-533

Concentrations and d34S and d13C values were determined on SO4, HCO3, CO2, and CH4 in interstitial water and gas samples from the uppermost 400 m of sediment on the Blake Outer Ridge. These measurements provide the basis for detailed interpretation of diagenetic processes associated with anaerobic respiration of electrons generated by organic- matter decomposition. The sediments are anaerobic at very shallow depths (<1 m) below the seafloor. Sulfate reduction is confined to the uppermost 15 m of sediment and results in a significant outflux of oxidized carbon from the sediments. At the base of the sulfate reduction zone, upward-diffusing CH4 is being oxidized, apparently in conjunction with SO4 reduction. CH4 generation by CO2 reduction is the most important metabolic process below the 15-m depth. CO2 removal is more rapid than CO2 input over the depth interval from 15 to 100 m, and results in a slight decrease in HCO3 concentration accompanied by a 40 per mil positive shift in d13C. The differences among coexisting CH4, CO2, and HCO3 are consistent with kinetic fractionation between CH4 and dissolved CO2, and equilibrium fractionation between CO2 and HCO3. At depths greater than 100 m, the rate of input of CO2 (d13C = -25 per mil) exceeds by 2 times the rate of removal of CO2 by conversion to CH4 (d13C of -60 to -65 per mil). This results in an increase of dissolved HCO3 concentration while maintaining d13C of HCO3 relatively constant at +10 per mil. Non-steady-state deposition has resulted in significantly higher organic carbon contents and unusually high (70 meq/l) pore-water alkalinities below 150 m. These high alkalinities are believed to be related more to spontaneous decarboxylation reactions than to biological processes. The general decrease in HCO3 concentration with constant d13C over the depth interval of 200 to 400 m probably reflects increased precipitation of authigenic carbonate. Input-output carbon isotope-mass balance calculations, and carbonate system equilibria in conjunction with observed CO2-CH4 ratios in the gas phase, independently suggest that CH4 concentrations on the order of 100 mmol/kg are present in the pore waters of Blake Outer Ridge sediments. This quantity of CH4 is believed to be insufficient to saturate pore waters and stabilize the CH4*6H2O gas hydrate. Results of these calculations are in conflict with the physical recovery of gas hydrate from 238 m, and with the indirect evidence (seismic reflectors, sediment frothing, slightly decreasing salinity and chlorinity with depth, and pressure core barrel observations) of gas-hydrate occurrence in these sediments. Resolution of this apparent conflict would be possible if CH4 generation were restricted to relatively thin (1-10 m) depth intervals, and did not occur uniformly at all depths throughout the sediment column, or if another methanogenic process (e.g., acetate fermentation) were a major contributor of gas.

Composition of sediment pore water and isotopes in carbonate and benthic foraminifera at DSDP Holes 68-501, 69-505 and Site 69-504

Two sites on the southern flank of the Costa Rica Rift were drilled on DSDP Legs 68 and 69, one on crust 3.9 m.y. old and the other on crust 5.9 m.y. old. The basement of the younger site is effectively cooled by the circulation of seawater. The basement of the older site has been sealed by sediment, and an interval in the uppermost 560 meters of basement recently reheated to temperatures of 60 to 120°C. Although the thickness of the sediments at the two sites is similar (150-240 m versus 270 m), the much rougher basement topography at the younger Site 505 produces occasional basement outcrops, through which 80 to 90% of the total heat loss apparently occurs by advection of warm seawater. This seawater has been heated only slightly, however; the temperature at the base of the sediments is only 9°C. Changes in its composition due to reaction with the basement basalts are negligible, as indicated by profiles of sediment pore water chemistry. Bacterial sulfate reduction in the sediments produces a decrease in SO4 (and Ca) and an increase in alkalinity (and Sr and NH3) as depth increases to an intermediate level, but at deeper levels these trends reverse, and all of these species plus Mg, K, Na, and chlorinity approach seawater values near basement. Si, however, is higher, and Li may be lower. At the older site, Site 501/504, where heat loss is entirely by conduction, the temperature at the sediment/basement contact is 59°C. Sediment pore water chemistry is heavily affected by reaction with the basaltic basement, as indicated by large decreases in d18O, Mg, alkalinity, Na, and K and an increase in Ca with increasing depth. The size of the changes in d18O, Mg, alkalinity, Ca, Sr, and SO4 varies laterally over 500 meters, indicating lateral gradients in pore water chemistry that are nearly as large as the vertical gradients. The lateral gradients are believed to result from similar lateral gradients in the composition of the basement formation water, which propagate upward through the sediments by diffusion. A model of the d18O profile suggests that the basement at Site 501/504 was sealed off from advection about 1 m.y. ago, so that reaction rates began to dominate the basement pore water chemistry. A limestone-chert diagenetic front began to move upward through the lower sediments less than 200,000 yr. ago.

Stable carbon isotope ratios and the formation of dolomite in ODP Leg 175 sites

We examine the link between organic matter degradation, anaerobic methane oxidation (AMO), and sulfate depletion and explore how these processes potentially influence dolomitization. We determined rates and depths of AMO and dolomite formation for a variety of organic-rich sites along the west African Margin using data from Ocean Drilling Program (ODP) Leg 175. Rates of AMO are calculated from the diffusive fluxes of CH4 and SO4, and rates of dolomite formation are calculated from the diffusive flux of Mg. We find that the rates of dolomite formation are relatively constant regardless of the depth at which it is forming, indicating that the diffusive fluxes of Mg and Ca are not limiting. Based upon the calculated log IAP values, log K(sp) values for dolomite were found to narrowly range between -16.1 and -16.4. Dolomite formation is controlled in part by competition between AMO and methanogenesis, which controls the speciation of dissolved CO2. AMO increases the concentration of CO3[2-] through sulfate reduction, favoring dolomite formation, while methanogenesis increases the pCO2 of the pore waters, inhibiting dolomite formation. By regulating the pCO2 and alkalinity, methanogenesis and AMO can regulate the formation of dolomite in organic-rich marine sediments. In addition to providing a mechanistic link between AMO and dolomite formation, our findings provide a method by which the stability constant of dolomite can be calculated in modern sediments and allow prediction of regions and depth domains in which dolomite may be forming.

Isotopic composition of ground ice, ebullition gases and thermokarst lake water, North America

Thermokarst lakes are thought to have been an important source of methane (CH4) during the last deglaciation when atmospheric CH4 concentrations increased rapidly. Here we demonstrate that meltwater from permafrost ice serves as an H source to CH4 production in thermokarst lakes, allowing for region-specific reconstructions of dD-CH4 emissions from Siberian and North American lakes. dD CH4 reflects regionally varying dD values of precipitation incorporated into ground ice at the time of its formation. Late Pleistocene-aged permafrost ground ice was the dominant H source to CH4 production in primary thermokarst lakes, whereas Holocene-aged permafrost ground ice contributed H to CH4 production in later generation lakes. We found that Alaskan thermokarst lake dD-CH4 was higher (-334 ± 17 per mil) than Siberian lake dD-CH4 (-381 ± 18 per mil). Weighted mean dD CH4 values for Beringian lakes ranged from -385 per mil to -382 per mil over the deglacial period. Bottom-up estimates suggest that Beringian thermokarst lakes contributed 15 ± 4 Tg CH4 /yr to the atmosphere during the Younger Dryas and 25 ± 5 Tg CH4 /yr during the Preboreal period. These estimates are supported by independent, top-down isotope mass balance calculations based on ice core dD-CH4 and d13C-CH4 records. Both approaches suggest that thermokarst lakes and boreal wetlands together were important sources of deglacial CH4.

(Supplemental Table 2) Plant characteristics, GHG fluxes, and soil chemical and microbial properties in experimental treatments in Alexandra Fiord

We evaluated above- and belowground ecosystem changes in a 16 year, combined fertilization and warming experiment in a High Arctic tundra deciduous shrub heath (Alexandra Fiord, Ellesmere Island, NU, Canada). Soil emissions of the three key greenhouse gases (GHGs) (carbon dioxide, methane, and nitrous oxide) were measured in mid-July 2009 using soil respiration chambers attached to a FTIR system. Soil chemical and biochemical properties including Q10 values for CO2, CH4, and N2O, Bacteria and Archaea assemblage composition, and the diversity and prevalence of key nitrogen cycling genes including bacterial amoA, crenarchaeal amoA, and nosZ were measured. Warming and fertilization caused strong increases in plant community cover and height but had limited effects on GHG fluxes and no substantial effect on soil chemistry or biochemistry. Similarly, there was a surprising lack of directional shifts in the soil microbial community as a whole or any change at all in microbial functional groups associated with CH4 consumption or N2O cycling in any treatment. Thus, it appears that while warming and increased nutrient availability have strongly affected the plant community over the last 16 years, the belowground ecosystem has not yet responded. This resistance of the soil ecosystem has resulted in limited changes in GHG fluxes in response to the experimental treatments.