Search USGS⌕ Search

Geology topics

Susan L. Brantley

Publications and source records attributed to Susan L. Brantley.

17 recordsLinked to original sources
USGS2021

Signatures of hydrologic function across the critical zone observatory network

Despite a multitude of small catchment studies, we lack a deep understanding of how variations in critical zone architecture lead to variations in hydrologic states and fluxes. This study characterizes hydrologic dynamics of 15 catchments of the U.S. Critical Zone Observatory (CZO) network where we hypothesized that our understanding of subsurface structure would illuminate patterns of hydrologic partitioning. The CZOs collect data sets that characterize the physical, chemical, and biological architecture of the subsurface, while also monitoring hydrologic fluxes such as streamflow, precipitation, and evapotranspiration. For the first time, we collate time series of hydrologic variables across the CZO network and begin the process of examining hydrologic signatures across sites. We find that catchments with low baseflow indices and high runoff sensitivity to storage receive most of their precipitation as rain and contain clay‐rich regolith profiles, prominent argillic horizons, and/or anthropogenic modifications. In contrast, sites with high baseflow indices and low runoff sensitivity to storage receive the majority of precipitation as snow and have more permeable regolith profiles. The seasonal variability of water balance components is a key control on the dynamic range of hydraulically connected water in the critical zone. These findings lead us to posit that water balance partitioning and streamflow hydraulics are linked through the coevolution of critical zone architecture but that much work remains to parse these controls out quantitatively.

Water Resources Research↗
USGS2017

Expanding the role of reactive transport models in critical zone processes

Models test our understanding of processes and can reach beyond the spatial and temporal scales of measurements. Multi-component Reactive Transport Models (RTMs), initially developed more than three decades ago, have been used extensively to explore the interactions of geothermal, hydrologic, geochemical, and geobiological processes in subsurface systems. Driven by extensive data sets now available from intensive measurement efforts, there is a pressing need to couple RTMs with other community models to explore non-linear interactions among the atmosphere, hydrosphere, biosphere, and geosphere. Here we briefly review the history of RTM development, summarize the current state of RTM approaches, and identify new research directions, opportunities, and infrastructure needs to broaden the use of RTMs. In particular, we envision the expanded use of RTMs in advancing process understanding in the Critical Zone, the veneer of the Earth that extends from the top of vegetation to the bottom of groundwater. We argue that, although parsimonious models are essential at larger scales, process-based models offer tools to explore the highly nonlinear coupling that characterizes natural systems. We present seven testable hypotheses that emphasize the unique capabilities of process-based RTMs for (1) elucidating chemical weathering and its physical and biogeochemical drivers; (2) understanding the interactions among roots, micro-organisms, carbon, water, and minerals in the rhizosphere; (3) assessing the effects of heterogeneity across spatial and temporal scales; and (4) integrating the vast quantity of novel data, including “omics” data (genomics, transcriptomics, proteomics, metabolomics), elemental concentration and speciation data, and isotope data into our understanding of complex earth surface systems. With strong support from data-driven sciences, we are now in an exciting era where integration of RTM framework into other community models will facilitate process understanding across disciplines and across scales.

Earth-Science Reviews↗
USGS2016

Quantifying chemical weathering rates along a precipitation gradient on Basse-Terre Island, French Guadeloupe: new insight from U-series isotopes in weathering rinds

