SEARCH · USGS Science

Results for “Geochimica et Cosmochimica Acta”

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

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26Linked to original sources

USGS2015

Global patterns and environmental controls of perchlorate and nitrate co-occurrence in arid and semi-arid environments

Natural perchlorate (ClO 4 − ) is of increasing interest due to its wide-spread occurrence on Earth and Mars, yet little information exists on the relative abundance of ClO 4 − compared to other major anions, its stability, or long-term variations in production that may impact the observed distributions. Our objectives were to evaluate the occurrence and fate of ClO 4 − in groundwater and soils/caliche in arid and semi-arid environments (southwestern United States, southern Africa, United Arab Emirates, China, Antarctica, and Chile) and the relationship of ClO 4 − to the more well-studied atmospherically deposited anions NO 3 − and Cl − as a means to understand the prevalent processes that affect the accumulation of these species over various time scales. ClO 4 − is globally distributed in soil and groundwater in arid and semi-arid regions on Earth at concentrations ranging from 10 −1 to 10 6 μg/kg. Generally, the ClO 4 − concentration in these regions increases with aridity index, but also depends on the duration of arid conditions. In many arid and semi-arid areas, NO 3 − and ClO 4 − co-occur at molar ratios (NO 3 − /ClO 4 − ) that vary between ∼10 4 and 10 5 . We hypothesize that atmospheric deposition ratios are largely preserved in hyper-arid areas that support little or no biological activity (e.g. plants or bacteria), but can be altered in areas with more active biological processes including N 2 fixation, N mineralization, nitrification, denitrification, and microbial ClO 4 − reduction, as indicated in part by NO 3 − isotope data. In contrast, much larger ranges of Cl − /ClO 4 − and Cl − /NO 3 − ratios indicate Cl − varies independently from both ClO 4 − and NO 3 − . The general lack of correlation between Cl − and ClO 4 − or NO 3 − implies that Cl − is not a good indicator of co-deposition and should be used with care when interpreting oxyanion cycling in arid systems. The Atacama Desert appears to be unique compared to all other terrestrial locations having a NO 3 − /ClO 4 − molar ratio ∼10 3 . The relative enrichment in ClO 4 − compared to Cl − or NO 3 − and unique isotopic composition of Atacama ClO 4 − may reflect either additional in-situ production mechanism(s) or higher relative atmospheric production rates in that specific region or in the geological past. Elevated concentrations of ClO 4 − reported on the surface of Mars, and its enrichment with respect to Cl − and NO 3 − , could reveal important clues regarding the climatic, hydrologic, and potentially biologic evolution of that planet. Given the highly conserved ratio of NO 3 − /ClO 4 − in non-biologically active areas on Earth, it may be possible to use alterations of this ratio as a biomarker on Mars and for interpreting major anion cycles and processes on both Mars and Earth, particularly with respect to the less-conserved NO 3 − pool terrestrially.

Geochimica et Cosmochimica Acta↗

USGS2015

Compound-specific sulfur isotope analysis of thiadiamondoids of oils from the Smackover Formation, USA

