Rates of silicate dissolution in deep-sea sediment: In situ measurement using U/U of pore fluids

Bulk dissolution rates for sediment from ODP Site 984A in the North Atlantic are determined using the ²³⁴U/²³⁸U activity ratios of pore water, bulk sediment, and leachates. Site 984A is one of only several sites where closely spaced pore water samples were obtained from the upper 60 meters of the core; the sedimentation rate is high (11-15 cm/ka), hence the sediments in the upper 60 meters are less than 500 ka old. The sediment is clayey silt and composed mostly of detritus derived from Iceland with a significant component of biogenic carbonate (up to 30%).The pore water ²³⁴U/²³⁸U activity ratios are higher than seawater values, in the range of 1.2 to 1.6, while the bulk sediment ²³⁴U/²³⁸U activity ratios are close to 1.0. The ²³⁴U/²³⁸U of the pore water reflects a balance between the mineral dissolution rate and the supply rate of excess ²³⁴U to the pore fluid by α-recoil injection of ²³⁴Th. The fraction of ²³⁸U decays that result in α-recoil injection of ²³⁴U to pore fluid is estimated to be 0.10 to 0.20 based on the ²³⁴U/²³⁸U of insoluble residue fractions. The calculated bulk dissolution rates, in units of g/g/yr are in the range of 4 × 10⁻⁷ to 2 × 10⁻⁶ yr⁻¹. There is significant down-hole variability in pore water ²³⁴U/²³⁸U activity ratios (and hence dissolution rates) on a scale of ca. 10 m. The inferred bulk dissolution rate constants are 100 to 10⁴ times slower than laboratory-determined rates, 100 times faster than rates inferred for older sediments based on Sr isotopes, and similar to weathering rates determined for terrestrial soils of similar age. The results of this study suggest that U isotopes can be used to measure in situ dissolution rates in fine-grained clastic materials.The rate estimates for sediments from ODP Site 984 confirm the strong dependence of reactivity on the age of the solid material: the bulk dissolution rate (R_d) of soils and deep-sea sediments can be approximately described by the expression R_d ≈ 0.1 Age⁻¹ for ages spanning 1000 to 5 × 10⁸ yr. The age of the material, which encompasses the grain size, surface area, and other chemical factors that contribute to the rate of dissolution, appears to be a much stronger determinant of dissolution rate than any single physical or chemical property of the system.     

Calcium isotope evidence for suppression of carbonate dissolution in carbonate-bearing organic-rich sediments

Pore fluid calcium isotope, calcium concentration and strontium concentration data are used to measure the rates of diagenetic dissolution and precipitation of calcite in deep-sea sediments containing abundant clay and organic material. This type of study of deep-sea sediment diagenesis provides unique information about the ultra-slow chemical reactions that occur in natural marine sediments that affect global geochemical cycles and the preservation of paleo-environmental information in carbonate fossils. For this study, calcium isotope ratios (δ^44/40Ca) of pore fluid calcium from Ocean Drilling Program (ODP) Sites 984 (North Atlantic) and 1082 (off the coast of West Africa) were measured to augment available pore fluid measurements of calcium and strontium concentration. Both study sites have high sedimentation rates and support quantitative sulfate reduction, methanogenesis and anaerobic methane oxidation. The pattern of change of δ^44/40Ca of pore fluid calcium versus depth at Sites 984 and 1082 differs markedly from that of previously studied deep-sea Sites like 590B and 807, which are composed of nearly pure carbonate sediment. In the 984 and 1082 pore fluids, δ^44/40Ca remains elevated near seawater values deep in the sediments, rather than shifting rapidly toward the δ^44/40Ca of carbonate solids. This observation indicates that the rate of calcite dissolution is far lower than at previously studied carbonate-rich sites. The data are fit using a numerical model, as well as more approximate analytical models, to estimate the rates of carbonate dissolution and precipitation and the relationship of these rates to the abundance of clay and organic material. Our models give mutually consistent results and indicate that calcite dissolution rates at Sites 984 and 1082 are roughly two orders of magnitude lower than at previously studied carbonate-rich sites, and the rate correlates with the abundance of clay. Our calculated rates are conservative for these sites (the actual rates could be significantly slower) because other processes that impact the calcium isotope composition of sedimentary pore fluid have not been included. The results provide direct geochemical evidence for the anecdotal observation that the best-preserved carbonate fossils are often found in clay or organic-rich sedimentary horizons. The results also suggest that the presence of clay minerals has a strong passivating effect on the surfaces of biogenic carbonate minerals, slowing dissolution dramatically even in relation to the already-slow rates typical of carbonate-rich sediments.     