Inside soil and saprolite, rock fragments can form weathering clasts (alteration rinds surrounding an unweathered core) and these weathering rinds provide an excellent field system for investigating the initiation of weathering and long term weathering rates. Recently, uranium-series (U-series) disequilibria have shown great potential for determining rind formation rates and quantifying factors controlling weathering advance rates in weathering rinds. To further investigate whether the U-series isotope technique can document differences in long term weathering rates as a function of precipitation, we conducted a new weathering rind study on tropical volcanic Basse-Terre Island in the Lesser Antilles Archipelago. In this study, for the first time we characterized weathering reactions and quantified weathering advance rates in multiple weathering rinds across a steep precipitation gradient. Electron microprobe (EMP) point measurements, bulk major element contents, and U-series isotope compositions were determined in two weathering clasts from the Deshaies watershed with mean annual precipitation (MAP) = 1800 mm and temperature (MAT) = 23 °C. On these clasts, five core-rind transects were measured for locations with different curvature (high, medium, and low) of the rind-core boundary. Results reveal that during rind formation the fraction of elemental loss decreases in the order: Ca ≈ Na > K ≈ Mg > Si ≈ Al > Zr ≈ Ti ≈ Fe. Such observations are consistent with the sequence of reactions after the initiation of weathering: specifically, glass matrix and primary minerals (plagioclase, pyroxene) weather to produce Fe oxyhydroxides, gibbsite and minor kaolinite. Uranium shows addition profiles in the rind due to the infiltration of U-containing soil pore water into the rind as dissolved U phases. U is then incorporated into the rind as Fe-Al oxides precipitate. Such processes lead to significant U-series isotope disequilibria in the rinds. This is the first time that multiple weathering clasts from the same watershed were analyzed for U-series isotope disequlibrian and show consistent results. The U-series disequilibria allowed for the determination of rind formation ages and weathering advance rates with a U-series mass balance model. The weathering advance rates generally decreased with decreasing curvature: ∼0.17 ± 0.10 mm/kyr for high curvature, ∼0.12 ± 0.05 mm/kyr for medium curvature, and ∼0.11 ± 0.04, 0.08 ± 0.03, 0.06 ± 0.03 mm/kyr for low curvature locations. The observed positive correlation between the curvature and the weathering rates is well supported by predictions of weathering models, i.e., that the curvature of the rind-core boundary controls the porosity creation and weathering advance rates at the clast scale. At the watershed scale, the new weathering advance rates derived on the low curvature transects for the relatively dry Deshaies watershed (average rate of 0.08 mm/kyr; MAP = 1800 mm and MAT = 23 °C) are ∼60% slower than the rind formation rates previously determined in the much wetter Bras David watershed (∼0.18 mm/kyr, low curvature transect; MAP = 3400 mm and MAT = 23 °C) also on Basse-Terre Island. Thus, a doubling of MAP roughly correlates with a doubling of weathering advance rate. The new rind study highlights the effect of precipitation on weathering rates over a time scale of ∼100 kyr. Weathering rinds are thus a suitable system for investigating long-term chemical weathering across environmental gradients, complementing short-term riverine solute fluxes.

Basse-Terre Island↗
USGS2013

Probing the deep critical zone beneath the Luquillo Experimental Forest, Puerto Rico

Recent work has suggested that weathering processes occurring in the subsurface produce the majority of silicate weathering products discharged to the world's oceans, thereby exerting a primary control on global temperature via the well-known positive feedback between silicate weathering and CO 2 . In addition, chemical and physical weathering processes deep within the critical zone create aquifers and control groundwater chemistry, watershed geometry and regolith formation rates. Despite this, most weathering studies are restricted to the shallow critical zone (e.g. soils, outcrops). Here we investigate the chemical weathering, fracturing and geomorphology of the deep critical zone in the Bisley watershed in the Luquillo Critical Zone Observatory, Puerto Rico, from two boreholes drilled to 37.2 and 27.0 m depth, from which continuous core samples were taken. Corestones exposed aboveground were also sampled. Weathered rinds developed on exposed corestones and along fracture surfaces on subsurface rocks slough off of exposed corestones once rinds attain a thickness up to ~1 cm, preventing the corestones from rounding due to diffusion limitation. Such corestones at the land surface are assumed to be what remains after exhumation of similar, fractured bedrock pieces that were observed in the drilled cores between thick layers of regolith. Some of these subsurface corestones are massive and others are highly fractured, whereas aboveground corestones are generally massive with little to no apparent fracturing. Subsurface corestones are larger and less fractured in the borehole drilled on a road where it crosses a ridge compared with the borehole drilled where the road crosses the stream channel. Both borehole profiles indicate that the weathering zone extends to well below the stream channel in this upland catchment; hence weathering depth is not controlled by the stream level within the catchment and not all of the water in the watershed is discharged to the stream