Thiadiamondoids (TDs) are diamond-like compounds with a sulfide bond located within the cage structure. These compounds were suggested as a molecular proxy for the occurrence and extent of thermochemical sulfate reduction (TSR). Compound-specific sulfur-isotope analysis of TDs may create a multi-parameter system, based on molecular and δ34S values that may be sensitive over a wider range of TSR and thermal maturation stages. In this study, we analyzed a suite of 12 Upper Jurassic oil and condensate samples generated from source rocks in the Smackover Formation to perform a systematic study of the sulfur isotope distribution in thiadiamondoids (one and two cages). For comparison we measured the δ34S composition of benzothiophenes (BTs) and dibenzothiophenes (DBTs). We also conducted pyrolysis experiments with petroleum and model compounds to have an insight into the formation mechanisms of TDs. The δ34S of the TDs varied significantly (ca 30&permil;) between the different oils depending on the degree of TSR alteration. The results showed that within the same oil, the one-cage TDs were relatively uniform, with 34S enriched values similar to those of the coexisting BTs. The two-cage TDs had more variable δ34S values that range from the δ34S values of BTs to those of the DBTs, but with general 34S depletion relative to one cage TDs. Hydrous pyrolysis experiments (360 °C, 40 h) with either CaSO4 or elemental S (equivalent S molar concentrations) and adamantane as a model compound demonstrate the formation of one cage TDs in relatively low yields (<0.2%). Higher concentrations of TDs were observed in the elemental sulfur experiments, most likely because of the higher rates of reaction with adamantane under these experimental conditions. These results show that the formation of TDs is not exclusive to TSR reactions, and that they can also form by reaction with reduced S species apart from sulfate reduction, though at low yields. Oxygenated compounds, most notably 2-thiaadamantanone and phenol, were also formed during these pyrolysis experiments. This may represent the first stage in the formation of sulfurized compounds and the oxidation of organic matter under TSR conditions. Pyrolysis experiments with elemental S and a TD-enriched oil showed that the δ34S values of the TDs did not change, whereas the BTs did change significantly. It is therefore concluded that TDs do not exchange S atoms with coexisting inorganic reduced sulfur species. They can only change their δ34S values via addition of newly generated TDs that form predominantly during TSR. We therefore suggest that TDs will preserve their δ34S values even under high-temperature reservoir conditions and will reflect the original sulfates δ34S value. The combination of TDs, BTs, and DBTs δ34S values and concentrations allowed for a more reliable detection of the occurrence and extent of TSR than either proxy alone. It showed that except for two oils, all of the oils that were measured in this study were affected by TSR or TSR-sourced H2S, to some degree. It is still not known if some of the oils with the lower concentrations of TDs and enriched δ34S values (close to sulfate minerals) were affected by TSR or by a secondary charge of 34S-enriched H2S.

Geochimica et Cosmochimica Acta↗

USGS2015

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°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↗

USGS2016

A thermodynamical model for the surface tension of silicate melts in contact with H 2 O gas

Surface tension plays an important role in the nucleation of H 2 O gas bubbles in magmatic melts and in the time-dependent rheology of bubble-bearing magmas. Despite several experimental studies, a physics based model of the surface tension of magmatic melts in contact with H 2 O is lacking. This paper employs gradient theory to develop a thermodynamical model of equilibrium surface tension of silicate melts in contact with H 2 O gas at low to moderate pressures. In the last decades, this approach has been successfully applied in studies of industrial mixtures but never to magmatic systems. We calibrate and verify the model against literature experimental data, obtained by the pendant drop method, and by inverting bubble nucleation experiments using the Classical Nucleation Theory (CNT). Our model reproduces the systematic decrease in surface tension with increased H 2 O pressure observed in the experiments. On the other hand, the effect of temperature is confirmed by the experiments only at high pressure. At atmospheric pressure, the model shows a decrease of surface tension with temperature. This is in contrast with a number of experimental observations and could be related to microstructural effects that cannot be reproduced by our model. Finally, our analysis indicates that the surface tension measured inverting the CNT may be lower than the value measured by the pendant drop method, most likely because of changes in surface tension controlled by the supersaturation.

Geochimica et Cosmochimica Acta↗

USGS2016

Mercury transformation and release differs with depth and time in a contaminated riparian soil during simulated flooding

Riparian soils are an important environment in the transport of mercury in rivers and wetlands, but the biogeochemical factors controlling mercury dynamics under transient redox conditions in these soils are not well understood. Mercury release and transformations in the Oa and underlying A horizons of a contaminated riparian soil were characterized in microcosms and an intact soil core under saturation conditions. Pore water dynamics of total mercury (HgT), methylmercury (MeHg), and dissolved gaseous mercury (Hg0(aq)) along with selected anions, major elements, and trace metals were characterized across redox transitions during 36 d of flooding in microcosms. Next, HgT dynamics were characterized over successive flooding (17 d), drying (28 d), and flooding (36 d) periods in the intact core. The observed mercury dynamics exhibit depth and temporal variability. At the onset of flooding in microcosms (1–3 d), mercury in the Oa horizon soil, present as a combination of ionic mercury (Hg(II)) bound to thiol groups in the soil organic matter (SOM) and nanoparticulate metacinnabar (b-HgS), was mobilized with organic matter of high molecular weight. Subsequently, under anoxic conditions, pore water HgT declined coincident with sulfate (3–11 d) and the proportion of nanoparticulate b-HgS in the Oa horizon soil increased slightly. Redox oscillations in the intact Oa horizon soil exhausted the mobile mercury pool associated with organic matter. In contrast, mercury in the A horizon soil, present predominantly as nanoparticulate b-HgS, was mobilized primarily as Hg0(aq) under strongly reducing conditions (5–18 d). The concentration of Hg0(aq) under dark reducing conditions correlated positively with byproducts of dissimilatory metal reduction (P(Fe,Mn)). Mercury dynamics in intact A horizon soil were consistent over two periods of flooding, indicating that nanoparticulate b-HgS was an accessible pool of mobile mercury over recurrent reducing conditions. The concentration of MeHg increased with flooding time in both the Oa and A horizon pore waters. Temporal changes in pore water constituents (iron, manganese, sulfate, inorganic carbon, headspace methane) all implicate microbial control of redox transitions. The mobilization of mercury in multiple forms, including HgT associated with organic matter, MeHg, and Hg0(aq), to pore waters during periodic soil flooding may contribute to mercury releases to adjacent surface waters and the recycling of the legacy mercury to the atmosphere.