The role of reaction affinity and secondary minerals in regulating chemical weathering rates at the Santa Cruz Soil Chronosequence, California

In order to explore the reasons for the apparent discrepancy between laboratory and field weathering rates and to determine the extent to which weathering rates are controlled by the approach to thermodynamic equilibrium, secondary mineral precipitation, and flow rates, a multicomponent reactive transport model (CrunchFlow) was used to interpret soil profile development and mineral precipitation and dissolution rates at the 226 ka Marine Terrace Chronosequence near Santa Cruz, CA. Aqueous compositions, fluid chemistry, transport, and mineral abundances are well characterized [White A. F., Schulz M. S., Vivit D. V., Blum A., Stonestrom D. A. and Anderson S. P. (2008) Chemical weathering of a Marine Terrace Chronosequence, Santa Cruz, California. I: interpreting the long-term controls on chemical weathering based on spatial and temporal element and mineral distributions. Geochim. Cosmochim. Acta72 (1), 36-68] and were used to constrain the reaction rates for the weathering and precipitating minerals in the reactive transport modeling. When primary mineral weathering rates are calculated with either of two experimentally determined rate constants, the nonlinear, parallel rate law formulation of Hellmann and Tisserand [Hellmann R. and Tisserand D. (2006) Dissolution kinetics as a function of the Gibbs free energy of reaction: An experimental study based on albite feldspar. Geochim. Cosmochim. Acta70 (2), 364-383] or the aluminum inhibition model proposed by Oelkers et al. [Oelkers E. H., Schott J. and Devidal J. L. (1994) The effect of aluminum, pH, and chemical affinity on the rates of aluminosilicate dissolution reactions. Geochim. Cosmochim. Acta58 (9), 2011-2024], modeling results are consistent with field-scale observations when independently constrained clay precipitation rates are accounted for. Experimental and field rates, therefore, can be reconciled at the Santa Cruz site.Additionally, observed maximum clay abundances in the argillic horizons occur at the depth and time where the reaction fronts of the primary minerals overlap. The modeling indicates that the argillic horizon at Santa Cruz can be explained almost entirely by weathering of primary minerals and in situ clay precipitation accompanied by undersaturation of kaolinite at the top of the profile. The rate constant for kaolinite precipitation was also determined based on model simulations of mineral abundances and dissolved Al, SiO₂(aq) and pH in pore waters. Changes in the rate of kaolinite precipitation or the flow rate do not affect the gradient of the primary mineral weathering profiles, but instead control the rate of propagation of the primary mineral weathering fronts and thus total mass removed from the weathering profile. Our analysis suggests that secondary clay precipitation is as important as aqueous transport in governing the amount of dissolution that occurs within a profile because clay minerals exert a strong control over the reaction affinity of the dissolving primary minerals. The modeling also indicates that the weathering advance rate and the total mass of mineral dissolved is controlled by the thermodynamic saturation of the primary dissolving phases plagioclase and K-feldspar, as is evident from the difference in propagation rates of the reaction fronts for the two minerals despite their very similar kinetic rate laws.     

U-Sr isotopic speedometer: Fluid flow and chemical weathering rates in aquifers

Both chemical weathering rates and fluid flow are difficult to measure in natural systems. However, these parameters are critical for understanding the hydrochemical evolution of aquifers, predicting the fate and transport of contaminants, and for water resources/water quality considerations. ⁸⁷Sr/⁸⁶Sr and (²³⁴U/²³⁸U) activity ratios are sensitive indicators of water-rock interaction, and thus provide a means of quantifying both flow and reactivity. The ⁸⁷Sr/⁸⁶Sr values in ground waters are controlled by the ratio of the dissolution rate to the flow rate. Similarly, the (²³⁴U/²³⁸U) ratio of natural ground waters is a balance between the flow rate and the dissolution of solids, and α-recoil loss of ²³⁴U from the solids. By coupling these two isotope systems it is possible to constrain both the long-term (ca. 100’s to 1000’s of years) flow rate and bulk dissolution rate along the flow path. Previous estimates of the ratio of the dissolution rate to the infiltration flux from Sr isotopes (⁸⁷Sr/⁸⁶Sr) are combined with a model for (²³⁴U/²³⁸U) to constrain the infiltration flux and dissolution rate for a 70-m deep vadose zone core from Hanford, Washington. The coupled model for both (²³⁴U/²³⁸U) ratios and the ⁸⁷Sr/⁸⁶Sr data suggests an infiltration flux of 5 ± 2 mm/yr, and bulk silicate dissolution rates between 10^−15.7 and 10^−16.5 mol/m²/s. The process of α-recoil enrichment, while primarily responsible for the observed variation in (²³⁴U/²³⁸U) of natural systems, is difficult to quantify. However, the rate of this process in natural systems affects the interpretation of most U-series data. Models for quantifying the α-recoil loss fraction based on geometric predictions, surface area constraints, and chemical methods are also presented. The agreement between the chemical and theoretical methods, such as direct measurement of (²³⁴U/²³⁸U) of the small grain size fraction and geometric calculations for that size fraction, is quite good.     

