Dehydration of dioctahedral aluminous phyllosilicates: thermodynamic modelling and implications for thermobarometric estimates

We propose a solid-solution model for dioctahedral aluminous phyllosilicates accounting for the main compositional variations, including hydration, observed in natural smectites, interlayered illite/smectite, illites, and phengites from diagenetic to high-grade metamorphic conditions. The suggested formalism involves dehydrated micas and hydrated pyrophyllite-like thermodynamic end-members. With these end-members, the equilibrium conditions of quartz + water + K-bearing mica-like phyllosilicates of fixed 2:1 composition are represented by a line in P–T space along which the interlayer water content varies. The relevant thermodynamic properties required for the calculation of equilibrium conditions were derived using a set of 250 natural data of known maximal temperature and pressure conditions, which covers a range between 25°C and few MPa to 800°C and 5 GPa. The temperatures calculated at fixed pressure with our model are in fair agreement with those reported in the literature for the 250 natural data. At low temperature and pressure, the amount of interlayer water in K-deficient phengite and illite is predicted to reach 100% of the apparent vacancies, which is consistent with previous values reported in the literature. Although the amount of interlayer water is predicted to decrease with pressure and temperature, it is calculated to be significant in K-deficient phengite from LT–HP pelites metamorphosed at about 350°C, 10 kbar. The presence of molecular water in the interlayer site of such phengites has been confirmed by FTIR mapping. Its implications for P–T estimates are discussed.     

Temperature — composition relationships of authigenic micaceous minerals in the Los Azufres geothermal system

Smectite, illite, celadonite and chlorite are the major products of alteration of rhyolites and andesites, in the upper part of the Los Azufres geothermal system. Changes in mineral assemblages and composition of phases are observed as a function of depth and host rock lithology. Two different sequences characterize the rhyolites and the andesites from the surface to a depth of about 1500 m: kaolinitesmectite (±interlayered illite/smectite)illitemuscovite (rhyolites), and kaolinitesmectite (±interlayered illite/smectite)illite+celadoniteillite+chloritechlorite (andesites).Illite, and chlorite at depth, are largely dominant. Similar substitutions and correlations among chemical constituents characterize illites from rhyolites or andesites, but their compositions in the two host rock lithologies exhibit slight but significant differences, especially in their Fe and Mg contents which are the highest in illites from andesites. Illite exhibits progressive changes in composition with depth: a strong increase in the K content in the interlayer, together with an increase of the Fe content in the octahedral site. These changes correspond to a slight increase in the molar fraction of Fe-(Mg) celadonite end-members, and mostly to a dramatic decrease of pyrophyllite with increasing temperature.Temperature of the mineralogical and compositional changes was estimated from fluid inclusions studies, combined with other geothermometric approaches (chemical geothermometers and direct measurements). Variation of X-pyrophyllite with temperature is proposed as a geothermometer for different host rock lithologies. Transitions between the stability fields of illite±interstratified illite-smectite and illite+chlorite is around 200±30° C, and between illite+chlorite and chlorite around 290±20° C.     

Estimation of land subsidence caused by loss of smectite-interlayer water in shallow aquifer systems

Traditionally, land subsidence that results from groundwater over-pumping has often been described by the theory of consolidation. The mechanism of land subsidence due to the dehydration of clay minerals is not well addressed. A model of the “hydration state of smectite”, and a “solid solution model of smectite dehydration”, incorporating a thermodynamic solid solution model and laboratory results concerning clay-water systems of swelling pressure, hydration state and basal spacing in smectite interlayer, are employed to examine the effect of the release of water from the smectite interlayer on land subsidence in the coastal area of the Chou-Shui River alluvial fan and the Yun Lin offshore industrial infrastructure complex in Taiwan. The results indicate that 9.56–22.80% of the total cumulative land subsidence to a depth of 300 m is consistent with smectite dehydration following the over-pumping of groundwater. This dehydration-related land subsidence occurred to a depth of 0–60 m, with subsidence due to smectite dehydration accounting for 6.20–13.32% of the primary consolidation. Additionally, the total amount of subsidence resulting from both smectite dehydration and primary consolidation is consistent with the subsidence observed in the field. This study reveals that smectite dehydration appears to be important in assessing and predicting land subsidence in shallow aquifer systems.     

Geochemical and mineralogical impacts of H2SO4 on clays between pH 5.0 and −3.0