Tennessee↗

USGS2016

Novel insights from NMR spectroscopy into seasonal changes in the composition of dissolved organic matter exported to the Bering Sea by the Yukon River

Seasonal (spring freshet, summer–autumn, and winter) variability in the chemical composition of dissolved organic matter (DOM) from the Yukon River was determined using advanced one- and two-dimensional (2D) solid-state NMR spectroscopy, coupled with isotopic measurements and UV–visible spectroscopy. Analyses were performed on two major DOM fractions, the hydrophobic organic acid (HPOA) and transphilic organic acid (TPIA) fractions obtained using XAD resins. Together these two fractions comprised 64–74% of the total DOM. Carboxyl-rich alicyclic molecules (CRAM) accounted for the majority of carbon atoms in the HPOA (63–77%) and TPIA (54–78%) samples, and more so in winter and summer than in spring samples. 2D and selective NMR data revealed association of abundant nonprotonated O-alkyl and quaternary alkyl C (OC np , OC np O and C q , 13–17% of HPOA and 15–20% of TPIA) and isolated O–CH structures with CRAM, which were not recognized in previous studies. Spectral editing and 2D NMR allowed for the discrimination of carbohydrate-like O-alkyl C from non-carbohydrate O-alkyl C. Whereas two spring freshet TPIA samples contained carbohydrate clusters such as carboxylated carbohydrates (16% and 26%), TPIA samples from other seasons or HPOA samples mostly had small amounts (<8%) of sugar rings dispersed in a nonpolar alkyl environment. Though nonprotonated aromatic C represented the largest fraction of aromatic C in all HPOA/TPIA isolates, only a small fraction (&sim;5% in HPOA and 3% in TPIA) was possibly associated with dissolved black carbon. Our results imply a relatively stable portion of DOM exported by the Yukon River across different seasons, due to the predominance of CRAM and their associated nonprotonated C–O and O–C–O structures, and elevated reactivity (bio- and photo-lability) of spring DOM due to the presence of terrestrial inputs enriched in carbohydrates and aromatic structures.

Alaska, Yukon↗

USGS2016

Deposition, accumulation, and alteration of Cl−, NO3−, ClO4− and ClO3− salts in a hyper-arid polar environment: Mass balance and isotopic constraints

The salt fraction in permafrost soils/sediments of the McMurdo Dry Valleys (MDV) of Antarctica can be used as a proxy for cold desert geochemical processes and paleoclimate reconstruction. Previous analyses of the salt fraction in MDV permafrost soils have largely been conducted in coastal regions where permafrost soils are variably affected by aqueous processes and mixed inputs from marine and stratospheric sources. We expand upon this work by evaluating permafrost soil/sediments in University Valley, located in the ultraxerous zone where both liquid water transport and marine influences are minimal. We determined the abundances of Cl − , NO 3 − , ClO 4 − and ClO 3 − in dry and ice-cemented soil/sediments, snow and glacier ice, and also characterized Cl − and NO 3 − isotopically. The data are not consistent with salt deposition in a sublimation till, nor with nuclear weapon testing fall-out, and instead point to a dominantly stratospheric source and to varying degrees of post depositional transformation depending on the substrate, from minimal alteration in bare soils to significant alteration (photodegradation and/or volatilization) in snow and glacier ice. Ionic abundances in the dry permafrost layer indicate limited vertical transport under the current climate conditions, likely due to percolation of snowmelt. Subtle changes in ClO 4 − /NO 3 − ratios and NO 3 − isotopic composition with depth and location may reflect both transport related fractionation and depositional history. Low molar ratios of ClO 3 − /ClO 4 − in surface soils compared to deposition and other arid systems suggest significant post depositional loss of ClO 3 − , possibly due to reduction by iron minerals, which may have important implications for oxy-chlorine species on Mars. Salt accumulation varies with distance along the valley and apparent accumulation times based on multiple methods range from ∼10 to 30 kyr near the glacier to 70–200 kyr near the valley mouth. The relatively young age of the salts and relatively low and homogeneous anion concentrations in the ice-cemented sediments point to either a mechanism of recent salt removal, or to relatively modern permafrost soils (<1 million years). Together, our results show that near surface salts in University Valley serve as an end-member of stratospheric sources not subject to biological processes or extensive remobilization.

