Toward an Understanding of Biogenic-silica Dissolution in Seawater – An Initial Rate Approach Applied between 40 and 90 °C

Summary The kinetics of phytoplankton frustule dissolution has generally been studied as the appearance of silicic acid in a batch reactor. Unfortunately, this approach, though often illuminating, has not so far been successful because of the difficulty of parameterising the full reaction curve. This current study shows how the initial rate approach to chemical kinetics offers a way around this bottleneck, thereby allowing much chemical kinetics information about frustule dissolution to be collected. The technique is shown to be flexible and suited to short reaction times which facilitate detailed quantitative kinetics investigation, indeed, as would be expected in a solution phase, kinetics study. The technique is exemplified by a dissolution study of uncleaned frustules of Cyclotella crypticaat 40 °C and above. The frustules were found to yield the same dissolution rate after 5 weeks dark storage, at 4 °C. Meanwhile, log dissolution rate was found to vary linearly with pH, with gradient 0.38 ± 0.01 (r ²=0.990). Linearity was upheld even at pHs as high as 14. Finally, a robust Arrhenius plot was established between 40 and 90 °C yielding an activation energy for dissolution of 84 ± 3 kJ mol ⁻¹. Follow through with the Eyring equation yielded an activation enthalpy, ΔH ^‡, and an activation entropy, ΔS ^‡, of 81 and 85 J mol ⁻¹ K ⁻¹, respectively. The discussion brings salient aspects of existing knowledge about diatom frustule dissolution kinetics into the wider context of silicate mineral dissolution.     

Calorimetric study of perovskite solid solutions in the CaSiO3–CaGeO3 system

 Enthalpies of drop solution (ΔH _drop-sol) of CaGeO₃, Ca(Si_0.1Ge_0.9)O₃, Ca(Si_0.2Ge_0.8)O₃, Ca(Si_0.3Ge_0.7)O₃ perovskite solid solutions and CaSiO₃ wollastonite were measured by high-temperature calorimetry using molten 2PbO · B₂O₃ solvent at 974 K. The obtained values were extrapolated linearly to the CaSiO₃ end member to give ΔH _drop-sol of CaSiO₃ perovskite of 0.2 ± 4.4 kJ mol⁻¹. The difference in ΔH _drop-sol between CaSiO₃, wollastonite, and perovskite gives a transformation enthalpy (wo → pv) of 104.4 ± 4.4 kJ mol⁻¹. The formation enthalpy of CaSiO₃ perovskite was determined as 14.8 ± 4.4 kJ mol⁻¹ from lime + quartz or −22.2 ± 4.5 kJ mol⁻¹ from lime + stishovite. A comparison of lattice energies among A²⁺B⁴⁺O₃ perovskites suggests that amorphization during decompression may be due to the destabilizing effect on CaSiO₃ perovskite from a large nonelectrostatic energy (repulsion energy) at atmospheric pressure. By using the formation enthalpy for CaSiO₃ perovskite, phase boundaries between β-Ca₂SiO₄ + CaSi₂O₅ and CaSiO₃ perovskite were calculated thermodynamically utilizing two different reference points [where ΔG(P,T )=0] as the measured phase boundary. The calculations suggest that the phase equilibrium boundary occurs between 11.5 and 12.5 GPa around 1500 K. Its slope is still not well constrained. Received: 20 September 2000 / Accepted: 17 January 2001     

Stability of clinohumite in the system MgO-SiO2-H2O

The low-pressure stability of clinohumite has been investigated in phase-equilibrium experiments on the reaction forsterite + brucite = clinohumite. The reaction was bracketed between 2.45 and 2.84 GPa at 650 °C, extending to between 1.37 and 1.57 GPa at 850 °C. At temperatures above the reaction brucite = periclase + vapour, the reaction clinohumite = forsterite + vapour was bracketed between 1.27 and 1.52 GPa at 900 °C, rising to between 1.90 and 2.00 GPa at 1000 °C. The position of the reaction forsterite + brucite = clinohumite is ∼0.5 GPa below the position determined in previous work, the difference arising either from pressure uncertainties in both studies, from enhanced reaction to clinohumite in this study due to the presence of excess brucite in the starting material, or from different concentrations of defects in the two samples. The brackets on the reaction were combined with other measured and estimated thermodynamic data for clinohumite to determine its enthalpy of formation and entropy, in a revised version of the thermodynamic dataset of Holland and Powell (1998). The values obtained were ΔH _f =−9607.29±3.05 kJ mol⁻¹, S=445 J mol⁻¹ K⁻¹. These data were used to calculate positions of other reactions involving clinohumite. The calculations suggest a larger stability field for clinohumite than implied by the results of previous experimental studies, indicating a need for more high-pressure phase-equilibrium studies to provide better thermodynamic data. Received: 30 April 1999 / Accepted: 8 November 1999     

