Heat capacities of aqueous solutions of sodium hydroxide and water ionization up to 300 °C at 10 MPa

A commercial (Setaram C80) calorimeter has been modified to measure the heat capacities of highly caustic solutions at temperatures up to 300 °C and pressures up to 20 MPa. The improvements have allowed more accurate determination of the isobaric volumetric heat capacities of chemically aggressive liquids at high temperatures. Test measurements with aqueous solutions of sodium chloride showed a reproducibility of about ±0.1%, with an accuracy of ∼0.3% or better, over the whole temperature range. Heat capacities of aqueous solutions of sodium hydroxide at concentrations from 0.5 to 8 mol/kg were measured at temperatures from 50 to 300 °C and a pressure of 10 MPa. Apparent molar isobaric heat capacities of NaOH(aq) were calculated using densities determined previously for the same solutions by vibrating-tube densimetry. Standard state (infinite dilution) partial molar isobaric heat capacities of NaOH(aq) were obtained by extrapolation using an extended Redlich-Meyer equation. Values of the standard heat capacity change for the ionization of water up to 300 °C were derived by combining the present results with the literature data for HCl(aq) and NaCl(aq).     

Empirical prediction of the Pitzer’s interaction parameters for cationic Al species with both SiO2(aq) and CO2(aq): Implications for the geochemical modelling of very saline solutions

The Pitzer’s interaction parameters, λ_N–M, involving the Mth cationic Al species Al³⁺ or AlOH²⁺ or AlO⁺ and the Nth neutral species SiO_2(aq) (at temperatures of 25–300 °C) or CO_2(aq) (at temperatures of 25–150 °C), have been evaluated through empirical linear relationships between λ_N–M and the surface electrostatic field of the ionic species of interest. These relationships have been obtained starting from the known λ_N–M for both SiO_2(aq) and CO_2(aq) with the main dissolved cations. The Pitzer’s interaction parameter thus estimated for the pair CO_2(aq)–Al³⁺ at 25 °C, 0.327, is 20–40% higher than the corresponding values obtained from CO₂ solubilities in concentrated solutions of AlCl₃, 0.272 ± 0.010 (2σ), and Al₂(SO₄)₃, 0.232 ± 0.002 (2σ), partly corroborating the empirical approach adopted in this study. To test the Pitzer’s interaction coefficients for cationic Al species with aqueous SiO₂, the log K values of the kaolinite dissolution reaction have been computed starting from available experimental data at 23–25 °C and ionic strengths of 0.0001–0.12 mol/kg adopting, alternatively, the Pitzer’s equations and the Debye–Hückel equation. A satisfactory agreement has been found between the log K values obtained through these two approaches, with maximum deviations of 0.11–0.12 log units. This good convergence of results is encouraging as it represents a necessary condition to prove the reliability of the Pitzer’s interaction coefficients estimated in this work. These results are a first step to take into account specific interactions among solutes in concentrated electrolyte solutions, such as those hosted in sedimentary basins or geothermal waters, for instance through the Pitzer’s equations. However, experimental or field data at higher ionic strengths are absolutely necessary to validate the reliability of the Pitzer’s interaction coefficients determined in this study.     

Comment on “Hydrolysis of neptunium(V) at variable temperatures (10-85 °C)” by L. Rao, T.G. Srinivasan, A.Yu. Garnov, P. Zanonato, P. Di Bernardo, and A. Bismondo

In a recent study [Rao, L., Srinivasan, T.G., Garnov, A.Yu., Zanonato, P., Di Bernardo, P., Bismondo, A., 2004. Hydrolysis of neptunium(V) at variable temperatures (10-85 °C). Geochim. Cosmochim. Acta68, 4821-4830.] the hydrolysis of Np(V) was investigated at 10-85 °C by absorption spectroscopy, potentiometry, and microcalorimetry along the titration of Np(V) solutions with tetramethylammonium hydroxide up to pH 10. However, there is strong evidence that the precautions to avoid competing reactions with carbonate were not sufficient and that the measured effects are not caused by the formation of Np(V) hydroxide complexes but primarily by the formation of Np(V) carbonate complexes. The reported equilibrium constants, enthalpies, entropies, and heat capacities for the complexes NpO₂OH(aq) and are severely in error and must not be used for the geochemical modeling of neptunium. If the hydrolysis constants reported by Rao et al. [Rao, L., Srinivasan, T.G., Garnov, A.Yu., Zanonato, P., Di Bernardo, P., Bismondo, A., 2004. Hydrolysis of neptunium(V) at variable temperatures (10-85 °C). Geochim. Cosmochim. Acta68, 4821-4830] are used to calculate neptunium solubilities in alkaline solutions relevant for nuclear waste repositories, the Np(V) concentrations are overestimated by orders of magnitude.     