Geochimica et Cosmochimica Acta↗

USGS2016

Magma transport and olivine crystallization depths in Kīlauea’s East Rift Zone inferred from experimentally rehomogenized melt inclusions

Concentrations of H 2 O and CO 2 in olivine-hosted melt inclusions can be used to estimate crystallization depths for the olivine host. However, the original dissolved CO 2 concentration of melt inclusions at the time of trapping can be difficult to measure directly because in many cases substantial CO 2 is transferred to shrinkage bubbles that form during post-entrapment cooling and crystallization. To investigate this problem, we heated olivine from the 1959 Kīlauea Iki and 1960 Kapoho (Hawai‘i) eruptions in a 1-atm furnace to temperatures above the melt inclusion trapping temperature to redissolve the CO 2 in shrinkage bubbles. The measured CO 2 concentrations of the experimentally rehomogenized inclusions (⩽590 ppm for Kīlauea Iki [n=10]; ⩽880 ppm for Kapoho, with one inclusion at 1863 ppm [n=38]) overlap with values for naturally quenched inclusions from the same samples, but experimentally rehomogenized inclusions have higher within-sample median CO 2 values than naturally quenched inclusions, indicating at least partial dissolution of CO 2 from the vapor bubble during heating. Comparison of our data with predictions from modeling of vapor bubble formation and published Raman data on the density of CO 2 in the vapor bubbles suggests that 55-85% of the dissolved CO 2 in the melt inclusions at the time of trapping was lost to post-entrapment shrinkage bubbles. Our results combined with the Raman data demonstrate that olivine from the early part of the Kīlauea Iki eruption crystallized at <6 km depth, with the majority of olivine in the 1-3 km depth range. These depths are consistent with the interpretation that the Kīlauea Iki magma was supplied from Kīlauea’s summit magma reservoir (&sim;2-5 km depth). In contrast, olivine from Kapoho, which was the rift zone extension of the Kīlauea Iki eruption, crystallized over a much wider range of depths (&sim;1-16 km). The wider depth range requires magma transport during the Kapoho eruption from deep beneath the summit region and/or from deep beneath Kīlauea’s east rift zone. The deeply derived olivine crystals and their host magma mixed with stored, more evolved magma in the rift zone, and the mixture was later erupted at Kapoho.

Hawaii↗

USGS2016

Study of thermochemical sulfate reduction mechanism using compound specific sulfur isotope analysis

The sulfur isotopic fractionation associated with the formation of organic sulfur compounds (OSCs) during thermochemical sulfate reduction (TSR) was studied using gold-tube pyrolysis experiments to simulate TSR. The reactants used included n -hexadecane ( n -C 16 ) as a model organic compound with sulfate, sulfite, or elemental sulfur as the sulfur source. At the end of each experiment, the S-isotopic composition and concentration of remaining sulfate, H 2 S, benzothiophene, dibenzothiophene, and 2-phenylthiophene (PT) were measured. The observed S-isotopic fractionations between sulfate and BT, DBT, and H 2 S in experimental simulations of TSR correlate well with a multi-stage model of the overall TSR process. Large kinetic isotope fractionations occur during the first, uncatalyzed stage of TSR, 12.4&permil; for H 2 S and as much as 22.2&permil; for BT. The fractionations decrease as the H 2 S concentration increases and the reaction enters the second, catalyzed stage. Once all of the oxidizable hydrocarbons have been consumed, sulfate reduction ceases and equilibrium partitioning then dictates the fractionation between H 2 S and sulfate (&sim;17&permil;). Experiments involving sparingly soluble CaSO 4 show that during the second catalytic phase of TSR the rate of sulfate reduction exceeds that of sulfate dissolution. In this case, there is no apparent isotopic fractionation between source sulfate and generated H 2 S, as all of the available sulfate is effectively reduced at all reaction times. When CaSO 4 is replaced with fully soluble Na 2 SO 4 , sulfate dissolution is no longer rate limiting and significant S-isotopic fractionation is observed. This supports the notion that CaSO 4 dissolution can lead to the apparent lack of fractionation between H 2 S and sulfate produced by TSR in nature. The S-isotopic composition of individual OSCs record information related to geochemical reactions that cannot be discerned from the δ 34 S values obtained from bulk phases such as H 2 S, oil, and sulfate minerals, and provide important mechanistic details about the overall TSR process.