Earth Surface Processes and Landforms↗
USGS2012

Where fast weathering creates thin regolith and slow weathering creates thick regolith

Weathering disaggregates rock into regolith – the fractured or granular earth material that sustains life on the continental land surface. Here, we investigate what controls the depth of regolith formed on ridges of two rock compositions with similar initial porosities in Virginia (USA). A priori, we predicted that the regolith on diabase would be thicker than on granite because the dominant mineral (feldspar) in the diabase weathers faster than its granitic counterpart. However, weathering advanced 20 deeper into the granite than the diabase. The 20  -thicker regolith is attributed mainly to connected micron-sized pores, microfractures formed around oxidizing biotite at 20 m depth, and the lower iron (Fe) content in the felsic rock. Such porosity allows pervasive advection and deep oxidation in the granite. These observations may explain why regolith worldwide is thicker on felsic compared to mafic rock under similar conditions. To understand regolith formation will require better understanding of such deep oxidation reactions and how they impact fluid flow during weathering.

Virginia↗
USGS2011

Proposed initiative would study Earth's weathering engine

At the Earth's surface, a complex suite of chemical, biological, and physical processes combines to create the engine that transforms bedrock into soil (Figure 1). Earth's weathering engine provides nutrients to nourish ecosystems and human society mediates the transport of toxic components within the biosphere, creates water flow paths that carve and weaken bedrock, and contributes to the evolution of landscapes at all temporal and spatial scales. At the longest time scales, the weathering engine sequesters CO 2 , thereby influencing long-term climate change. Despite the importance of soil, our knowledge of the rate of soil formation is limited because the weathering zone forms a complex, ever-changing interface, and because scientific approaches and funding paradigms have not promoted integrated research agendas to investigate such complex interactions. No national initiative has promoted a systems approach to investigation of weathering science across the broad array of geology, soil science, ecology and hydrology. Such a program is certainly needed, and this article describes a platform on which to build the initiative to answer the following question: How does the Earth weathering engine break down rock to nourish ecosystems, carve terrestrial landscapes, and control carbon dioxide in the global atmosphere?

Eos Science News↗
USGS2011

Halloysite nanotubes and bacteria at the saprolite-bedrock interface, Rio Icacos watershed, Puerto Rico

Quartz diorite bedrock underlying the Luquillo Mountains of eastern Puerto Rico undergoes weathering at one of the fastest documented rates for granitic rocks in the world. Although tropical temperatures and precipitation promote rapid weathering in this location, increased bacterial densities in the regolith immediately above the bedrock suggest that microorganisms contribute to mineral weathering as well. Deep saprolite and saprock samples were obtained at the bedrock interface in an upland location (Guaba Ridge) in the Rio Icacos watershed for examination by environmental scanning electron microscopy (ESEM). In ESEM images, mineral nanotubes were observed to occur frequently in association with coccus- and rod-shaped structures resembling bacteria. These nanotubes (50–140-nm width and 150–2700-nm length) were identified as halloysite using transmission electron microscopy. Observations of multiple nanotubes on the surfaces of an individual cell are consistent with the cell's exterior functional groups interacting with Si in pore water to facilitate halloysite nucleation. We propose that one mechanism by which bacteria contribute to the rapid weathering of quartz diorite minerals in this regolith is by lowering the free energy for secondary mineral formation. The presence of bacterial surfaces may result in more rapid removal of Si from solution, thereby increasing the dissolution rates of primary minerals.

Puerto Rico↗
USGS2010

Phosphorus and iron cycling in deep saprolite, Luquillo Mountains, Puerto Rico