Effect of water on the heat capacity of polymerized aluminosilicate glasses and melts

The effect of water on heat capacity has been determined for four series of hydrated synthetic aluminosilicate glasses and supercooled liquids close to albite, phonolite, trachyte, and leucogranite compositions. Heat capacities were measured at atmospheric pressure by differential scanning calorimetry for water contents between 0 and 4.9 wt % from 300 K to about 100 K above the glass transition temperature (T_g). The partial molar heat capacity of water in polymerized aluminosilicate glasses, which can be considered as independent of composition, is (J/mol K). In liquids containing at least 1 wt % H₂O, the partial molar heat capacity of water is about 85 J/mol K. From speciation data, the effects of water as hydroxyl groups and as molecular water have tentatively been estimated, with partial molar heat capacities of 153 ± 18 and 41 ± 14 J/mol K, respectively. In all cases, water strongly increases the configurational heat capacity at T_g and exerts a marked depressing effect on T_g, in close agreement with the results of viscosity experiments on the same series of glasses. Consistent with the Adam and Gibbs theory of relaxation processes, the departure of the viscosity of hydrous melts from Arrhenian variations correlates with the magnitude of configurational heat capacities.     

Solubility products of amorphous ferric arsenate and crystalline scorodite (FeAsO4 · 2H2O) and their application to arsenic behavior in buried mine tailings

Published solubility data for amorphous ferric arsenate and scorodite have been reevaluated using the geochemical code PHREEQC with a modified thermodynamic database for the arsenic species. Solubility product calculations have emphasized measurements obtained under conditions of congruent dissolution of ferric arsenate (pH < 3), and have taken into account ion activity coefficients, and ferric hydroxide, ferric sulfate, and ferric arsenate complexes which have association constants of 10^4.04 (FeH₂AsO₄²⁺), 10^9.86 (FeHAsO₄⁺), and 10^18.9 (FeAsO₄). Derived solubility products of amorphous ferric arsenate and crystalline scorodite (as log K_sp) are −23.0 ± 0.3 and −25.83 ± 0.07, respectively, at 25 °C and 1 bar pressure. In an application of the solubility results, acid raffinate solutions (molar Fe/As = 3.6) from the JEB uranium mill at McClean Lake in northern Saskatchewan were neutralized with lime to pH 2-8. Poorly crystalline scorodite precipitated below pH 3, removing perhaps 98% of the As(V) from solution, with ferric oxyhydroxide (FO) phases precipitated starting between pH 2 and 3. Between pH 2.18 and 7.37, the apparent log K_sp of ferric arsenate decreased from −22.80 to −24.67, while that of FO (as Fe(OH)₃) increased from −39.49 to −33.5. Adsorption of As(V) by FO can also explain the decrease in the small amounts of As(V)(aq) that remain in solution above pH 2-3. The same general As(V) behavior is observed in the pore waters of neutralized tailings buried for 5 yr at depths of up to 32 m in the JEB tailings management facility (TMF), where arsenic in the pore water decreases to 1-2 mg/L with increasing age and depth. In the TMF, average apparent log K_sp values for ferric arsenate and ferric hydroxide are −25.74 ± 0.88 and −37.03 ± 0.58, respectively. In the laboratory tests and in the TMF, the increasing crystallinity of scorodite and the amorphous character of the coexisting FO phase increases the stability field of scorodite relative to that of the FO to near-neutral pH values. The kinetic inability of amorphous FO to crystallize probably results from the presence of high concentrations of sulfate and arsenate.     

华北石炭——二叠纪煤系粘土岩夹矸中铝的氢氧化物矿物研究

华北石炭——二叠纪煤系太原组和山西组不同煤层粘土岩夹矸中,广泛含有一水软铝石、一水硬铝石和三水铝石。大多数一水铝石矿物呈薄纹层状、似层状及层状聚集于富岩质及富有有机质粘土岩夹矸条带中,有的呈细脉状沿炭质条带及粘土岩夹矸层理的节理裂隙分布,并与玉髓、方解石、高岭石或迪开石和铵云母细脉共生。铝的氢氧化物矿物可以作为成岩古温度和成岩阶段划分的标志矿物。最后还讨论了该矿物的成因。     