Natural and constructed clay liners are routinely used to contain waste and wastewater. The impact of acidic solutions on the geochemistry and mineralogy of clays has been widely investigated in relation to acid mine drainage systems at pH > 1.0. The impact of H₂SO₄ leachate characterized by pH < 1.0 and potentially negative pH values on the geochemistry and mineralogy of clays is, however, not clear. Thus, laboratory batch experiments were conducted on three natural clay samples with different mass ratios of smectite, illite and kaolinite to investigate the impact of H₂SO₄ on the geochemistry and mineralogy of aluminosilicates from pH 5.0 to −3.0. Batch testing was conducted at seven pH treatments (5.0, 3.0, 1.0, 0.0, −1.0, −2.0 and −3.0) using standardized H₂SO₄ solutions for four exposure periods (14, 90, 180, and 365 d). Aqueous geochemical and XRD analyses showed: increased dissolution of aluminosilicates with decreasing pH and increasing exposure period, that smectite was more susceptible to dissolution than illite and kaolinite, precipitation of an amorphous silica phase occurred at pH ? 0.0, and anhydrite precipitated in Ca-rich clays at pH ? −1.0. In addition, global dissolution rates were calculated for the clays and showed good agreement to literature smectite, illite and kaolinite dissolution rates, which suggests global dissolution rates for complex clays could be determined from monomineralic studies. A stepwise conceptual model of the impact of H₂SO₄ on aluminosilicate geochemistry and mineralogy between pH 5.0 and −3.0 is proposed.     

Hydration thermodynamics of the SWy-1 montmorillonite saturated with alkali and alkaline-earth cations: A predictive model

The aim of the present work was to study the thermodynamic equilibria between water and a homo-ionic montmorillonite SWy-1 saturated by different cations. The choice of this smectite is justified by the large set of experimental data available from the literature for eight different interlayer cations: Na⁺, K⁺, Rb⁺, Cs⁺, Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺. In particular, studies by [Cases et al., 1992], [Cases et al., 1997] and [Bérend et al., 1995] are providing heat of adsorption data, pairs of desorption-adsorption isotherms, and information about the partition of adsorption-desorption water molecules between external surfaces and internal spaces. By calculating the effective amount of hydration water as the difference between the so-called gravimetric water and the surface covering water, a thermodynamic model was then developed, based on the concept of Ransom and Helgeson (1994) considering an asymmetric subregular binary solid solution between a fully hydrated and a anhydrous smectite. A set of six thermodynamic parameters ( and four Margules parameters) was extracted by a least square method from measurements of enthalpies of adsorption and paired adsorption-desorption isotherms for each interlayer cation. These six initial parameters were then used to determine a complete set of standard thermodynamic hydration parameters (, heat capacity, molar volume, and number of interlayer H₂O) and quantify, for each cation, the number of moles of hydration water molecules as a function of relative humidity and temperature. The validation of the standard state thermodynamic properties of hydration for each end member was carried out using three approaches: (1) a comparison with experimental isotherms obtained on hetero-ionic and homo-ionic SWy-1 smectite at different temperatures; (2) a comparison with the experimental integral enthalpy and entropy of hydration of the SWy-1 smectite; and (3) a comparison with experimental isotherms acquired on various smectites (Upton, MX80, Arizona) with different layer charges.Eventually, the present work demonstrates that, from a limited number of measurements, it is possible to provide the hydration thermodynamic parameters for hydrated smectites with different compositions and under different conditions of temperature and relative humidity, using the newly developed predictive model.     

Geochemical chaos: Periodic and nonperiodic growth of mixed-layer phyllosilicates

