Variation in calcite dissolution rates:: A fundamental problem?

A comparison of published calcite dissolution rates measured far from equilibrium at a pH of ∼ 6 and above shows well over an order of magnitude in variation. Recently published AFM step velocities extend this range further still. In an effort to understand the source of this variation, and to provide additional constraint from a new analytical approach, we have measured dissolution rates by vertical scanning interferometry. In areas of the calcite cleavage surface dominated by etch pits, our measured dissolution rate is 10^−10.95 mol/cm²/s (PCO₂ 10^−3.41 atm, pH 8.82), 5 to ∼100 times slower than published rates derived from bulk powder experiments, although similar to rates derived from AFM step velocities. On cleavage surfaces free of local etch pit development, dissolution is limited by a slow, “global” rate (10^−11.68 mol/cm²/s). Although these differences confirm the importance of etch pit (defect) distribution as a controlling mechanism in calcite dissolution, they also suggest that “bulk” calcite dissolution rates observed in powder experiments may derive substantial enhancement from grain boundaries having high step and kink density. We also observed significant rate inhibition by introduction of dissolved manganese. At 2.0 μM Mn, the rate diminished to 10^−12.4 mol/cm²/s, and the well formed rhombic etch pits that characterized dissolution in pure solution were absent. These results are in good agreement with the pattern of manganese inhibition in published AFM step velocities, assuming a step density on smooth terraces of ∼9 μm⁻¹.     

Direct determinations of the rates of rhyolite dissolution and clay formation over 52,000 years and comparison with laboratory measurements

Four porous, glass-dominated rhyolites from Kozushima Island, different in age and extent of weathering, were studied. Because the four rhyolites are homogeneously weathered to considerable depth, and because their initial chemical compositions were equal, the different rock characteristics can provide information about rates of rhyolite dissolution and clay mineral formation over ∼52,000 yr. Because glass surfaces retreat without surface roughening, surface area (measured by Brunauer-Emmett-Teller method; BET) was assumed to be approximately constant over time. The field dissolution rate, as inferred from the rate of loss of Si, was ∼6 × 10⁻¹⁹ mol cm⁻² s⁻¹. The estimated clay mineral formation rate was ∼1 × 10⁻¹⁹ mol cm⁻² s⁻¹. About 20% of dissolved Si precipitated as clays. In order to investigate the factors affecting the field dissolution rate, dissolution experiments that used powdered and block rhyolite samples were conducted. Under relevant field conditions (20°C and pH 6∼7), the rates were ∼5 × 10⁻¹⁷ and ∼5 × 10⁻¹⁸ mol cm⁻² s⁻¹ for powdered rhyolite and blocks, respectively. The dissolution rates obtained in this study decrease in the order powder > block > field. Because all surface areas were directly measured by BET, the differences are not attributable to the errors in surface area. The most plausible explanations of the slower rates are the lower degree of flushing and resultant high-solution saturation states in the pores (both in the field and in the rhyolite blocks used in experiments) plus the formation of alteration/hydrated layers at the glass surface.     

Kinetics of reaction between O2 and Mn(II) species in aqueous solutions

The objective of this research is to assess critically the experimental rate data for O₂ oxidation of dissolved Mn(II) species at 25°C and to interpret the rates in terms of the solution species of Mn(II) in natural waters. A species kinetic rate expression for parallel paths expresses the total rate of Mn(II) oxidation as Σk_i a_ij, where k_i is the rate constant of species i and a_ij is the species concentration fraction in solution j. Among the species considered in the rate expression are Mn(II) hydrolysis products, carbonate complexes, ammonia complexes, and halide and sulfate complexes, in addition to the free aqueous ion. Experiments in three different laboratory buffers and in seawater yield an apparent rate constant for Mn(II) disappearance, k_app,j ranging from 8.6 × 10⁻⁵ to 2.5 × 10⁻² (M⁻¹s⁻¹), between pH 8.03 and 9.30, respectively. Observed values of k_app exceed predictions based on Marcus outer-sphere electron transfer theory by more than four orders of magnitude, lending strong support to the proposal that Mn(II) + O₂ electron transfer follows an inner-sphere path. A multiple linear regression analysis fit of the observed rates to the species kinetic rate expression yields the following oxidation rate constants (M⁻¹s⁻¹) for the most reactive species: MnOH⁺, 1.66 × 10⁻²; Mn(OH)₂, 2.09 × 10¹; and Mn(CO₃)₂²⁻, 8.13 × 10⁻². The species kinetic rate expression accounts for the influence of pH and carbonate on oxidation rates of Mn(II), through complex formation and acid-base equilibria of both reactive and unreactive species. At pH ∼8, the greater fraction of the total rate is carried by MnOH⁺. At pH greater than ∼8.4, the species Mn(OH)₂ and Mn(CO₃)₂²⁻ make the greater contributions to the total rate.     