Rapid weathering and erosion rates in mountainous tropical watersheds lead to highly variable soil and saprolite thicknesses which in turn impact nutrient fluxes and biological populations. In the Luquillo Mountains of Puerto Rico, a 5-m thick saprolite contains high microorganism densities at the surface and at depth overlying bedrock. We test the hypotheses that the organisms at depth are limited by the availability of two nutrients, P and Fe. Many tropical soils are P-limited, rather than N-limited, and dissolution of apatite is the dominant source of P. We document patterns of apatite weathering and of bioavailable Fe derived from the weathering of primary minerals hornblende and biotite in cores augered to 7.5 m on a ridgetop as compared to spheroidally weathering bedrock sampled in a nearby roadcut. Iron isotopic compositions of 0.5 N HCl extracts of soil and saprolite range from about δ 56 Fe = 0 to − 0.1‰ throughout the saprolite except at the surface and at 5 m depth where δ 56 Fe = − 0.26 to − 0.64‰. The enrichment of light isotopes in HCl-extractable Fe in the soil and at the saprolite–bedrock interface is consistent with active Fe cycling and consistent with the locations of high cell densities and Fe(II)-oxidizing bacteria, identified previously. To evaluate the potential P-limitation of Fe-cycling bacteria in the profile, solid-state concentrations of P were measured as a function of depth in the soil, saprolite, and weathering bedrock. Weathering apatite crystals were examined in thin sections and an apatite dissolution rate of 6.8 × 10 − 14 mol m − 2 s − 1 was calculated. While surface communities depend on recycled nutrients and atmospheric inputs, deep communities survive primarily on nutrients released by the weathering bedrock and thus are tightly coupled to processes related to saprolite formation including mineral weathering. While low available P may limit microbial activity within the middle saprolite, fluxes of P from apatite weathering should be sufficient to support robust growth of microorganisms in the deep saprolite.

Puerto Rico↗
USGS2010

Effects of temperature on silicate weathering: Solute fluxes and chemical weathering in a temperate rain forest watershed, Jamieson Creek, British Columbia

Chemical weathering of silicate minerals has long been known as a sink for atmospheric CO 2 , and feedbacks between weathering and climate are believed to affect global climate. While warmer temperatures are believed to increase rates of weathering, weathering in cool climates can be accelerated by increased mineral exposure due to mechanical weathering by ice. In this study, chemical weathering of silicate minerals is investigated in a small temperate watershed. The Jamieson Creek watershed is covered by mature coniferous forest and receives high annual precipitation (4000 mm), mostly in the form of rainfall, and is underlain by quartz diorite bedrock and glacial till. Analysis of pore water concentration gradients indicates that weathering in hydraulically unsaturated ablation till is dominated by dissolution of plagioclase and hornblende. However, a watershed scale solute mass balance indicates high relative fluxes of K and Ca, indicating preferential leaching of these solutes possibly from the relatively unweathered lodgement till. Weathering rates for plagioclase and hornblende calculated from a watershed scale solute mass balance are similar in magnitude to rates determined using pore water concentration gradients. When compared to the Rio Icacos basin in Puerto Rico, a pristine tropical watershed with similar annual precipitation and bedrock, but with dissimilar regolith properties, fluxes of weathering products in stream discharge from the warmer site are 1.8 to 16.2-fold higher, respectively, and regolith profile-averaged plagioclase weathering rates are 3.8 to 9.0-fold higher. This suggests that the Arrhenius effect, which predicts a 3.5- to 9-fold increase in the dissolution rate of plagioclase as temperature is increased from 3.4° to 22 °C, may explain the greater weathering fluxes and rates at the Rio Icacos site. However, more modest differences in K and Ca fluxes between the two sites are attributed to accelerated leaching of those solutes from glacial till at Jamieson Creek. Our findings suggest that under conditions of high rainfall and favorable topography, weathering rates of silicate minerals in warm tropical systems will tend to be higher than in cool temperate systems, even if the temperate system is has been perturbed by an episode of glaciation that deposits regolith high in fresh mineral surface area.

British Columbia↗
USGS2010

Shifting microbial community structure across a marine terrace grassland chronosequence, Santa Cruz, California