Interstratification—periodic or nonperiodic stacking of two different silicate layers along a c^∗-axis—is common in phyllosilicates. Published evidence indicates that some interstratified minerals precipitate directly from aqueous solutions. In this paper, we have demonstrated, based on chaos theory, that both periodic and nonperiodic interstratification can autonomously arise from simple kinetics of mineral growth from a solution. Growth of a mixed-layer mineral is assumed to proceed layer by layer, and each layer starts with the formation of a base (Si, Al)-O tetrahedral sheet, whose structural configuration in a-b dimensions determines the type of new layer that forms. The sequence of layer stacking can be described by a one-dimensional map (i.e., a difference equation), which accounts for two competing factors: (1) the affinity of each end-member structural component for attaching to the surface of the preceding layer, and (2) the strain energy created by stacking next to each other two silicate layers with different structural configurations. Chaotic (or nonperiodic) interstratification emerges when the contacting solution becomes slightly supersaturated with respect to both structural components. The transition from one interstratification pattern to another reflects a change in chemical environment during mineral crystallization. Our model can successfully predict the occurrence of mixed-layer phyllosilicates and the associated layer stacking sequences observed in both hydrothermal alteration and sediment diagenesis. The model suggests that the diagenetic transition of smectite → nonperiodic illite/smectite → ordered illite/smectite → illite may reflect relative changes in the saturation degree of pore water with respect to two end-member phases as a result of increasing burial temperatures.     

A regular solution site-mixing model for illites

The dependence on the composition of the thermodynamic stability of an illite can be treated in terms of a regular solution site-mixing model. Four end-member micas (muscovite, pyrophyllite, phlogopite, and annite) were mixed to simulate an illite in this study. In the model, random mixing of cations was assumed over each given class of cation sites. Mixing over cation sites between different classes of cation sites was not allowed. The resulting free energy and chemical potential equations contain four site interaction parameters: three for octahedral site interactions and one for interlayer and tetrahedral site interactions. These parameters cannot presently be evaluated because of a lack of experimental data on Fe³⁺-free illites. The model does imply that the octahedral site interaction parameters must be significantly more positive than the interlayer and tetrahedral site interaction parameter to account for the dominant dioctahedral nature of most natural illites. This constraint is necessary to balance out the increase in stability due to the configurational entropy of an illite having a major trioctahedral component.The model can be extended to cover a wider range of illite compositions by the inclusion of an end-member mica containing Fe³⁺ ions in the octahedral sites. At present the thermodynamic properties of such an end-member are unknown.     

Effects of layer-charge distribution on the thermodynamic and microscopic properties of Cs-smectite

To explore the effects of layer-charge distribution on the thermodynamic and microscopic properties of Cs-smectites, classical molecular dynamic simulations are performed to derive the swelling curves, distributions and mobility of interlayer species, and Cs binding structures. Three representative smectites with distinct layer-charge distributions are used as model clay frameworks and interlayer water content is set within a wide range from 0 to 380 mg_water/g_clay. All the three smectites swell in a similar way, presenting the characteristic swelling plateaus and similar trends of swelling energetic profiles. The full-monolayer hydrate, corresponding to the global minima of the immersion energy, is the most stable hydrated state of Cs-smectites. The calculated diffusion coefficients of interlayer species disclose the confining effects in all smectites: both water molecules and ions diffuse slower than corresponding bulk cases and they are much more mobile in the direction parallel to the clay surfaces than perpendicular to them. The formed inner-sphere complex structures are very similar in different smectites: ions bind on the H-sites or T-sites and water molecules form cage-like caps covering the ions. Layer-charge distribution is found to have significant influences on the mobility of interlayer species and preference of ion binding sites. A general sequence is proposed to elucidate the preferences of various hexagonal sites (H-sites) and triangular sites (T-sites), that is, tetrahedrally substituted H-sites > nonsubstituted H-sites > tetrahedrally substituted T-sites > nonsubstituted T-sites, but the influence of octahedral substitutions on the preference of the neighboring sites is not obvious. Analysis of mobility indicates that H-sites are more stable Cs-fixation positions than T-sites and smectite with higher fraction of octahedral charges seems to be the most effective barrier material no matter how water content varies although all smectites can immobilize Cs ions in relatively dry conditions. These findings will not only facilitate basic research in geochemistry and material sciences, but also promote the barrier material designing.     

Physico/chemical stability of smectite clays

There is convincing evidence from field data that smectite clay undergoes conversion primarily to illite and chlorite if it is fully water-saturated and heated. The conversion may take place through mixed-layer formation with increasing illite/smectite ratio at higher temperatures and pressures. This process requires dehydration of the interlamellar space, for which either an external pressure or drying are needed. An alternative mechanism that takes place without dehydration, is dissolution of smectite and neoformation of illite. Both processes imply reorganization of the smectite crystal lattice for which the activation energy is fairly high, meaning that the conversion is negligible at temperatures lower than about 60°C. At elevated temperatures the conversion rate is controlled by the access to potassium for either mechanism.

