Art F. White

EFFECTS OF CLIMATE ON CHEMICAL WEATHERING IN WATERSHEDS

Abstract Climatic effects on chemical weathering are evaluated by correlating variations On solute concentrations and fluxes with temperature, precipitation, runoff, and evapotranspiration (ET) for a worldwide distribution of sixty-eight watersheds underlain by granitoid rock types. Stream solute concentrations are strongly correlated with proportional ET loss, and evaporative concentration makes stream solute concentrations an inappropriate surrogate for chemical weathering. Chemical fluxes are unaffected by ET, and SiO2 and Na weathering fluxes exhibit systematic increases with precipitation, runoff, and temperature. However, warm and wet watersheds produce anomalously rapid weathering rates. A proposed model that provides an improved prediction of weathering rates over climatic extremes Os the product of linear precipitation and Arrhenius temperature functions. The resulting apparent activation energies based on SiO2 and Na fluxes are 59.4 and 62.5 kJ · mol-1, respectively. The coupling between temperature and precipitation emphasizes the importance of tropical regions On global silicate weathering fluxes, and suggests it is not representative to use continental averages for temperature and precipitation On the weathering rate functions of global carbon cycling and climatic change models. Fluxes of K, Ca, and Mg exhibit no climatic correlation, implying that other processes, such as ion exchange, nutrient cycling, and variations On lithology, obscure any climatic signal. The correlation between yearly variations On precipitation and solute fluxes within individual watersheds Os stronger than the correlation between precipitation and solute fluxes of watersheds with different climatic regimes. This underscores the significance of transport-induced variability On controlling stream chemistry, and the importance of distinguishing between short-term and long-term climatic trends. No correlation exists between chemical fluxes and topographic relief or the extent of recent glaciation, implying that physical erosion rates do not have a critical influence on chemical weathering rates.

CHEMICAL WEATHERING IN A TROPICAL WATERSHED, LUQUILLO MOUNTAINS, PUERTO RICO : I. LONG-TERM VERSUS SHORT-TERM WEATHERING FLUXES

Abstract The pristine Rio Icacos watershed in the Luquillo Mountains in eastern Puerto Rico has the fastest documented weathering rate of silicate rocks on the Earth’s surface. A regolith propagation rate of 58 m Ma−1, calculated from iso-volumetric saprolite formation from quartz diorite, is comparable to the estimated denudation rate (25–50 Ma−1) but is an order of magnitude faster than the global average weathering rate (6 Ma−1). Weathering occurs in two distinct environments; plagioclase and hornblende react at the saprock interface and biotite and quartz weather in the overlying thick saprolitic regolith. These environments produce distinctly different water chemistries, with K, Mg, and Si increasing linearly with depth in saprolite porewaters and with stream waters dominated by Ca, Na, and Si. Such differences are atypical of less intense weathering in temperate watersheds. Porewater chemistry in the shallow regolith is controlled by closed-system recycling of inorganic nutrients such as K. Long-term elemental fluxes through the regolith (e.g., Si = 1.7 × 10−8 moles m−2 s−1) are calculated from mass losses based on changes in porosity and chemistry between the regolith and bedrock and from the age of the regolith surface (200 Ma). Mass losses attributed to solute fluxes are determined using a step-wise infiltration model which calculates mineral inputs to the shallow and deep saprolite porewaters and to stream water. Pressure heads decrease with depth in the shallow regolith (−2.03 m H2O m−1), indicating that both increasing capillary tension and graviometric potential control porewater infiltration. Interpolation of experimental hydraulic conductivities produces an infiltration rate of 1 m yr−1 at average field moisture saturation which is comparable with LiBr tracer tests and with base discharge from the watershed. Short term weathering fluxes calculated from solute chemistries and infiltration rates (e.g., Si = 1.4 × 10−8 moles m−2 s−1) are compared to watershed flux rates (e.g., Si = 2.7 × 10−8 moles m−2 s−1). Consistency between three independently determined sets of weathering fluxes imply that possible changes in precipitation, temperature, and vegetation over the last several hundred thousand years have not significantly impacted weathering rates in the Luquillo Mountains of Puerto Rico. This has important ramifications for tropical environments and global climate change.

Chemical weathering rates of a soil chronosequence on granitic alluvium: I. Quantification of mineralogical and surface area changes and calculation of primary silicate reaction rates