Fixation of radionuclides in the 238U decay series in the vicinity of mineralized zones: 1. The Austatom Uranium Prospect,Northern Territory,Australia

The minimum age of a zone of secondary uranium mineralization, located at the Austatom Prospect in the Alligator Rivers region of Australia, is estimated to be 3.6 × 10⁵y. This is derived from a geochronological model based on retarded leaching of ²³⁴U with respect to ²³⁸U and on ratios within the ore of these members of the ²³⁸U decay series. Although kaolinite is a dominant mineral in the weathered schist-host-rocks, retarded dissolution of ²³⁴U occurs only in the presence of the clay minerals illite and montmorillonite. In their absence the reverse occurs. A model is proposed to explain the results. Ratios of ²³⁰Th to ²³⁸U indicate that the mineralization has probably remained stationary within the weathered schist for at least 1 to 2 × 10⁵y. Future use of clay minerals as buffers in radioactive waste repositories is supported by the excellent long-term retention obtained for oxidized uranium, probably due in part to isomorphic substitution into the clay crystal lattice.     

Transport of U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of ²³⁸U, ²³⁴U and ²³⁰Th was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim was to assess the in situ, long-term migration behaviour of ²³⁴U in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both ²³⁸U and ²³⁰Th immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of ²³⁴U from the rock matrix to the pore water by the decay of ²³⁸U will be established. Diffusive redistribution separates ²³⁴U from its immobile parent ²³⁸U resulting in bulk rock ²³⁴U/²³⁸U activity disequilibria. These may provide a means of estimating the mobility of ²³⁴U in the rock if the diffusion rate of ²³⁴U is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by β-scanning, which shows that the (dominant) clay is richer in both elements.     

A coupled non-isothermal reactive transport model for long-term geochemical evolution of a HLW repository in clay

A numerical model of coupled saturated/unsaturated water flow, heat transfer and multi-component reactive solute transport is presented to evaluate the long-term geochemical evolution in bentonite, concrete and clay formation for a potential geological radioactive waste repository. Changes in formation porosity caused by mineral dissolution/precipitation reactions are taken into account. Simulations were carried out with a general-purpose multicomponent reactive transport code, CORE^2D V4. Numerical results show that pH in the bentonite porewater can vary from neutral to up to 13 over a time scale of 1 Ma although dissolution of silica minerals and precipitation of secondary calcium silicate hydrate (CSH) minerals in bentonite buffer the effect of the hyperalkaline plume. Mineral precipitation reduces the volume of pore space in bentonite close to the bentonite–concrete interface due to the precipitation of CSH minerals. Model results indicate that bentonite porosity decreases less than 25%. The hyperalkaline plume from the concrete only extends to a distance of 0.7 m in the clay formation over the time range of 1 Ma.     

Preferential solution of234U from recoil tracks and234U/238U radioactive disequilibrium in natural waters

A mathematical model to calculate the²³⁴U/²³⁸U activity ratio (AR) in an aqueous phase in contact with rock/soil is presented. The model relies on the supply of²³⁸U by dissolution and that of²³⁴U by dissolution and preferential release from radiation damaged regions (recoil tracks). The model predicts that values of²³⁴U/²³⁸U AR>1 in the aqueous phase can be obtained only from weathering “virgin” surfaces. Thus, to account for the observed steady-state supply of²³⁴U excess to the oceans by the preferential leaching model, ‘virgin’ rock/soil surfaces would have to be continually exposed and weathered. The²³⁸U concentration and²³⁴U/²³⁸U AR in continental waters allow us to estimate the exposure rates of “virgin” rock/soil surfaces.     

Modeling non-steady state radioisotope transport in the vadose zone - A case study using uranium isotopes at Peña Blanca, Mexico