An ongoing detailed investigation of this subject has led to a tentative model for the smectite-to-illite conversion in natural sediments and in canister-embedding clay in high-level radioactive waste (HLW) repositories.     

Experimental investigation of main controls to methane adsorption in clay-rich rocks

In this study a series of CH₄ adsorption experiments on clay-rich rocks were conducted at 35 °C, 50 °C and 65 °C and at CH₄ pressure up to 15 MPa under dry conditions. The clay-dominated rock samples used are fresh samples from quarries and mines. Samples are individually dominated by montmorillonite, kaolinite, illite, chlorite, and interstratified illite/smectite. The experimental results show that clay mineral type greatly affects CH₄ sorption capacity under the experimental conditions. In terms of relative CH₄ sorption capacity: montmorillonite ? illite/smectite mixed layer > kaolinite > chlorite > illite. Physisorption is the dominant process for CH₄ absorption on clay minerals, as a result, there is a linear correlation between CH₄ sorption capacity and BET surface area in these clay-mineral dominated rocks. The abundance of micro-mesopores in the size range of a few to a few 10 s of nanometers in montmorillonite clay and illite–smectite interstratified clay results in large BET surface area values for these mineral species.     

粘土矿物在油气初次运移中作用的探讨

本文在大量实际分析资料基础上,总结出我国含油气盆地中蒙脱石的转化类型,并讨论了它们在油气初次运移中的作用。在蒙脱石向伊利石连续转化型的盆地中,矿物转化对初次运移有明显作用;在蒙脱石向伊利石不连续转化或少量转化类型盆地中,由于矿物转化与生油期不相匹配或可转化的蒙脱石含量太少,对初次运移的作用不明显;在蒙脱石无转化和缺少蒙脱石的盆地中,基本不存在蒙脱石的转化过程,因而对初次运移不起作用。     

The effect of clay dehydration on land subsidence in the Yun-Lin coastal area, Taiwan

 The smectite dehydration theory developed by Ransom and Helgeson was applied for simulation of land subsidence in the Yun-Lin coastal area, Taiwan. The volumetric reduction of smectite clay at equilibrium state was computed by assuming that the dehydration of interlayer water in smectite clay can be described with a regular solid solution reaction. By using the in situ stratigraphic data collected from the subsidence monitoring wells in the simulated area, the amounts of land subsidence caused by smectite dehydration in three scenarios with pressure variation were calculated. The results indicate that significant amounts of land subsidence can be attributed to smectite dehydration. This finding reveals that smectite dehydration is of importance for assessment and prediction of land subsidence. Additionally, the results also indicate the overburden weight has a larger effect on clay dehydration than the effective stress change resulting from over-pumping, although both of them induce relatively minor variations on land subsidence. Received: 23 February 2000 · Accepted: 23 March 2000     

Regional retrograde alteration of sub-greenschist facies chlorite to smectite

Dioctahedral smectite is present as a retrograde alteration product of chlorite in Permian-Triassic red slates of the Malaguide Complex in Sierra de Espuña (Betic Cordillera). Mineral assemblages and textures, illite crystallinity indices, and fluid inclusion data indicate sub-greenschist facies conditions that reached at least 180°C in the higher-grade tectonic unit of the Malaguide Complex, preceding formation of smectite. Smectite, having K as the dominant interlayer cation, occurs ubiquitously intercalated with trioctahedral chlorite as thin packets of layers and as individual layers that commonly change to chlorite along layers. Although some chlorite is typically homogeneous and trioctahedral, much chlorite shows signs of alteration and has compositions corresponding to different degrees of smectite contaimination. The incompatibility of metamorphic grade with the occurrence of smectite, the general association of chlorite and smectite, and the textural relations collectively show that dioctahedral smectite is derived through replacement of trioctahedral chlorite. Such replacement occurs on a regional basis and demonstrates that caution must be used in interpreting the occurrence of smectite in pelites as being due to prograde processes. Alteration of trioctahedral chlorite under oxidizing conditions due to introduction of phreatic water after uplift of the Betic Cordillera is proposed as the cause of formation of smectite.     

Dating incipient metamorphism using 40Ar/39Ar geochronology and XRD modeling: a case study from the Swiss Alps

