Iron isotope fractionation in subterranean estuaries

Dissolved Fe concentrations in subterranean estuaries, like their river-seawater counterparts, are strongly controlled by non-conservative behavior during mixing of groundwater and seawater in coastal aquifers. Previous studies at a subterranean estuary of Waquoit Bay on Cape Cod, USA demonstrate extensive precipitation of groundwater-borne dissolved ferrous iron and subsequent accumulation of iron oxides onto subsurface sands. Waquoit Bay is thus an excellent natural laboratory to assess the mechanisms of Fe-isotope fractionation in redox-stratified environments and determine potential Fe-isotope signatures of groundwater sources to coastal seawater. Here, we report Fe isotope compositions of iron-coated sands and porewaters beneath the intertidal zone of Waquoit Bay. The distribution of pore water Fe shows two distinct sources of Fe: one residing in the upward rising plume of Fe-rich groundwater and the second in the salt-wedge zone of pore water. The groundwater source has high Fe(II) concentration consistent with anoxic conditions and yield δ⁵⁶Fe values between 0.3 and −1.3‰. In contrast, sediment porewaters occurring in the mixing zone of the subterranean estuary have very low δ⁵⁶Fe values down to −5‰. These low δ⁵⁶Fe values reflect Fe-redox cycling and result from the preferential retention of heavy Fe-isotopes onto newly formed Fe-oxyhydroxides. Analysis of Fe-oxides precipitated onto subsurface sands in two cores from the subterranean estuary revealed strong δ⁵⁶Fe and Fe concentration gradients over less than 2m, yielding an overall range of δ⁵⁶Fe values between −2 and 1.5‰. The relationship between Fe concentration and δ⁵⁶Fe of Fe-rich sands can be modeled by the progressive precipitation of Fe-oxides along fluid flow through the subterranean estuary. These results demonstrate that large-scale Fe isotope fractionation (up to 5‰) can occur in subterranean estuaries, which could lead to coastal seawater characterized by very low δ⁵⁶Fe values relative to river values.     

Anaerobic oxidation of methane (AOM) in marine sediments from the Skagerrak (Denmark): II. Reaction-transport modeling

A steady-state reaction-transport model is applied to sediments retrieved by gravity core from two stations (S10 and S13) in the Skagerrak to determine the main kinetic and thermodynamic controls on anaerobic oxidation of methane (AOM). The model considers an extended biomass-implicit reaction network for organic carbon degradation, which includes extracellular hydrolysis of macromolecular organic matter, fermentation, sulfate reduction, methanogenesis, AOM, acetogenesis and acetotrophy. Catabolic reaction rates are determined using a modified Monod rate expression that explicitly accounts for limitation by the in situ catabolic energy yields. The fraction of total sulfate reduction due to AOM in the sulfate-methane transition zone (SMTZ) at each site is calculated. The model provides an explanation for the methane tailing phenomenon which is observed here and in other marine sediments, whereby methane diffuses up from the SMTZ to the top of the core without being consumed. The tailing is due to bioenergetic limitation of AOM in the sulfate reduction zone, because the methane concentration is too low to engender favorable thermodynamic drive. AOM is also bioenergetically inhibited below the SMTZ at both sites because of high hydrogen concentrations (∼3-6 nM). The model results imply there is no straightforward relationship between pore water concentrations and the minimum catabolic energy needed to support life because of the highly coupled nature of the reaction network. Best model fits are obtained with a minimum energy for AOM of ∼11 kJ mol⁻¹, which is within the range reported in the literature for anaerobic processes.     

Seasonal δS variations in two high elevation snow pits measured by S-S double spike thermal ionization mass spectrometry