Speciation and thermodynamic properties for cobalt chloride complexes in hydrothermal fluids at 35-440 °C and 600 bar: An in-situ XAS study

Aqueous Co(II) chloride complexes play a crucial role in cobalt transport and deposition in ore-forming hydrothermal systems, ore processing plants, and in the corrosion of special Co-bearing alloys. Reactive transport modelling of cobalt in hydrothermal fluids relies on the availability of thermodynamic properties for Co complexes over a wide range of temperature, pressure and salinity. Synchrotron X-ray absorption spectroscopy was used to determine the speciation of cobalt(II) in 0-6 m chloride solutions at temperatures between 35 and 440 °C at a constant pressure of 600 bar. Qualitative analysis of XANES spectra shows that octahedral species predominate in solution at 35 °C, while tetrahedral species become increasingly important with increasing temperature. Ab initio XANES calculations and EXAFS analyses suggest that in high temperature solutions the main species at high salinity (Cl:Co >> 2) is CoCl₄²⁻, while a lower order tetrahedral complex, most likely CoCl₂(H₂O)_2(aq), predominates at low salinity (Cl:Co ratios ∼2). EXAFS analyses further revealed the bonding distances for the octahedral Co(H₂O)₆²⁺ (^octCo-O = 2.075(19) Å), tetrahedral CoCl₄²⁻ (^tetCo-Cl = 2.252(19) Å) and tetrahedral CoCl₂(H₂O)_2(aq) (^tetCo-O = 2.038(54) Å and ^tetCo-Cl = 2.210(56) Å). An analysis of the Co(II) speciation in sodium bromide solutions shows a similar trend, with tetrahedral bromide complexes becoming predominant at higher temperature/salinity than in the chloride system. EXAFS analysis confirms that the limiting complex at high bromide concentration at high temperature is CoBr₄²⁻. Finally, XANES spectra were used to derive the thermodynamic properties for the CoCl₄²⁻ and CoCl₂(H₂O)_2(aq) complexes, enabling thermodynamic modelling of cobalt transport in hydrothermal fluids. Solubility calculations show that tetrahedral CoCl₄²⁻ is responsible for transport of cobalt in hydrothermal solutions with moderate chloride concentration (∼2 m NaCl) at temperatures of 250 °C and higher, and both cooling and dilution processes can cause deposition of cobalt from hydrothermal fluids.     

The structure of monomeric and dimeric uranyl adsorption complexes on gibbsite: A combined DFT and EXAFS study

We investigated the structure of uranyl sorption complexes on gibbsite (pH 5.6-9.7) by two independent methods, density functional theory (DFT) calculations and extended X-ray absorption fine structure (EXAFS) spectroscopy at the U-L_III edge. To model the gibbsite surface with DFT, we tested two Al (hydr)oxide clusters, a dimer and a hexamer. Based on polarization, structure, and relaxation energies during geometry optimization, the hexamer cluster was found to be the more appropriate model. An additional advantage of the hexamer model is that it represents both edges and basal faces of gibbsite. The DFT calculations of (monomeric) uranyl sorption complexes show an energetic preference for the corner-sharing versus the edge-sharing configuration on gibbsite edges. The energy difference is so small, however, that possibly both surface species may coexist. In contrast to the edge sites, sorption to basal sites was energetically not favorable. EXAFS spectroscopy revealed in all investigated samples the same interatomic distances of the uranyl coordination environment (R_U-Oax ≈ 1.80 Å, R_U-Oeq ≈ 2.40 Å), and towards the gibbsite surface (R_U-O ≈ 2.87 Å, R_U-Al ≈ 3.38 Å). In addition, two U-U distances were observed, 3.92 Å at pH 9.7 and 4.30 Å at pH 5.6, both with coordination numbers of ∼1. The short U-U distance is close to that of the aqueous uranyl hydroxo dimer, UO₂(OH)₂, reported as 3.875 Å in the literature, but significantly longer than that of aqueous trimers (3.81-3.82 Å), suggesting sorption of uranyl dimers at alkaline pH. The longer U-U distance (4.30 Å) at acidic pH, however, is not in line with known aqueous uranyl polymer complexes. Based on the EXAFS findings we further refined dimeric surface complexes with DFT. We propose two structural models: in the acidic region, the observed long U-U distance can be explained with a distortion of the uranyl dimer to form both a corner-sharing and an edge-sharing linkage to neighboring Al octahedra, leading to R_U-U = 4.150 Å. In the alkaline region, a corner-sharing uranyl dimer complex is the most favorable. The U-O path at ∼2.87 Å in the EXAFS spectra arises from the oxygen atom linking two Al cations in corner-sharing arrangement. The adsorption structures obtained by DFT calculations are in good agreement with the structural parameters from EXAFS analysis: U-Al (3.394 Å), U-U (3.949 Å), and U-O (2.823 Å) for the alkaline pH model, and U-Al (3.279 Å), U-U (4.150 Å), and U-O (2.743 Å) for the acidic pH model. This work shows that by combining EXAFS and DFT, consistent structural models for uranyl sorption complexes can be obtained, which are relevant to predict the migration behavior of uranium at nuclear facilities.     