Transient and quasi-steady-state dissolution of biotite at 22-25°C in high pH, sodium, nitrate, and aluminate solutions

Biotite dissolution under conditions of high pH and high aluminum, sodium, and nitrate concentrations analogous to those found in tank wastes at the Hanford Site was investigated using continuously stirred flow-through reactors at 22 to 25 °C. Experiments were designed to simulate tank leaks into the Hanford vadose zone where Fe(II) from biotite is the dominant reducing agent available to immobilize certain contaminants. Both non-steady-state and steady-state dissolution kinetics were quantified; interest in non-steady-state kinetics derives from the inherently transitory nature of tank leaks. Biotite was conditioned in pH 8 solutions to simulate the alkaline environment of the Hanford sediment, and then reacted in pH 10-14 solutions, some including 0.055 M Al(NO₃)₃ and/or 2 M or 6 M NaNO₃. Initial dissolution transients (intervals of rapid release rates that decay to slower steady-state rates) showed fast preferential release of K followed by near-stoichiometric release of Si, Al, and Mg, and slower release of Fe. Each increase in pH resulted in a second transient with the greatest amounts of Si, Al, and K released at pH 14, followed by pHs 13, 12, 11, and 10. Fe release also was highest at pH 14, but unchanging at pHs 10-13 within experimental error. Transient releases at high pH are attributed to dissolution of amphoteric secondary phases such as ferrihydrite that are inferred from saturation calculations and solid analyses to form during the conditioning interval. Transient release of Si was inhibited by the presence of 0.055 M Al(NO₃)₃; the effects of Al(NO₃)₃ and NaNO₃ on the release rates of Al, Fe, Mg, and K were variable and generally outweighed by the effect of pH. Quasi-steady-state release rates were slowest at pH 11-12 (10^−12.2 mol biotite m⁻² s⁻¹ for Si) and increased in either direction in pH away from this minimum (to 10^−11.5 at pHs 8 and 14 for Si). Fe release rates at high pH were sufficient to account for observed Cr(VI) reduction at Hanford. The net release rates of the major framework cations, from which the biotite dissolution rate is inferred, may reflect the precipitation of secondary phases or the alteration of biotite to vermiculite. The most extensive solid-phase alterations were observed in Na-enriched solutions.     

The origin and mechanisms of salinization of the lower Jordan river

The chemical and isotopic (⁸⁷Sr/⁸⁶Sr, δ¹¹B, δ³⁴S_sulfate, δ¹⁸O_water, δ¹⁵N_nitrate) compositions of water from the Lower Jordan River and its major tributaries between the Sea of Galilee and the Dead Sea were determined in order to reveal the origin of the salinity of the Jordan River. We identified three separate hydrological zones along the flow of the river:

(1)

A northern section (20 km downstream of its source) where the base flow composed of diverted saline and wastewaters is modified due to discharge of shallow sulfate-rich groundwater, characterized by low ⁸⁷Sr/⁸⁶Sr (0.7072), δ³⁴S_sulfate (−2‰), high δ¹¹B (∼36‰), δ¹⁵N_nitrate (∼15‰) and high δ¹⁸O_water (−2 to-3‰) values. The shallow groundwater is derived from agricultural drainage water mixed with natural saline groundwater and discharges to both the Jordan and Yarmouk rivers. The contribution of the groundwater component in the Jordan River flow, deduced from mixing relationships of solutes and strontium isotopes, varies from 20 to 50% of the total flow.

(2)

A central zone (20-50 km downstream from its source) where salt variations are minimal and the rise of ⁸⁷Sr/⁸⁶Sr and SO₄/Cl ratios reflects predominance of eastern surface water flows.

(3)

A southern section (50-100 km downstream of its source) where the total dissolved solids of the Jordan River increase, particularly during the spring (70-80 km) and summer (80-100 km) to values as high as 11.1 g/L. Variations in the chemical and isotopic compositions of river water along the southern section suggest that the Zarqa River (⁸⁷Sr/⁸⁶Sr∼0.70865; δ¹¹B∼25‰) has a negligible affect on the Jordan River. Instead, the river quality is influenced primarily by groundwater discharge composed of sulfate-rich saline groundwater (Cl^-=31-180 mM; SO₄/Cl∼0.2-0.5; Br/Cl∼2-3×10^-3; ⁸⁷Sr/⁸⁶Sr∼0.70805; δ¹¹B∼30‰; δ¹⁵N_nitrate ∼17‰, δ³⁴S_sulfate=4-10‰), and Ca-chloride Rift valley brines (Cl^-=846-1500 mM; Br/Cl∼6-8×10^-3; ⁸⁷Sr/⁸⁶Sr∼0.7080; δ¹¹B>40‰; δ³⁴S_sulfate=4-10‰). Mixing calculations indicate that the groundwater discharged to the river is composed of varying proportions of brines and sulfate-rich saline groundwater. Solute mass balance calculations point to a ∼10% contribution of saline groundwater (Cl⁻=282 to 564 mM) to the river. A high nitrate level (up to 2.5 mM) in the groundwater suggests that drainage of wastewater derived irrigation water is an important source for the groundwater. This irrigation water appears to leach Pleistocene sediments of the Jordan Valley resulting in elevated sulfate contents and altered strontium and boron isotopic compositions of the groundwater that in turn impacts the water quality of the lower Jordan River.

     