Thermodynamic data of the high-pressure phase Mg5Al5Si6O21(OH)7 (Mg-sursassite)

 Calorimetric and P–V–T data for the high-pressure phase Mg₅Al₅Si₆O₂₁(OH)₇ (Mg-sursassite) have been obtained. The enthalpy of drop solution of three different samples was measured by high-temperature oxide melt calorimetry in two laboratories (UC Davis, California, and Ruhr University Bochum, Germany) using lead borate (2PbO·B₂O₃) at T=700 ^∘C as solvent. The resulting values were used to calculate the enthalpy of formation from different thermodynamic datasets; they range from −221.1 to −259.4 kJ mol⁻¹ (formation from the oxides) respectively −13892.2 to −13927.9 kJ mol⁻¹ (formation from the elements). The heat capacity of Mg₅Al₅Si₆O₂₁(OH)₇ has been measured from T=50 ^∘C to T=500 ^∘C by differential scanning calorimetry in step-scanning mode. A Berman and Brown (1985)-type four-term equation represents the heat capacity over the entire temperature range to within the experimental uncertainty: C _P (Mg-sursassite) =(1571.104 −10560.89×T ^−0.5−26217890.0 ×T ⁻²+1798861000.0×T ⁻³) J K⁻¹ mol⁻¹ (T in K). The P V T behaviour of Mg-sursassite has been determined under high pressures and high temperatures up to 8 GPa and 800 ^∘C using a MAX 80 cubic anvil high-pressure apparatus. The samples were mixed with Vaseline to ensure hydrostatic pressure-transmitting conditions, NaCl served as an internal standard for pressure calibration. By fitting a Birch-Murnaghan EOS to the data, the bulk modulus was determined as 116.0±1.3 GPa, (K ^′=4), V _T,0 =446.49 ³ exp[∫(0.33±0.05) × 10⁻⁴ + (0.65±0.85)×10⁻⁸ T dT], (K _T/T) _P  = −0.011± 0.004 GPa K⁻¹. The thermodynamic data obtained for Mg-sursassite are consistent with phase equilibrium data reported recently (Fockenberg 1998); the best agreement was obtained with Δ_f H ⁰ ₂₉₈ (Mg-sursassite) = −13901.33 kJ mol⁻¹, and S ⁰ ₂₉₈ (Mg-sursassite) = 614.61 J K⁻¹ mol⁻¹. Received: 21 September 2000 / Accepted: 26 February 2001     

Pb diffusion in rutile

Diffusion of Pb was measured in natural and synthetic rutile under dry, 1 atmosphere conditions, using mixtures of Pb titanate or Pb sulfide and TiO₂ as the sources of diffusant. Pb depth profiles were then measured with Rutherford Backscattering Spectrometry (RBS). Over the temperature range 700–1100 °C, the following Arrhenius relation was obtained for the synthetic rutile: D=3.9 × 10⁻¹⁰exp(−250 ± 12 kJ mol⁻¹/RT) m²s⁻¹. Results for diffusion in natural and synthetic rutile were quite similar, despite significant differences in trace element compositions. Mean closure temperatures calculated from the diffusion parameters are around 600 °C for rutile grains of ∼100 μm size. This is about 100 °C higher than rutile closure temperature determinations from past field-based studies, suggesting that rutile is more resistant to Pb loss through volume diffusion than previously thought. Received: 28 June 1999 / Accepted: 29 December 1999     