Thermodynamic studies of pyrrhotite-pyrite equilibria in the Ag-Fe-S system by solid-state galvanic cell technique at 518-723 K and total pressure of 1 atm

The reaction FeS₂(cr) + 2Ag(cr) = ‘FeS’(cr) + Ag₂S(cr) was studied by measuring the temperature dependence of the electromotive force (EMF) of the all-solid-state galvanic cell with common gas space:

(-)Pt|Ag|AgI|Ag₂S,‘FeS’,FeS₂|Pt(+)     

Composition and structure of an iron-bearing, layered double hydroxide (LDH) - Green rust sodium sulphate

Mixed-valent Fe(II),Fe(III)-layered hydroxide, known as green rust, was synthesized from slightly basic, sodium sulphate solutions in an oxygen-free glove box. Solution conditions were monitored with pH and Eh electrodes and optimized to ensure a pure sulphate green-rust phase. The solid was characterised using Mössbauer spectroscopy, X-ray diffraction, scanning electron microscopy and atomic force microscopy. The composition of the solution from which the green rust precipitated was established by mass and absorption spectroscopy. The sulphate form of green rust is composed of brucite-like layers with Fe(II) and Fe(III) in an ordered distribution. The interlayers contain sulphate, water and sodium in an arrangement characteristic for the nikischerite group. The crystal structure is highly disordered by stacking faults. The composition, formula and crystallographic parameters are: NaFe(II)₆Fe(III)₃(SO₄)₂(OH)₁₈·12H₂O, space group P-3, a = 9.528(6) Å, c = 10.968(8) Å and Z = 1. Green rust sodium sulphate, GR_Na,SO4, crystallizes in thin, hexagonal plates. Particles range from less than 50 nm to 2 μm in diameter and are 40 nm thick or less. The material is redox active and reaction rates are fast. Extremely small particle size and high surface area contribute to rapid oxidation, transforming green rust to an Fe(III)-phase within minutes.     

A calorimetric investigation of spessartine: Vibrational and magnetic heat capacity

The heat capacity (C_p) of two synthetic spessartine samples (Sps) was measured on 20-30 mg-size samples in the temperature range 2-864 K by relaxation calorimetry (RC) and differential scanning calorimetry (DSC). The polycrystalline spessartine samples were synthesized in two different laboratories at high pressures and temperatures from glass and oxide-mixture starting materials and characterized by X-ray powder diffraction and electron-microprobe analysis. The low-temperature heat capacity data show a prominent lambda transition with a peak at 6.2 K, which is interpreted to be the result of a paramagnetic-antiferromagnetic phase transition. The DSC data around ambient T agree excellently with the RC data and can be represented by the C_p polynomial for T > 250 K:

     

Iron and manganese diagenesis in deep sea volcanogenic sediments and the origins of pore water colloids