Current models using U- and Th-series disequilibria to study radioisotope transport in groundwater systems mostly consider a steady-state situation. These models have limited applicability to the vadose zone (UZ) where the concentration and migratory behavior of radioisotopes in fluid are often transitory. We present here, as a first attempt of its kind, a model simulating the non-steady state, intermittent fluid transport in vadose layers. It provides quantitative constraints on in-situ migration of dissolved and colloidal radioisotopes in terms of retardation factor and rock-water interaction (or water transit) time. For uranium, the simulation predicts that intermittent flushing in the UZ leads to a linear relationship between reciprocal U concentration and ²³⁴U/²³⁸U ratio in percolating waters, with the intercept and slope bearing information on the rates of dissolution and α-recoil of U isotopes, respectively. The general validity of the model appears to be borne out by the measurement of uranium isotopes in UZ waters collected at various times over a period during 1995-2006 from a site in the Peña Blanca mining district, Mexico, where the Nopal I uranium deposit is located. Enhanced ²³⁴U/²³⁸U ratios in vadose-zone waters resulting from lengthened non-flushing time as prescribed by the model provide an interpretative basis for using ²³⁴U/²³⁸U in cave calcites to reconstruct the regional changes in hydrology and climate. We also provide a theoretical account of the model’s potential applications using radium isotopes.     

Weathering rates from top to bottom in a carbonate environment

This study investigates U-series, Sr isotopes, major and trace elements in a chalk aquifer system located in Eastern France. Soil and rock samples were collected along depth profiles down to 45 m in four localities as an attempt to investigate the weathering processes in the soil, the unsaturated zone and the saturated zone of the aquifer. Interstitial water was extracted from soils and rocks by a centrifugation technique. U-series offer a powerful tool to calculate weathering rates because the relative mobility of the U- and Th-isotopes can be precisely measured and it does not require the determination of a reference state as in other approaches. As expected, the data show very large mobile element depletion in the soil with large ²³⁰Th excess relative to ²³⁸U, while the rocks show more limited but not insignificant mobile element depletion. The U-series data have been used to constrain weathering rates based on a 1-D reactive transport model. Weathering rates in the near surface are about 10–100 times faster than at depth. However, when integrated over the depth of the cores, including the unsaturated and the saturated zones, this underground weathering represents more than 30% of the total weathering flux, assuming congruent dissolution of carbonates. The (²³⁴U/²³⁸U) ratios in interstitial water are consistent with solid samples showing ²³⁴U depletion near the surface and an excess ²³⁴U at depth. A leaching experiment performed on chalk shows that the excess ²³⁴U in natural waters percolating through carbonate rocks results both from preferential ²³⁴U leaching and direct recoil in the interstitial water. A new approach was used to derive the recoil ejection factor based on BET measurements and the fractal dimension of chalk surface. Consideration of preferential leaching and recoil allows a more accurate modeling of weathering rates.     

Experimental evidence for U-U fractionation during granite weathering with implications for U/U in natural waters

The daughter to parent (²³⁴U/²³⁸U) activity ratio in natural waters is often out of secular radioactive equilibrium. The major reason for this disequilibrium is related to the energetic α-decay of ²³⁸U and differential release of ²³⁴U relative to ²³⁸U. This disequilibrium originates from (1) preferential release of more loosely bound ²³⁴U from damaged mineral lattice sites or; (2) direct recoil of ²³⁴Th into surrounding media from near mineral surface boundaries, however, it is unclear which of the two mechanisms is most important in nature. To better quantify the effects of preferential release of ²³⁴U, two continuous laboratory granite leaching experiments conducted over 1100 h were performed. The leachates were characterized by declining U concentrations with time and (²³⁴U/²³⁸U) initially greater than unity (up to 1.15), which changed to below unity during leaching (∼0.95). The early elevated (²³⁴U/²³⁸U) suggests that additional ²³⁴U is released into solution by preferential release of ²³⁴U from mineral phases. However, the excess ²³⁴U constitutes a finite pool of easy leachable ²³⁴U and the (²³⁴U/²³⁸U) values become lower than unity when this pool is used up. A model based on first-order kinetics, dissolution rates and preferential release of ²³⁴U from damaged lattice sites was developed and is able to quantitatively predict the observed pattern of (²³⁴U/²³⁸U) values and U concentrations for the two granite leaching experiments. Extending the modeling to longer time scales more comparable to natural systems shows that the production of waters with high (²³⁴U/²³⁸U) ratios can be achieved in two distinct regimes (1) slow weathering where the rate of directly recoiled ²³⁴U near mineral surfaces into waters is high; (2) fast weathering where the role of incipient chemical weathering and preferential release of loosely bound ²³⁴U are important. The model is able to explain apparent opposite correlations between physical erosion rates and (²³⁴U/²³⁸U) in waters and it provides a new framework that will be useful for examining weathering regimes, their timescales and their coupling with physical erosion.     

