An absolute paleotemperature record from 10 to 6 Ka inferred from fluid inclusion D/H ratios of a stalagmite from Vancouver Island, British Columbia, Canada

By using continuous helium flow during the crushing of calcite speleothem samples, we are able to recover liberated inclusion waters without isotopic fractionation. A paleotemperature record for the Jacklah Jill Cave locality, Vancouver Island, BC, was obtained from a 30-cm tall stalagmite that grew 10.3-6.3 Ka ago, using δ¹⁸O values of the crushed calcite and of the inclusion water as inferred from its δD. It is found that the locality experienced mean annual temperature variations up to 11 °C over a 4-Ka period in the early Holocene. At the beginning of the period, local temperature quickly increased from a minimum of ∼1 °C to around 10 °C, but this early climate optimum, about 3 °C warmer than today, only lasted for ∼1200 years. About 8.6 Ka ago, temperature had declined to ∼7 °C, approximately the same as the modern cave temperature. Since then, the study area has experienced only minor temperature fluctuations, but there was a brief fall to ∼4 °C at around 7 Ka ago, which might be caused by a short lived expansion of local alpine glaciers. The long-term T-dependence of δD was 1.47‰/°C, identical to the value in modern precipitation.     

Silicon isotope fractionation in bamboo and its significance to the biogeochemical cycle of silicon

A systematic investigation on silica contents and silicon isotope compositions of bamboos was undertaken. Seven bamboo plants and related soils were collected from seven locations in China. The roots, stem, branch and leaves for each plant were sampled and their silica contents and silicon isotope compositions were determined. The silica contents and silicon isotope compositions of bulk and water-soluble fraction of soils were also measured. The silica contents of studied bamboo organs vary from 0.30% to 9.95%. Within bamboo plant the silica contents show an increasing trend from stem, through branch, to leaves. In bamboo roots the silica is exclusively in the endodermis cells, but in stem, branch and leaves, the silica is accumulated mainly in epidermal cells. The silicon isotope compositions of bamboos exhibit significant variation, from −2.3‰ to 1.8‰, and large and systematic silicon isotope fractionation was observed within each bamboo. The δ³⁰Si values decrease from roots to stem, but then increase from stem, through branch, to leaves. The ranges of δ³⁰Si values within each bamboo vary from 1.0‰ to 3.3‰. Considering the total range of silicon isotope composition in terrestrial samples is only 7‰, the observed silicon isotope variation in single bamboo is significant and remarkable. This kind of silicon isotope variation might be caused by isotope fractionation in a Rayleigh process when SiO₂ precipitated in stem, branches and leaves gradually from plant fluid. In this process the Si isotope fractionation factor between dissolved Si and precipitated Si in bamboo (α_pre-sol) is estimated to be 0.9981. However, other factors should be considered to explain the decrease of δ³⁰Si value from roots to stem, including larger ratio of dissolved H₄SiO₄ to precipitated SiO₂ in roots than in stem. There is a positive correlation between the δ³⁰Si values of water-soluble fractions in soils and those of bulk bamboos, indicating that the dissolved silicon in pore water and phytoliths in soil is the direct sources of silicon taken up by bamboo roots. A biochemical silicon isotope fractionation exists in process of silicon uptake by bamboo roots. Its silicon isotope fractionation factor (α_bam-wa) is estimated to be 0.9988. Considering the distribution patterns of SiO₂ contents and δ³⁰Si values among different bamboo organs, evapotranspiration may be the driving force for an upward flow of a silicon-bearing fluid and silica precipitation. Passive silicon uptake and transportation may be important for bamboo, although the role of active uptake of silicic acid by roots may not be neglected. The samples with relatively high δ³⁰Si values all grew in soils showing high content of organic materials. In contrast, the samples with relatively low δ³⁰Si values all grew in soil showing low content of organic materials. The silicon isotope composition of bamboo may reflect the local soil type and growth conditions. Our study suggests that bamboos may play an important role in global silicon cycle.     

Atomistic simulation of the mechanisms of noble gas incorporation in minerals

Atomistic simulations have been carried out to investigate the mechanisms of noble gas incorporation in minerals using both the traditional two-region approach and the “supercell” method. The traditional two-region approach has been used to calculate defect energies for Ne, Ar, Kr and Xe incorporation in MgO, CaO, diopside and forsterite in the static limit and at one atmosphere pressure. The possibilities of noble gas incorporation via both substitution and interstitial mechanisms are studied. The favored mechanism varies from mineral to mineral and from noble gas to noble gas. In all minerals studied, the variation of the solution energies of noble gas substitution with atomic radius appears approximately parabolic, analogous to those for 1⁺, 2⁺, 3⁺ and 4⁺ trace element incorporation on crystal lattice sites. Noble gas solution energies thus also fall on a curve, similar to those previously observed for cations with different charges, but with much lower curvature.The “supercell” method has been used to investigate the pressure dependence of noble gas incorporation in the same systems. Results indicate a large variation of the solubility of the larger noble gases, Kr and Xe with pressure. In addition, explicit simulation of incorporation at the (0 0 1) surface of MgO shows that the solubility of the heavier noble gases may be considerably enhanced by the presence of interfaces.     

Low-oxygen and chemical kinetic constraints on the geochemical niche of neutrophilic iron(II) oxidizing microorganisms