δ³⁴S and sulfate concentrations were determined in snow pit samples using a thermal ionization mass spectrometric technique capable of 0.2‰ accuracy and requires ≈5 μg (0.16 μmol) natural S. The technique utilizes a ³³S-³⁶S double spike for instrumental mass fractionation correction, and has been applied to snow pit samples collected from the Inilchek Glacier, Kyrgyzstan and from Summit, Greenland. These δ³⁴S determinations provide the first high-resolution seasonal data for these sites, and are used to estimate seasonal sulfate sources. Deuterium (δD) and oxygen (δ¹⁸O) isotope data show that the Inilchek and Summit snow pit samples represent precipitation over ≈20 months.The δ³⁴S values for the Inilchek ranged from +2.6 ± 0.4‰ to +7.6 ± 0.4‰ on sample sizes ranging from 0.3 to 1.8 μmol S. δ³⁴S values for Greenland ranged from +3.6 ± 0.7‰ to +13.3 ± 5‰ for sample sizes ranging from 0.05 to 0.29 μmol S. The concentration ranged from 92.6 ± 0.4 to 1049 ± 4 ng/g for the Inilchek and 18 ± 9 to 93 ± 6 ng/g for the Greenland snow pit. Anthropogenic sulfate dominates throughout the sampled time interval for both sites based on mass balance considerations. Additionally, both sites exhibit a seasonal signature in both δ³⁴S and concentration. The thermal ionization mass spectrometric technique has three advantages compared to gas source isotopic methods: (1) sample size requirements of this technique are 10-fold less permitting access to the higher resolution S isotope record of low concentration snow and ice, (2) the double spike technique permits δ³⁴S and S concentration to be determined simultaneously, and (3) the double spike is an internal standard.     

The Fe-C system at 5 GPa and implications for Earth’s core

Earth’s core may contain C, and it has been suggested that C in the core could stabilize the formation of a solid inner core composed of Fe₃C. We experimentally examined the Fe-C system at a pressure of 5 GPa and determined the Fe-C phase diagram at this pressure. In addition, we measured solid metal/liquid metal partition coefficients for 17 trace elements and examined the partitioning behavior between Fe₃C and liquid metal for 14 trace elements. Solid metal/liquid metal partition coefficients are similar to those found in one atmosphere studies, indicating that the effect of pressure to 5 GPa is negligible. All measured Fe₃C/liquid metal partition coefficients investigated are less than one, such that all trace elements prefer the C-rich liquid to Fe₃C. Fe₃C/liquid metal partition coefficients tend to decrease with decreasing atomic radii within a given period. Of particular interest, our 5 GPa Fe-C phase diagram does not show any evidence that the Fe-Fe₃C eutectic composition shifts to lower C contents with increasing pressure, which is central to the previous reasoning that the inner core may be composed of Fe₃C.     

Rare earth elements (REE) and yttrium in stream waters, stream sediments, and Fe-Mn oxyhydroxides: Fractionation, speciation, and controls over REE + Y patterns in the surface environment

We have collected ∼500 stream waters and associated bed-load sediments over an ∼400 km² region of Eastern Canada and analyzed these samples for Fe, Mn, and the rare earth elements (REE + Y). In addition to analyzing the stream sediments by total digestion (multi-acid dissolution with metaborate fusion), we also leached the sediments with 0.25 M hydroxylamine hydrochloride (in 0.05 M HCl), to determine the REE + Y associated with amorphous Fe- and Mn-oxyhydroxide phases. We are thus able to partition the REE into “dissolved” (<0.45 μm), labile (hydroxylamine) and detrital sediment fractions to investigate REE fractionation, and in particular, with respect to the development of Ce and Eu anomalies in oxygenated surface environments. Surface waters are typically LREE depleted ([La/Sm]_NASC ranges from 0.16 to 5.84, average = 0.604, n = 410; where the REE are normalized to the North America Shale Composite), have strongly negative Ce anomalies ([Ce/Ce^∗]_NASC ranges from 0.02 to 1.25, average = 0.277, n = 354), and commonly have positive Eu anomalies ([Eu/Eu^∗]_NASC ranges from 0.295 to 1.77, average = 0.764, n = 84). In contrast, the total sediment have flatter REE + Y patterns relative to NASC ([La/Sm]_NASC ranges from 0.352 to 1.12, average = 0.778, n = 451) and are slightly middle REE enriched ([Gd/Yb]_NASC ranges from 0.55 to 3.75, average = 1.42). Most total sediments have negative Ce and Eu anomalies ([Ce/Ce^∗]_NASC ranges from 0.097 to 2.12, average = 0.799 and [Eu/Eu^∗]_NASC ranges from 0.39 to 1.43, average = 0.802). The partial extraction sediments are commonly less LREE depleted than the total sediments ([La/Sm]_NASC ranges from 0.24 to 3.31, average = 0.901, n = 4537), more MREE enriched ([Gd/Yb]_NASC ranges from 0.765 to 6.28, average = 1.97) and Ce and Eu anomalies (negative and positive) are more pronounced.The partial extraction recovered, on average ∼20% of the Fe in the total sediment, ∼80% of the Mn, and 21-29% of the REEs (Ce = 19% and Y = 32%). Comparison between REEs in water, partial extraction and total sediment analyses indicates that REEs + Y in the stream sediments have two primary sources, the host lithologies (i.e., mechanical dispersion) and hydromorphically transported (the labile fraction). Furthermore, Eu appears to be more mobile than the other REE, whereas Ce is preferentially removed from solution and accumulates in the stream sediments in a less labile form than the other REEs + Y. Despite poor statistical correlations between the REEs + Y and Mn in either the total sediment or partial extractions, based on apparent distribution coefficients and the pH of the stream waters, we suggest that either sediment organic matter and/or possibly δ-MnO₂/FeOOH are likely the predominant sinks for Ce, and to a lesser extent the other REE, in the stream sediments.     