Geochimica et Cosmochimica Acta↗

USGS2016

Isotopic constraints on the genesis and evolution of basanitic lavas at Haleakala, Island of Maui, Hawaii

To understand the dynamics of solid mantle upwelling and melting in the Hawaiian plume, we present new major and trace element data, Nd, Sr, Hf, and Pb isotopic compositions, and 238 U– 230 Th– 226 Ra and 235 U– 231 Pa– 227 Ac activities for 13 Haleakala Crater nepheline normative basanites with ages ranging from ∼900 to 4100 yr B.P. These basanites of the Hana Volcanics exhibit an enrichment in incompatible trace elements and a more depleted isotopic signature than similarly aged Hawaiian shield lavas from Kilauea and Mauna Loa. Here we posit that as the Pacific lithosphere beneath the active shield volcanoes moves away from the center of the Hawaiian plume, increased incorporation of an intrinsic depleted component with relatively low 206 Pb/ 204 Pb produces the source of the basanites of the Hana Volcanics. Haleakala Crater basanites have average ( 230 Th/ 238 U) of 1.23 ( n = 13), average age-corrected ( 226 Ra/ 230 Th) of 1.25 ( n = 13), and average ( 231 Pa/ 235 U) of 1.67 ( n = 4), significantly higher than Kilauea and Mauna Loa tholeiites. U-series modeling shows that solid mantle upwelling velocity for Haleakala Crater basanites ranges from ∼0.7 to 1.0 cm/yr, compared to ∼10 to 20 cm/yr for tholeiites and ∼1 to 2 cm/yr for alkali basalts. These modeling results indicate that solid mantle upwelling rates and porosity of the melting zone are lower for Hana Volcanics basanites than for shield-stage tholeiites from Kilauea and Mauna Loa and alkali basalts from Hualalai. The melting rate, which is directly proportional to both the solid mantle upwelling rate and the degree of melting, is therefore greatest in the center of the Hawaiian plume and lower on its periphery. Our results indicate that solid mantle upwelling velocity is at least 10 times higher at the center of the plume than at its periphery under Haleakala.

Hawaii↗

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↗

USGS2017

Anoxic nitrate reduction coupled with iron oxidation and attenuation of dissolved arsenic and phosphate in a sand and gravel aquifer

Nitrate has become an increasingly abundant potential electron acceptor for Fe(II) oxidation in groundwater, but this redox couple has not been well characterized within aquifer settings. To investigate this reaction and some of its implications for redox-sensitive groundwater contaminants, we conducted an in situ field study in a wastewater-contaminated aquifer on Cape Cod. Long-term (15 year) geochemical monitoring within the contaminant plume indicated interacting zones with variable nitrate-, Fe(II)-, phosphate-, As(V)-, and As(III)-containing groundwater. Nitrate and phosphate were derived predominantly from wastewater disposal, whereas Fe(II), As(III), and As(V) were mobilized from the aquifer sediments. Multiple natural gradient, anoxic tracer tests were conducted in which nitrate and bromide were injected into nitrate-free, Fe(II)-containing groundwater. Prior to injection, aqueous Fe(II) concentrations were approximately 175 μM, but sorbed Fe(II) accounted for greater than 90% of the total reactive Fe(II) in the aquifer. Nitrate reduction was stimulated within 1 m of transport for 100 μM and 1000 μM nitrate additions, initially producing stoichiometric quantities of nitrous oxide (>300 μM N). In subsequent injections at the same site, nitrate was reduced even more rapidly and produced less nitrous oxide, especially over longer transport distances. Fe(II) and nitrate concentrations decreased together and were accompanied by Fe(III) oxyhydroxide precipitation and decreases in dissolved phosphate, As(III), and As(V) concentrations. Nitrate N and O isotope fractionation effects during nitrate reduction were approximately equal (ε 15 N/ε 18 O = 1.11) and were similar to those reported for laboratory studies of biological nitrate reduction, including denitrification, but unlike some reported effects on nitrate by denitrification in aquifers. All constituents affected by the in situ tracer experiments returned to pre-injection levels after several weeks. Additionally, Fe(II)-oxidizing, nitrate-reducing microbial enrichment cultures were obtained from aquifer sediments. Growth experiments with the cultures sequentially produced nitrite and nitrous oxide from nitrate while simultaneously oxidizing Fe(II). Field and culture results suggest that nitrogen oxide reduction and Fe(II) oxidation in the aquifer are a complex interaction of coupled biotic and abiotic reactions. Overall, the results of this study demonstrate that anoxic nitrate-dependent iron oxidation can occur in groundwater; that it could control iron speciation; and that the process can impact the mobility of other chemical species (e.g., phosphate and arsenic) not directly involved in the oxidation–reduction reaction.