Changes in the biomass and structure of soil microbial communities have the potential to impact ecosystems via interactions with plants and weathering minerals. Previous studies of forested long-term (1000s – 100,000s of years) chronosequences suggest that surface microbial communities change with soil age. However, significant gaps remain in our understanding of long-term soil microbial community dynamics, especially for non-forested ecosystems and in subsurface soil horizons. We investigated soil chemistry, aboveground plant productivity, and soil microbial communities across a grassland chronosequence (65,000–226,000 yrs old) located near Santa Cruz, CA. Aboveground net primary productivity (ANPP) initially increased to a maximum and then decreased for the older soils. We used polar lipid fatty acids (PLFA) to investigate microbial communities including both surface (<0.1 m) and subsurface (≥0.2 m) soil horizons. PLFAs characteristic of Gram-positive bacteria and actinobacteria increased as a fraction of the microbial community with depth while the fungal fraction decreased relative to the surface. Differences among microbial communities from each chronosequence soil were found primarily in the subsurface where older subsurface soils had smaller microbial community biomass, a higher proportion of fungi, and a different community structure than the younger subsurface soil. Subsurface microbial community shifts in biomass and community structure correlated with, and were likely driven by, decreasing soil P availability and Ca concentrations, respectively. Trends in soil chemistry as a function of soil age led to the separation of the biological (≤1 m depth) and geochemical (>1 m) cycles in the old, slowly eroding landscape we investigated, indicating that this separation, commonly observed in tropical and subtropical ecosystems, can also occur in temperate climates. This study is the first to investigate subsurface microbial communities in a long-term chronosequence. Our results highlight connections between soil chemistry and both the aboveground and belowground parts of an ecosystem.

California↗
USGS2009

Approaches to modeling weathered regolith

Sustainable soils are a requirement for maintaining human civilizations ( Carter and Dale 1974 ; Lal 1989 ). However, as the “most complicated biomaterial on the planet” ( Young and Crawford 2004 ), soils represent one of the most difficult systems to understand and model with respect to chemical, physical, and biological coupling over time (Fig. 1 ). Despite the complexity of these interactions, certain patterns in soil properties and development are universally observed and have been used in soil science as a means for classification. Elemental, mineralogical, or isotopic concentrations in soils plotted versus depth beneath the land surface comprise such patterns. Soil depth profiles are often reported for solid soil materials, and, less frequently, for solutes in soil pore waters. These profiles cross a large range in spatial scales that traditionally have been studied by different disciplines. For example, shallow, biologically active horizons are commonly defined as the soil zone in agronomic studies whereas the mobile layer of the regolith is referred to as soil in geomorphological studies. In contrast, many geochemical studies target chemical weathering to tens or even hundreds of meters in depth, sometimes extending the definition of “soils” to include the entire regolith down to parent bedrock or alluvium. Soil profiles also exhibit a large range in temporal scales ( Amundson 2004 ; Brantley 2008b ). Solid-state profiles document chemical and mineralogical changes integrated over the time scales of evolution of regolith from protolith. This “geologic time” can vary from tens to hundreds of years for weathered material developed on moraines deposited by active glaciers ( Anderson et al. 1997 ), to millions or possibly hundreds of millions of years of regolith evolution as documented in laterites and bauxites on stable cratons ( Nahon 1986 ). In contrast, solute profiles reflect much shorter time scales corresponding to the residence time of the soil water which commonly ranges from days to decades ( Stonestrom et al. 1998 ). Factors impacting soil minerals can therefore be related to geologically old processes while those impacting pore waters are related to contemporary processes. We first discuss a geochemical frame work for modeling soil profiles, including a simple scheme that depends on the extent of enrichment or depletion. Such profiles are comprised of reaction fronts affected by chemical, hydrologic, geologic and biologic processes that control soil evolution. We then present a hierarchy of models that have been used to interpret both solid state and solute compositions in regolith. The more simple approaches to model depletion in soils, using analytical models, are first described. The most elementary of these is a linear model that calculates rate constants from the slopes of either solid or solute weathering gradients: these rate constants represent lumped parameters that describe weathering in terms of an integrated reaction rate. Two other analytical models are then presented that have been used to fit solid state elemental profiles with exponential and sigmoidal functions. All of these analytical approaches are derived for models of soils as containing a limited number of components, phases, and species. At a more complex level, numerical models are then presented to elucidate how kinetic and transport parameters as well as chemical, hydrologic, and physical soil data can be incorporated. We consider two forms of these models, first relatively simple spreadsheet calculators and then more sophisticated multi-component, multi-phase reactive-transport numerical codes. Our treatment of reactive transport modeling is relatively cursory, in recognition of the treatment in the chapter by Steefel and Maher (2009 , this volume). Because these models incorporate more phases, components, and species than the other approaches and explicitly model the more fundamental reaction mechanisms involved, they generally have a greater need for parameterization. In our conclusion section, we discuss how this hierarchy of approaches can yield generalizations about soils that are often complementary.