High-temperature heat capacity and thermodynamic properties for end-member titanite (CaTiSiO5)

 The heat capacity of end-member titanite and (CaTiSiO₅) glass has been measured in the range 328–938 K using differential scanning calorimetry. The data show a weak λ-shaped anomaly at 483 ± 5 K, presumably associated with the well-known low-pressure P2₁/a ⇆ A2/a transition, in good agreement with previous studies. A value of 0.196 ± 0.007 kJ mol⁻¹ for the enthalpy of the P2₁/a ⇆ A2/a transition was determined by integration of the area under the curve for a temperature interval of 438–528 K, bracketing the anomaly. The heat capacity data for end-member titanite and (CaTiSiO₅) glass can be reproduced within <1% using the derived empirical equations (temperature in K, pressure in bars):

Enthalpy of formation of CaSi2O5, a quenched high-pressure phase with pentacoordinate silicon

 The monoclinic titanite-like high-pressure form of calcium disilicate has been synthesized and quenched to ambient conditions to form the triclinic low-pressure phase containing silicon in four-, five- and sixfold coordination. The enthalpy of formation of the quench product has been measured by high-temperature oxide melt calorimetry. The value obtained from samples from a series of several synthesis experiments is ΔH _f = (−26.32 ± 4.27) kJ mol⁻¹ for the formation from the component oxides, or ΔH _f  = (−2482.81 ± 4.59) kJ mol⁻¹ for the formation from the elements. The result is identical within experimental error to available estimates, although the previously predicted energy difference between the monoclinic and triclinic phases could not be verified. Received: 16 February 2000 / Accepted: 14 July 2000     

The normal grain growth behaviour of nominally pure calcitic aggregates

The normal grain growth behaviour of four different, but all nominally pure, calcite powders (99%+ analytic grade calcite, 99.7% chalk, 99.97% crushed Iceland Spar, 99.95%+ chelometric grade calcite) has been investigated as a function of temperature (550, 600, 650, 700 °C) and confining pressure (100, 190 MPa) under both “dry” and hydrostatic (P _fluid = P _total) conditions. The initial particle size of both the analytic grade and chelometric grade calcite was about 5 μm, and that of the chalk was about 3 μm, while the experiments on the Iceland Spar were conducted on powders of three different initial particle sizes (3.4, 7.5, 38.5 μm). On each material, at each pressure/temperature condition 6 to 15 experiments, equally spaced in log time from 15 minutes to 50 days, were conducted. Under dry conditions all four materials recrystallized to aggregates which contained less than 2% porosity and which had a grain size of between 4 and 20 μm (depending on the initial particle size). Subsequently the aggregates coarsened by normal grain growth, with the kinetics of the growth process being controlled by the rate at which the grain boundaries could drag the residual pores with them as they migrated. Under nominally identical conditions both the mechanism and rates of pore drag differed greatly for the different materials, implying that this process is highly sensitive to trace solute impurity concentrations. This sensitivity renders the task of providing a systematic account of dry calcite grain growth kinetics highly problematic. Under hydrostatic conditions all the powders followed the same normal grain growth kinetics in which the growth process was rate-controlled by diffusion through the pore fluid on the grain boundaries. An activation enthalpy of 162.6 kJ mol⁻¹ and an activation volume of 34.35 cm³ mol⁻¹ was obtained for this process. Received: 23 May 1996 / Accepted: 8 July 1997     

Oxygen diffusion in anhydrous sanidine feldspar

Oxygen exchange experiments have been performed between single crystals of sanidine feldspar and oxygen gas enriched in ¹⁸O, at temperatures in the range 869–1053 °C, total pressure 1 atmosphere, for times up to 28 days. Oxygen isotope diffusion profiles in a direction perpendicular to (001) were determined with an ion microprobe. The experimental data obey a single Arrhenius relationship of the form D = 8.4 × 10⁻¹¹ exp. (−245 ± 15 kJ mol⁻¹/RT) m²s⁻¹. The results indicate that oxygen diffusion in anhydrous sanidine feldspar is marginally slower than oxygen diffusion in anhydrous anorthite. Comparison with published atomistic simulation studies suggests that oxygen transport in feldspar is by an interstitial mechanism. Received: 17 October 1997 / Accepted: 6 July 1998     