Volcanogenic sediments are typically rich in Fe and Mn-bearing minerals that undergo substantial alteration during early marine diagenesis, however their impact on the global biogeochemical cycling of Fe and Mn has not been widely addressed. This study compares the near surface (0-20 cm below sea floor [cmbsf]) aqueous (<0.02 μm) and aqueous + colloidal here in after ‘dissolved’ (<0.2 μm) pore water Fe and Mn distributions, and ancillary O_2(aq), and solid-phase reactive Fe distributions, between two volcanogenic sediment settings: [1] a deep sea tephra-rich deposit neighbouring the volcanically active island of Montserrat and [2] mixed biosiliceous-volcanogenic sediments from abyssal depths near the volcanically inactive Crozet Islands archipelago. Shallow penetration of O_2(aq) into Montserrat sediments was observed (<1 cmbsf), and inferred to partially reflect oxidation of fine grained Fe(II) minerals, whereas penetration of O_2(aq) into abyssal Crozet sediments was >5 cmbsf and largely controlled by the oxidation of organic matter. Dissolved Fe and Mn distributions in Montserrat pore waters were lowest in the surface oxic-layer (0.3 μM Fe; 32 μM Mn), with maxima (20 μM Fe; 200 μM Mn) in the upper 1-15 cmbsf. Unlike Montserrat, Fe and Mn in Crozet pore waters were ubiquitously partitioned between 0.2 μm and 0.02 μm filtrations, indicating that the pore water distributions of Fe and Mn in the (traditionally termed) ‘dissolved’ size fraction are dominated by colloids, with respective mean abundances of 80% and 61%. Plausible mechanisms for the origin and composition of pore water colloids are discussed, and include prolonged exposure of Crozet surface sediments to early diagenesis compared to Montserrat, favouring nano-particulate goethite formation, and the elevated dissolved Si concentrations, which are shown to encourage fine-grained smectite formation. In addition, organic matter may stabilise authigenic Fe and Mn in the Crozet pore waters. We conclude that volcanogenic sediment diagenesis leads to a flux of colloidal material to the overlying bottom water, which may impact significantly on deep ocean biogeochemistry. Diffusive flux estimates from Montserrat suggest that diagenesis within tephra deposits of active island volcanism may also be an important source of dissolved Mn to the bottom waters, and therefore a source for the widespread hydrogenous MnO_x deposits found in the Caribbean region.     

Influence of physiology and climate on δD of leaf wax n-alkanes from C3 and C4 grasses

We measured hydrogen isotope compositions (δD) of high-molecular-weight n-alkanes (C₂₇-C₃₃) from grasses grown in greenhouses and collected from the US Great Plains. In both cases, n-alkanes from C₄ grasses are enriched in D by more than 20‰ relative to those from C₃ grasses. The apparent enrichment factor (ε_C29-GW) between C₂₉n-alkane and greenhouse water is −165 ± 12‰ for C₃ grasses and −140 ± 15‰ for C₄ grasses. For samples from the Great Plains, δD values of C₂₉n-alkanes range from −280 to −136‰, with values for C₄ grasses ca. 21‰ more positive than those for C₃ grasses from the same site. Differences in C₃ and C₄ grass n-alkane δD values are consistent with the shorter interveinal distance in C₄ grass leaves, and greater back-diffusion of enriched water from stomata to veins, than in C₃ grass leaves. Great Plains’ grass n-alkane isotopic ratios largely reflect precipitation δD values. However, the offset or apparent fractionation between n-alkanes and precipitation is not uniform and varies with annual precipitation and relative humidity, suggesting climatic controls on lipid δD values. The dryer sites exhibit smaller absolute apparent fractionation indicative of D-enrichment of source waters through transpiration and/or soil evaporation. To explore the relationship between climate and n-alkane δD values, we develop three models. (1) The ‘direct analog’ model estimates δD_C29 values simply by applying the apparent enrichment factors, ε_C29-GW, observed in greenhouse grasses to precipitation δD values from the Great Plains. (2) The ‘leaf-water’ model uses a Craig-Gordon model to estimate transpirational D-enrichment for both greenhouse and field sites. The transpiration-corrected enrichment factors between C₂₉ and bulk leaf-water, ε_C29-GW, calculated from the greenhouse samples (−181‰ for C₃ and −157‰ for C₄) are applied to estimate δD_C29 values relative to modeled bulk leaf-water δD values. (3) The ‘soil- and leaf-water’ model estimates the combined effects of soil evaporation, modeled by analogy with a flow-through lake, and transpiration on δD_C29 values. Predictions improve with the addition of the explicit consideration of transpiration and soil evaporation, indicating that they are both important processes in determining plant lipid δD values. D-enrichment caused by these evaporative processes is controlled by relative humidity, suggesting that important climatic information is recorded in leaf wax n-alkane δD values. Calibration studies such as this one provide a baseline for future studies of plant-water-deuterium systematics and form the foundation for interpretation of plant wax hydrogen isotope ratios as a paleo-aridity proxy.     

Characteristics of alkenones synthesized by a bloom of emiliania huxleyi in the Bering Sea