Six samples of a single carbonate-rich unit of the Swiss Préalpes, progressively metamorphosed from diagenesis to deep anchizone, yield ⁴⁰Ar/³⁹Ar spectra with variably developed staircase patterns, consistent with mixtures of detrital mica and neocrystallized mixed-layer illite/smectite. The lowest temperature heating steps for different size fractions (2–6?μm and 6–20?μm) converge to ～40?Ma providing an imprecise, maximum age of regional metamorphism. A method is described for distinguishing and quantifying the amount of pre-existing detrital mica versus neoformed illite layer in the illite/smectite formed during Tertiary Alpine metamorphism by comparison of X-ray diffraction patterns with Newmod^© simulations. In the least metamorphosed samples the illite/smectite contains ～65% neoformed illite, and this illite accounts for approximately 17% of all dioctahedral phyllosilicate minerals in the rock (e.g., detrital mica and illite/smectite). In contrast, the illite/smectite from the more strongly metamorphosed samples contains >97% neoformed illite, which accounts for ～70% to >90% of all dioctahedral phyllosilicate minerals. Phyllosilicate morphologies viewed by scanning electron microscopy are consistent with these estimates. A process of dissolution/reprecipitation is inferred as a mechanism for the growth of the neoformed phyllosilicates. A plot of neoformed illite content versus ⁴⁰Ar/³⁹Ar total fusion age yields a near-linear curve with an extrapolated age of 27?Ma for 100% neoformed dioctahedral phyllosilicates. This age is interpreted as the time of incipient metamorphism and is consistent with independent biostratigraphic constraints. Model ⁴⁰Ar/³⁹Ar age spectra constructed with the XRD simulation results correspond well to the experimental data and illustrate the changes in degassing properties of progressively metamorphosed mixtures of detrital mica and neoformed illite.     

Relation between the dissolution rates of single minerals and reservoir rocks in acidified pore waters

A series of kinetic experiments has been carried out to investigate the rates of dissolution (release of Al and Si) of common sandstone minerals in response to acidification of pore waters (pH = 3), using an experimental procedure designed to maximise the proportion of solid to fluid, and to minimise possible damage from agitation. The results have then been compared with those from experiments using disaggregated sandstones from two North Sea reservoirs. Experiments were carried out at 25 °C and 80 °C and in 0.01, 0.1 and 1 M NaCl solutions, with a pH of 3. Hydrochloric acid was used as the source of acidity and rate constants were determined based on both release of Al and Si. Mineral dissolution rates were closely comparable to literature values, despite the different experimental technique, except in the case of smectite where particle aggregation appears to have inhibited reaction. The dissolution rates calculated for reservoir sandstones based on their modal mineralogy and surface areas agree within a factor of 2 with the measured vales. Based on the reaction rates measured here, reservoir rocks rich in feldspar, illite and/or smectite are likely to react most rapidly with acidified pore waters.     

Spatial variation in fluid flow and geochemical fluxes across the sediment-seawater interface at the Carlos Ribeiro mud volcano (Gulf of Cadiz)

Submarine mud volcanism is an important pathway for transfer of deep-sourced fluids enriched in hydrocarbons and other elements into the ocean. Numerous mud volcanoes (MVs) have been discovered along oceanic plate margins, and integrated elemental fluxes are potentially significant for oceanic chemical budgets. Here, we present the first detailed study of the spatial variation in fluid and chemical fluxes at the Carlos Ribeiro MV in the Gulf of Cadiz. To this end, we combine analyses of the chemical composition of pore fluids with a 1-D transport-reaction model to quantify fluid fluxes, and fluxes of boron, lithium and methane, across the sediment-seawater interface. The pore fluids are significantly depleted in chloride, but enriched in lithium, boron and hydrocarbons, relative to seawater. Pore water profiles of sulphate, hydrogen sulphide and total alkalinity indicate that anaerobic oxidation of methane occurs at 34-180 cm depth below seafloor. Clay mineral dehydration, and in particular the transformation of smectite to illite, produces pore fluids that are depleted in chloride and potassium. Profiles of boron, lithium and potassium are closely related, which suggests that lithium and boron are released from the sediments during this transformation. Pore fluids are expelled into the water column by advection; fluid flow velocities are 4 cm yr⁻¹ at the apex of the MV but they rapidly decrease to 0.4 cm yr⁻¹ at the periphery. The associated fluxes of boron, lithium and methane vary between 7-301, 0.5-6 and 0-806 mmol m⁻² yr⁻¹, respectively. We demonstrate that fluxes of Li and B due to mud volcanism may be important on a global scale, however, release of methane into the overlying water column is suppressed by microbial methanotrophy.     