The influence of Fe-rich coatings on the dissolution of anorthite at pH 2.6

Surface coatings are very common on mineral grains in soils but most laboratory dissolution experiments are carried out on pristine, uncoated mineral grains. An experiment designed to unambiguously isolate the effect of surface coatings on mineral dissolution from any influence of solution saturation state is reported. Two aliquots of 53 to 63 μm anorthite feldspar powder were used. One was dissolved in pH 2.6 HCl, the other in pH 2.6 FeCl₃ solution, both for ∼6000 h in flow-through reactors. An amorphous Fe-rich, Al-, Ca- and Si-free orange precipitate coated the anorthite dissolved in the FeCl₃ solution. BET surface area of the anorthite increased from 0.16 to 1.65 m² g⁻¹ in the HCl experiment and to 3.89 m² g⁻¹ in the FeCl₃ experiment. The increase in surface area in the HCl experiment was due to the formation of etch pits on the anorthite grain surface whilst the additional increase in the FeCl₃ experiment was due to the micro- and meso-porous nature of the orange precipitate. This precipitate did not inhibit or slow the dissolution of the anorthite. Steady state dissolution rates for the anorthite dissolved in the HCl and FeCl₃ were ∼2.5 and 3.2 × 10⁻¹⁰ mol_feldspar m⁻² s⁻¹ respectively. These rates are not significantly different after the cumulative uncertainty of 17% in their value due to uncertainty in the inputs parameters used in their calculation is taken into account. Results from this experiment support previous theoretical and inference-based conclusions that porous coatings should not inhibit mineral dissolution.     

Effects of pH and temperature on dimerization rate of glycine: Evaluation of favorable environmental conditions for chemical evolution of life

To evaluate favorable environmental conditions for the chemical evolution of life, we studied the effects of pH and temperature on the dimerization rate of glycine (Gly: NH₂-CH₂-COOH), one of the simplest amino acids. Gly dimerizes to form glycylglycine (GlyGly), and GlyGly further reacts to form diketopiperazine (DKP). Gly solutions with pH ranging from 3.1 to 10.9 were heated for 1-14 days at 140 °C, and changes in concentrations of Gly, GlyGly, and DKP were evaluated. At pH 9.8, the experiments were conducted at 120, 140, 160, and 180 °C. The dimerization rate of Gly was nearly constant at pH 3-7 and increased with increasing pH from 7 to 9.8 and then decreased with further increases in pH. To elucidate the reason for this pH dependency, we evaluated the role of the three dissociation states of Gly (cationic state: Gly⁺, zwitterionic state: Gly^±, and anionic state: Gly⁻). For pH >6, the dominant forms are Gly^± and Gly⁻, and the molar fraction of Gly^± decreases and that of Gly⁻ increases with increasing pH. The dimerization rate was determined for each dissociation state. The reaction between Gly^± and Gly⁻ was found to be the fastest; the rate constant of the reaction between Gly^± and Gly⁻ was 10 times the size of that between Gly⁻ and Gly⁻ and 98 times that between Gly^± and Gly^±. The dimerization rate became greatest at pH 9.8 because the molar fractions of Gly^± and Gly⁻ are approximately equal at this pH. The dimerization rate increased with temperature, and an activation energy of 88 kJ mol⁻¹ was obtained. Based on these results and previous reports on the stability of amino acids under hydrothermal conditions, we determined that Gly dimerizes most efficiently under alkaline pH (∼9.8) at about 150 °C.     

Does reactive surface area depend on grain size? Results from pH 3, 25 °C far-from-equilibrium flow-through dissolution experiments on anorthite and biotite

