Confirmation of mass-independent Ni isotopic variability in iron meteorites

We report high-precision analyses of internally-normalised Ni isotope ratios in 12 bulk iron meteorites. Our measurements of ⁶⁰Ni/⁶¹Ni, ⁶²Ni/⁶¹Ni and ⁶⁴Ni/⁶¹Ni normalised to ⁵⁸Ni/⁶¹Ni and expressed in parts per ten thousand (?) relative to NIST SRM 986 as and , vary by 0.146, 0.228 and 0.687, respectively. The precision on a typical analysis is 0.03?, 0.05? and 0.08? for , and , respectively, which is comparable to our sample reproducibility. We show that this ‘mass-independent’ Ni isotope variability cannot be ascribed to interferences, inaccurate correction of instrumental or natural mass-dependent fractionation, fractionation controlled by nuclear field shift effects, nor the influence of cosmic ray spallation. These results thus document the presence of mass-independent Ni isotopic heterogeneity in bulk meteoritic samples, as previously proposed by Regelous et al. (2008) (EPSL 272, 330-338), but our new analyses are more precise and include determination of ⁶⁴Ni. Intriguingly, we find that terrestrial materials do not yield homogenous internally-normalised Ni isotope compositions, which, as pointed out by Young et al. (2002) (GCA 66, 1095-1104), may be the expected result of using the exponential (kinetic) law and atomic masses to normalise all fractionation processes. The certified Ni isotope reference material NIST SRM 986 defines zero in this study, while appropriate ratios for the bulk silicate Earth are given by the peridotites JP-1 and DTS-2 and, relative to NIST SRM 986, yield deviations in , and of −0.006?, 0.036? and 0.119?, respectively. There is a strong positive correlation between and in iron meteorites analyses, with a slope of 3.03 ± 0.71. The variations of Ni isotope anomalies in iron meteorites are consistent with heterogeneous distribution of a nucleosynthetic component from a type Ia supernova into the proto-solar nebula.     

The ionic strength dependence of lead (II) carbonate complexation in perchlorate media

Lead speciation in many aqueous geochemical systems is dominated by carbonate complexation. However, direct observations of Pb²⁺ complexation by carbonate ions are few in number. This work represents the first investigation of the equilibrium over a range of ionic strength. Through spectrophotometric observations of formation at 25 °C in NaHCO₃-NaClO₄ solutions, formation constants of the form were determined between 0.001 and 5.0 molal ionic strength. Formation constant results were well represented by the equation:

     

Fluoride geochemistry of thermal waters in Yellowstone National Park: I. Aqueous fluoride speciation

Thermal water samples from Yellowstone National Park (YNP) have a wide range of pH (1-10), temperature, and high concentrations of fluoride (up to 50 mg/l). High fluoride concentrations are found in waters with field pH higher than 6 (except those in Crater Hills) and temperatures higher than 50 °C based on data from more than 750 water samples covering most thermal areas in YNP from 1975 to 2008. In this study, more than 140 water samples from YNP collected in 2006-2009 were analyzed for free-fluoride activity by ion-selective electrode (ISE) method as an independent check on the reliability of fluoride speciation calculations. The free to total fluoride concentration ratio ranged from <1% at low pH values to >99% at high pH. The wide range in fluoride activity can be explained by strong complexing with H⁺ and Al³⁺ under acidic conditions and lack of complexing under basic conditions. Differences between the free-fluoride activities calculated with the WATEQ4F code and those measured by ISE were within 0.3-30% for more than 90% of samples at or above 10⁻⁶ molar, providing corroboration for chemical speciation models for a wide range of pH and chemistry of YNP thermal waters. Calculated speciation results show that free fluoride, F⁻, and major complexes (, AlF²⁺, and ) account for more than 95% of total fluoride. Occasionally, some complex species like , FeF²⁺, , MgF⁺ and may comprise 1-10% when the concentrations of the appropriate components are high. According to the simulation results by PHREEQC and calculated results, the ratio of main fluoride species to total fluoride varies as a function of pH and the concentrations and ratios of F and Al.     

Effects of magnesium ions on near-equilibrium calcite dissolution: Step kinetics and morphology