Fe–Mg partitioning between ringwoodite and magnesiowüstite and the effect of pressure, temperature and oxygen fugacity

 The partitioning of Mg and Fe between magnesiowüstite and ringwoodite solid solutions has been measured between 15 and 23 GPa and 1200–1600 ^∘C using both Fe and Re capsule materials to vary the oxidation conditions. The partitioning results show a clear dependence on the capsule material used due to the variation in Fe³⁺ concentrations as a consequence of the different oxidation environments. Using results from experiments performed in Fe capsules, where metallic Fe was also added to the starting materials, the difference in the interaction parameters for the two solid solutions (W _FeMg ^mw−W _FeMg ^ring) is calculated to be 8.5±1 kJ mol⁻¹. Similar experiments performed in Re metal capsules result in a value for W _FeMg ^mw−W _FeMg ^ring that is apparently 4 kJ higher, if all Fe is assumed to be FeO. Electron energy-loss near-edge structure (ELNES) spectroscopic analyses, however, show Fe³⁺ concentrations to be approximately three times higher in magnesiowüstite produced in Re capsules than in Fe capsules and that Fe³⁺ partitions preferentially into magnesiowüstite, with K _{D Fe3+} ^ring/mw estimated between 0.1 and 0.6. Using an existing activity composition model for magnesiowüstite, a least–squares fit to the partitioning data collected in Fe capsules results in a value for the ringwoodite interaction parameter (W _FeMg ^ring) of 3.5±1 kJ mol⁻¹. The equivalent regular interaction parameter for magnesiowüstite (W _FeMg ^mw) is 12.1±1.8 kJ mol. These determinations take into account the Fe³⁺ concentrations that occur in both phases in the presence of metallic Fe. The free energy change in J mol⁻¹ for the Fe exchange reaction can be described, over the range of experimental conditions, by 912 + 4.15 (T−298)+18.9P with T in K, P in kbar. The estimated volume change for this reaction is smaller than that predicted using current compilations of equation of state data and is much closer to the volume change at ambient conditions. These results are therefore a useful test of high pressure and temperature equation of state data. Using thermodynamic data consistent with this study the reaction of ringwoodite to form magnesiowüstite and stishovite is calculated from the data collected using Fe capsules. Comparison of these results with previous studies shows that the presence of Fe³⁺ in phases produced in multianvil experiments using Re capsules can have a marked effect on apparent phase relations and determined thermodynamic properties. Received: 13 September 2000 / Accepted: 25 March 2001     

Iron Diffusion in Single-Crystal Diopside

 Iron tracer diffusion experiments in diopside have been performed using natural and synthetic single crystals of diopside, and stable iron tracers enriched in ⁵⁴Fe, at temperatures in the range 950–1100 °C, total pressure 1 atm, for times up to 29 days. Iron isotope diffusion profiles were determined with an ion microprobe. For experiments performed at log pO₂ = −13, in directions parallel to the c axis and the b axis of two natural, low iron (Fe ∼ 1.8 at %) diopsides, the data obey a single Arrhenius relationship of the form D = 6.22_−5.9 ^+49.6×10⁻¹⁵ exp(−161.5 ± 35.0 kJ mol⁻¹/RT) m² s⁻¹. A single datum for iron diffusion in iron-free, single-crystal diopside at 1050 °C, is approximately 1 order of magnitude slower than in the natural crystals. The pO₂ dependence of iron diffusion in natural crystals at 1050 °C (power exponent = 0.229 ± 0.036) indicates a vacancy mechanism; this is consistent with the results of unpublished atomistic simulation studies. There is no evidence of anisotropy for iron diffusion in diopside. Received: 16 March 1999 / Accepted: 10 April 2000     

Thermal decomposition of Ferrian chamosite: an infrared emission spectroscopic study

The thermal behaviour of ripidolite, an iron-rich chlorite, has been studied in situ by infrared emission spectroscopy up to 800 °C. The more di,trioctahedral nature due to significant amounts of Fe³⁺ is reflected, in addition to the two bands around 3420 and 3560 cm⁻¹, by an extra band around 3345 cm⁻¹. This extra band is absent in pure dioctahedral chlorites without Fe³⁺. These bands have been assigned to (AlAl)O-OH, (SiAl)O-OH and (SiSi)O-OH stretching modes with increasing frequencies. The bands disappear upon dehydroxylation around 650 °C. A similar behaviour is observed for the corresponding libration modes around 716, 759 and 802 cm⁻¹. The stretching and bending modes of the inner-OH of the octahedral sheet in the 2:1 clay-like layer are observed around 3645, 943 and 904 cm⁻¹. Although the bands decrease in intensity, they remain present up to 800 °C as dehydroxylation of the octahedral sheet is not yet complete at this temperature. The presence of two bending modes is explained as being due to a differentiation between Mg-OH and Fe-OH modes. At 650 °C a new sharp band is observed around 502 cm⁻¹ assigned to a (Fe,Mg)-O-Al bending mode caused by the formation of a spinel-like interlayer phase after dehydroxylation. Received: 4 June 1999 / Accepted: 6 August 1999     