Laboratory determined mineral weathering rates need to be normalised to allow their extrapolation to natural systems. The principle normalisation terms used in the literature are mass, and geometric- and BET specific surface area (SSA). The purpose of this study was to determine how dissolution rates normalised to these terms vary with grain size. Different size fractions of anorthite and biotite ranging from 180-150 to 20-10 μm were dissolved in pH 3, HCl at 25 °C in flow through reactors under far from equilibrium conditions. Steady state dissolution rates after 5376 h (anorthite) and 4992 h (biotite) were calculated from Si concentrations and were normalised to initial- and final- mass and geometric-, geometric edge- (biotite), and BET SSA. For anorthite, rates normalised to initial- and final-BET SSA ranged from 0.33 to 2.77 × 10⁻¹⁰ mol_feldspar m⁻² s⁻¹, rates normalised to initial- and final-geometric SSA ranged from 5.74 to 8.88 × 10⁻¹⁰ mol_feldspar m⁻² s⁻¹ and rates normalised to initial- and final-mass ranged from 0.11 to 1.65 mol_feldspar g⁻¹ s⁻¹. For biotite, rates normalised to initial- and final-BET SSA ranged from 1.02 to 2.03 × 10⁻¹² mol_biotite m⁻² s⁻¹, rates normalised to initial- and final-geometric SSA ranged from 3.26 to 16.21 × 10⁻¹² mol_biotite m⁻² s⁻¹, rates normalised to initial- and final-geometric edge SSA ranged from 59.46 to 111.32 × 10⁻¹² mol_biotite m⁻² s⁻¹ and rates normalised to initial- and final-mass ranged from 0.81 to 6.93 × 10⁻¹² mol_biotite g⁻¹ s⁻¹. For all normalising terms rates varied significantly (p ? 0.05) with grain size. The normalising terms which gave least variation in dissolution rate between grain sizes for anorthite were initial BET SSA and initial- and final-geometric SSA. This is consistent with: (1) dissolution being dominated by the slower dissolving but area dominant non-etched surfaces of the grains and, (2) the walls of etch pits and other dissolution features being relatively unreactive. These steady state normalised dissolution rates are likely to be constant with time. Normalisation to final BET SSA did not give constant ratios across grain size due to a non-uniform distribution of dissolution features. After dissolution coarser grains had a greater density of dissolution features with BET-measurable but unreactive wall surface area than the finer grains. The normalising term which gave the least variation in dissolution rates between grain sizes for biotite was initial BET SSA. Initial- and final-geometric edge SSA and final BET SSA gave the next least varied rates. The basal surfaces dissolved sufficiently rapidly to influence bulk dissolution rate and prevent geometric edge SSA normalised dissolution rates showing the least variation. Simple modelling indicated that biotite grain edges dissolved 71-132 times faster than basal surfaces. In this experiment, initial BET SSA best integrated the different areas and reactivities of the edge and basal surfaces of biotite. Steady state dissolution rates are likely to vary with time as dissolution alters the ratio of edge to basal surface area. Therefore they would be more properly termed pseudo-steady state rates, only appearing constant because the time period over which they were measured (1512 h) was less than the time period over which they would change significantly.     

Early diagenetic remineralization of sedimentary organic C in the Gulf of Papua deltaic complex (Papua New Guinea): Net loss of terrestrial C and diagenetic fractionation of C isotopes

Oceania supplies ∼40% of the global riverine flux of organic carbon, approximately half of which is injected onto broad continental shelves and processed in shallow deltaic systems. The Gulf of Papua, on the south coast of the large island of New Guinea, is one such deltaic clinoform complex. It receives ∼4 Mt yr⁻¹ particulate terrestrial organic carbon with initial particle C_org loading ∼0.7 mg m⁻². C_org loading is reduced to ∼0.3 mg m⁻² in the topset-upper foreset zones of the delta despite additional inputs of mangrove and planktonic detritus, and high net sediment accumulation rates of 1-4 cm yr⁻¹. Carbon isotopic analyses (δ¹³C, Δ¹⁴C) of ΣCO₂ and C_org demonstrate rapid (<100 yr) remineralization of both terrestrial (δ¹³C <−28.6) and marine C_org (δ¹³C ∼−20.5) ranging in average age from modern (bomb) (Δ¹⁴C ∼60) to ∼1000 yr (Δ¹⁴C ∼−140). Efficient and rapid remineralization in the topset-upper foreset zone is promoted by frequent physical reworking, bioturbation, exposure, and reoxidation of deposits. The seafloor in these regions, particularly <20 m, apparently functions as a periodically mixed, suboxic batch reactor dominated by microbial biomass. Although terrestrial sources can be the primary metabolic substrates at inshore sites, relatively young marine C_org often preferentially dominates pore water ΣCO₂ relative to bulk C_org in the upper foreset. Thus a small quantity of young, rapidly recycled marine organic material is often superimposed on a generally older, less reactive terrestrial background. Whereas the pore water ΣCO₂ reflects both rapidly cycled marine and terrestrial sources, terrestrial material dominates the slower overall net loss of C_org from particles in the topset-upper foreset zone (i.e. recycled marine C_org leaves little residue). Preferential utilization of C_org subpools and diagenetic fractionation of C isotopes supports the reactive continuum model as a conceptual basis for net decomposition kinetics. Early diagenetic fractionation of C isotopes relative to the bulk sedimentary C_org composition can produce changes in ¹⁴C activity independent of radioactive decay. In the Gulf of Papua topset-upper foreset, Δ¹⁴C of pore water ΣCO₂ averaged ∼ 300‰ greater than C_org sediment between ∼1-3 m depth in deposits. Diagenetic fractionation and decomposition aging of sedimentary C_org compromises simple application of ¹⁴C for determination of sediment accumulation rates in diagenetically reactive deposits.     