Massachusetts↗

USGS2017

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↗

USGS2017

Rare earth element behavior during groundwater – seawater mixing along the Kona Coast of Hawaii

Groundwater and seawater samples were collected from nearshore wells and offshore along the Kona Coast of the Big Island of Hawaii to investigate rare earth element (REE) behavior in local subterranean estuaries. Previous investigations showed that submarine groundwater discharge (SGD) is the predominant flux of terrestrial waters to the coastal ocean along the arid Kona Coast of Hawaii. Groundwater and seawater samples were filtered through 0.45 μm and 0.02 μm pore-size filters to evaluate the importance of colloidal and soluble (i.e., truly dissolved ionic species and/or low molecular weight [LMW] colloids) fractions of the REEs in the local subterranean estuaries. Mixing experiments using groundwater collected immediately down gradient from a wastewater treatment facility (WWTF) proximal to the Kaloko-Hanokohau National Historic Park, and more “pristine” groundwater from a well constructed in a lava tube at Kiholo Bay, were mixed with local seawater to study the effect of solution composition (i.e., pH, salinity) on the concentrations and fractionation behavior of the REEs as groundwater mixes with seawater in Kona Coast subterranean estuaries. The mixed waters were also filtered through 0.45 or 0.02 μm filters to ascertain the behavior of colloidal and soluble fractions of the REEs across the salinity gradient in each mixing experiment. Concentrations of the REEs were statistically identical (two-tailed Student t -test, 95% confidence) between the sequentially filtered sample aliquots, indicating that the REEs occur as dissolved ionic species and/or LMW colloids in Kona Coast groundwaters. The mixing experiments revealed that the REEs are released to solution from suspended particles or colloids when Kona Coast groundwater waters mix with local seawater. The order of release that accompanies increasing pH and salinity follows light REE (LREE) > middle REE (MREE) > heavy REE (HREE). Release of REEs in the mixing experiments is driven by decreases in the free metal ion activity in solution and the concomitant increase in the amount of each REE that occurs in solution as dicarbonato complexes [i.e., Ln(CO 3 ) 2 - ] as pH increases across the salinity gradient. Input-normalized REE patterns of Kona Coast groundwater and coastal seawater are nearly identical and relatively flat compared to North Pacific seawater, indicating that SGD is the chief source of these trace elements to the ocean along the Kona Coast. Additionally, REE concentrations of the coastal seawater are between 10 and 50 times higher than previously reported open-ocean seawater values from the North Pacific, further demonstrating the importance of SGD fluxes of REEs to these coastal waters. Taken together, these observations indicate that large-scale removal of REEs, which characterizes the behavior of REEs in the low salinity reaches of many surface estuaries, is not a feature of the subterranean estuary along the Kona Coast. A large positive gadolinium (Gd) anomaly characterizes groundwater from the vicinity of the WWTF. The positive Gd anomaly can be traced to the coastal ocean, providing further evidence of the impact of SGD on the coastal waters. Estimates of the SGD fluxes of the REEs to the coastal ocean along the Kona Coast (i.e., 1.3 – 2.6 mmol Nd day -1 ) are similar to recent estimates of SGD fluxes of REEs along Florida’s east coast and to Rhode Island Sound, all of which points to the importance of SGD as significant flux of REEs to the coastal ocean.