Dissolution kinetics at the aqueous solution-calcite interface at 50 °C were investigated using in situ atomic force microscopy (AFM) to reveal the influence of magnesium concentration and solution saturation state on calcite dissolution kinetics and surface morphology. Under near-equilibrium conditions, dissolved Mg²⁺ displayed negligible inhibitory effects on calcite dissolution even at concentrations of . Upon the introduction of , the solution saturation state with respect to calcite, , acted as a “switch” for magnesium inhibition whereby no significant changes in step kinetics were observed at Ω_calcite<0.2, whereas a sudden inhibition from Mg²⁺ was activated at Ω_calcite?0.2. The presence of the Ω-switch in dissolution kinetics indicates the presence of critical undersaturation in accordance with thermodynamic principles. The etch pits formed in solutions with exhibited a unique distorted rhombic profile, different from those formed in Mg-free solutions and in de-ionized water. Such unique etch pit morphology may be associated with the anisotropy in net detachment rates of counter-propagating kink sites upon the addition of Mg²⁺.     

Trace element partitioning between ilmenite, armalcolite and anhydrous silicate melt: Implications for the formation of lunar high-Ti mare basalts

We performed a series of experiments at high pressures and temperatures to determine the partitioning of a wide range of trace elements between ilmenite (Ilm), armalcolite (Arm) and anhydrous lunar silicate melt, to constrain geochemical models of the formation of titanium-rich melts in the Moon. Experiments were performed in graphite-lined platinum capsules at pressures and temperatures ranging from 1.1 to 2.3 GPa and 1300-1400 °C using a synthetic Ti-enriched Apollo ‘black glass’ composition in the CaO-FeO-MgO-Al₂O₃-TiO₂-SiO₂ system. Ilmenite-melt and armalcolite-melt partition coefficients (D) show highly incompatible values for the rare earth elements (REE) with the light REE more incompatible compared to the heavy REE ( 0.0020 ± 0.0010 to 0.069 ± 0.010 for ilmenite; 0.0048 ± 0.0023 to 0.041 ± 0.008 for armalcolite). D values for the high field strength elements vary from highly incompatible for Th, U and to a lesser extent W (for ilmenite: 0.0013 ± 0.0008, 0.0035 ± 0.0015 and 0.039 ± 0.005, and for armalcolite 0.008 ± 0.003, 0.0048 ± 0.0022 and 0.062 ± 0.03), to mildly incompatible for Nb, Ta, Zr, and Hf (e.g. 0.28 ± 0.05 and : 0.76 ± 0.07). Both minerals fractionate the high field strength elements with D_Ta/D_Nb and D_Hf/D_Zr between 1.3 and 1.6 for ilmenite and 1.3 and 1.4 for armalcolite. Armalcolite is slightly more efficient at fractionating Hf from W during lunar magma ocean crystallisation, with D_Hf/D_W = 12-13 compared to 6.7-7.5 for ilmenite. The transition metals vary from mildly incompatible to compatible, with the highest compatibilities for Cr in ilmenite (D ∼ 7.5) and V in armalcolite (D ∼ 8.1). D values show no clear variation with pressure in the small range covered.Crystal lattice strain modelling of D values for di-, tri- and tetravalent trace elements shows that in ilmenite, divalent elements prefer to substitute for Fe while armalcolite data suggest REE replacing Mg. Tetravalent cations appear to preferentially substitute for Ti in both minerals, with the exception of Th and U that likely substitute for the larger Fe or Mg cations. Crystal lattice strain modelling is also used to identify and correct for very small (∼0.3 wt.%) melt contamination of trace element concentration determinations in crystals.Our results are used to model the Lu-Hf-Ti concentrations of lunar high-Ti mare basalts. The combination of their subchondritic Lu/Hf ratios and high TiO₂ contents requires preferential dissolution of ilmenite or armalcolite from late-stage, lunar magma ocean cumulates into low-Ti partial melts of deeper pyroxene-rich cumulates.     

A new value for the stable oxygen isotope fractionation between dissolved sulfate ion and water

Although the stable oxygen isotope fractionation between dissolved sulfate ion and H₂O (hereafter ) is of physico-chemical and biogeochemical significance, no experimental value has been established until present. The primary reason being that uncatalyzed oxygen exchange between and H₂O is extremely slow, taking ∼10⁵ years at room temperature. For lack of a better approach, values of 16‰ and 31‰ at 25 °C have been assumed in the past, based on theoretical ‘gas-phase’ calculations and extrapolation of laboratory results obtained at temperatures >75 °C that actually pertain to the bisulfate system. Here I use novel quantum-chemistry calculations, which take into account detailed solute-water interactions to establish a new value for of 23‰ at 25 °C. The results of the corresponding calculations for the bisulfate ion are in agreement with observations. The new theoretical values show that sediment -data, which reflect oxygen isotope equilibration between sulfate and ambient water during microbial sulfate reduction, are consistent with the abiotic equilibrium between and water.     