Mineral weathering rates are determined for a series of soils ranging in age from 0.2–3000 Ky developed on alluvial terraces near Merced in the Central Valley of California. Mineralogical and elemental abundances exhibit time-dependent trends documenting the chemical evolution of granitic sand to residual kaolinite and quartz. Mineral losses with time occur in the order: hornblende > plagioclase > K-feldspar. Maximum volume decreases of >50% occur in the older soils. BET surface areas of the bulk soils increase with age, as do specific surface areas of aluminosilicate mineral fractions such as plagioclase, which increases from 0.4–1.5 m2 g−1 over 600 Ky. Quartz surface areas are lower and change less with time (0.11–0.23 m2 g−1). BET surface areas correspond to increasing external surface roughness (λ = 10–600) and relatively constant internal surface area (≈ 1.3 m2 g−1). SEM observations confirm both surface pitting and development of internal porosity. A numerical model describes aluminosilicate dissolution rates as a function of changes in residual mineral abundance, grain size distributions, and mineral surface areas with time. A simple geometric treatment, assuming spherical grains and no surface roughness, predicts average dissolution rates (plagioclase, 10−17.4; K-feldspar, 10−17.8; and hornblende, 10−17.5 mol cm−1 s−1) that are constant with time and comparable to previous estimates of soil weathering. Average rates, based on BET surface area measurements and variable surface roughnesses, are much slower (plagioclase, 10−19.9; K-feldspar, 10−20.5; and hornblende 10−20.1 mol cm−2 s−1). Rates for individual soil horizons decrease by a factor of 101.5 over 3000 Ky indicating that the surface reactivities of minerals decrease as the physical surface areas increase. Rate constants based on BET estimates for the Merced soils are factors of 103–104 slower than reported experimental dissolution rates determined from freshly prepared silicates with low surface roughness (λ < 10). This study demonstrates that the utility of experimental rate constants to predict weathering in soils is limited without consideration of variable surface areas and processes that control the evolution of surface reactivity with time.

Demonstration of significant abiotic iron isotope fractionation in nature

Field and laboratory studies reveal that the mineral ferrihydrite, formed as a result of abiotic oxidation of aqueous ferrous to ferric Fe, contains Fe that is isotopically heavy relative to coexisting aqueous Fe. Because the electron transfer step of the oxidation process at pH >5 is essentially irreversible and should favor the lighter Fe isotopes in the ferric iron product, this result suggests that relatively heavy Fe isotopes are preferentially partitioned into the readily oxidized Fe(II)(OH) x (aq) species or their transition complexes prior to oxidation. The apparent Fe isotope fractionation factor, α ferrihydrite- water , depends primarily on the relative abundances of the Fe(II) (aq) species. This study demonstrates that abiotic processes can fractionate the Fe isotopes to the same extent as biotic processes, and thus Fe isotopes on their own do not provide an effective biosignature.

Kinetics of Water-Rock Interaction

Analysis of Rates of Geochemical Reactions.- Transition State Theory and Molecular Orbital Calculations Applied to Rates and Reaction Mechanisms in Geochemical Kinetics.- The Mineral-Water Interface.- Kinetics of Sorption-Desorption.- Kinetics of Mineral Dissolution.- Data Fitting Techniques with Applications to Mineral Dissolution Kinetics.- Nucleation, Growth, and Aggregation of Mineral Phases: Mechanisms and Kinetic Controls.- Microbiological Controls on Geochemical Kinetics 1: Fundamentals and Case Study on Microbial Fe(III) Oxide Reduction.- Microbiological Controls on Geochemical Kinetics 2: Case Study on Microbial Oxidation of Metal Sulfide Minerals and Future Prospects.- Quantitative Approaches to Characterizing Natural Chemical Weathering Rates.- Geochemical Kinetics and Transport.- Isotope Geochemistry as a Tool for Deciphering Kinetics of Water-Rock Interaction.- Kinetics of Global Geochemical Cycles.

The role of disseminated calcite in the chemical weathering of granitoid rocks

Accessory calcite, present at concentrations between 300 and 3000 mg kg−1, occurs in fresh granitoid rocks sampled from the Merced watershed in Yosemite National Park, CA, USA; Loch Vale in Rocky Mountain National Park CO USA; the Panola watershed, GA USA; and the Rio Icacos, Puerto Rico. Calcite occurs as fillings in microfractures, as disseminated grains within the silicate matrix, and as replacement of calcic cores in plagioclase. Flow-through column experiments, using de-ionized water saturated with 0.05 atm. CO2, produced effluents from the fresh granitoid rocks that were dominated by Ca and bicarbonate and thermodynamically saturated with calcite. During reactions up to 1.7 yr, calcite dissolution progressively decreased and was superceded by steady state dissolution of silicates, principally biotite. Mass balance calculations indicate that most calcite had been removed during this time and accounted for 57–98% of the total Ca released from these rocks. Experimental effluents from surfically weathered granitoids from the same watersheds were consistently dominated by silicate dissolution. The lack of excess Ca and alkalinity indicated that calcite had been previously removed by natural weathering. The extent of Ca enrichment in watershed discharge fluxes corresponds to the amounts of calcite exposed in granitoid rocks. High Ca/Na ratios relative to plagioclase stoichiometries indicate excess Ca in the Yosemite, Loch Vale, and other alpine watersheds in the Sierra Nevada and Rocky Mountains of the western United States. This Ca enrichment correlates with strong preferential weathering of calcite relative to plagioclase in exfoliated granitoids in glaciated terrains. In contrast, Ca/Na flux ratios are comparable to or less than the Ca/Na ratios for plagioclase in the subtropical Panola and tropical Rio Icacos watersheds, in which deeply weathered regoliths exhibit concurrent losses of calcite and much larger masses of plagioclase during transport-limited weathering. These results indicate that the weathering of accessory calcite may strongly influence Ca and alkalinity fluxes from silicate rocks during and following periods of glaciation and tectonism but is much less important for older stable geomorphic surfaces.