We investigated the characteristics of the alkenones produced by a bloom of Emiliania huxleyi in the eastern Bering Sea in 2000. Alkenones were detected in surface waters between 57°N and 63°N, where phosphate concentrations were low and the ammonium/nitrate ratio was high. The total alkenone content (C_37:2, C_37:3, and C_37:4) ranged from 22.0 to 349 μg g⁻¹ in suspended particles and from 0.109 to 1.42 μg g⁻¹ in surface sediments. This suggests that a large proportion of the particulate alkenones synthesized in the surface water rapidly degraded within the water column and/or at the water-sediment interface of the Bering Shelf. The change in the stable carbon isotopic composition (δ¹³C) of C_37:3 alkenone could not be explained only by variation in [CO₂(aq)] in the surface water but also depended on the growth rate of E. huxleyi. The alkenone unsaturation index (U^K′₃₇) was converted into an alkenone “temperature” with three equations [Prahl et al 1988], [Sikes et al 1997] and [Müller et al 1998]; Sikes et al.’s (1997) equation gave the best correlation with the observed sea surface temperature (SST) in the eastern Bering Sea. However, some temperatures estimated by Sikes et al.’s (1997) equation from the U^K′₃₇ varied from the observed SST, possibly because of the rapidly changing rate of alkenone synthesis in the logarithmic growth stage or the low rate of alkenone synthesis when nutrients were limiting. Temperatures estimated from U^K′₃₇ in the surface sediments (6.8-8.2°C) matched the observed SST in September (7-8°C) but differed from the annual average SST of 4 to 5°C, suggesting that most of the alkenone in the eastern Bering Sea was synthesized during limited periods, for instance, in September. The relative amounts of C_37:4 alkenone as proportions of the total alkenones (referred to as C_37:4%) were high, ranging from 18.3 to 41.4%. Low-salinity water (<32 psu) within the study area would have contributed to the high C_37:4% because a negative linear relationship between C_37:4% and salinity was found in this study.     

Revised values for the Gibbs free energy of formation of [Al(OH)4 aq−], diaspore,boehmite and bayerite at 298.15 K and 1 bar,the thermodynamic properties of kaolinite to 800 K and 1 bar,and the heats of solution of several gibbsite samples

Solution calorimetric measurements compared with solubility determinations from the literature for the same samples of gibbsite have provided a direct thermochemical cycle through which the Gibbs free energy of formation of [Al(OH)_{4 aq}^?] can be determined. The Gibbs free energy of formation of [Al(OH)_{4 aq}^?] at 298.15 K is ?1305 ± 1 kJ/mol. These heat-of-solution results show no significant difference in the thermodynamic properties of gibbsite particles in the range from 50 to 0.05 μm.The Gibbs free energies of formation at 298.15 K and 1 bar pressure of diaspore, boehmite and bayerite are ?9210 ± 5.0, ?918.4 ± 2.1 and ?1153 ± 2 kJ/mol based upon the Gibbs free energy of [A1(OH)_{4 aq}^?] calculated in this paper and the acceptance of ?1582.2 ± 1.3 and ?1154.9 ± 1.2 kJ/mol for the Gibbs free energy of formation of corundum and gibbsite, respectively.Values for the Gibbs free energy formation of [Al(OH)_{2 aq}⁺] and [AlO_{2 aq}^?] were also calculated as ?914.2 ± 2.1 and ?830.9 ± 2.1 kJ/mol, respectively. The use of [AlC_{2 aq}^?] as a chemical species is discouraged.A revised Gibbs free energy of formation for [H₄SiO_4aq⁰] was recalculated from calorimetric data yielding a value of ?1307.5 ± 1.7 kJ/mol which is in good agreement with the results obtained from several solubility studies.Smoothed values for the thermodynamic functions C_P⁰, ($\text{H}{}_{\mathrm{T}}{}^0\text{- H}{}_{298}{}^0\text{)}\text{T}$, $\text{(}\text{G}{}_{\mathrm{T}}{}^0\text{- H}{}_{298}{}^0\text{)}\text{T}$, S_T⁰ - S₀⁰, $\text{ΔH}{}_{\mathrm{?,298}}{}^0$ kaolinite are listed at integral temperatures between 298.15 and 800 K. The heat capacity of kaolinite at temperatures between 250 and 800 K may be calculated from the following equation: C_P⁰ = 1430.26 ? 0.78850 T + 3.0340 × 10^?4T² ?1.85158 × 10^?4T²$\text{1}\text{2}\text{+ 8.3341 × 10}{}^6\text{T}{}^{\mathrm{?2}}$.The thermodynamic properties of most of the geologically important Al-bearing phases have been referenced to the same reference state for Al, namely gibbsite.     

Dissolution of illite in saline-acidic solutions at 25 °C