Effect of groundwater pH and ionic strength on strontium sorption in aquifer sediments: Implications for Sr mobility at contaminated nuclear sites

Strontium-90 is a beta emitting radionuclide produced during nuclear fission, and is a problem contaminant at many nuclear facilities. Transport of ⁹⁰Sr in groundwaters is primarily controlled by sorption reactions with aquifer sediments. The extent of sorption is controlled by the geochemistry of the groundwater and sediment mineralogy. Here, batch sorption experiments were used to examine the sorption behaviour of ⁹⁰Sr in sediment–water systems representative of the UK Sellafield nuclear site based on groundwater and contaminant fluid compositions. In experiments with low ionic strength groundwaters (<0.01 mol L⁻¹), pH variation is the main control on sorption. The sorption edge for ⁹⁰Sr was observed between pH 4 and 6 with maximum sorption occurring (K_d ∼ 10³ L kg⁻¹) at pH 6–8. At ionic strengths above 10 mmol L⁻¹, and at pH values between 6 and 8, cation exchange processes reduced ⁹⁰Sr uptake to the sediment. This exchange process explains the lower ⁹⁰Sr sorption (K_d ∼ 40 L kg⁻¹) in the presence of artificial Magnox tank liquor (IS = 29 mmol L⁻¹). Strontium K-edge EXAFS spectra collected from sediments incubated with Sr²⁺ in either HCO₃-buffered groundwater or artificial Magnox tank liquor, revealed a coordination environment of ∼9 O atoms at 2.58–2.61 Å after 10 days. This is equivalent to the Sr²⁺ hydration sphere for the aqueous ion and indicates that Sr occurs primarily in outer sphere sorption complexes. No change was observed in the Sr sorption environment with EXAFS analysis after 365 days incubation. Sequential extractions performed on sediments after 365 days also found that ∼80% of solid associated ⁹⁰Sr was exchangeable with 1 M MgCl₂ in all experiments. These results suggest that over long periods, ⁹⁰Sr in contaminated sediments will remain primarily in weakly bound surface complexes. Therefore, if groundwater ionic strength increases (e.g. by saline intrusion related to sea level rise or by design during site remediation) then substantial remobilisation of ⁹⁰Sr is to be expected.     

Kinetic and equilibrium Fe isotope fractionation between aqueous Fe(II) and Fe(III)

Equilibrium and kinetic Fe isotope fractionation between aqueous ferrous and ferric species measured over a range of chloride concentrations (0, 11, 110 mM Cl⁻) and at two temperatures (0 and 22°C) indicate that Fe isotope fractionation is a function of temperature, but independent of chloride contents over the range studied. Using ⁵⁷Fe-enriched tracer experiments the kinetics of isotopic exchange can be fit by a second-order rate equation, or a first-order equation with respect to both ferrous and ferric iron. The exchange is rapid at 22°C, ∼60-80% complete within 5 seconds, whereas at 0°C, exchange rates are about an order of magnitude slower. Isotopic exchange rates vary with chloride contents, where ferrous-ferric isotope exchange rates were ∼25 to 40% slower in the 11 mM HCl solution compared to the 0 mM Cl⁻ (∼10 mM HNO₃) solutions; isotope exchange rates are comparable in the 0 and 110 mM Cl⁻ solutions.The average measured equilibrium isotope fractionations, Δ_{Fe(III)-Fe(II)}, in 0, 11, and 111 mM Cl⁻ solutions at 22°C are identical within experimental error at +2.76±0.09, +2.87±0.22, and +2.76±0.06 ‰, respectively. This is very similar to the value measured by Johnson et al. (2002a) in dilute HCl solutions. At 0°C, the average measured Δ_{Fe(III)-Fe(II)} fractionations are +3.25±0.38, +3.51±0.14 and +3.56±0.16 ‰ for 0, 11, and 111 mM Cl⁻ solutions. Assessment of the effects of partial re-equilibration on isotope fractionation during species separation suggests that the measured isotope fractionations are on average too low by ∼0.20 ‰ and ∼0.13 ‰ for the 22°C and 0°C experiments, respectively. Using corrected fractionation factors, we can define the temperature dependence of the isotope fractionation from 0°C to 22°C as: where the isotopic fractionation is independent of Cl⁻ contents over the range used in these experiments. These results confirm that the Fe(III)-Fe(II) fractionation is approximately half that predicted from spectroscopic data, and suggests that, at least in moderate Cl⁻ contents, the isotopic fractionation is relatively insensitive to Fe-Cl speciation.     