Argon and nitrogen isotopes of trapped air in the GISP2 ice core during the Holocene epoch (0-11,500 B.P.): Methodology and implications for gas loss processes

Argon and nitrogen isotopes of air in polar ice cores provide constraints on past temperature and firn thickness, with relevance to past climate. We developed a method to simultaneously measure nitrogen and argon isotopes in trapped air from the same sample of polar ice. This method reduces the time required for analysis, allowing large numbers of measurements. We applied this method to the entire Holocene sequence of the GISP2 ice core (82.37-1692.22 m) with a 10-20 year sampling interval (670 depths). δ⁴⁰Ar and δ¹⁵N show elevated values in the oldest part of the dataset, consistent with a thicker firn layer and increased temperature gradient in the firn due to the legacy of the abrupt warming at the end of the Younger Dryas interval and the gradual warming during the Preboreal interval (11.5-10.0 ka). The Preboreal Oscillation and the 8.2k event are clearly recorded. The data show remarkable stability after the 8.2k event.Available data suggests that post-coring gas loss involves two distinct types of fractionation. First, smaller molecules with less than a certain threshold size leak through the ice lattice with little isotopic fractionation. Second, gas composition changes via gas loss through microcracks, which induces isotopic fractionation. These two gas loss processes can explain most trends in our data and in other ice core records.     

Thermodynamic model for arsenic speciation in sulfidic waters: A novel use of ab initio computations

Recent discoveries demonstrate that the chemistry of arsenic in sulfidic waters is much more complex that previously believed. One implication is that all earlier thermodynamic data on stabilities of As thioanions require revision. Previously used experimental approaches for determining As thioanion stabilities may be inadequate to deal with the full range of complexity. Here we use computational as well as empirical information to construct a provisional model for equilibrium As thioanion distributions in sulfidic waters. Whereas previous authors have argued for either As(III) or As(V) thioanions, the new model predicts that both are important and can occur simultaneously under commonly encountered pH and ΣS^−II conditions. At the order of magnitude level, the model reasonably predicts the solubility of As₂S₃ in sulfidic solutions, provides tentative peak assignments for published Raman spectroscopic data and plausibly accounts for how sulfide modifies the bacterial toxicity of As. The model yields a thermodynamic justification for how sulfide, which is usually regarded as a reducing agent, can counter-intuitively drive oxidation of As(III) to As(V), as has been observed both in the laboratory and in the field. Despite its uncertain accuracy, the model serves as a useful source of new, testable hypotheses about As geochemistry and highlights crucial experimental data needs.     

Calcite precipitation instability under laminar, open-channel flow