A short range interaction model applied to illite/smectite mixed-layer minerals

We present a solid solution model, able to take into account short-range interactions up to the third nearest neighbor. This model has been applied to illite/smectite mixed-layer clay minerals and tests the relative thermodynamic stability of each illite/smectite stacking sequence. The non-ideal energy terms have been calculated, showing a strong decrease of the interaction energy as the distance between elements grows. From ΔG _m calculations the model predicts the following succession in I/S stacking sequences as X _i tends to 1:R1, R2 and then R3. This succession agrees well with some mineralogical observations in natural series. Received: July 12, 1996 / Revised, accepted: March 21, 1997     

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.     

A method for activity calculations in saline and mixed solvent solutions at elevated temperature and pressure: A framework for geological phase equilibria calculations

Quantitative thermodynamic calculations that involve aqueous fluids have proved difficult because of the complexity of the interactions that occur within the fluids. Existing thermodynamic models are difficult to apply to mixed solvent or highly saline solutions at P > 0.3 GPa and T > 300 °C. This work constructs a method for activity-composition calculations in saline, mixed solvent, supercritical aqueous solutions. Mixing is formulated on a mole-fraction scale in terms of a set of independent end-members that describe composition and speciation within the solution. The ideal mixing term takes speciation into account and avoids problems with the common ion effect. Non-ideal interactions are represented by an activity coefficient term that combines a limited form of Debye-Hückel and a van Laar formulation. This approach, referred to as the DH-ASF model, is thermodynamically valid over a wide range of P, T and fluid composition. The value of the model lies in its broad applicability, and small number of calibration parameters. Experimental data from the literature for the systems NaCl-H₂O, KCl-H₂O, H₂O-SiO₂-CO₂, H₂O-NaCl-CO₂, H₂O-NaCl-SiO₂ and for H₂O-albite melts have been used to calibrate the DH-ASF model. Calculations were performed using Thermocalc, computer software that calculates equilibria for mineral-based chemical systems.¹ The model represents the data to within experimental error in most cases. Conditions modelled include pressures between 0.2 and 1.4 GPa, temperatures between 500 and 900 °C, and x_H2O from 0.1 to 1. Calibrated parameters are consistent with expectations based on the conceptual model for the fluid, and are relatively insensitive to changes in pressure and temperature for most examples. The DH-ASF model is thermodynamically valid for a range of P-T conditions that includes pressures from 0.1 to 2 GPa and temperatures from 200 to 1000 °C. A lack of experimental data restricts calibration of the model for many end-members. However, it may be possible to neglect parameters associated with end-members present in small amount. In this case, or with new experimental data for calibrations, the DH-ASF model allows previously inaccessible geological systems and processes to be modelled.     

Mechanisms of siderophore sorption to smectite and siderophore-enhanced release of structural Fe

Sorption of the trihydroxamate siderophores desferrioxamine-B and -D (DFOB and DFOD, respectively) and of the monohydroxamate ligand acetohydroxamic acid (aHA) to smectite were examined in batch sorption studies (pH 5.5, 0.1 M ionic strength) coupled with X-ray diffraction (XRD) and Fourier transform infrared (FTIR) spectroscopy. Both DFOB and DFOD, which have similar molecular structures but different charge properties (cationic versus neutral, respectively) showed a high affinity for smectite. In contrast, the smaller aHA molecule did not sorb appreciably. XRD analysis indicated that DFOB and DFOD each absorbed in the interlamellar region of the clay to give d-spacings of 13.4-13.7 Å at equilibrium solution concentrations <250 μM. FTIR spectra of sorbed DFOB and DFOD indicated that the conformation of each species was distinct from its conformation in the crystalline or dissolved states. At elevated initial solution concentrations of 500-1500 μM, DFOB formed a bilayer in the clay interlayer. Changes in the FTIR spectra of the DFOB-loaded clay samples at these higher surface loadings were consistent with the presence of a metal-siderophore complex in the interlayer. DFOB and DFOD both enhanced Fe and Al release from smectite, but aHA did not. Possible dissolution mechanisms are discussed in light of the FTIR and batch dissolution results.