New indices and calibrations derived from the distribution of crenarchaeal isoprenoid tetraether lipids: Implications for past sea surface temperature reconstructions

Several studies have shown that there is a strong relationship between the distribution of crenarchaeotal isoprenoid glycerol dibiphytanyl glycerol tetraethers (GDGTs) and sea surface temperature (SST). Based on this, a ratio of certain GDGTs, called TEX₈₆ (TetraEther indeX of tetraethers consisting of 86 carbon atoms), was developed as a SST proxy. In this study, we determined the distribution of crenarchaeotal isoprenoid GDGTs in 116 core-top sediments mostly from (sub)polar oceans and combined these data with previously published core-top data. Using this extended global core-top dataset (n = 426), we re-assessed the relationship of crenarchaeal isoprenoid GDGTs with SST. We excluded data from the Red Sea from the global core-top dataset to define new indices and calibration models, as the Red Sea with its elevated salinity appeared to behave differently compared to other parts of the oceans. We tested our new indices and calibration models on three different paleo datasets, representing different temperature ranges. Our results indicate that the crenarchaeol regio-isomer plays a more important role for temperature adaptation in (sub)tropical oceans than in (sub)polar oceans, suggesting that there may be differences in membrane adaptation of the resident crenarchaeotal communities at different temperatures. We, therefore, suggest to apply two different calibration models. For the whole calibration temperature range (−3 to 30 °C), a modified version of TEX₈₆ with a logarithmic function which does not include the crenarchaeol regio-isomer, called , is shown to correlate best with SST: (r² = 0.86, n=396, p <0.0001). Application of on sediments from the subpolar Southern Ocean results in realistic absolute SST estimates and a similar SST trend compared to a diatom SST record from the same core. , which is defined as the logarithmic function of TEX₈₆, yields the best correlation with SST, when the data from the (sub)polar oceans are removed: (r² = 00.87, n = 255, p < 0.0001). Furthermore, gives the best correlation for mescosm data with temperatures ranging between 10 and 46 °C. For Quaternary sediments from the tropical Arabian Sea, both and yield similar trends and SST estimates. However, the extrapolation of calibration on a sediment record from a greenhouse world ocean predicts more reliable absolute SST estimates and relative SST changes in agreement with estimates based on the δ¹⁸O of planktonic foraminifera. Based on the comparison of and derived SSTs using the core top data, we recommend applying above 15 °C and below 15 °C. In cases where paleorecords encompass temperatures both below and above 15 °C, we suggest to use .     

Elemental and isotopic abundances of lithium in chondrule constituents in the Allende meteorite

Elemental and isotopic abundances of lithium in chondrule constituents in the Allende CV3 meteorite were determined using secondary ion mass spectrometry. Olivines and mesostasis dominated by a feldspathic phase are depleted in Li ( and , respectively). In contrast, low-Ca pyroxenes and mesostasis dominated by a Na-rich phase are enriched in Li ( and , respectively) and the interchondrule matrix is generally enriched in Li ( on average). The Li isotopic abundance of olivine ranges from to 21. The spatial distributions of elemental and isotopic abundances of Li in olivines within individual chondrules exhibit no systematic pattern. This suggests that the distribution of Li in olivine was not disturbed during aqueous alteration or thermal metamorphism on the Allende meteorite parent body. Although mesostasis is the last crystallizing phase from a chondrule melt and is expected to be enriched in Li, in the Allende meteorite it is generally depleted in Li. We suggest that during aqueous alteration on the CV asteroid, Li in mesostasis was leached out by aqueous fluids. The Li-enriched Na-rich mesostasis was probably produced later by infiltration of Na-rich fluids. It seems likely that aqueous fluids sequestered alkali elements from the Allende-chondrite region in the CV parent asteroid, although significant amounts of Li are preserved in ferrous olivine in the interchondrule matrix.     

Organic geochemistry of Cuban oils—I. The northern geological province

