USGS Science⌕ Search

SEARCH · USGS Science

Results for “Chemistry”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,765 records · Page 98Linked to original sources

Mineralogical, chemical and K-Ar isotopic changes in Kreyenhagen Shale whole rocks and <2 µm clay fractions during natural burial and hydrous-pyrolysis experimental maturation

Progressive maturation of the Eocene Kreyenhagen Shale from the San Joaquin Basin of California was studied by combining mineralogical and chemical analyses with K&ndash;Ar dating of whole rocks and <2 &mu;m clay fractions from naturally buried samples and laboratory induced maturation by hydrous pyrolysis of an immature outcrop sample. The K&ndash;Ar age decreases from 89.9 &plusmn; 3.9 and 72.4 &plusmn; 4.2 Ma for the outcrop whole rock and its <2 &mu;m fraction, respectively, to 29.7 &plusmn; 1.5 and 21.0 &plusmn; 0.7 Ma for the equivalent materials buried to 5167 m. The natural maturation does not produce K&ndash;Ar ages in the historical sense, but rather K/Ar ratios of relative K and radiogenic 40 Ar amounts resulting from a combined crystallization of authigenic and alteration of initial detrital K-bearing minerals of the rocks. The Al/K ratio of the naturally matured rocks is essentially constant for the entire depth sequence, indicating that there is no detectable variation in the crystallo-chemical organization of the K-bearing alumino-silicates with depth. No supply of K from outside of the rock volumes occurred, which indicates a closed-system behavior for it. Conversely, the content of the total organic carbon (TOC) content decreases significantly with burial, based on the progressive increasing Al/TOC ratio of the whole rocks. The initial varied mineralogy and chemistry of the rocks and their <2 &mu;m fractions resulting from differences in detrital sources and depositional settings give scattered results that homogenize progressively during burial due to increased authigenesis, and concomitant increased alteration of the detrital material. Hydrous pyrolysis was intended to alleviate the problem of mineral and chemical variations in initially deposited rocks of naturally matured sequences. However, experiments on aliquots from thermally immature Kreyenhagen Shale outcrop sample did not mimic the results from naturally buried samples. Experiments conducted for 72 h at temperatures from 270 to 365 &deg;C did not induce significant changes at temperatures above 310 &deg;C in the mineralogical composition and K&ndash;Ar ages of the rock and <2 &mu;m fraction. The K&ndash;Ar ages of the <2 &mu;m fraction range from 72.4 &plusmn; 4.2 Ma in the outcrop sample to 62.4 &plusmn; 3.4 Ma in the sample heated the most at 365 &deg;C for 216 h. This slight decrease in age outlines some loss of radiogenic 40 Ar, together with losses of organic matter as oil, gas, and aqueous organic species. Large amounts of smectite layers in the illite&ndash;smectite mixed layers of the pyrolyzed outcrop <2 &mu;m fraction remain during thermal experiments, especially above 310 &deg;C. With no illitization detected above 310 &deg;C, smectite appears to have inhibited rather than promoted generation of expelled oil from decomposition of bitumen. This hindrance is interpreted to result from bitumen impregnating the smectite interlayer sites and rock matrix. Bitumen remains in the <2 &mu;m fraction despite leaching with H 2 O 2 . Its presence in the smectite interlayers is apparent by the inability of the clay fraction to fully expand or collapse once bitumen generation from the thermal decomposition of the kerogen is completed, and by almost invariable K&ndash;Ar ages confirming for the lack of any K supply and/or radiogenic 40 Ar removal. This suggests that once bitumen impregnates the porosity of a progressively maturing source rock, the pore system is no longer wetted by water and smectite to illite conversion ceases. Experimental attempts to evaluate the smectite conversion to illite should preferentially use low-TOC rocks to avoid inhibition of the reaction by bitumen impregnation.

California↗

Normalization of stable isotope data for carbonate minerals: implementation of IUPAC guideline