Solubility of nanocrystalline scorodite and amorphous ferric arsenate: Implications for stabilization of arsenic in mine wastes

Solubility experiments were performed on nanocrystalline scorodite and amorphous ferric arsenate. Nanocrystalline scorodite occurs as stubby prismatic crystals measuring about 50 nm and having a specific surface area of 39.88 ± 0.07 m²/g whereas ferric arsenate is amorphous and occurs as aggregated clusters measuring about 50–100 nm with a specific surface area of 17.95 ± 0.19 m²/g. Similar to its crystalline counterpart, nanocrystalline scorodite has a solubility of about 0.25 mg/L at around pH 3–4 but has increased solubilities at low and high pH (i.e. <2 and >6). Nanocrystalline scorodite dissolves incongruently at about pH > 2.5 whereas ferric arsenate dissolution is incongruent at all the pH ranges tested (pH 2–5). It appears that the solubility of scorodite is not influenced by particle size. The dissolution rate of nanocrystalline scorodite is 2.64 × 10⁻¹⁰ mol m⁻² s⁻¹ at pH 1 and 3.25 × 10⁻¹¹ mol m⁻² s⁻¹ at pH 2. These rates are 3–4 orders of magnitude slower than the oxidative dissolution of pyrite and 5 orders of magnitude slower than that of arsenopyrite. Ferric arsenate dissolution rates range from 6.14 × 10⁻⁹ mol m⁻² s⁻¹ at pH 2 to 1.66 × 10⁻⁹ mol m⁻² s⁻¹ at pH 5. Among the common As minerals, scorodite has the lowest solubility and dissolution rate. Whereas ferric arsenate is not a suitable compound for As control in mine effluents, nanocrystalline scorodite that can be easily precipitated at ambient pressure and temperature conditions would be satisfactory in meeting the regulatory guidelines at pH 3–4.     

Chemical erosion in the eastern Himalaya: Major ion composition of the Brahmaputra and δC of dissolved inorganic carbon

Major ion composition of waters, δ¹³C of its DIC (dissolved inorganic carbon), and the clay mineral composition of bank sediments in the Brahmaputra River System (draining India and Bangladesh) have been measured to understand chemical weathering and erosion and the factors controlling these processes in the eastern Himalaya. The time-series samples, collected biweekly at Guwahati, from the Brahmaputra mainstream, were also analyzed for the major ion composition. Clay mineralogy and chemical index of alteration (CIA) of sediments suggest that weathering intensity is relatively poor in comparison to that in the Ganga basin. This is attributed to higher runoff and associated physical erosion occurring in the Brahmaputra basin. The results of this study show, for the first time, spatial and temporal variations in chemical and silicate erosion rates in the Brahmaputra basin. The subbasins of the Brahmaputra watershed exhibit chemical erosion rates varying by about an order of magnitude. The Eastern Syntaxis basin dominates the erosion with a rate of ∼300 t km⁻² y⁻¹, one of the highest among the world river basins and comparable to those reported for some of the basaltic terrains. In contrast, the flat, cold, and relatively more arid Tibetan basin undergoes much slower chemical erosion (∼40 t km⁻² y⁻¹). The abundance of total dissolved solids (TDS, 102-203 mg/L) in the time-series samples collected over a period of one year shows variations in accordance with the annual discharge, except one of them, cause for which is attributable to flash floods. Na* (Na corrected for cyclic component) shows a strong positive correlation with Si, indicating their common source: silicate weathering. Estimates of silicate cations (Na_sil+K_sil+Ca_sil+Mg_sil) suggest that about half of the dissolved cations in the Brahmaputra are derived from silicates, a proportion higher than that for the Ganga system. The CO₂ consumption rate due to silicate weathering in the Brahmaputra watershed is ∼6 × 10⁵ moles km⁻² y⁻¹; whereas that in the Eastern Syntaxis subbasin is ∼19 × 10⁵ moles km⁻² y⁻¹, similar to the estimates for some of the basaltic terrains. This study suggests that the Eastern Syntaxis basin of the Brahmaputra is one of most intensely chemically eroding regions of the globe; and that runoff and physical erosion are the controlling factors of chemical erosion in the eastern Himalaya.     