The effect of temperature on experimental and natural chemical weathering rates of granitoid rocks

The effects of climatic temperature variations (5–35°C) on chemical weathering are investigated both experimentally using flow-through columns containing fresh and weathered granitoid rocks and for natural granitoid weathering in watersheds based on annual solute discharge. Although experimental Na and Si effluent concentrations are significantly higher in the fresh relative to the weathered granitoids, the proportional increases in concentration with increasing temperature are similar. Si and Na exhibit comparable average apparent activation energies (Ea) of 56 and 61 kJ/mol, respectively, which are similar to those reported for experimental feldspar dissolution measured over larger temperature ranges. A coupled temperature–precipitation model, using an expanded database for solute discharge fluxes from a global distribution of 86 granitoid watersheds, produces an apparent activation energy for Si (51 kJ/mol), which is also comparable to those derived from the experimental study. This correlation reinforces evidence that temperature does significantly impact natural silicate weathering rates. Effluent K concentrations in the column study are elevated with respect to other cations compared to watershed discharge due to the rapid oxidation/dissolution of biotite. K concentrations are less sensitive to temperature, resulting in a lower average Ea value (27 kJ/mol) indicative of K loss from lower energy interlayer sites in biotite. At lower temperatures, initial cation release from biotite is significantly faster than cation release from plagioclase. This agrees with reported higher K/Na ratios in cold glacial watersheds relative to warmer temperate environments. Increased release of less radiogenic Sr from plagioclase relative to biotite at increasing temperature produces corresponding decreases in 87Sr/86Sr ratios in the column effluents. A simple mixing calculation using effluent K/Na ratios, Sr concentrations and 87Sr/86Sr ratios for biotite and plagioclase approximates stoichiometric cation ratios from biotite/plagioclase dissolution at warmer temperatures (35°C), but progressively overestimates the relative proportion of biotite with decreasing temperature. Ca, Mg, and Sr concentrations closely correlate, exhibit no consistent trends with temperature, and are controlled by trace amounts of calcite or exchange within weathered biotite. The inability of the watershed model to differentiate a climate signal for such species correlates with the lower temperature dependence observed in the experimental studies.

REDUCTION OF AQUEOUS TRANSITION METAL SPECIES ON THE SURFACES OF FE( II)-CONTAINING OXIDES

Experimental studies demonstrate that structural Fe(II) in magnetite and ilmenite heterogeneously reduce aqueous ferric, cupric, vanadate, and chromate ions at the oxide surfaces over a pH range of 1–7 at 25°C. For an aqueous transition metal m, such reactions are 3[Fe2+Fe23+]O4(magnetite)+2/nmz→4[Fe23+]O3(magnetite)+Fe2++2/nmz−n and 3[Fe2+Ti]O3(ilmenite)+2/nmz→Fe23+Ti3O9(pseudorutile)+Fe2++2/nmz−n, where z is the valance state and n is the charge transfer number. The half cell potential range for solid state oxidation [Fe(II)] → [Fe(III)] is −0.34 to −0.65 V, making structural Fe(II) a stronger reducing agent than aqueous Fe2+ (−0.77 V). Reduction rates for aqueous metal species are linear with time (up to 36 h), decrease with pH, and have rate constants between 0.1 and 3.3 × 10−10 mol m−2 s−1. Iron is released to solution both from the above reactions and from dissolution of the oxide surface. In the presence of chromate, Fe2+ is oxidized homogeneously in solution to Fe3+. X-ray photoelectron spectroscopy (XPS) denotes a Fe(III) oxide surface containing reduced Cr(III) and V(IV) species. Magnetite and ilmenite electrode potentials are insensitive to increases in divalent transition metals including Zn(II), Co(II), Mn(II), and Ni(II) and reduced V(IV) and Cr(III) but exhibit a log-linear concentration-potential response to Fe(III) and Cu(II). Complex positive electrode responses occur with increasing Cr(VI) and V(V) concentrations. Potential dynamic scans indicate that the high oxidation potential of dichromate is capable of suppressing the cathodic reductive dissolution of magnetite. Oxide electrode potentials are determined by the Fe(II)/Fe(III) composition of the oxide surface and respond to aqueous ion potentials which accelerate this oxidation process. Natural magnetite sands weathered under anoxic conditions are electrochemically reactive as demonstrated by rapid chromate reduction and the release of aqueous Fe(III) to experimental solution. In contrast, magnetite weathered under oxidizing vadose conditions show minimum reactivity toward chromate ions. The ability of Fe(II) oxides to reduce transition metals in soils and groundwaters will be strongly dependent on the redox environment.