Neutrophilic iron oxidizing bacteria (FeOB) must actively compete with rapid abiotic processes governing Fe(II) oxidation and as a result have adapted to primarily inhabit low-O₂ environments where they can more successfully compete with abiotic Fe(II) oxidation. The spatial distribution of these microorganisms can be observed through the chemical gradients they affect, as measured using in situ voltammetric analysis for dissolved Fe(II), Fe(III), O₂, and FeS_(aq). Field and laboratory determination of the chemical environments inhabited by the FeOB were coupled with detailed kinetic competition studies for abiotic and biotic oxidation processes using a pure culture of FeOB to quantify the geochemical niche these organisms inhabit. In gradient culture tubes, the maximum oxygen levels, which were associated with growth bands of Sideroxydans lithotrophicus (ES-1, a novel FeOB), were 15-50 μM. Kinetic measurements made on S. lithotrophicus compared biotic/abiotic (killed control) Fe oxidation rates. The biotic rate can be a significant and measurable fraction of the total Fe oxidation rate below O₂ concentrations of approximately 50 μM, but biotic Fe(II) oxidation (via the biotic/abiotic rate comparison) becomes difficult to detect at higher O₂ levels. These results are further supported by observations of conditions supporting FeOB communities in field settings. Variablity in cell densities and cellular activity as well as variations in hydrous ferrous oxide mineral quantities significantly affect the laboratory kinetic rates. The microbial habitat (or geochemical niche) where FeOB are active is thus largely controlled by the competition between abiotic and biotic kinetics, which are dependent on Fe(II) concentration, P_O2, temperature and pH in addition to the surface area of hydrous ferric oxide minerals and the cell density/activity of FeOB. Additional field and lab culture observations suggest a potentially important role for the iron-sulfide aqueous molecular cluster, FeS_(aq), in the overall cycling of iron associated with the environments these microorganisms inhabit.     

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.     

The persistence of lead from past gasoline emissions and mining drainage in a large riparian system: Evidence from lead isotopes in the Sacramento River, California

Lead concentrations and isotope ratios measured in river water colloids and streambed sediment samples along 426 km of the Sacramento River, California reveal that the influence of lead from the historical mining of massive sulfide deposits in the West Shasta Cu-mining district (at the headwaters of the Sacramento River) is confined to a 60 km stretch of river immediately downstream of that mining region, whereas inputs from past leaded gasoline emissions and historical hydraulic Au-mining in the Sierra Nevadan foothills are the dominant lead sources in the remaining 370 km of the river. Binary mixing calculations suggest that more than 50% of the lead in the Sacramento River outside of the region of influence of the West Shasta Cu-mining district is derived from past depositions of leaded gasoline emissions. This predominance is the first direct documentation of the geographic extent of gasoline lead persistence throughout a large riparian system (>160,000 km²) and corroborates previous observations based on samples taken at the mouth of the Sacramento River. In addition, new analyses of sediment samples from the hydraulic gold mines of the Sierra Nevada foothills confirm the present-day fluxes into the Sacramento River of contaminant metals derived from historical hydraulic Au-mining that occurred during the latter half of the 19th and early part of the 20th centuries. These fluxes occur predominantly during periods of elevated river discharge associated with heavy winter precipitation in northern California. In the broadest context, the study demonstrates the potential for altered precipitation patterns resulting from climate change to affect the mobility and transport of soil-bound contaminants in the surface environment.     

Geochemistry, petrology and ages of the lunar meteorites Kalahari 008 and 009: New constraints on early lunar evolution

Kalahari 008 and 009 are two lunar meteorites that were found close to each other in Botswana. Kalahari 008 is a typical lunar anorthositic breccia; Kalahari 009 a monomict breccia with basaltic composition and mineralogy. Based on minor and trace elements Kalahari 009 is classified as VLT (very-low-Ti) mare basalt with extremely low contents of incompatible elements, including the REE. The Lu-Hf data define an age of 4286 ± 95 Ma indicating that Kalahari 009 is one of the oldest known basalt samples from the Moon. It provides evidence for lunar basalt volcanism prior to 4.1 Ga (pre-Nectarian) and may represent the first sample from a cryptomare. The very radiogenic initial ¹⁷⁶Hf/¹⁷⁷Hf (εHf = +12.9 ± 4.6), the low REE, Th and Ti concentrations indicate that Kalahari 009 formed from re-melting of mantle material that had undergone strong incompatible trace element depletion early in lunar history. This unusually depleted composition points toward a hitherto unsampled basalt source region for the lunar interior that may represent a new depleted endmember source for low-Ti mare basalt volcanism. Apparently, the Moon became chemically very heterogeneous at an early stage in its history and different cumulate sources are responsible for the diverse mare basalt types.Evidence that Kalahari 008 and 009 may be paired includes the similar fayalite content of their olivine, the identical initial Hf isotope composition, the exceptionally low exposure ages of both rocks and the fact that they were found close to each other. Since cryptomaria are covered by highland ejecta, it is possible that these rocks are from the boundary area, where basalt deposits are covered by highland ejecta. The concentrations of cosmogenic radionuclides and trapped noble gases are unusually low in both rocks, although Kalahari 008 contains slightly higher concentrations. A likely reason for this difference is that Kalahari 008 is a polymict breccia containing a briefly exposed regolith, while Kalahari 009 is a monomict brecciated rock that may never have been at the surface of the Moon.Altogether, the compositions of Kalahari 008 and 009 permit new insight into early lunar evolution, as both meteorites sample lunar reservoirs hitherto unsampled by spacecraft missions. The very low Th and REE content of Kalahari 009 as well as the depletion in Sm and the lack of a KREEP-like signature in Kalahari 008 point to a possible source far from the influence of the Procellarum-KREEP Terrane, possibly the lunar farside.     

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.     

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.