Structural thermal behavior of paragonite and its dehydroxylate: a high-temperature single-crystal study

 The lattice constants of paragonite-2M₁, NaAl₂(AlSi₃)O₁₀(OH)₂, were determined to 800 °C by the single-crystal diffraction method. Mean thermal expansion coefficients, in the range 25–600 °C, were: α_a = 1.51(8) × 10⁻⁵, α_b = 1.94(6) × 10⁻⁵, α_c = 2.15(7) ×  10⁻⁵ °C⁻¹, and α_V = 5.9(2) × 10⁻⁵ °C⁻¹. At T higher than 600 °C, cell parameters showed a change in expansion rate due to a dehydroxylation process. The structural refinements of natural paragonite, carried out at 25, 210, 450 and 600 °C, before dehydroxylation, showed that the larger thermal expansion along the c parameter was mainly due to interlayer thickness dilatation. In the 25–600 °C range, Si,Al tetrahedra remained quite unchanged, whereas the other polyhedra expanded linearly with expansion rate proportional to their volume. The polyhedron around the interlayer cation Na became more regular with temperature. Tetrahedral rotation angle α changed from 16.2 to 12.9°. The structure of the new phase, nominally NaAl₂ (AlSi₃)O₁₁, obtained as a consequence of dehydroxylation, had a cell volume 4.2% larger than that of paragonite. It was refined at room temperature and its expansion coefficients determined in the range 25–800 °C. The most significant structural difference from paragonite was the presence of Al in fivefold coordination, according to a distorted trigonal bipyramid. Results confirm the structural effects of the dehydration mechanism of micas and dioctahedral 2:1 layer silicates. By combining thermal expansion and compressibility data, the following approximate equation of state in the PTV space was obtained for paragonite: V/V ₀ = 1 + 5.9(2) × 10⁻⁵ T(°C) − 0.00153(4) P(kbar). Received: 12 July 1999 / Revised, accepted: 7 December 1999     

Phase equilibrium and calorimetric study of the transition of MnTiO3 from the ilmenite to the lithium niobate structure and implications for the stability field of perovskite

The phase boundary between MnTiO₃ I (ilmenite structure) and MnTiO₃ II (lithium niobate structure) has been determined by analysis of quench products from reversal experiments in a cubic anvil apparatus at 1073–1673 K and 43–75 kbar using mixtures of MnTiO₃ I and II as starting materials. Tight brackets of the boundary give P(kbar)=121.2−0.045 T(K). Thermodynamic analysis of this boundary gives ΔH^o=5300±1000 J·mol⁻¹, ΔS^o = 1.98 ±1J·K⁻¹· mol⁻¹. The enthalpy of transformation obtained directly by transposed-temperature-drop calorimetry is 8359 ±2575 J·mol⁻¹. Possible topologies of the phase relations among the ilmenite, lithium niobate, and perovskite polymorphs are constrained using the above data and the observed (reversible with hysteresis) transformation of II to III at 298 K and 20–30 kbar (Ross et al. 1989). The observed II–III transition is likely to lie on a metastable extension of the II–III boundary into the ilmenite field. However the reversed I–II boundary, with its negative dP/ dT does represent stable equilibrium between ilmenite and lithium niobate, as opposed to the lithium niobate being a quench product of perovskite. We suggest a topology in which the perovskite occurs stably at low T and high P with a triple point (I, II, III) at or below 1073 K near 70 kbar. The I–II boundary would have a negative P-T slope while the II–III and I–III boundaries would be positive, implying that entropy decreases in the order lithium niobate, ilmenite, perovskite. The inferred positive slope of the ilmenite-perovskite transition in MnTiO₃ is different from the negative slopes in silicates and germanates. These thermochemical parameters are discussed in terms of crystal structure and lattice vibrations.     