Reviews in Mineralogy and Geochemistry↗
USGS2008

Preface

No abstract available.

Book chapter↗
USGS2004

Mineral dissolution in the Cape Cod aquifer, Massachusetts, USA: I . Reaction stoichiometry and impact of accessory feldspar and glauconite on strontium isotopes, solute concentrations, and REY distribution

To compare relative reaction rates of mineral dissolution in a mineralogically simple groundwater aquifer, we studied the controls on solute concentrations, Sr isotopes, and rare earth element and yttrium (REY) systematics in the Cape Cod aquifer. This aquifer comprises mostly carbonate-free Pleistocene sediments that are about 90% quartz with minor K-feldspar, plagioclase, glauconite, and Fe-oxides. Silica concentrations and pH in the groundwater increase systematically with increasing depth, while Sr isotopic ratios decrease. No clear relationship between 87 Sr/ 86 Sr and Sr concentration is observed. At all depths, the 87 Sr/ 86 Sr ratio of the groundwater is considerably lower than the Sr isotopic ratio of the bulk sediment or its K-feldspar component, but similar to that of a plagioclase-rich accessory separate obtained from the sediment. The Si- 87 Sr/ 86 Sr-depth relationships are consistent with dissolution of accessory plagioclase. In addition, solutes such as Sr, Ca, and particularly K show concentration spikes superimposed on their respective general trends. The K-Sr- 87 Sr/ 86 Sr systematics suggests that accessory glauconite is another major solute source to Cape Cod groundwater. Although the authigenic glauconite in the Cape Cod sediment is rich in Rb, it is low in in-grown radiogenic 87 Sr because of its young Pleistocene age. The low 87 Sr/ 86 Sr ratios are consistent with equilibration of glauconite with seawater. The impact of glauconite is inferred to vary due to its variable abundance in the sediments. In the Cape Cod groundwater, the variation of REY concentrations with sampling depth resembles that of K and Rb, but differs from that of Ca and Sr. Shale-normalized REY patterns are light REY depleted, show negative Ce anomalies and super-chondritic Y/Ho ratios, but no Eu anomalies. REY input from feldspar, therefore, is insignificant compared to input from a K-Rb-bearing phase, inferred to be glauconite. These results emphasize that interpretation of groundwater chemistry, even in relatively simple aquifers, may be complicated by solute contributions from “exotic” accessory minerals such as glauconite. To detect such peculiarities, groundwater studies should combine the study of elemental concentration and isotopic composition of several solutes that show different geochemical behavior.

Massachusetts↗
USGS2003

The effect of time on the weathering of silicate minerals: Why do weathering rates differ in the laboratory and field?