Does variscite control phosphate availability in acidic natural waters? An experimental study of variscite dissolution rates

The dissolution rates of natural, well crystallized variscite (AlPO₄·2H₂O) were determined from the evolution of aqueous Al and P concentrations in closed and open-system mixed-flow reactors at 25 °C and pH from 1.5 to 9.0. Measured dissolution rates decrease with increasing pH, from 6 × 10⁻¹⁶ mol/cm²/s at pH 1.5 to 5 × 10⁻¹⁷ mol/cm²/s at pH 5.89, and then increase with increasing pH to 4 × 10⁻¹⁶ mol/cm²/s at pH 9.0. Geochemical modeling calculations, performed using measured dissolution rates, indicate that it would take no more than a few weeks or months to equilibrate a mildly acidic, Al and P-free solution with variscite. Hence, variscite can buffer aqueous phosphate concentrations in mildly acidic near surface environments. This conclusion is confirmed by consideration of the compositions of natural waters.     

Kinetics of iron (II) oxidation determined in the field

This paper presents the results of extensive field trials measuring rates of Fe(II) oxidation at a number of Fe-bearing mine drainage discharges in the UK. Batch experiments were carried out with samples taken at regular intervals and Fe(II) concentration determined spectrophotometrically using 2′2-bipyridyl as the complexing agent. Initial concentrations for Fe(II) were 5.65-76.5 mg/L. Temperature, pH and dissolved O₂ (DO) were logged every 10 s, with pH at the start of the experiments in the range 5.64-6.95 and alkalinity ranging from 73 to 741 mg/L CaCO₃ equivalent. A numerical model based on a fourth order Runge-Kutta method was developed to calculate values for k₁, the rate constant for homogeneous oxidation, from the experimental data. The measured values of pH, temperature, [Fe(II)] and DO were input into the model with resulting values for k₁ found to be in the range 2.7 × 10¹⁴-2.7 × 10¹⁶ M⁻² atm⁻¹ min⁻¹. These values for k₁ are 1-3 orders of magnitude higher than previously reported for laboratory studies at a similar pH. Comparison of the observed Fe(II) oxidation rates to data published by other authors show a good correlation with heterogenous oxidation rates and may indicate the importance of autocatalysis in these systems. These higher than expected rates of Fe oxidation could have a significant impact on the design of treatment schemes for the remediation of mine drainage and other Fe-bearing ground waters in the future.     

Atmospheric deposition of Be and Be in New Zealand rain (1996-98)

Beryllium isotope concentrations were determined in monthly rainfall collections at three sites across New Zealand (36 to 45° S), from October 1996 to November 1998. At the northern sites of Leigh (near Auckland) and Gracefield (near Wellington), ⁷Be and ¹⁰Be concentrations are relatively constant at 1.2 to 1.4 × 10⁷ atoms kg⁻¹ rain and 2.1 to 2.6 × 10⁷ atoms kg⁻¹ rain, respectively. These concentrations correspond to annual flux rates of ∼15 and ∼27 × 10⁹ atoms m⁻² y⁻¹, respectively. At the southern site of Dunedin, concentrations are similar to those at the northern sites, but flux rates are significantly lower at ∼ 9 and ∼19 × 10⁹ atoms m⁻² y⁻¹, respectively, because of lower average rainfall east of the Southern Alps mountain range. These results are broadly similar to those reported by Brown et al. (1989) and Knies et al. (1994) for rain from midlatitude sites in the USA sampled from 1986 to 1994. Mean ⁷Be/¹⁰Be ratios for New Zealand (0.47 to 0.61) are, however, significantly lower than for the USA (0.69 to 0.78), due in part to the addition of ¹⁰Be from re-suspended dust. Subtraction of the dust component increases the New Zealand ⁷Be/¹⁰Be ratios to 0.70 (Leigh), 0.65 (Gracefield) and 0.50 (Dunedin). The adjusted results provide evidence for transfer of older stratospheric air to the troposphere in late-spring-summer, an effect which is strongest in the north. The overall reduction of ⁷Be/¹⁰Be from north to south implies an increase in residence time from ∼ 80 to ∼100 d for Be isotopes in the atmosphere above New Zealand.     

Aeration to degas CO2, increase pH,and increase iron oxidation rates for efficient treatment of net alkaline mine drainage