Kinetic Modeling of Diagenesis of Eogene Lacustrine Sandstone Reservoirs in the Jianghan Basin,Southeastern China

In the Tuoshi oilfield,located in the Cenozoic Jianghan Basin of southeastern China ,there have been found hydrocarbon reservoirs hosted in lacustrine sandstones of the Eogene Xingouzui Formation.The main diagenetic features identified in these sandstones include the dissolution of detrital K-feldspar and albite grains,the precipitation of quartz as overgrowths and /or cements ,and the precipitation and /or transformation of clay minerals.These diagenetic features were interpreted to have occurred in early,intermediate and late stages,based on the burial depth.The kinetics of fluid-mineral reactions and the concentrations of aqueous species au each stage of diagenesis were simulated numerically for these lacustrine sandstones,using a quasi-sta-tionary state approximation that incorporates simultaneous chemical reactions in a time-space continuum.During the early diagenetic stage,pore fluid was weakly acidic,which resulted in dissolution of K-feldspar and albite and,therefore,led to the release of K^ ,Na^ ,Al^3 and SiO2(aq) into the diagenetic fluid.The increased K^ ,Na^ ,Al^3 and SiO2(aq) concentrations in the diagenetic fluid caused the precipitation of quartz,kaolinite and illite.At the beginning of the intermediate diagenetic stage the concentration of H^ was built up due to the decomposition of organic matter,which was responsible for further dissolution of K-feldspar and albite and pre-cipitation of quartz,kaolinite,and illite.During the late diagenetic stage,the pore fluid was weakly alkaline,K-feldspar became stable and was precipitated with quartz and clay minerals.When the burial depth was greater than 3000 m,the pore fluids became supersaturated with respect to allbite,but undersaturated with respect to quartz,resulting in the precipitation of albite and the dissolution of quartz.The diagenetic reactions forecasted in the numerical modeling closely matched the diagenet-ic features identified by petrographic examination, and therefore,can help us to gain a better understanding of the diagenetic processes and associated porosity evolution in sandstone reservoirs.     

Processes controlling 234U and 238U isotope fractionation and helium in the groundwater of the St. Lawrence Lowlands,Quebec: The potential role of natural rock fracturing

The goal of this study is to explain the origin of ²³⁴U–²³⁸U fractionation in groundwater from sedimentary aquifers of the St. Lawrence Lowlands (Quebec, Canada), and its relationship with ³He/⁴He ratios, to gain insight regarding the evolution of groundwater in the region. (²³⁴U/²³⁸U) activity ratios, or (²³⁴U/²³⁸U)_act, were measured in 23 groundwater samples from shallow Quaternary unconsolidated sediments and from the deeper fractured regional aquifer of the Becancour River watershed. The lowest (²³⁴U/²³⁸U)_act, 1.14 ± 0.01, was measured in Ca–HCO₃-type freshwater from the Quaternary Shallower Aquifer, where bulk dissolution of the carbonate allows U to migrate into water with little ²³⁴U–²³⁸U isotopic fractionation. The (²³⁴U/²³⁸U)_act increases to 6.07 ± 0.14 in Na–HCO₃–Cl-type groundwater. Preferential migration of ²³⁴U into water by α-recoil is the underlying process responsible for this isotopic fractionation. An inverse relationship between (²³⁴U/²³⁸U)_act and ³He/⁴He ratios has been observed. This relationship reflects the mixing of newly recharged water, with (²³⁴U/²³⁸U)_act close to the secular equilibrium and containing atmospheric/tritiogenic helium, and mildly-mineralized older water (¹⁴C ages of 6.6 kyrs), with (²³⁴U/²³⁸U)_act of ≥6.07 and large amounts of radiogenic ⁴He, in excess of the steady-state amount produced in situ. The simultaneous fractionation of (²³⁴U/²³⁸U)_act and the addition of excess ⁴He could be locally controlled by stress-induced rock fracturing. This process increases the surface area of the aquifer matrix exposed to pore water, from which produced ⁴He and ²³⁴U can be released by α-recoil and diffusion. This process would also facilitate the release of radiogenic helium at rates greater than those supported by steady-state U–Th production in the rock. Consequently, sources internal to the aquifers could cause the radiogenic ⁴He excesses measured in groundwater.     

Uranium isotopes as a tracer of sources of dissolved solutes in the Han River,South Korea