The dissolution rate of illite, a common clay mineral in Australian soils, was studied in saline-acidic solutions under far from equilibrium conditions. The clay fraction of Na-saturated Silver Hill illite (K_1.38Na_0.05)(Al_2.87Mg_0.46Fe³⁺_0.39Fe²⁺_0.28Ti_0.07)[Si_7.02Al_0.98]O₂₀(OH)₄ was used for this study. The dissolution rates were measured using flow-through reactors at 25 ± 1 °C, solution pH range of 1.0-4.25 (H₂SO₄) and at two ionic strengths (0.01 and 0.25 M) maintained using NaCl solution. Illite dissolution rates were calculated from the steady state release rates of Al and Si. The dissolution stoichiometry was determined from Al/Si, K/Si, Mg/Si and Fe/Si ratios. The release rates of cations were highly incongruent during the initial stage of experiments, with a preferential release of Al and K over Si in majority of the experiments. An Al/Si ratio >1 was observed at pH 2 and 3 while a ratio close to the stoichiometric composition was observed at pH 1 and 4 at the higher ionic strength. A relatively higher K⁺ release rate was observed at I = 0.25 in 2-4 pH range than at I = 0.01, possibly due to ion exchange reaction between Na⁺ from the solution and K⁺ from interlayer sites of illite. The steady state release rates of K, Fe and Mg were higher than Si over the entire pH range investigated in the study. From the point of view of the dominant structural cations (Si and Al), stoichiometric dissolution of illite occurred at pH 1-4 in the higher ionic strength experiments and at pH ?3 for the lower ionic strength experiments. The experiment at pH 4.25 and at the lower ionic strength exhibited lower R_Al (dissolution rate calculated from steady state Al release) than R_Si (dissolution rate calculated from steady state Si release), possibly due to the adsorption of dissolved Al as the output solutions were undersaturated with respect to gibbsite. The dissolution of illite appears to proceed with the removal of interlayer K followed by the dissolution of octahedral cations (Fe, Mg and Al), the dissolution of Si is the limiting step in the illite dissolution process. A dissolution rate law showing the dependence of illite dissolution rate on proton concentration in the acid-sulfate solutions was derived from the steady state dissolution rates and can be used in predicting the impact of illite dissolution in saline acid-sulfate environments. The fractional reaction orders of 0.32 (I = 0.25) and 0.36 (I = 0.01) obtained in the study for illite dissolution are similar to the values reported for smectite. The dissolution rate of illite is mainly controlled by solution pH and no effect of ionic strength was observed on the dissolution rates.     

Melting of enriched Archean subcontinental lithospheric mantle: Evidence from the ca. 1760 Ma volcanic rocks of the Xiong’er Group,southern margin of the North China Craton

The Xiong’er Group is an important geologic unit in the southern margin of the North China Craton. It is dominated by the volcanic rocks, dated at 1763 ± 15 Ma, that have SiO₂ contents ranging from 52.10 wt% to 73.51 wt%. These volcanic rocks are sub-alkaline and can be classified into three subgroups: basaltic andesites, andesites and rhyolites. They unexceptionally show enrichment of light rare earth elements (LREE) and share similar trace element patterns. Depletions in Nb, Ta, Sr, P and Ti relative to the adjacent elements are evident for all the samples. The volcanic rocks are evolved with low MgO contents (0.29–5.88 wt%) and accordingly low Mg^# values of 11–53. The Nd isotopes are enriched and show a weak variation with ?_Nd(t) = −7.12 to −9.63. Zircon Hf isotopes are also enriched with ?_Hf(t) = −12.02 ± 0.45. The volcanic rocks of the Xiong’er Group are interpreted to represent fractional crystallization of a common mantle source. The volcanic rocks might have been generated by high-degree partial melting of a lithospheric mantle that was originally modified by the oceanic subduction in the Late Archean. This brings a correlation with the subduction-modified lithospheric mantle in an extensional setting during breakup of the Columbia supercontinent in the late Paleoproterozoic, rather than in an arc setting. The elevated SiO₂ contents and evolved radiogenic isotope features indicate the possible incorporation into their source of lower crustal materials that have similar Nd isotopic characteristics to the subcontinental lithospheric mantle. The existence of extensive Xiong’er volcanic rocks (60,000 km²) indicates an early large-scale subduction-related metasomatism in the area and probably suggest a flat subduction model for the plate-margin magmatism in the Late Archean.     

Coupled phase and aqueous species equilibrium of the H2O-CO2-NaCl-CaCO3 system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite

A model is developed for the calculation of coupled phase and aqueous species equilibrium in the H₂O-CO₂-NaCl-CaCO₃ system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite. The vapor-liquid-solid (calcite, halite) equilibrium together with the chemical equilibrium of H⁺, Na⁺, Ca²⁺, , Ca(OH)⁺, OH⁻, Cl⁻, , , CO_2(aq) and CaCO_3(aq) in the aqueous liquid phase as a function of temperature, pressure, NaCl concentrations, CO_2(aq) concentrations can be calculated, with accuracy close to those of experiments in the stated T-P-m range, hence calcite solubility, CO₂ gas solubility, alkalinity and pH values can be accurately calculated. The merit and advantage of this model is its predictability, the model was generally not constructed by fitting experimental data.One of the focuses of this study is to predict calcite solubility, with accuracy consistent with the works in previous experimental studies. The resulted model reproduces the following: (1) as temperature increases, the calcite solubility decreases. For example, when temperature increases from 273 to 373 K, calcite solubility decreases by about 50%; (2) with the increase of pressure, calcite solubility increases. For example, at 373 K changing pressure from 10 to 500 bar may increase calcite solubility by as much as 30%; (3) dissolved CO₂ can increase calcite solubility substantially; (4) increasing concentration of NaCl up to 2 m will increase calcite solubility, but further increasing NaCl solubility beyond 2 m will decrease its solubility.The functionality of pH value, alkalinity, CO₂ gas solubility, and the concentrations of many aqueous species with temperature, pressure and NaCl_(aq) concentrations can be found from the application of this model. Online calculation is made available on www.geochem-model.org/models/h2o_co2_nacl_caco3/calc.php.     

Oxygen isotope fractionation between hydroxide minerals and water

The modified increment method has been applied to the calculation of oxygen isotope fractionation factors for hydroxide minerals. The results suggest the following sequence of ¹⁸O-enrichment in the common hydroxides: limonite > gibbsite > goethite > brucite > diaspore. The hydroxides are significantly enriched in ¹⁸O relative to the corresponding oxides. The sequence of ¹⁸O-enrichment in the hydroxides and oxides of trivalent cations is as follows: M(OH)₃ > MO(OH) > M₂O₃. There are also considerable fractionations within the polymorphos of Al(OH)₃. The internally consistent fractionation factors for hydroxide–water systems are obtained for the temperature range of 0 to 1200 °C, which are comparable with the data derived from synthesis experiments and natural samples at surficial temperatures. Temperature dependence of oxygen isotope fractionations between goethite, gibbsite, boehmite and diaspore and water are significant enough for the purpose of geothermometry. Thus the hydroxide–water pairs hold great promise of serving as reliable paleothermometers in surficial geological environments. Received: 22 January 1997 / Revised, accepted: 2 June 1997     

Water diffusion in dacitic melt

H₂O diffusion in dacitic melt was investigated at 0.48-0.95 GPa and 786-893 K in a piston-cylinder apparatus. The diffusion couple design was used, in which a nominally dry dacitic glass makes one half and is juxtaposed with a hydrous dacitic glass containing up to ∼8 wt.% total water (H₂O_t). H₂O concentration profiles were measured on quenched glasses with infrared microspectroscopy. The H₂O diffusivity in dacite increases rapidly with water content under experimental conditions, similar to previous measurements at the same temperature but at pressure <0.15 GPa. However, compared with the low-pressure data, H₂O diffusion at high pressure is systematically slower. H₂O diffusion profiles in dacite can be modeled by assuming molecular H₂O (H₂O_m) is the diffusing species. Total H₂O diffusivity D_H2Ot within 786-1798 K, 0-1 GPa, and 0-8 wt.% H₂O_t can be expressed as: where D_H2Ot is in m²/s, T is temperature in K, P is pressure in GPa, K = exp(1.49 − 2634/T) is the equilibrium constant of speciation reaction (H₂O_m+O?2OH) in the melt, X = C/18.015/[C/18.015 + (100 − C)/33.82], C is wt.% of H₂O_t, and 18.015 and 33.82 g/mol correspond to the molar masses of H₂O and anhydrous dacite on a single oxygen basis. Compared to H₂O diffusion in rhyolite, diffusivity in dacite is lower at intermediate temperatures but higher at superliquidus temperatures. This general H₂O diffusivity expression can be applied to a broad range of geological conditions, including both magma chamber processes and volcanic eruption dynamics from conduit to the surface.