Passive treatment systems for mine drainage use no energy other than gravity, but they require greater area than active treatment systems. Researchers are considering “hybrid” systems that have passive and active components for increased efficiency, especially where space limitations render passive-only technology ineffective. Flow-through reactor field experiments were conducted at two large net-alkaline anthracite mine discharges in central Pennsylvania. Assuming an Fe removal rate of 20 g m⁻² day⁻¹ and Fe loading from field data, 3.6 × 10³ and 3.0 × 10⁴ m² oxidation ponds would be required for the passive treatment of Site 21 and Packer 5 discharges, respectively. However, only a small area is available at each site. This paper demonstrates aeration to drive off CO₂, increase pH, and increase Fe(II) oxidation rates, enabling treatment within a small area compared to passive treatment methods, and introduces a geochemical model to accurately predict these rates as well as semi-passive treatment system sizing parameters. Both net-alkaline discharges were suboxic with a pH of ≈5.7, Fe(II) concentration of ≈16 mg L⁻¹, and low Mn and Al concentrations. Flow rates were ≈4000 L min⁻¹ at Site 21 and 15,000 L min⁻¹ at Packer 5. Three-h aeration experiments with flow rates scaled to a 14-L reactor resulted in pH increases from 5.7 to greater than 7, temperature increases from 12 to 22 °C, dissolved O₂ increases to saturation with respect to the atmosphere, and Fe(II) concentration decreases from 16 to <0.05 mg L⁻¹. A 17,000-L pilot-scale reactor at Site 21 produced similar results although aeration was not as complete as in the smaller reactor. Two non-aerated experiments at Site 21 with 13 and 25-h run times resulted in pH changes of ?0.2 and Fe(II) concentration decreases of less than 3 mg L⁻¹.     

Chemical weathering in the Krishna Basin and Western Ghats of the Deccan Traps, India: Rates of basalt weathering and their controls

Rates of chemical and silicate weathering of the Deccan Trap basalts, India, have been determined through major ion measurements in the headwaters of the Krishna and the Bhima rivers, their tributaries, and the west flowing streams of the Western Ghats, all of which flow almost entirely through the Deccan basalts.Samples (n = 63) for this study were collected from 23 rivers during two consecutive monsoon seasons of 2001 and 2002. The Total dissolved solid (TDS) in the samples range from 27 to 640 mg l⁻¹. The rivers draining the Western Ghats that flow through patches of cation deficient lateritic soils have lower TDS (average: 74 mg l⁻¹), whereas the Bhima (except at origin) and its tributaries that seem to receive Na, Cl, and SO₄ from saline soils and anthropogenic inputs have values in excess of 170 mg l⁻¹. Many of the rivers sampled are supersaturated with respect to calcite. The chemical weathering rates (CWR) of “selected” basins, which exclude rivers supersaturated in calcite and which have high Cl and SO₄, are in range of ∼3 to ∼60 t km⁻² y⁻¹. This yields an area-weighted average CWR of ∼16 t km⁻² y⁻¹ for the Deccan Traps. This is a factor of ∼2 lower than that reported for the Narmada-Tapti-Wainganga (NTW) systems draining the more northern regions of the Deccan. The difference can be because of (i) natural variations in CWR among the different basins of the Deccan, (ii) “selection” of river basin for CWR calculation in this study, and (iii) possible contribution of major ions from sources, in addition to basalts, to rivers of the northern Deccan Traps.Silicate weathering rates (SWR) in the selected basins calculated using dissolved Mg as an index varies between ∼3 to ∼60 t km⁻² y⁻¹, nearly identical to their CWR. The Ca/Mg and Na/Mg in these rivers, after correcting for rain input, are quite similar to those in average basalts of the region, suggesting near congruent release of Ca, Mg, and Na from basalts to rivers. Comparison of calculated and measured silicate-Ca in these rivers indicates that at most ∼30% of Ca can be of nonsilicate origin, a likely source being carbonates in basalts and sediments.The chemical and silicate weathering rates of the west flowing rivers of the Deccan are ∼4 times higher than the east flowing rivers. This difference is due to the correspondingly higher rainfall and runoff in the western region and thus reemphasises the dominant role of runoff in regulating weathering rates. The silicon weathering rate (SWR) in the Krishna Basin is ∼15 t km⁻² y⁻¹, within a factor of ∼2 to those in the Yamuna, Bhagirathi, and Alaknanda basins of the Himalaya, suggesting that under favourable conditions (intense physical weathering, high runoff) granites and the other silicates in the Himalaya weather at rates similar to those of Deccan basalts. The CO₂ consumption rate for the Deccan is deduced to be ∼3.6 × 10⁵ moles km⁻² y⁻¹ based on the SWR. The rate, though, is two to three times lower than reported for the NTW rivers system; it still reinforces the earlier findings that, in general, basalts weather more rapidly than other silicates and that they significantly influence the atmospheric CO₂ budget on long-term scales.