Origin of the Zálesí U–Ni–Co–As–Ag/Bi deposit,Bohemian Massif,Czech Republic: fluid inclusion and stable isotope constraints

The Zálesí vein-type deposit is hosted by Early Paleozoic high-grade metamorphic rocks on the northern margin of the Bohemian Massif. The mineralization is composed of three main stages: uraninite, arsenide, and sulfide. The mineral assemblages formed at low temperatures (~80 to 130°C, locally even lower) and low pressures (<100 bars). The salinity of the aqueous hydrothermal fluids (0 to 27 wt.% salts) and their chemical composition vary significantly. Early fluids of the oldest uraninite stage contain a small admixture of a clathrate-forming gas, possibly CO₂. Salinity correlates with oxygen isotope signature of the fluid and suggests mixing of brines [δ ¹⁸O around +2‰ relative to standard mean ocean water (SMOW)] with meteoric waters (δ ¹⁸O around −4‰ SMOW). The fluid is characterized by highly variable halogen ratios (molar Br/Cl = 0.8 × 10⁻³ to 5.3 × 10⁻³; molar I/Cl = 5.7 × 10⁻⁶ to 891 × 10⁻⁶) indicating a dominantly external origin for the brines, i.e., from evaporated seawater, which mixed with iodine-enriched halite dissolution brine. The cationic composition of these fluids indicates extensive interaction of the initial brines with their country rocks, likely associated with leaching of sulfur, carbon, and metals. The brines possibly originated from Permian–Triassic evaporites in the neighboring Polish Basin, infiltrated into the basement during post-Variscan extension and were finally expelled along faults giving rise to the vein-type mineralization. Cenozoic reactivation by low-salinity, low-δ ¹⁸O (around −10‰ SMOW) fluids of mainly meteoric origin resulted in partial replacement of primary uraninite by coffinite-like mineral aggregates.     

Fluid chemistry and processes at the Porgera gold deposit, Papua New Guinea

The Porgera gold deposit in Papua New Guinea is a world-class example of an alkalic-type epithermal gold system (stage II), which overprints a precursor stage of magmatic-hydrothermal gold mineralization (stage I). Gas and ion chromatographic analyses of fluid inclusions contained in vein minerals from both mineralization stages have been carried out in order to constrain the compositions of the fluids involved in, and the processes attending, ore deposition. These data indicate the presence of three end-member liquids, the most dilute of which was present throughout the mineralization history and is interpreted to represent evolved groundwater of meteoric origin. Its composition is estimated to have been approximately 500 mM Na⁺, 10 mM K⁺, 5 mM Li⁺, 250 mM Cl⁻, 0.15 mM Br⁻, and 0.01 mM I⁻, plus significant concentrations of dissolved gases. More saline liquids were also present during the two main stages of ore formation, and although their compositions differ, both are interpreted to have been derived at least in part from magmatic fluids, and to have been the media by which gold was introduced into the system. Stage I minerals contain fluid inclusions which decrease in salinity towards this dilute end-member composition through the vein paragenesis, reflecting progressive dilution at depth of the magmatic fluid source by groundwaters. Ore deposition is thought to have been caused largely by simple cooling and/or wallrock reactions, although limited in situ fluid mixing may also have occurred. The most saline fluids, present in early quartz and pyrite, contain at least 810 mM Na⁺, 530 mM Ca²⁺, 130 mM K⁺, 12 mM Li⁺, 87 mM SO₄ ²⁻, 960 mM Cl⁻, 1.1 mM Br⁻, and 0.05 mM I⁻, plus significant but variable concentrations of dissolved gases. Fluid inclusions from stage II hydraulic breccia veins reveal the presence of two distinct liquids with contrasting salinities, which were present at different times during vein formation. A higher salinity liquid appears to have predominated during mineralization, whereas lower salinity groundwaters filled the structures during intervening periods. The ore-forming fluid may have been forcibly injected into the veins from depth during fracturing and depressurization events, displacing the resident groundwaters in the process. The original composition of this fluid is estimated to have been at least 1770 mM Na⁺, 59 mM K⁺, 180 mM Li⁺, 210 mM SO₄ ²⁻, 680 mM Cl⁻, 1.4 mM Br⁻, and 0.09 mM I⁻, plus 1.5 mol% CO₂, 0.19 mol% CH₄, and 0.04 mol% N₂. Gas chromatographic analyses of fluid inclusions from stage II samples show a decrease in total gas content between early unmineralized veins and post-mineralization vuggy quartz (suitable samples could not be obtained from the ore stage itself). Post-mineralization samples plot along an experimental gas-saturation curve in the CO₂-CH₄-H₂O-NaCl system, obtained at conditions similar to those attending stage II ore deposition at Porgera (200–300 bar, ˜165 °C). These results are interpreted to indicate a period of depressurization-induced phase separation during hydraulic fracturing, which resulted in rich ore deposition. Volatile gases such as CH₄ and N₂, in addition to CO₂ in solution, are shown to have a significant negative effect on total gas solubility. This effect may be of critical importance in lowering the temperature and increasing the depth (pressure) at which phase separation can occur in epithermal systems. Received: 28 November 1995 / Accepted: 17 July 1996     