Hawaii↗

USGS2017

Long-term flow-through column experiments and their relevance to natural granitoid weathering rates

Four pairs of fresh and partly-weathered granitoids, obtained from well-characterized watersheds—Merced River, CA, USA; Panola, GA, USA; Loch Vale, CO, USA, and Rio Icacos, Puerto Rico—were reacted in columns under ambient laboratory conditions for 13.8 yrs, the longest running experimental weathering study to date. Low total column mass losses (<1 wt. %), correlated with the absence of pitting or surface roughening of primary silicate grains. BET surface area (S BET ) increased, primarily due to Fe-oxyhydroxide precipitation. Surface areas returned to within factors of 2 to 3 of their original values after dithionite extraction. Miscible displacement experiments indicated homogeneous plug flow with negligible immobile water, commonly cited for column experiments. Fresh granitoid effluent solute concentrations initially declined rapidly, followed by much slower decreases over the next decade. Weathered granitoid effluent concentrations increased modestly over the same time period, indicating losses of natural Fe-oxide and/or clay coatings and the increased exposure of primary mineral surfaces. Corresponding (fresh and weathered) elemental effluent concentrations trended toward convergence during the last decade of reaction. NETPATH/PHREEQC code simulations indicated non-stoichiometric dissolution involving Ca release from disseminated calcite and excess K release from interlayer biotite. Effluent 87 Sr/ 85 Sr ratios reflected a progressive weathering sequence beginning and ending with 87 Sr/ 85 Sr values of plagioclase with an additional calcite input and a radiogenic biotite excursion proportional to the granitoid ages. Effluents became thermodynamically saturated with goethite and gibbsite, slightly under-saturated with kaolinite and strongly under-saturated with plagioclase, consistent with kinetically-limited weathering in which solutes such as Na varied with column flow rates. Effluent Na concentrations showed no clear trend with time during the last decade of reaction (fresh granitoids) or increased slowly with time (weathered granitoids). Analysis of cumulative Na release indicated that plagioclase dissolution achieved steady state in 3 of the 4 fresh granitoids during the last decade of reaction. Surface-area normalized plagioclase dissolution rates exhibited a narrow range (0.95 to 1.26 10 -13 moles m -2 s -1 ), in spite of significant stoichiometric differences (An 0.21 to An 0.50 ). Rates were an order of magnitude slower than previously reported in shorter duration experiments but generally 2 to 3 orders of magnitude faster than corresponding natural analogs. CrunchFlow simulations indicated that more than a hundredfold decrease in column flow rates would be required to produce near-saturation reaction affinities that would start to slow plagioclase weathering to real-world levels. Extending simulations to approximate long term weathering in naturally weathered profiles required additional decreases in the intrinsic plagioclase dissolution and kaolinite precipitation rates and relatively large decreases in the fluid flow rate, implying that exposure to reactive mineral surfaces is significantly limited in the natural environment compared to column experiments.

Geochimica et Cosmochimica Acta↗

USGS2017

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↗

USGS2017

Oxygen and U-Th isotopes and the timescales of hydrothermal exchange and melting in granitoid wall rocks at Mount Mazama, Crater Lake, Oregon