Carbonate minerals provide a rich source of geochemical information because their &delta; 13 C and &delta; 18 O values provide information about surface and subsurface Earth processes. However, a significant problem is that the same &delta; 18 O value is not reported for the identical carbonate sample when analyzed in different isotope laboratories in spite of the fact that the International Union of Pure and Applied Chemistry (IUPAC) has provided reporting guidelines for two decades. This issue arises because (1) the &delta; 18 O measurements are performed on CO 2 evolved by reaction of carbonates with phosphoric acid, (2) the acid-liberated CO 2 is isotopically fractionated (enriched in 18 O) because it contains only two-thirds of the oxygen from the solid carbonate, (3) this oxygen isotopic fractionation factor is a function of mineralogy, temperature, concentration of the phosphoric acid, and &delta; 18 O value of water in the phosphoric acid, (4) researchers may use any one of an assortment of oxygen isotopic fractionation factors that have been published for various minerals at various reaction temperatures, and (5) it sometimes is not clear how one should calculate &delta; 18 O VPDB values on a scale normalized such that the &delta; 18 O value of SLAP reference water is &minus;55.5 &permil; relative to VSMOW reference water. To enable researchers worldwide to publish the same &delta; 18 O value (within experimental uncertainty) for the same carbonate sample, we have re-evaluated reported acid fractionation factors for calcite at 25, 50, and 75 &deg;C and propose a revised relation for the temperature dependence of oxygen isotopic acid fractionation factor, &alpha; CO 2 ( ACID ) -calcite &alpha;CO2(ACID)-calcite , of 1000ln &alpha; CO 2 ( ACID ) -calcite=3.48(10 3 /T)-1.47 1000ln&alpha;CO2(ACID)-calcite=3.48(103/T)-1.47 Turn MathJax on where T is temperature in kelvin. At 25 &deg;C, &alpha; CO 2 ( ACID ) -calcite=1.01025 &alpha;CO2(ACID)-calcite=1.01025 , the most commonly accepted value for this quantity. We propose a normalization protocol in which (1) the internationally distributed carbonate isotopic reference materials NBS 18 and NBS 19 are interspersed among carbonate samples analyzed by treatment with phosphoric acid, (2) the &delta; 18 O values of the calcite reference materials and the carbonate samples are calculated, respectively, by using the &alpha; CO 2 ( ACID ) -calcite &alpha;CO2(ACID)-calcite relation above and oxygen-isotope acid fractionation factors appropriate for the sample mineralogy and reaction temperature, (3) the &delta; 18 O values of solid carbonate samples are determined on the VPDB scale ( &delta; 18 O VPDB ) with IUPAC-recommended scale expansion such that the &delta; 18 O of SLAP reference water is &minus;55.5 &permil; relative to VSMOW reference water by normalizing &delta; 18 O values of carbonate samples with 2014-IUPAC-recommended &delta; 18 O values of NBS 18 and NBS 19, and (4) &delta; 18 O values on the VPDB scale are converted to &delta; 18 O values on the VSMOW-SLAP scale by using IUPAC recommendations. To ease calculations in the protocol, a software application titled &ldquo;Carbon and Oxygen Isotopic Normalization Tool for Carbonates&rdquo; is available that relies upon IUPAC-recommended &delta; 13 C and &delta; 18 O values of carbonate isotopic reference materials ( http://isotopes.usgs.gov/research/topics/carbonatesnormalizationtool.html ).

Geochimica et Cosmochimica Acta↗

Hydrothermal alteration and diagenesis of terrestrial lacustrine pillow basalts: Coordination of hyperspectral imaging with laboratory measurements

We investigate an outcrop of &sim;187 Ma lacustrine pillow basalts of the Talcott Formation exposed in Meriden, Connecticut, USA, focusing on coordinated analyses of one pillow lava to characterize the aqueous history of these basalts in the Hartford Basin. This work uses a suite of multidisciplinary measurements, including hyperspectral imaging, other spectroscopic techniques, and chemical and mineralogical analyses, from the microscopic scale up to the scale of an outcrop. The phases identified in the sample are albite, large iron oxides, and titanite throughout; calcite in vesicles; calcic clinopyroxene, aegirine, and Fe/Mg-bearing clay in the rind; and fine-grained hematite and pyroxenes in the interior. Using imaging spectroscopy, the chemistry and mineralogy results extend to the hand sample and larger outcrop. From all of the analyses, we suggest that the pillow basalts were altered initially after emplacement, either by heated lake water or magmatic fluids, at temperatures of at least 400-600&deg;C, and the calcic clinopyroxenes and aegirine identified in the rind are a preserved record of that alteration. As the hydrothermal system cooled to slightly lower temperatures, clays formed in the rind, and, during this alteration, the sample oxidized to form hematite in the matrix of the interior and Fe 3+ in the pyroxenes in the rind. During the waning stages of the hydrothermal system, calcite precipitated in vesicles within the rind. Later, diagenetic processes albitized the sample, with albite replacing plagioclase, lining vesicles, and accreting onto the exterior of the sample. This albitization or Na-metasomatism occurred when the lake within the Hartford Basin evaporated during a drier past climatic era, resulting in Na-rich brines. As Ca-rich plagioclase altered to albite, Ca was released into solution, eventually precipitating as calcite in previously-unfilled vesicles, dominantly in the interior of the pillow. Coordinated analyses of this sample permit identification of the alteration phases and help synthesize the aqueous history of pillow lavas of the Talcott formation. These results are also relevant to Mars, where volcanically-resurfaced open basin lakes have been found, and this Hartford Basin outcrop may be a valuable analog for any potential volcano-lacustrine interactions. The results can also help to inform the utility and optimization of potentially complementary, synergistic, and uniquely-suited techniques for characterization of hydrothermally-altered terrains.

Connecticut↗

Lithological influences on contemporary and long-term regolith weathering at the Luquillo Critical Zone Observatory

Lithologic differences give rise to the differential weatherability of the Earth’s surface and globally variable silicate weathering fluxes, which provide an important negative feedback on climate over geologic timescales. To isolate the influence of lithology on weathering rates and mechanisms, we compare two nearby catchments in the Luquillo Critical Zone Observatory in Puerto Rico, which have similar climate history, relief and vegetation, but differ in bedrock lithology. Regolith and pore water samples with depth were collected from two ridgetops and at three sites along a slope transect in the volcaniclastic Bisley catchment and compared to existing data from the granitic Río Icacos catchment. The depth variations of solid-state and pore water chemistry and quantitative mineralogy were used to calculate mass transfer (tau) and weathering solute profiles, which in turn were used to determine weathering mechanisms and to estimate weathering rates. Regolith formed on both lithologies is highly leached of most labile elements, although Mg and K are less depleted in the granitic than in the volcaniclastic profiles, reflecting residual biotite in the granitic regolith not present in the volcaniclastics. Profiles of both lithologies that terminate at bedrock corestones are less weathered at depth, near the rock-regolith interfaces. Mg fluxes in the volcaniclastics derive primarily from dissolution of chlorite near the rock-regolith interface and from dissolution of illite and secondary phases in the upper regolith, whereas in the granitic profile, Mg and K fluxes derive from biotite dissolution. Long-term mineral dissolution rates and weathering fluxes were determined by integrating mass losses over the thickness of solid-state weathering fronts, and are therefore averages over the timescale of regolith development. Resulting long-term dissolution rates for minerals in the volcaniclastic regolith include chlorite: 8.9 × 10 −14 mol m −2 s −1 , illite: 2.1 × 10 −14 mol m −2 s −1 and kaolinite: 4.0 × 10 −14 mol m −2 s −1 . Long-term weathering fluxes are several orders of magnitude lower in the granitic regolith than in the volcaniclastic, despite higher abundances of several elements in the granitic regolith. Contemporary weathering fluxes were determined from net (rain-corrected) solute profiles and thus represent rates over the residence time of water in the regolith. Contemporary weathering fluxes within the granitic regolith are similar to the long-term fluxes. In contrast, the long-term fluxes are faster than the contemporary fluxes in the volcaniclastic regolith. Contemporary fluxes in the granitic regolith are generally also slightly faster than in the volcaniclastic. The differences in weathering fluxes over space and time between these two watersheds indicate significant lithologic control of chemical weathering mechanisms and rates.

Bisley watersheds↗

Relationship between water and aragonite barium concentrations in aquaria reared juvenile corals

Coral barium to calcium (Ba/Ca) ratios have been used to reconstruct records of upwelling, river and groundwater discharge, and sediment and dust input to the coastal ocean. However, this proxy has not yet been explicitly tested to determine if Ba inclusion in the coral skeleton is directly proportional to seawater Ba concentration and to further determine how additional factors such as temperature and calcification rate control coral Ba/Ca ratios. We measured the inclusion of Ba within aquaria reared juvenile corals ( Favia fragum ) at three temperatures (∼27.7, 24.6 and 22.5 °C) and three seawater Ba concentrations (73, 230 and 450 nmol kg −1 ). Coral polyps were settled on tiles conditioned with encrusting coralline algae, which complicated chemical analysis of the coral skeletal material grown during the aquaria experiments. We utilized Sr/Ca ratios of encrusting coralline algae (as low as 3.4 mmol mol −1 ) to correct coral Ba/Ca for this contamination, which was determined to be 26 ± 11% using a two end member mixing model. Notably, there was a large range in Ba/Ca across all treatments, however, we found that Ba inclusion was linear across the full concentration range. The temperature sensitivity of the distribution coefficient is within the range of previously reported values. Finally, calcification rate, which displayed large variability, was not correlated to the distribution coefficient. The observed temperature dependence predicts a change in coral Ba/Ca ratios of 1.1 μmol mol −1 from 20 to 28 °C for typical coastal ocean Ba concentrations of 50 nmol kg −1 . Given the linear uptake of Ba by corals observed in this study, coral proxy records that demonstrate peaks of 10–25 μmol mol −1 would require coastal seawater Ba of between 60 and 145 nmol kg −1 . Further validation of the coral Ba/Ca proxy requires evaluation of changes in seawater chemistry associated with the environmental perturbation recorded by the coral as well as verification of these results for Porites species, which are widely used in paleo reconstructions.

Geochimica et Cosmochimica Acta↗

Root-driven weathering impacts on mineral-organic associations in deep soils over pedogenic time scales

Plant roots are critical weathering agents in deep soils, yet the impact of resulting mineral transformations on the vast deep soil carbon (C) reservoir are largely unknown. Root-driven weathering of primary minerals may cause the formation of reactive secondary minerals, which protect mineral-organic associations (MOAs) for centuries or millennia. Conversely, root-driven weathering may also transform secondary minerals, potentially enhancing the bioavailability of C previously protected in MOAs. Here we examined the impact of root-driven weathering on MOAs and their capacity to store C over pedogenic time scales. To accomplish this, we examined deep horizons (100-160 cm) that experienced root-driven weathering in four soils of increasing ages (65-226 kyr) of the Santa Cruz Marine Terrace chronosequence. Specifically, we compared discrete rhizosphere zones subject to root-driven weathering, with adjacent zones that experienced no root growth. Using a combination of radiocarbon, mass spectrometry, 57Fe Mössbauer spectroscopy, high-resolution mass spectrometry, and X-ray spectromicroscopy approaches, we characterized transformations of MOAs in relation to changes in C content, Δ14C values, and chemistry across the chronosequence. We found that the onset of root-driven weathering (65-90kyr) increased the amount of C associated with poorly crystalline iron (Fe) and aluminum (Al) phases, particularly highly disordered nano-particulate goethite (np-goethite). This increase coincided with greater C concentrations, lower Δ14C values, and greater abundance of what is likely microbially-derived C. Continued root-driven weathering (137-226kyr) did not significantly change the amount of C associated with crystalline Fe and Al phases, but resulted in a decline in the amount of C associated with poorly crystalline Fe and Al phases. This decline coincided with a decrease in C concentrations, an increase in 14C values, and a shift toward plant-derived C. In contrast, soil not affected by root-driven weathering showed comparatively low amounts of C bound to poorly crystalline Fe and Al phases regardless of soil age and, correspondingly, lower C concentrations. Our results demonstrate that root-driven formation and disruption of MOAs are direct controls on both C accrual and loss in deep soil. This finding suggests that root impacts on soil C storage are dependent on soil weathering stage, a consideration that is critical for future predictions of the vulnerability of deep soil C to global change.

California↗

Position-specific distribution of hydrogen isotopes in natural propane: Effects of thermal cracking, equilibration and biodegradation

Intramolecular isotope distributions, including isotope clumping and position specific fractionation, can provide proxies for the formation temperature and formation and destruction pathways of molecules. In this study, we explore the position-specific hydrogen isotope distribution in propane. We analyzed propane samples from 10 different petroleum systems with high-resolution molecular mass spectrometry. Our results show that the hydrogen isotope fractionation between central and terminal positions of natural propanes ranges from −102‰ to +205‰, a much larger range than that expected for thermodynamic equilibrium at their source and reservoir temperatures (36–63‰). Based on these findings, we propose that the hydrogen isotope structure of catagenic propane is largely controlled by irreversible processes, expressing kinetic isotope effects (KIEs). Kinetic control on hydrogen isotope composition of the products of thermal cracking is supported by a hydrous pyrolysis experiment using the Woodford Shale as substrate, in which we observed isotopic disequilibrium in the early stage of pyrolysis. We make a more general prediction of KIE signatures associated with kerogen cracking by simulating this chemistry in a kinetic Monte Carlo model for different types of kerogens. In contrast, unconventional shale fluids or hot conventional reservoirs contain propane with an isotopic structure close to equilibrium, presumably reflecting internal and/or heterogeneous exchange during high temperature storage (ca. 100–150 °C). In relatively cold (<100 °C) conventional gas accumulations, propane can discharge from its source to a colder reservoir, rapidly enough to preserve disequilibrium signatures even if the source rock thermal maturity is high. These findings imply that long times at elevated temperatures are required to equilibrate the hydrogen isotopic structure of propane in natural gas host rocks and reservoirs. We further defined the kinetics of propane equilibration through hydrogen isotope exchange experiments under hydrous conditions; these experiments show that hydrogen in propane is exchangeable over laboratory timescales when exposed to clay minerals such as kaolinite. This implies rather rapid transfer of propane from sources to cold reservoirs in some of the conventional petroleum systems. Propane is also susceptible to microbial degradation in both oxic and anoxic environments. Biodegradation of propane in the Hadrian and Diana Hoover oil fields (Gulf of Mexico) results in strong increases in central-terminal hydrogen isotope fractionation. This reflects preferential attack on the central position, consistent with previous studies.

Geochimica et Cosmochimica Acta↗

Development of soil radiocarbon profiles in a reactive transport framework

Today, there is a greater appreciation for the importance of the physical protection of carbon (C) through interactions with mineral surfaces, isolation from microbes, and the important role of transport in shaping soil properties and controlling moisture limitations on decomposition. As our paradigm for soil organic carbon (SOC) preservation changes, so too should our representation of the underlying processes in soil models. Reactive transport models (RTMs) provide a framework capable of assessing the interactive influence of soil chemistry and transport processes on the accumulation and turnover of SOC. In this study, we present new developments in the isotopically enabled RTM “CrunchTope,” which is capable of explicitly tracking the three isotopes of carbon ( 12 C, 13 C, and 14 C) and their fractionation between multiple coexisting and interacting solid, liquid and gas phases. This modeling framework opens the door to new applications of depth-resolved RTMs models in application to SOC and deeper subsurface carbon reservoirs. Here, we demonstrate SOC accumulation and radiocarbon aging for long-timescale models of soil development in CrunchTope. Our goal is to assess advantages and limitations of such an approach and to identify the type and complexity of reaction networks that are required to adequately apply this model to SOC dynamics. We assess the behavior of this model relative to a high-resolution dataset of SOC content, stable isotope composition, and radiocarbon ages as well as physical and hydrologic data measured from a chronosequence of soils located near Santa Cruz, California. Starting from a previously published model using a simplified reaction network with a single class of carbon, we sequentially incorporate multiple C reservoirs subject to both reactivity and transport pathways. Our results indicate that multiple SOC pools with different mean ages of C do not inherently emerge as a result of including reactions which are conventionally expected to provide a diversity of transit times, i.e., sorption and complexation of SOC on mineral surfaces. Instead, transit times emerge as a result of the timescales of the reactions represented in the reaction network. For mineral associated C, the RTM framework imposes dynamic equilibrium with the fluid phase dissolved organic C, such that no distinction in radiocarbon ages is achieved between these pools. Aged C can be produced by including a solid-phase C reservoir, with a rate-limited solubilization coefficient. Aging of SOC in this way is more akin to selective preservation than to mineral protection and, while such a mechanism may be at play in many soils, mineral protection is thought to be at least as important. As such, our results indicate that additional parameterization is required to reproduce the heterogeneity of carbon transit times that result from organo-mineral interactions. These efforts show the promise of a modeling approach where the varied transit time of soil C emerges from the dynamic physical and hydrologic properties of the model rather than from the a priori assignment of operationally defined pools.

Geochimica et Cosmochimica Acta↗

Hydrogen isotope behavior during rhyolite glass hydration under hydrothermal conditions

The diffusion of molecular water (H 2 O m ) from the environment into volcanic glass can hydrate the glass up to several wt% at low temperature over long timescales. During this process, the water imprints its hydrogen isotope composition (δD H2O ) to the glass (δD gl ) offset by a glass-H 2 O fractionation factor (ΔD gl-H2O = δD gl – δD H2O ) which is approximately −33‰ at Earth surface temperatures. Glasses hydrate much more rapidly at higher, sub-magmatic temperatures as they interact with H 2 O during eruption, transport, and emplacement. To aid in the interpretation of δD gl in natural samples, we present hydrogen isotope results from vapor hydration experiments conducted at 175–375 °C for durations of hours to months using natural volcanic glasses. The results can be divided into two thermal regimes: above 250 °C and below 250 °C. Lower temperature experiments yield raw ΔD gl-H2O values in the range of −33 ± 11‰. Experiments at 225 °C using both positive and negative initial ΔD gl-H2O values converge on this range of values, suggesting this range represents the approximate equilibrium fractionation for H isotopes between glass and H 2 O vapor (10 3 lnα gl-H2O ) below 250 °C. Variation in ΔD gl-H2O (−33 ± 11‰) between different experiments and glasses may arise from incomplete hydration, analytical uncertainty, differences in glass chemistry, and/or subordinate kinetic isotope effects. Experiments above 250 °C yield unexpectedly low δD gl values with ΔD gl-H2O values of ≤–85‰. While alteration alone is incapable of explaining the data, these run products have more extensive surface alteration and are not interpreted to reflect equilibrium fractionation between glass and H 2 O vapor. Fourier transform infrared spectroscopy (FTIR) shows that glass can hydrate with as much as 5.9 wt% H 2 O m and 1.0 wt% hydroxl (OH − ) in the highest P-T experiment at 375 °C and 21.1 MPa. Therefore, we employ a 1D isotope diffusion–reaction model of glass hydration to evaluate the roles of equilibrium fractionation, isotope diffusion, water speciation reactions internal to the glass, and changing boundary conditions (e.g. alteration and dissolution). At lower temperatures, the best fitting model results to experimental data for low silica rhyolite (LSR) glasses require only an equilibrium fractionation factor and yield 10 3 lnα gl-H2O values of −33‰ ± 5‰ and −25‰ ± 5‰ at 175 °C and 225 °C, respectively. At higher temperatures, ΔD gl-H2O is dominated by boundary layer effects during glass hydration and glass surface alteration. The modeled bulk δD gl value is highly responsive to changes in the δD gl boundary condition regardless of the magnitude of other kinetic effects. Observed glass dissolution and surficial secondary mineral formation are likely to impose a disequilibrium boundary layer that drives extreme δD gl fractionation with progressive glass hydration. These results indicate that the observed ΔD gl-H2O of ∼−33 ± 11‰ can be cautiously applied as an equilibrium 10 3 lnα gl-H2O value to natural silicic glasses hydrated below 250 °C to identify hydration sources. This approximate ΔD gl-H2O may be applicable to even higher temperature glasses hydrated on short timescales (of seconds to minutes) in phreatomagmatic or submarine eruptions before H 2 O in the glass is primarily affected by boundary layer effects associated with alteration on the glass surface.

Geochimica et Cosmochimica Acta↗

The dynamic influence of subsurface geological processes on the assembly and diversification of thermophilic microbial communities in continental hydrothermal systems

An accepted paradigm of hydrothermal systems is the process of phase separation, or boiling, of a deep, homogeneous hydrothermal fluid as it ascends through the subsurface resulting in gas rich and gas poor fluids. While phase separation helps to explain first-order patterns in the chemistry and biology of a hot spring’s surficial expression, we know little about the subsurface architecture beneath “phase-separated” pools and the timescales over which phase separation processes occur. Essentially, we have a two-dimensional understanding of a four-dimensional process. By combining geophysical, geochemical, isotopic, and microbiological measurements of two adjacent phase-separated hot springs in Norris Geyser Basin, Yellowstone National Park, we provide a four-dimensional assessment of phase separation processes and their biological manifestation. We uniquely show that Yellowstone’s hydrothermal waters originate from a deep sedimentary aquifer and that both meteoric recharge and shallow reactive transport processes are required to establish the geobiological feedbacks that drive bimodal distributions in the geochemical and microbial composition of hot springs. Specifically, over periods of tens of years, gas-enriched fluids containing volcanic sulfide mix with meteoric waters resulting in microbially-mediated production of sulfuric acid by thermoacidophilic Archaea in the near subsurface. In contrast, over periods of hundreds of years, anoxic residual liquid rises to the surface where it is infused with atmospheric gas fostering Archaea and Bacteria that are largely dependent on oxygen. As such, our results provide formative insight into the causative links between subsurface geological processes, the development of geochemical fluids, and the assembly and diversification of thermophilic microbial communities in hydrothermal systems.

Wyoming↗

Quantifying form resistance is essential for estimating summer low and bankfull flow from stream survey channel morphology

Reliable estimates of low flow and flood discharge at ungaged locations are required for evaluating stream flow alteration, designing culverts and stream crossings, and interpreting regional surveys of habitat and biotic condition. Very few stream gaging stations are located on small, remote streams, which typically have complex channel morphology. Adequate gaging is also lacking on larger streams that are remote, smaller than those typically gaged, or have channel morphology not conducive to installation of gages. Complex channels typically contain large scale hydraulic roughness elements that dominate flow patterns (i.e., form roughness), making it difficult to measure channel cross-section area and water velocity, or to measure channel volume even where discharge is known. In channels with large channel form roughness, it is equally difficult to estimate discharge using commonly applied equations based on slope and channel dimensions or basin area. We employed a novel approach that explicitly accounts for hydraulic resistance from large wood and riffle-pool morphology (form roughness) in calculating low flow and bankfull discharge from stream and river physical habitat data collected from 4229 stream and river sites in the conterminous US (CONUS) sampled in 2008–9 and 2013–14 as part of the US Environmental Protection Agency's N ational R ivers and S treams A ssessment (NRSA). Hydraulic resistance derived from form roughness clearly dominated resistance derived from bed particles (particle resistance) during summer low flows in wadeable streams across the spectrum of channel slopes and substrate sizes smaller than boulders. Under bankfull conditions, the influence of form resistance relative to particle resistance was diminished, but form resistance still dominated except in large low gradient rivers lacking complex channels, and in streams or rivers with boulder-size bed particles. We validated our hydraulic resistance estimates by comparing measured discharges with calculated discharges that used those hydraulic resistance estimates along with measured NRSA channel morphology data. Morphology-based summer discharge (low flow) estimates and direct field measurements of discharge in 2333 wadeable CONUS streams showed reasonable agreement (median difference < 2.5×) for discharges ranging from 3.6 × 10 −5 to 123 m 3 /s and drainage areas of 0.12 to 171,000 km 2 . In a subset of 759 of NRSA's larger wadeable stream and non-wadeable river sites where nearby U.S. Geological Survey (USGS) gage data were available and adequate, our morphology-based summer discharge estimates agreed fairly well (median difference < 2.0×) with USGS 20-yr average August mean flows ranging from 0.003 to 16,000 m 3 /s. Similarly, morphology-based estimates of bankfull flow ranging from 0.3 to 100,000 m 3 /s agreed reasonably well with the 1.5-yr recurrence interval flood in these gaged sites (median deviation <2.2×). These findings demonstrate the importance of quantifying flow resistance from large-scale form roughness features in natural channels and provide a novel approach for estimating discharge from widely available survey data. This will allow examination of discharge and its ecological influence across the full range of stream and river sizes sampled by NRSA or other synoptic surveys where comprehensive measures of biota, physical habitat, and chemistry are also made. Although these morphology-based estimates exhibit some variability, they are adequate for examining regional patterns in discharge and flow alteration and their association with instream biota and anthropogenic disturbances, providing summer low and bankfull flow information where reliable estimates are lacking.

Geomorphology↗

Controls on the Karaha-Telaga Bodas geothermal reservoir, Indonesia

Karaha-Telaga Bodas is a partially vapor-dominated, fracture-controlled geothermal system located adjacent to Galunggung Volcano in western Java, Indonesia. The geothermal system consists of: (1) a caprock, ranging from several hundred to 1600 m in thickness, and characterized by a steep, conductive temperature gradient and low permeability; (2) an underlying vapor-dominated zone that extends below sea level; and (3) a deep liquid-dominated zone with measured temperatures up to 353 ??C. Heat is provided by a tabular granodiorite stock encountered at about 3 km depth. A structural analysis of the geothermal system shows that the effective base of the reservoir is controlled either by the boundary between brittle and ductile deformational regimes or by the closure and collapse of fractures within volcanic rocks located above the brittle/ductile transition. The base of the caprock is determined by the distribution of initially low-permeability lithologies above the reservoir; the extent of pervasive clay alteration that has significantly reduced primary rock permeabilities; the distribution of secondary minerals deposited by descending waters; and, locally, by a downward change from a strike-slip to an extensional stress regime. Fluid-producing zones are controlled by both matrix and fracture permeabilities. High matrix permeabilities are associated with lacustrine, pyroclastic, and epiclastic deposits. Productive fractures are those showing the greatest tendency to slip and dilate under the present-day stress conditions. Although the reservoir appears to be in pressure communication across its length, fluid, and gas chemistries vary laterally, suggesting the presence of isolated convection cells. ?? 2006 CNR.

Geothermics↗

The distribution and composition of REE-bearing minerals in placers of the Atlantic and Gulf coastal plains, USA

Rare earth element (REE) resources are currently of great interest because of their importance as raw materials for high-technology manufacturing. The REE-phosphates monazite (light REE enriched) and xenotime (heavy REE enriched) resist weathering and can accumulate in placer deposits as part of the heavy mineral assemblage. The Atlantic and Gulf coastal plains of the southeastern United States are known to host heavy mineral deposits with economic concentrations of zircon, ilmenite and rutile. This study provides a perspective on the distribution and composition of REE phosphate minerals in the region. The elemental chemistry and mineralogy of sands and associated heavy-mineral assemblages from new and archived sediment samples across the coastal plains are examined, along with phase-specific compositions of monazite, xenotime and zircon. Both monazite and xenotime are present across the coastal plains. The phase-specific compositions allow monazite content to be estimated using La as a geochemical proxy. Similarly, both Y and Yb are geochemical proxies for xenotime, but their additional presence in zircon and monazite require a correction to prevent overestimation of xenotime content. Applying this correction, maps of monazite and xenotime content across the coastal plains were generated using sample coverage from the National Geochemical Database (NGS) and National Uranium Resource Evaluation (NURE). The NGS and NURE approach of sampling stream sediments in small watersheds links samples to nearby lithologies. The results show an approximately 40 km-wide band of primarily Cretaceous, marine sediments bordering the Piedmont province from North Carolina to Alabama in which monazite and xenotime content are relatively high (up to 4.4 wt. % in < 150 &mu;m bulk sediment). Strong correlations between concentrations of the two phases were found, with estimated monazite:xenotime ratios ranging approximately 6:1 to 12:1 depending upon the dataset analyzed. From a resource perspective, xenotime correlation with monazite indicates a heavy REE potential in coastal plain placers, and exploration may be warranted within the identified coastal plain band along the boundary of the Piedmont region.

Journal of Geochemical Exploration↗

Ore mineralization in the Mofete and San Vito geothermal fields, Campi Flegrei volcanic complex, Naples, Italy

The Mofete and San Vito geothermal fields, located west of Naples, Italy, are part of the Campi Flegrei volcanic complex. In the 1970s, exploratory wells were drilled to a depth of ~3000 m in an attempt to locate high-enthalpy fluids for potential power production. Drill core samples from Mofete wells (MF1, MF2, and MF5) and from San Vito wells (SV1 and SV3) contain authigenic ore mineralization. Pyrite, pyrrhotite, and galena are abundant. Less common are chalcopyrite, sphalerite, arsenopyrite, and scheelite; rare are millerite, violarite, native bismuth, tellurobismuthite, cassiterite, molybdenite, and acanthite. Mineral chemistry was determined by electron microprobe wavelength dispersive spectroscopy aided by a scanning electron microscope equipped with energy-dispersive spectroscopy. The mineral assemblage suggests a low sulfidation environment and the absence of pyrrhotite in the MF1 well and upper part of the SV1 well indicates variable sulfur activity. Both molybdenite and scheelite were identified in samples SV1–2860 and SV3–2353 and scheelite in the SV3 well is zoned with variable Mo 6+ content; low Mo 6+ zones show blue cathodoluminescence, whereas, zones with high Mo 6+ content are yellow to brown. Zoned scheelite and the occurrence of both Mo-bearing minerals attest to the variability of ƒO 2 and ƒS 2 in the geothermal fluid.

Naples↗

The dilemma of the Jiaodong gold deposits: Are they unique?

The ca. 126–120 Ma Au deposits of the Jiaodong Peninsula, eastern China, define the country's largest gold province with an overall endowment estimated as >3000 t Au. The vein and disseminated ores are hosted by NE- to NNE-trending brittle normal faults that parallel the margins of ca. 165–150 Ma, deeply emplaced, lower crustal melt granites. The deposits are sited along the faults for many tens of kilometers and the larger orebodies are associated with dilatational jogs. Country rocks to the granites are Precambrian high-grade metamorphic rocks located on both sides of a Triassic suture between the North and South China blocks. During early Mesozoic convergent deformation, the ore-hosting structures developed as ductile thrust faults that were subsequently reactivated during Early Cretaceous “Yanshanian” intracontinental extensional deformation and associated gold formation. Classification of the gold deposits remains problematic. Many features resemble those typical of orogenic Au including the linear structural distribution of the deposits, mineralization style, ore and alteration assemblages, and ore fluid chemistry. However, Phanerozoic orogenic Au deposits are formed by prograde metamorphism of accreted oceanic rocks in Cordilleran-style orogens. The Jiaodong deposits, in contrast, formed within two Precambrian blocks approximately 2 billion years after devolatilization of the country rocks, and thus require a model that involves alternative fluid and metal sources for the ores. A widespread suite of ca. 130–123 Ma granodiorites overlaps temporally with the ores, but shows a poor spatial association with the deposits. Furthermore, the deposit distribution and mineralization style is atypical of ores formed from nearby magmas. The ore concentration requires fluid focusing during some type of sub-crustal thermal event, which could be broadly related to a combination of coeval lithospheric thinning, asthenospheric upwelling, paleo-Pacific plate subduction, and seismicity along the continental-scale Tan-Lu fault. Possible ore genesis scenarios include those where ore fluids were produced directly by the metamorphism of oceanic lithosphere and overlying sediment on the subducting paleo-Pacific slab, or by devolatilization of an enriched mantle wedge above the slab. Both the sulfur and gold could be sourced from either the oceanic sediments or the serpentinized mantle. A better understanding of the architecture of the paleo-Pacific slab during Early Cretaceous below the eastern margin of China is essential to determination of the validity of possible models.

Shandong Province↗

Micrometer-scale U–Pb age domains in eucrite zircons, impact re-setting, and the thermal history of the HED parent body

Meteoritic zircons are rare, but some are documented to occur in asteroidal meteorites, including those of the howardite&ndash;eucrite&ndash;diogenite (HED) achondrite clan (Rubin, A. [1997]. Meteorit. Planet. Sci. 32, 231&ndash;247). The HEDs are widely considered to originate from the Asteroid 4 Vesta. Vesta and the other large main belt asteroids record an early bombardment history. To explore this record, we describe sub-micrometer distributions of trace elements (U, Th) and 235,238U&ndash;207,206Pb ages from four zircons (>7&ndash;40 &mu;m &empty;) separated from bulk samples of the brecciated eucrite Millbillillie. Ultra-high resolution (&sim;100 nm) ion microprobe depth profiles reveal different zircon age domains correlative to mineral chemistry and to possible impact scenarios. Our new U&ndash;Pb zircon geochronology shows that Vesta&rsquo;s crust solidified within a few million years of Solar System formation (4561 &plusmn; 13 Ma), in good agreement with previous work (e.g. Carlson, R.W., Lugmair, G.W. [2000]. Timescales of planetesimal formation and differentiation based on extinct and extant radioisotopes. In: Canup, R., Righter, K. (Eds.), Origin of the Earth and Moon. University of Arizona Press, Tucson, pp. 25&ndash;44). Younger zircon age domains (ca. 4530 Ma) also record crustal processes, but these are interpreted to be exogenous because they are well after the effective extinction of 26Al (t1/2 = 0.72 Myr). An origin via impact-resetting was evaluated with a suite of analytical impact models. Output shows that if a single impactor was responsible for the ca. 4530 Ma zircon ages, it had to have been ⩾10 km in diameter and at high enough velocity (>5 km s&minus;1) to account for the thermal field required to re-set U&ndash;Pb ages. Such an impact would have penetrated at least 10 km into Vesta&rsquo;s crust. Later events at ca. 4200 Ma are documented in HED apatite 235,238U&ndash;207,206Pb ages (Zhou, Q. et al. [2011]. Early basaltic volcanism and Late Heavy Bombardment on Vesta: U&ndash;Pb ages of small zircons and phosphates in eucrites. Lunar Planet. Sci. 42. Abstract #2575) and 40&ndash;39Ar age spectra (Bogard, D.D. [2011]. Chem. Erde 71, 207&ndash;226). Yet younger ages, including those coincident with the Late Heavy Bombardment (LHB; ca. 3900 Ma), are absent from Millbillillie zircon. This is attributable to primordial changes to the velocity distributions of impactors in the asteroid belt, and differences in mineral closure temperatures (Tc zircon ≫ apatite).

Icarus↗

Sequestration of non-pure carbon dioxide streams in iron oxyhydroxide-containing saline repositories

Iron oxyhydroxide, goethite (&alpha;-FeOOH), was evaluated as a potential formation mineral reactant for trapping CO 2 in a mineral phase such as siderite (FeCO 3 ), when a mixture of CO 2 -SO 2 flue gas is injected into a saline aquifer. Two thermodynamic simulations were conducted, equilibrating a CO 2 -SO 2 fluid mixture with a NaCl-brine and Fe-rich rocks at 150 &deg;C and 300 bar. The modeling studies evaluated mineral and fluid composition at equilibrium and the influence of pH buffering in the system. Results show siderite precipitates both in the buffered and unbuffered system; however, the presence of an alkaline pH buffer enhances the stability of the carbonate. Based on the model, an experiment was designed to compare with thermodynamic predictions. A CO 2 -SO 2 gas mixture was reacted in 150 ml of NaCl-NaOH brine containing 10 g of goethite at 150 &deg;C and 300 bar for 24 days. Mineralogical and brine chemistry confirmed siderite as the predominant reaction product in the system. Seventy-six mg of CO 2 are sequestered in siderite per 10 g of goethite.

International Journal of Greenhouse Gas Control↗

Abundance and distribution of selected elements in soils, stream sediments, and selected forage plants from desert tortoise habitats in the Mojave and Colorado deserts, USA

A baseline and background chemical survey was conducted in southeastern California, USA, to identify potential sources of toxicants in natural and anthropogenically-altered habitats of the threatened desert tortoise (Gopherus agassizii). Soil, stream sediment, and plant samples were collected from six tortoise habitat study areas in the Mojave and Colorado deserts and analysed for up to 66 different elements. The chemical analyses provided new information on the abundances and distributions of selected elements in this region. Soil, stream-sediment, and plant analyses showed distinct variations in bulk chemistries from locality to locality. Variations were, in general, consistent with the many types of exposed rock units in the region, their highly variable bulk mineralogies, and chemical contents. Of elements in soils that might have been toxic to tortoises, only As seemed to be anomalous region-wide. Some soil and plant anomalies were clearly anthropogenic. In the Rand and Atolia mining districts, soil anomalies for As, Au, Cd, Hg, Sb, and(or) W and plant anomalies for As, Sb, and(or) W extend as far as ???15 km outward from the present area of mining; soils containing anomalous Hg were found at least 6 km away from old piles of tailings. The anomalous concentrations of As and Hg may have been the source of elevated levels of these elements found in ill tortoises from the region. In the Goldstone mining district, soil anomalies extended several km from the mining area. These areas probably represented anthropogenic surface contamination of dust redistributed by wind, vehicles, and rainfall. One of two study areas transected by a paved road (Chemehuevi Valley) showed weakly elevated levels of Pb, which extended as far as ???22 m from the pavement edge and were probably related to vehicle exhaust. No soil or plant samples from historically used military areas (Goldstone, Goffs, Chemehuevi Valley, Chuckwalla Bench) contained anomalous concentrations of the elements As, Cu, Cr, Fe, Pb, or Zn that could be ascribed to military maneuvers, vehicles, or ordnance. For future studies, the distribution and abundance of elements in the tortoise forage plants need to be evaluated for the respective roles of dust and systemic uptake. Additional chemical data from tortoise necropsies and nutritional studies are needed to determine the effects of potentially toxic elements in tortoise habitats on their health. ?? 2006 Elsevier Ltd. All rights reserved.

Journal of Arid Environments↗