The uranium (U) content and ²³⁴U/²³⁸U activity ratio were determined for water samples collected from Korea's Han River in spring, summer, and winter 2006 to provide data that might constrain the origin of U isotope fractionation in river water and the link between U isotope systematics in river waters and the lithological nature of the corresponding bedrock. The large difference in the major dissolved loads between the two major branches of the Han River, the North Han River (NHR) and South Han River (SHR), is reflected in the contrasting U content and ²³⁴U/²³⁸U activity ratio between the tributaries: low U content (0.08–0.75 nM; average, 0.34 nM) and small ²³⁴U/²³⁸U activity ratio (1.03–1.22; average, 1.09) in the NHR; and high U content (0.65–1.98 nM; average, 1.44 nM) and large ²³⁴U/²³⁸U activity ratio (1.05–1.45; average, 1.24) in the SHR. The large spatial differences in U content and ²³⁴U/²³⁸U activity ratio are closely related to both lithological differences between the two tributaries and groundwater input. The low U content and small ²³⁴U/²³⁸U activity ratio in the NHR arise mainly from a combination of surface and meteoric weathering of the dominant silicate rocks in this branch and congruent dissolution of already weathered (secular equilibrium) materials. In contrast, the high U content and large ²³⁴U/²³⁸U activity ratio in the SHR are ascribed to the dissolution of carbonates and black shales along with significant inputs of deep groundwater.     

U/Th systematics and ages of authigenic carbonates from Hydrate Ridge, Cascadia Margin: recorders of fluid flow variations

Uranium (U) concentrations and activity ratios (δ²³⁴U) of authigenic carbonates are sensitive recorders of different fluid compositions at submarine seeps of hydrocarbon-rich fluids (“cold seeps”) at Hydrate Ridge, off the coast of Oregon, USA. The low U concentrations (mean: 1.3 ± 0.4 μg/g) and high δ²³⁴U values (165-317‰) of gas hydrate carbonates reflect the influence of sedimentary pore water indicating that these carbonates were formed under reducing conditions below or at the seafloor. Their ²³⁰Th/²³⁴U ages span a time interval from 0.8 to 6.4 ka and cluster around 1.2 and 4.7 ka. In contrast, chemoherm carbonates precipitate from marine bottom water marked by relatively high U concentrations (mean: 5.2 ± 0.8 μg/g) and a mean δ²³⁴U ratio of 166 ± 3‰. Their U isotopes reflect the δ²³⁴U ratios of the bottom water being enriched in ²³⁴U relative to normal seawater. Simple mass balance calculations based on U concentrations and their corresponding δ²³⁴U ratios reveal a contribution of about 11% of sedimentary pore water to the bottom water. From the U pore water flux and the reconstructed U pore water concentration a mean flow rate of about 147 ± 68 cm/a can be estimated. ²³⁰Th/²³⁴U ages of chemoherm carbonates range from 7.3 to 267.6 ka. ²³⁰Th/²³⁴U ages of two chemoherms (Alvin and SE-Knoll chemoherm) correspond to time intervals of low sealevel stands in marine isotope stages (MIS) 2, 4, 5, 6, 7 and 8. This observation indicates that fluid flow at cold seep sites sensitively reflects pressure changes of the hydraulic head in the sediments. The δ¹⁸O_PDB ratios of the chemoherm carbonates support the hypothesis of precipitation during glacial times. Deviations of the chemoherm δ¹⁸O values from the marine δ¹⁸O record can be interpreted as to reflect temporally and spatially varying bottom water and/or vent fluid temperatures during carbonate precipitation between 2.6 and 8.6°C.     

The dissolution kinetics of a granite and its minerals—Implications for comparison between laboratory and field dissolution rates

The present study compares the dissolution rates of plagioclase, microcline and biotite/chlorite from a bulk granite to the dissolution rates of the same minerals in mineral-rich fractions that were separated from the granite sample. The dissolution rate of plagioclase is enhanced with time as a result of exposure of its surface sites due to the removal of an iron oxide coating. Removal of the iron coating was slower in the experiment with the bulk granite than in the mineral-rich fractions due to a higher Fe concentration from biotite dissolution. As a result, the increase in plagioclase dissolution rate was initially slower in the experiment with the bulk granite. The measured steady state dissolution rates of both plagioclase (6.2 ± 1.2 × 10⁻¹¹ mol g⁻¹ s⁻¹) and microcline (1.6 ± 0.3 × 10⁻¹¹ mol g⁻¹ s⁻¹) were the same in experiments conducted with the plagioclase-rich fraction, the alkali feldspar-rich fraction and the bulk granite.Based on the observed release rates of the major elements, we suggest that the biotite/chlorite-rich fraction dissolved non-congruently under near-equilibrium conditions. In contrast, the biotite and chlorite within the bulk granite sample dissolved congruently under far from equilibrium conditions. These differences result from variations in the degree of saturation of the solutions with respect to both the dissolving biotite/chlorite and to nontronite, which probably was precipitating during dissolution of the biotite and chlorite-rich fraction. Following drying of the bulk granite, the dissolution rate of biotite was significantly enhanced, whereas the dissolution rate of plagioclase decreased.The presence of coatings, wetting and drying cycles and near equilibrium conditions all significantly affect mineral dissolution rates in the field in comparison to the dissolution rate of fully wetted clean minerals under far from equilibrium laboratory conditions. To bridge the gap between the field and the laboratory mineral dissolution rates, these effects on dissolution rate should be further studied.     