The isoprenoid, hopanoid and steroid compositions of 15 oils from the most productive oil fields in Cuba were studied to determine source-rock depositional environments and organic matter sources. The oils, which are from the northern geological province of Cuba and can be defined by the position with respect to the overthrust belt, can be grouped into two families: those from the Remedios (1) and those from the Placetas (2) tectonostratigraphic units.Remedios oils contain 17α(H)-diahopane, high relative amounts of 18α(H)-22,29,30-trisnorneohopane and diasteranes, which is indicative of generation from clay-rich source rocks. The crude oils of the Placetas Unit exhibit a sterane/hopane composition consistent with a carbonate origin. Nevertheless, these polycyclic hydrocarbons exhibit significant changes in composition, indicating that several organic matter sources, e.g. a carbonate/evaporitic origin of Varadero and Varadero Sur oils, have contributed to the oils from this Unit. The Remedios oils are more mature [evaluated from the and sterane indices] than the Placetas oils.A wide range of biodegradation levels are encountered in these oils (from 0 to 7–8 using the scale derived by Volkman et al. (1983). The high relative abundance of 25-norhopanes is a distinctive feature of Remedios oils. The presence of these compounds in lightly biodegraded or nondegraded oils corresponds to a mixing of paleobiodegraded oil with more recently sourced nondegraded oil in the reservoir. The most biodegraded oil, Cantel, exhibits and sterane ratios, and hopane ratios that have been altered by microbial attack.     

The solubility of platinum and gold in NaCl brines at 1.5 kbar, 600 to 800°C: A laser ablation ICP-MS pilot study of synthetic fluid inclusions

The concentration and distribution of Pt and Au in a fluid-melt system has been investigated by reacting the metals with S-free, single-phase aqueous brines (20, 50, 70 wt% eq. NaCl) ± peraluminous melt at a confining pressure of 1.5 kbar and temperatures of 600 to 800 °C, trapping the fluid in synthetic fluid inclusions (quartz-hosted) and vesicles (silicate melt-hosted), and quantifying the metal content of the trapped fluid and glass by laser ablation ICP-MS. HCl concentration was buffered using the assemblage albite-andalusite-quartz and fO₂ was buffered using the assemblage Ni-NiO. Over the range of experimental conditions, measured concentrations of Pt and Au in the brines (, ) are on on the order of 1-10³ ppm. Concentrations of Pt and Au in the melt (, ) are ∼35-100 ppb and ∼400-1200 ppb, respectively. Nernst partition coefficients (, ) are on the order of 10²-10³ and vary as a function of (non-Henry’s Law behavior). Trapped fluids show a significant range of metal concentrations within populations of inclusions from single experiments (∼ 1 log unit variability for Au; ∼2-3 log unit variability for Pt). Variability in metal concentration within single inclusion groups is attributed to premature brine entrapment (prior to metal-fluid-melt equilibrium being reached); this allows us to make only minimum estimates of metal solubility using metal concentrations from primary inclusions. The data show two trends: (i) maximum and average values of and in inclusions decrease ∼2 orders of magnitude as fluid salinity () increases from ∼4 to 40 molal (20 to 70 wt % eq. NaCl) at a constant temperature; (ii) maximum and average values of increase approximately 1 order of magnitude for every 100°C increase temperature at a fixed . The observed behavior may be described by the general expression:

     

Sphalerite oxidation pathways detected by oxygen and sulfur isotope studies

Sphalerite oxidation is a common process under acid-mine drainage (AMD) conditions and results in the release of , Zn and potentially toxic trace metals, which can pollute rivers and oceans. However, there are only a few studies on the mechanisms of aerobic sphalerite oxidation. Oxygen and S isotope investigations of the produced may contribute to the understanding of sphalerite oxidation mechanisms so helping to interpret field data from AMD sites. Therefore, batch oxidation experiments with an Fe-rich sphalerite were performed under aerobic abiotic conditions at different initial pH values (2 and 6) for different lengths of time (2–100 days). The O and S isotope composition of the produced indicated changing oxidation pathways during the experiments. During the first 20 days of the experiments at both initial pH values, molecular O₂ was the exclusive O source of . Furthermore, the lack of S isotope enrichment processes between and sphalerite indicated that O₂ was the electron acceptor from sphalerite S. As the oxidation proceeded, a sufficient amount of released Fe(II) was oxidized to Fe(III) by O₂. Therefore, electrons could be transferred from sphalerite S sites to adsorbed hydrous Fe(III) and O from the hydration sphere of Fe was incorporated into the produced as indicated by decreasing δ¹⁸O_SO4 values which became more similar to the δ¹⁸O_H2O values. The enrichment of ³²S in relative to the sphalerite may also result from sphalerite oxidation by Fe(III).The incorporation of O₂ into during the oxidation of sphalerite was associated with an O isotope enrichment factor ε_SO4–O2 of ca. −22‰. The O isotope enrichment factor ε_SO4–H2O was determined to be ?4.1‰. A comparison with O and S studies of other sulfides suggests that there is no general oxidation mechanism for acid-soluble sulfides.     

Molybdenum, tungsten and manganese partitioning in the system pyrrhotite-Fe-S-O melt-rhyolite melt: Impact of sulfide segregation on arc magma evolution

Experiments were performed to determine the partitioning of molybdenum, tungsten and manganese among a rhyolitic melt (melt), pyrrhotite (po), and an immiscible Fe-S-O melt (Fe-S-O). Sulfide phases such as these may be isolated from a silicate melt along with other crystallizing phases during the evolution of arc magma, and partition coefficients are required to model the effect of this process on molybdenum and tungsten budgets.We developed an experimental design to take advantage of properties of the phases under study. Careful control of temperature allowed pyrrhotite and magnetite to be stable along with an Fe-S-O melt, and this phase assemblage allowed the composition of run-product pyrrhotite to be used to calculate both fS₂ and fO₂ for the experiments. At run temperature, (1042 ± 2 °C), a rhyolitic melt can be formed at low pressure, under nominally dry conditions, which removed the need for confining pressure as well as externally imposed fugacities. The silica-saturated melt allowed the charges to be contained in sealed evacuated silica tubes without danger of reaction, and with closed system behavior for molybdenum and tungsten.Experiments were run for durations up to 2000 min. Molybdenite (mb) and wolframite (wo) were added to the experiments as sources for molybdenum and tungsten, respectively. Manganese was added to the system as both a component of the starting rhyolitic pumice, and of Mn-bearing wolframite. Oxygen fugacity in these experiments was fixed at the Ni-NiO oxygen fugacity buffer. Sulfur fugacity was 10⁻¹ bar. Run products were analyzed by EPMA and LA-ICP-MS. Analysis of the run products yielded ( standard deviation of the mean): , , , and . The partition coefficients for manganese in this system are and .Simple Rayleigh fractionation modeling suggests that oxidized felsic melts produced through fractional crystallization may have lost as much as 14% of their initial molybdenum, but only 2% of their initial tungsten, through the removal of an Fe-S-O melt along with crystalline phases. Modeling consistent with conditions of oxygen and sulfur fugacity influenced by assimilation of sulfide (with low concentrations of molybdenum and tungsten) from, for example, sedimentary rock, results in evolved magmas significantly depleted in molybdenum, but only moderately depleted in tungsten. The molybdenum:tungsten ratio can vary by two orders of magnitude. These systematics may help to explain some of the variability in metal ratios of intrusion-related hydrothermal ore deposits.     

First-principles estimates of equilibrium magnesium isotope fractionation in silicate, oxide, carbonate and hexaaquamagnesium(2+) crystals

Equilibrium mass-dependent magnesium isotope fractionation factors are estimated for a range of crystalline compounds including oxides, silicates, carbonates, and salts containing the solvation complex. Fractionation factors for the gas-phase species Mg and MgO are also presented. Fractionation factors are calculated with density functional perturbation theory (DFPT), using norm-conserving pseudopotentials. The results suggest that there will be substantial inter-mineral fractionation, particularly between tetrahedrally coordinated Mg²⁺ in spinel (MgAl₂O₄) and the more common octahedrally coordinated Mg²⁺-sites in silicate and carbonate minerals. Isotope fractionations calculated for Mg²⁺ in hexaaquamagnesium(2+) salts are in good agreement with previous fractionation models of based on large molecular clusters (Black et al., 2007), but show possibly more significant disagreement with a more recent study (Rustad et al., 2010). These models further suggest that solvated , in the form of , will have higher ²⁶Mg/²⁴Mg than coexisting magnesite and dolomite. Calculated fractionations are consistent with Mg-isotope fractionations observed in peridotite mineral separates and inorganic carbonate precipitates. Predicted large, temperature-sensitive spinel-silicate fractionations, in particular, may find use in determining equilibration temperatures of peridotites and other high-temperature rock types.     

Sorption of yttrium and rare earth elements by amorphous ferric hydroxide: Influence of solution complexation with carbonate

The influence of solution complexation on the sorption of yttrium and the rare earth elements (YREEs) by amorphous ferric hydroxide was investigated at 25 °C over a range of pH (4.0-7.1) and carbonate concentrations . Distribution coefficients, defined as , where [MS_i]_T is the total concentration of sorbed YREE, M_T is the total YREE concentration in solution, and [S_i] is the concentration of amorphous ferric hydroxide, initially increased in magnitude with increasing carbonate concentration, and then decreased. The initial increase of is due to sorption of YREE carbonate complexes , in addition to sorption of free YREE ions (M³⁺). The subsequent decrease of , which is more extensive for the heavy REEs, is due to the increasing intensity of YREE solution complexation by carbonate ions. The competition for YREEs between solution complexation and surface complexation was modeled via the equation:

     

Comparative carbonate and hydroxide complexation of cations in seawater

Comparative concentrations of carbonate and hydroxide complexes in natural solutions can be expressed in terms of reactions with bicarbonate that have no explicit pH dependence (). Stability constants for this reaction with n = 1 were determined using conventional formation constant data expressed in terms of hydroxide and carbonate. Available data indicate that stability constants appropriate to seawater at 25 °C expressed in the form are on the order of 10^4.2 for a wide range of cations (M^z+) with z = +1, +2 and +3. Φ₁ is sufficiently large that species appear to substantially dominate MOH^z−1 species in seawater. Evaluations of comparative stepwise carbonate and hydroxide stability constant behavior leading to the formation of n = 2 and n = 3 complexes suggest that carbonate complexes generally dominate hydroxide complexes in seawater, even for cations whose inorganic speciation schemes in seawater are currently presumed to be strongly dominated by hydrolyzed forms (). Calculated stability constants, and , indicate that the importance of carbonate complexation is sufficiently large that carbonate and hydroxide complexes would be generally comparable even if calculated Φ₂ and Φ₃ values are overestimated by two or more orders of magnitude. Inclusion of mixed ligand species in carbonate-hydroxide speciation models allows cation complexation intensities (M_T/[M^z+]) to be expressed in the following form:

     

Flow and nutrient dynamics in a subterranean estuary (Waquoit Bay, MA, USA): Field data and reactive transport modeling

A two-dimensional (2D) reactive transport model is used to investigate the controls on nutrient (, , PO₄) dynamics in a coastal aquifer. The model couples density-dependent flow to a reaction network which includes oxic degradation of organic matter, denitrification, iron oxide reduction, nitrification, Fe²⁺ oxidation and sorption of PO₄ onto iron oxides. Porewater measurements from a well transect at Waquoit Bay, MA, USA indicate the presence of a reducing plume with high Fe²⁺, , DOC (dissolved organic carbon) and PO₄ concentrations overlying a more oxidizing -rich plume. These two plumes travel nearly conservatively until they start to overlap in the intertidal coastal sediments prior to discharge into the bay. In this zone, the aeration of the surface beach sediments drives nitrification and allows the precipitation of iron oxide, which leads to the removal of PO₄ through sorption. Model simulations suggest that removal of through denitrification is inhibited by the limited overlap between the two freshwater plumes, as well as by the refractory nature of terrestrial DOC. Submarine groundwater discharge is a significant source of to the bay.     

Sorption of phosphate onto calcite; results from batch experiments and surface complexation modeling

The adsorption of phosphate onto calcite was studied in a series of batch experiments. To avoid the precipitation of phosphate-containing minerals the experiments were conducted using a short reaction time (3 h) and low concentrations of phosphate (?50 μM). Sorption of phosphate on calcite was studied in 11 different calcite-equilibrated solutions that varied in pH, P_CO2, ionic strength and activity of Ca²⁺, and . Our results show strong sorption of phosphate onto calcite. The kinetics of phosphate sorption onto calcite are fast; adsorption is complete within 2-3 h while desorption is complete in less than 0.5 h. The reversibility of the sorption process indicates that phosphate is not incorporated into the calcite crystal lattice under our experimental conditions. Precipitation of phosphate-containing phases does not seem to take place in systems with ?50 μM total phosphate, in spite of a high degree of super-saturation with respect to hydroxyapatite (SI_HAP ? 7.83). The amount of phosphate adsorbed varied with the solution composition, in particular, adsorption increases as the activity decreases (at constant pH) and as pH increases (at constant activity). The primary effect of ionic strength on phosphate sorption onto calcite is its influence on the activity of the different aqueous phosphate species. The experimental results were modeled satisfactorily using the constant capacitance model with >CaPO₄Ca⁰ and either >CaHPO₄Ca⁺ or > as the adsorbed surface species. Generally the model captures the variation in phosphate adsorption onto calcite as a function of solution composition, though it was necessary to include two types of sorption sites (strong and weak) in the model to reproduce the convex shape of the sorption isotherms.