We present a 2D numerical model for the growth of calcite from supersaturated aqueous solutions under laminar, open-channel flow conditions. The model couples solution chemistry, precipitation at solution/calcite interfaces, hydrodynamics, diffusion and degassing. The model output is compared with experimental results obtained using an oversaturated calcite solution produced by mixing CaCl₂ and Na₂CO₃. The precipitation rate is observed to increase when the supersaturated solution flows over an obstruction, leading to a growth instability that causes the formation of terraces. At relatively high flow rates, the most important mechanism for this behaviour seems to be hydrodynamic advection of dissolved species either towards or away from the calcite surface, depending on location relative to the obstruction, which deforms the concentration gradients. At lower flow rates, steepening of diffusion gradients around protrusions becomes important. Enhanced degassing over the obstruction due to shallowing and pressure drop is not important on small scales. Diffusion controlled transport close to the calcite surface can lead to a fingering-type growth instability, which generates porous textures. Our results are consistent with existing diffusive boundary layer theory, but for flow over non-smooth surfaces, simple calcite precipitation models that include empirical correlations between fluid flow rate and calcite precipitation rate are inaccurate.     

Models in geography revisited

Agent based models (ABMs) have many applications, and illustrate a rapidly developing field of enquiry, spanning both the physical-mathematical and human-social sciences. ABMs are seen as most appropriate in situations where decisions or actions are distributed around particular locations, and in which structure is viewed as emergent from the interaction between individuals. ABMs may be used either as representational devices, to reproduce the patterns observed or desired in the system of interest, or as foundational tools contributing to the development of social or economic theory. The role and status of models and modelling is itself an instantiation of a wider debate concerning representation and explanation. Today, a case can be made that the nature of explanation and the use of scientific interpretation reflect much less definite and exclusive positions and permit more diverse approaches than hitherto. The underlying proposition of this commentary is, therefore, that the time is right for a positive application of ABMs within the discipline of geography, and for a rediscovery and reappraisal of the richness and depth of insight in the model-building enterprise more generally. First, the context for ABM development and application is set with reference to the agency-structure debate. Second, some aspects of the heritage of models in geography is presented, based upon reviews of two benchmark publications bearing that title. Next, some of the most significant characteristics, uses, potentials and limitations of ABMs, are reviewed. Finally, some constructive ways forward are suggested, as informed by theory and method from the interpretative social sciences.     

Models of natural and human dynamics in forest landscapes: Cross-site and cross-cultural synthesis

We synthesize the study of coupled natural and human systems across sites and cultures through a process of simplification and abstraction based on multiple dimensions of human-nature connectedness: satisfaction of basic needs, psycho-cultural connectedness and regulation of use of natural resources. We thus provide both a place-based and general understanding of value-driven anthropogenic environmental change and response. Two questions guide this research: what are the crucial stakeholder values that drive land use decisions and thus land cover change? And how can knowledge of these values be used to make decisions and policies that sustain both the human and natural systems in a place? To explore these questions we build simulation models of four study sites, two in the State of Texas, United States, and two in Venezuela. All include protected areas, though they differ in the specifics of vegetation and land use. In the Texas sites, relatively affluent individuals are legally converting forests to residential, commercial, and industrial uses, while in Venezuela landless settlers are extra-legally converting forests for purposes of subsistence agriculture. Contemporary modeling techniques now facilitate simulations of stakeholder and ecosystem dynamics revealing emergent patterns. Such coupled human and natural systems are currently recognized as a form of biocomplexity. Our modeling framework is flexible enough to allow adaptation to each of the study sites, capturing the essential features of the respective natural and anthropogenic land use changes and stakeholder reactions. The interactions between human stakeholders are simulated using multi-agent models that act on forest landscape models, and receive feedback of the effects of these actions on ecological habitats and hydrological response. The multi-agent models employ a formal logic-based method for the Venezuelan sites and a decision analysis approach using multi-attribute utility functions for the Texas sites, differing more in style and emphasis than in substance. Our natural-systems models are generic and can be tailored according to site-specific conditions. Similar models of tree growth and patch transitions are used for all the study sites and the differing responses to environmental variables are specified for each local species and terrain conditions.