Uranium isotopic evidence for the origin of the Bahariya iron deposits, Egypt

The economic iron ore deposits of Egypt are located at Bahariya Oasis in the Lower Middle Eocene limestone. The main iron minerals are goethite, hematite, siderite, pyrite, and jarosite. Manganese minerals are pyrolusite and manganite. Gangue minerals are barite, glauconite, gibbsite, alunite, quartz, halite, kaolinite, illite, smectite, palygorskite, and halloysite. Geochemical comparison between the ore and the Nubia sandstone showed that the ore is depleted in the residual elements (Al, Ti, V, and Ni) and enriched in the mobile elements (Fe, Mn, Zn, Ba, and U) which indicates that the Bahariya iron ore is not a lateritic deposit despite the deep weathering in this area. On the other hand, the Nubia sandstone showed depletion in the mobile elements, which demonstrates the leaching process in the Nubia Aquifer. The presence of such indicator minerals as jarosite, alunite, glauconite, gibbsite, palygorskite, and halloysite indicate that the ore was deposited under strong acidic conditions in fresh water.Isotopic analyses of the uranium in the amorphous and crystalline phases of the ore, in the country rocks, and dissolved in the Nubia Aquifer water, all support the conclusion that U and Fe were precipitated together from warm ascending groundwater. U and Fe display strong co-variation in the ore, and the ²³⁴U/²³⁸U activity ratio of the newly precipitated U in the country rock and the leached component of U in the groundwater are identical. There is only slightly more uranium in the amorphous phase than in the crystalline and only a slightly lower ²³⁴U/²³⁸U activity ratio, suggesting that the iron in the two phases have a similar origin. Comparison of the excess ²³⁴U in the water and in the total ore leads to the conclusion that the precipitation of the U, and by inference the iron, occurred within the last million years. However, that both precipitation and leaching of U have occurred over the last 300,000 years is evidenced by the extreme ²³⁰Th/²³⁴U disequilibria observed in some of the samples. Some of the amorphous depositional events have been very recent, perhaps within the last 10,000 years.     

The transport of U- and Th-series nuclides in sandy confined aquifers

Abundances of ²³⁸U, ²³⁴U, ²³²Th, ²²⁶Ra, ²²⁸Ra, ²²⁴Ra, and ²²²Rn were measured in groundwaters of the Ojo Alamo aquifer in northwest New Mexico. This is an arid area with annual precipitation of ∼22 cm. The purpose was to investigate the transport of U-Th series nuclides and their daughter products in an old, slow-moving groundwater mass as a means of understanding water-rock interactions and to compare the results with a temperate zone aquifer. It was found that ²³²Th is approximately at saturation and supports the view of Tricca et al. (2001) that Th is precipitated irreversibly upon weathering, leaving surface coatings of ²³²Th and ²³⁰Th on aquifer grains. Uranium in the aquifer waters has very high [²³⁴U/²³⁸U] ∼ 9 and low ²³⁸U concentrations. These levels can be explained by low weathering rates in the aquifer (w_238U ∼ 2 × 10⁻¹⁸ to 2 × 10⁻¹⁷s⁻¹) using a continuous flow, water-rock interaction model. The Ra isotopes are roughly in secular equilibrium despite their very different mean lifetimes. The ²²²Rn and ²²⁸Ra isotopes in the aquifer correspond to ∼10% of the net production rate of the bulk rock. This is interpreted to reflect an earlier formed irreversible surface coating of Th that provides Ra and Rn to the aquifer waters. The surface waters that appear to be feeding the aquifer have low [²³⁴U/²³⁸U] and high ²³⁸U concentrations. The flow model shows that it is not possible to obtain the high [²³⁴U/²³⁸U] and low [²³⁸U] values in the aquifer from a source like the present vadose zone input. It follows that the old aquifer waters studied cannot be fed by the present vadose zone input unless they are greatly diluted with waters with very low U concentrations. If the present sampling of vadose zone sources is representative of the present input, then this requires that there was a major change in water input with much larger rainfall some several thousand years ago. This may represent a climatic change in the Southwest.     