The correlation between decreasing reaction rates of silicate minerals and increasing duration of chemical weathering was investigated for both experimental and field conditions. Column studies, using freshly prepared Panola Granite, produced ambient plagioclase weathering rates that decreased parabolically over 6 years to a final rate of 7.0×10 −14 mol m −2 s −1 . In contrast, the corresponding plagioclase reaction rate for partially kaolinized Panola Granite, after reaching steady-state weathering after 2 months of reactions, was significantly less (2.1×10 −15 mol m −2 s −1 ). Both rates were normalized to plagioclase content and BET surface area. Extrapolation of decreasing rates for the fresh plagioclase with time indicated that several thousand years of reaction would be required to replicate the rate of the naturally weathered plagioclase under identical experimental conditions. Both rates would remain orders of magnitude faster than field weathering rates previously measured for a weathering profile in the Panola Granite. Additional trends in weathering rates with time were established from a tabulation of previously reported experimental and field rates for plagioclase, K-feldspar, hornblende and biotite. Discrepancies in the literature, produced by normalization of weathering rates with respect to surface areas measured by gas absorption (BET) and geometric methods, were overcome by developing a time-dependent roughness factor. Regression curves through the corrected rates produced strong correlations with time that were similar for the four silicate minerals. The average silicate weathering rate R (mol m −2 s −1 ) was described by the power function R=3.1&#xD7;10 &#x2212;13 t &#x2212;0.61 "> R=3.1×10 -13 t -0.61 which was similar to the relationship describing the decrease in the fresh Panola plagioclase with time and suggesting control by transport-limited reaction. The time dependency of silicate weathering is discussed in terms of processes intrinsic to the silicate mineral and those extrinsic to the weathering environment. Intrinsic surface area, which increases with the duration of weathering, was shown to account for a third of the exponential decrease in the weathering rate shown by the above equation. Other factors, including progressive depletion of energetically reactive surfaces and accumulation of leached layers and secondary precipitates, must explain differences for fresh and weathered plagioclase reacting under identical experimental conditions. Extrinsic controls, including low permeability, high mineral/fluid ratios and increased solute concentrations, produce weathering reactions close to thermodynamic equilibrium under field conditions compared to highly unsaturated conditions during experimental reaction of fresh and weathered plagioclase. These differences explain the time-dependent difference in field and lab rates.

Chemical Geology↗
USGS2001

Fractionation of Fe isotopes by soil microbes and organic acids

Small natural variations in Fe isotopes have been attributed to biological cycling. However, without understanding the mechanism of fractionation, it is impossible to interpret such variations. Here we show that the &delta; 56 Fe of Fe dissolved from a silicate soil mineral by siderophore-producing bacteria is as much as 0.8% lighter than bulk Fe in the mineral. A smaller isotopic shift is observed for Fe released abiotically by two chelates, and the magnitude of the shift increases with affinity of the ligand for Fe, consistent with a kinetic isotope effect during hydrolysis of Fe at the mineral surface. Fe dissolved abiotically without chelates shows no isotopic shift. The &delta; 56 Fe of the exchange fraction on soil grains is also lighter by ~0.6%-1% than Fe from both hornblende and iron oxyhydroxides. The kinetic isotope effect is therefore preserved in open systems such as soils. when recorded in the rock record, Fe isotopic fractionation could document Fe transport by organic molecules or by microbes where such entities were present in the geologic past.

Geology↗
USGS1995

Groundwater transport of crater-lake brine at Poa´s Volcano, Costa Rica

Poa´s Volcano is an active stratovolcano in Costa Rica that has a lake in its active crater. The crater lake has high temperatures (50–90 °C), high acidity (pH ≈ 0.0), and a high dissolved-solids content (100 g/kg). The volcano has numerous freshwater springs on its flanks, but a few on the northwestern flank are highly acidic (pH = 1.6–2.5) and have high dissolved-solids concentrations (2–22 g/kg). This study analyzes the regional groundwater system at Poa´s and demonstrates the likelihood that the water discharging from the acidic springs in the Rio Agrio watershed originates at the acidic crater lake. Both heat and solute transport are analyzed on a regional scale through numerical simulations using the HST3D finite-difference model, which solves the coupled equations for fluid flow, heat transport, and solute transport. The code allows fluid viscosity and density to be functions of both temperature and solute concentration. The simulations use estimates for recharge to the mountain and a range of values and various distributions of permeability and porosity. Several sensitivity analyses are performed to test how the uncertainty in many of the model parameters affects the simulation results. These uncertainties yield an estimated range of travel times from the crater lake to the Rio Agrio springs of 1–30 years, which is in close agreement with the results of tritium analyses of the springs. Calculated groundwater fluxes into and out of the crater lake are both about several hundred kg/s. These fluxes must be accounted for in water budgets of the crater lake.

Poa's Volcano↗