The enthalpy of transformation of Ca(OH)2-I (portlandite) to Ca(OH)2-II (EuI2 structure) by low-temperature DSC

 Thermodynamic properties of high-pressure minerals that are not recoverable from synthesis experiments by conventional quenching methods (“unquenchable” phases) usually are calculated from equation of state data and phase diagram topologies. The present study shows that, with cryogenic methods of recovery and sample treatment, phases with a suitable decomposition rate can be made accessible to direct thermodynamic measurements. A set of samples of Ca(OH)₂-II has been synthesized in a multianvil device and subsequently recovered by cooling the high-pressure assembly with liquid nitrogen. Upon heating from liquid nitrogen to room temperature, the material transformed back to Ca(OH)₂-I. The heat effect of this backtransformation was measured by differential scanning calorimetry. A commercial differential scanning calorimeter (Netzsch DSC 404), modified to allow sample loading at liquid nitrogen temperature was used to heat the material from −150 to +200 °C at rates varying between 5 and 15 °C min⁻¹. The transformation started around −50 °C very gradually, and peaked at about 0 °C. To obtain a baseline correction, each sample was scanned under exactly the same conditions after the backtransformation was complete. Because of the relative sluggishness, onset and offset temperatures were not well defined as compared to fast (e.g., melting) reactions. To aid in integration, the resulting signals were successfully fitted using a generic asymmetric peak model. The enthalpy of backtransformation was determined to be ΔH =−10.37 ± 0.50 kJ mol⁻¹. From previous in situ X-ray diffraction experiments, the location of the direct transformation in P-T space has been constrained to 5.7 ± 0.4 GPa at 500 °C (Kunz et al. 1996). With the reaction volume known from the same study, and assuming that ΔC _p of the transformation remains negligible between the conditions of our measurements and 500 °C, our result gives an estimate of the entropy of transition and the P-T slope of the reaction curve. To a first approximation, the values ΔS = −16.00 ± 0.65 J(mol · K)⁻¹ and dP/dT = 0.0040 ± 0.0002 GPa/K have been determined. These results need to be refined by equation of state data for Ca(OH)₂-II. Received: 30 December 1999 / Accepted: 10 April 2000     

Diffusion of tetravalent cations in zircon

Diffusion rates for the three tetravalent cations U, Th and Hf have been measured in synthetic zircon. Diffusant sources included oxide powders and ground pre-synthesized silicates. Rutherford backscattering spectrometry (RBS) was used to measure depth profiles. Over the temperature range 1400–1650 °C, the following Arrhenius relations were obtained (diffusion coefficients in m²sec⁻¹): log D _Th = (1.936 ± 0.9820) + (− 792 ± 34 kJ mol⁻¹ /2.303 RT) log D _U = (0.212 ± 2.440) + (− 726 ± 83 kJ mol⁻¹ /2.303 RT) log D _Hf = (3.206 ± 1.592) + (− 812 ± 54 kJ mol⁻¹ /2.303 RT) The data show a systematic increase in diffusivity with decreasing ionic radius (i.e., faster diffusion rates for Hf than for U or Th), a trend also observed in our earlier study of rare earth diffusion in zircon. Diffusive fractionation may be a factor in the Lu-Hf system given the much slower diffusion rates of tetravalent cations when compared with the trivalent rare earths. The very slow diffusion rates measured for these tetravalent cations suggest that they are essentially immobile under most geologic conditions, permitting the preservation of fine-scale chemical zoning and isotopic signatures of inherited cores. Received: 12 July 1996 / Accepted: 2 December 1996     