Uranium-series disequilibrium in tuffs from Yucca Mountain,Nevada, as evidence of pore-fluid flow over the last million years

Samples of tuff from boreholes drilled into fault zones in the Exploratory Studies Facility (ESF) and relatively unfractured rock of the Cross Drift tunnels, at Yucca Mountain, Nevada, have been analysed by U-series methods. This work is part of a project to verify the finding of fast flow-paths through the tuff to ESF level, indicated by the presence of ‘bomb’ ³⁶Cl in pore fluids. Secular radioactive equilibrium in the U decay series, (i.e. when the radioactivity ratios ²³⁴U/²³⁸U, ²³⁰Th /²³⁴U and ²²⁶Ra/²³⁰Th all equal 1.00) might be expected if the tuff samples have not experienced radionuclide loss due to rock-water interaction occurring within the last million years. However, most fractured and unfractured samples were found to have a small deficiency of ²³⁴U (weighted mean ²³⁴U/²³⁸U=0.95±0.01) and a small excess of ²³⁰Th (weighted mean ²³⁰Th/²³⁴U 1.10±0.02). The ²²⁶Ra/²³⁰Th ratios are close to secular equilibrium (weighted mean=0.94±0.07). These data indicate that ²³⁴U has been removed from the rock samples in the last ∼350 ka, probably by pore fluids. Within the precision of the measurement, it would appear that ²²⁶Ra has not been mobilized and removed from the tuff, although there may be some localised ²²⁶Ra redistribution as suggested by a few ratio values that are significantly different from 1.0. Because both fractured and unfractured tuffs show approximately the same deficiency of ²³⁴U, this indicates that pore fluids are moving equally through fractured and unfractured rock. More importantly, fractured rock appears not to be a dominant pathway for groundwater flow (otherwise the ratio would be more strongly affected and the Th and Ra isotopic ratios would likely also show disequilibrium). Application of a simple mass-balance model suggests that surface infiltration rate is over an order of magnitude greater than the rate indicated by other infiltration models and that residence time of pore fluids at ESF level is about 400 a. Processes of U sorption, precipitation and re-solution are believed to be occurring and would account for these anomalous results but have not been included in the model. Despite the difficulties, the U-series data suggest that fractured rock, specifically the Sundance and Drill Hole Wash faults, are not preferred flow paths for groundwater flowing through the Topopah Spring tuff and, by implication, rapid-flow, within 50 a, from the surface to the level of the ESF is improbable.     

Chemical weathering in a tropical watershed,Luquillo Mountains,Puerto Rico III: quartz dissolution rates

The paucity of weathering rates for quartz in the natural environment stems both from the slow rate at which quartz dissolves and the difficulty in differentiating solute Si contributed by quartz from that derived from other silicate minerals. This study, a first effort in quantifying natural rates of quartz dissolution, takes advantage of extremely rapid tropical weathering, simple regolith mineralogy, and detailed information on hydrologic and chemical transport. Quartz abundances and grain sizes are relatively constant with depth in a thick saprolite. Limited quartz dissolution is indicated by solution rounding of primary angularity and by the formation of etch pits. A low correlation of surface area (0.14 and 0.42 m² g⁻¹) with grain size indicates that internal microfractures and pitting are the principal contributors to total surface area.Pore water silica concentration increases linearly with depth. On a molar basis, between one and three quarters of pore water silica is derived from quartz with the remainder contributed from biotite weathering. Average solute Si remains thermodynamically undersaturated with respect to recently revised estimates of quartz solubility (<180 μM) but exceeds estimated critical saturation concentrations controlling the initiation of etch pit formation (>17–81 μM). Etch pitting is more abundant on grains in the upper saprolite and is associated with pore waters lower in dissolved silica. Rate constants describing quartz dissolution increase with decreasing depth (from 10^−14.5–10^−15.1 mol m⁻² s⁻¹), which correlate with both greater thermodynamic undersaturation and increasing etch pit densities. Unlike for many aluminosilicates, the calculated natural weathering rates of quartz fall slightly below the rate constants previously reported for experimental studies (10^−12.4–10^−14.2 mol m⁻² s⁻¹). This agreement reflects the structural simplicity of quartz, dilute solutes, and near-hydrologic saturation.