Chemical weathering of a soil chronosequence on granitoid alluvium: II. Mineralogic and isotopic constraints on the behavior of strontium

The use of strontium isotopes to evaluate mineral weathering and identify sources of base cations in catchment waters requires an understanding of the behavior of Sr in the soil environment as a function of time. Our approach is to model the temporal evolution of 87Sr/86Sr of the cation exchange pool in a soil chronosequence developed on alluvium derived from central Sierra Nevada granitoids during the past 3 Ma. With increasing soil age, 87Sr/86Sr of ammonium-acetate extractable Sr initially decreases from values typical of K-feldspar to those of plagioclase and hornblende and then remains constant, even though plagioclase and hornblende are absent from the soils after approximately 1 Ma of weathering. The temporal variation of 87Sr/86Sr of exchangeable Sr is modeled by progressively equilibrating Sr derived from mineral weathering and atmospheric deposition with Sr on exchange sites as waters infiltrate a soil column. Observed decreases in quartz-normalized modal abundances of plagioclase, hornblende, and K-feldspar with time, and the distinct87Sr/86Sr values of these minerals can be used to calculate Sr flux from weathering reactions. Hydrobiotites in the soils have nearly constant modal abundances, chemistry, and 87Sr/86Sr over the chronosequence and provide negligible Sr input to weathering solutions. The model requires time and soil horizon-dependent changes in the amount of exchangeable Sr and the efficiency of Sr exchange, as well as a biologic cycling term. The model predicts that exchangeable Sr initially has 87Sr/86Sr identical to that of K-feldspar, and thus could be dominated by Sr leached from K-feldspar following deposition of the alluvium. The maximum value of 87Sr/86Sr observed in dilute stream waters associated with granitoids of the Yosemite region is likewise similar to that of the K-feldspars, suggesting that K-feldspar and not biotite may be the dominant source of radiogenic Sr in the streams. This study reveals that, when attempting to use Strontium isotopes to identify sources of base cations in catchment waters and biomass, both preferential leaching of Sr from minerals during incipient soil development and changing Sr exchange efficiency must be considered along with chemical contributions due to mineral dissolution.

Electrochemistry and dissolution kinetics of magnetite and ilmenite

Natural samples of magnetite and ilmenite were experimentally weathered in pH 1–7 anoxic solutions at temperatures of 2–65 °C. Reaction of magnetite is described as [Fe2+Fe23+]O4(magnetite) + 2H+ → γ[Fe23+]O3(maghemite) + Fe2+ + H2O. Dynamic polarization experiments using magnetite electrodes confirmed that this reaction is controlled by two electrochemical half cells, 3[Fe2+Fe23+]O4(magnetite) → 4γ[Fe23+]O3(maghemite) + Fe2+ + 2e− and [Fe2+Fe23+]O4(magnetite) + 8 H+ + 2e− → 3Fe2+ + 4H2O, which result in solid state Fe3+ reduction, formation of an oxidized layer and release of Fe(II) to solution. XPS data revealed that iron is present in the ferric state in the surfaces of reacted magnetite and ilmenite and that the TiFe ratio increased with reaction pH for ilmenite. Short-term (<36 h) release rates of Fe(II) were linear with time. Between pH 1 and 7, rates varied between 0.3 and 13 × 10−14 mol · cm−2 · s−1 for magnetite and 0.05 and 12.3 × 10−14 mol · cm−2 · s−1 for ilmenite. These rates are two orders of magnitude slower than electrochemical rates determined by Tafel and polarization resistance measurements. Discrepancies are due to both differences in geometric and BET surface area estimates and in the oxidation state of the mineral surface. In long-term closed-system experiments (<120 days), Fe(II) release slowed with time due to the passivation of the surfaces by increasing thicknesses of oxide surface layers. A shrinking core model, coupling surface reaction and diffusion transport, predicted that at neutral pH, the mean residence time for sand-size grains of magnetite and ilmenite will exceed 107 years. This agrees with long-term stability of these oxides in the geologic record.