The origin of the Tongkeng-Changpo tin deposit,Dachang metal district,Guangxi, China: clues from fluid inclusions and He isotope systematics

Tongkeng-Changpo is the largest tin deposit within the giant Dachang polymetallic tin ore field in Guangxi, southern China, which is part of a large skarn system associated with Cretaceous granitoids. The Tongkeng-Changpo mineralization consists of veins and stockworks in the upper levels and replacement stratiform orebodies (mantos) at lower levels. Based on textural relationships, three major mineralizing stages can be recognized: stage I with cassiterite, sulphides, stannite, tourmaline, and quartz; stage II with cassiterite, sulphides, sulphosalts, quartz, and calcite; and stage III with calcite as the main phase. The study of fluid inclusions has shown that there are two main fluid types: CO₂ and NaCl-H₂O. Homogenization temperatures are 270 to 365°C, 210 to 240°C, and 140 to 190°C for stages I, II, and III, respectively. Salinities range from 1 to 7 wt.% NaCl equiv. in the early ore stage and 3 to 10 wt.% NaCl equiv. in the late stages. Laser Raman Spectroscopy indicates that the inclusion fluids in stages I and II were of carbono-aqueous composition, with minor amounts of CH₄ and H₂S, whereas those in stage III were aqueous. Helium isotopic analyses of inclusion fluids indicate that the ³He/⁴He ratios in the ore veins are in between 1.2 to 2.9 Ra (Ra = 1.4 × 10⁻⁶, modern atmospheric ratio), and range from 1.6 to 2.5 Ra in the stratiform orebodies. This range of ³He/⁴He ratios is significantly higher than that of crustal fluids (0.01–0.05 Ra). The similar characteristics of fluid inclusions and their He isotopic composition, as well as age constraints, indicate that the ore veins and stratiform orebodies of the Tongkeng-Changpo deposit formed from the same hydrothermal system, likely related to granite intrusions of the Mesozoic Yanshanian tectono-thermal event. In addition, the high R/Ra ratios indicate a mantle contribution in the ore fluids.     

OH absorption coefficients of rutile and cassiterite deduced from nuclear reaction analysis and FTIR spectroscopy

Summary The OH content of four rutile and two cassiterite single-crystals was studied by nuclear reaction analysis (NRA) and by polarised FTIR microspectroscopy. The OH absorption bands of both minerals are centered around 3300 cm⁻¹ with different absorption features. The analytical H₂O content determined by NRA ranges from 70 to 820 wt.ppm. The integrated molar absorption coefficients deduced from the total integrated OH absorbances are equal to 38000 lċmol⁻¹ _H2Oċcm⁻² for rutile and 65000 lċmol⁻¹ _H2Oċcm⁻² for cassiterite. For both minerals the absorption coefficients are significantly smaller than those expected from the linear calibration curves given by Paterson (1982) and by Libowitzky and Rossman (1997).

---

Zusammenfassung OH-Absorptionskoeffizienten von Rutil und Cassiterit ermittelt durch Kernreaktions-Analyse und FTIR Spektroskopie Der OH-Gehalt von vier Rutil- und zwei Cassiterit-Einkristallen wurde mittels Kernreaktions-Analyse (NRA) und polarisierter FTIR Mikrospektroskopie untersucht. Die OH Absorptionsbanden beider Minerale sind um 3.300 cm⁻¹ zentriert, mit unterschiedlichen Absorptionserscheinungen. Der analytische H₂O-Gehalt, der mit NRA bestimmt wurde, schwankt von 70 bis 820 Gew.ppm. Die integrierten molaren Absorptionskoeffizienten, die auf den gesamten integrierten OH-Absorptionen basieren, betragen etwa 38.000 lċmol⁻¹ _H2Oċcm⁻² für Rutil and 65.000 lċmol⁻¹ _H2Oċcm⁻² für Cassiterit. Für beide Minerale sind die Absorptionskoeffizienten signifikant kleiner als die, die auf Grund der linearen Kalibrationskurven von Paterson (1982) und Libowitzky und Rossmann (1997) zu erwarten sind.

---

---

Received January 4, 2000; revised version accepted April 10, 2000