We report new whole rock U-Th and in-situ oxygen isotope compositions for partially melted (0–50 vol% melt), low-δ 18 O Pleistocene granitoid blocks ejected during the ∼7.7 ka caldera-forming eruption of Mt. Mazama (Crater Lake, Oregon). The blocks are interpreted to represent wall rocks of the climactic magma chamber that, prior to eruption, experienced variable amounts of exchange with meteoric hydrothermal fluids and subsequent partial melting. U-Th and oxygen isotope results allow us to examine the timescales of hydrothermal circulation and partial melting, and provide an “outside in” perspective on the buildup to the climactic eruption of Mt. Mazama. Oxygen isotope compositions measured in the cores and rims of individual quartz ( n = 126) and plagioclase ( n = 91) crystals, and for transects across ten quartz crystals, document zonation in quartz (Δ 18 O Core-Rim ≤ 0.1–5.5‰), but show homogeneity in plagioclase (Δ 18 O Core-Rim ≤ ±0.8‰). We propose that oxygen isotope zonation in quartz records hydrothermal exchange followed by high-temperature exchange in response to partial melting caused by injection of basaltic to andesitic recharge magma into the deeper portions of the chamber. Results of modeling of oxygen diffusion in quartz indicates that hydrothermal exchange in quartz occurred over a period of ∼1000–63,000 years. Models also suggest that the onset of melting of the granitoids occurred a minimum of ∼10–200 years prior to the Mazama climactic eruption, an inference which is broadly consistent with results for magnetite homogenization and for Zr diffusion in melt previously reported by others. Uranium-thorium isotope compositions of most granitoid blocks are in 238 U excess, and are in agreement with a 238 U enriched array previously measured for volcanic rocks at Mt. Mazama. Uranium excess in the granitoids is likely due to enrichment via hydrothermal circulation, given their low δ 18 O values. The sample with the highest U excess (≥5.8%) also has the most 18 O isotope depletion (average δ 18 O plag = −4.0‰). The granitoids are a probable assimilant and source of U excess in volcanic rocks from Mt. Mazama. Two granitoids have Th excess and low δ 18 O values, interpreted to record leaching of U during hydrothermal alteration. A U-Th isochron based on the U excess array of the granitoids and volcanic rocks indicates that hydrothermal circulation initiated ∼40–75 kyrs before the climactic eruption, potentially marking the initiation of a persistent upper-crustal magma chamber. The U-Th ages are consistent with the maximum timescales inferred for hydrothermal alteration based on oxygen isotope zoning in quartz.

Oregon↗

USGS2017

Mantle and crustal gases of the Colorado Plateau: Geochemistry, sources, and migration pathways

The Colorado Plateau hosts several large accumulations of naturally occurring, non-hydrocarbon gases, including CO 2 , N 2 , and the noble gases, making it a good field location to study the fluxes of these gases within the crust and to the atmosphere. In this study, we present a compilation of 1252 published gas-composition measurements. The data reveal at least three natural gas associations in the field area, which are dominated by hydrocarbons, CO 2, and N 2 + He + Ar, respectively. Most gas accumulations of the region exhibit compositions that are intermediate between the three end members. The first non-hydrocarbon gas association is characterized by very high-purity CO 2 , in excess of 75 mol% (hereafter, %). Many of these high-purity CO 2 fields have recently been well described and interpreted as magmatic in origin. The second non-hydrocarbon gas association is less well described on the Colorado Plateau. It exhibits He concentrations on the order of 1–10%, and centered log ratio biplots show that He occurs proportionally to both N 2 and Ar. Overall ratios of N 2 to He to Ar are ≈100:10:1 and correlation in concentrations of these gases suggests that they have been sourced from the same reservoir and/or by a common process. To complement the analysis of the gas-composition data, stable isotope and noble-gas isotope measurements are compiled or newly reported from 11 representative fields (previously published data from 4 fields and new data from 7 fields). Gas sampled from the Harley Dome gas field in Utah contains nearly pure N 2 + He + Ar. The various compositional and stable and noble gas isotopic data for this gas indicate that noble gas molecule/isotope ratios are near crustal radiogenic production values and also suggest a crustal N 2 source. Across the field area, most of the high-purity N 2 + He + Ar gas accumulations are associated with the mapped surface trace of structures or sutures in the Precambrian basement and are often accumulated in lower parts of the overlying Phanerozoic sedimentary cover. The high-purity gas association mostly occurs in areas interior to the plateau that are characterized by a narrow range of elevated, moderate heat flow values (53–74 mW/m 2 ) in the ancient (1.8–1.6 Ga) basement terranes of the region. Collectively, the geochemical and geological data suggest that (1) the N 2 + He + Ar gas association is sourced from a crustal reservoir, (2) the gas association migrates preferentially along structures in the Precambrian basement, and (3) the sourcing process relates to heating of the crust. Prospecting for noble-gas accumulations may target areas with elevated Cenozoic heat flow, ancient crust, and deep crustal structures that focus gas migration. High-purity CO 2 gas may also migrate through regional basement structures, however, there is not always a clear spatial association. Rather, CO 2 accumulations are more clearly associated with zones of high heat flow (>63 mW/m 2 ) that sit above hot upper mantle and are proximal to Cenozoic volcanic rocksnear the plateau margins. These observations are consistent with previous interpretations of a magmatic gas source, which were based on geochemical measurements.

Colorado Plateau↗