Bathymetric and geochemical investigation of Kawah Ijen Crater Lake, East Java, Indonesia

A bathymetric survey of Kawah Ijen crater lake was conducted by acoustic sounding in 1996 to compare the lake morphology with those measured in 1922, 1925 and 1938, and to calculate the present lake volume. Even though the lake experienced several hydrothermal eruptions, the maximum depth became shallower (182 m) than before (200 m), resulting in a reduced lake volume (3.0×10⁷ m³).Fifty-two major and minor constituents including rare earth elements and polythionates (PT) of the lake waters at various depths were determined by ICP-AES, ICP-MS and HPLC, respectively. These ions except for several volatile elements are taken up by lake fringe through congruent dissolution of pyroclastics of Kawah Ijen volcano. Most ions are homogeneously distributed throughout the lake, although PT showed a considerable vertical variation. Rare earth elements (REE) in the Kawah Ijen water as well as those from other hyper-acidic crater lakes show distribution patterns likely due to the three rock dissolution (preferential, congruent and residual) types, and their logarithmic concentrations linearly depend upon the pH values of the lake waters.Using the PT degradation kinetics data, production rates of PT, injection rates of SO₂ and H₂S into the lake were estimated to be 114, 86 and 30 tons/day, respectively. Also travel time of the spring water at the Banyupahit Riverhead from Kawah Ijen was estimated to be 600–1000 days through the consideration of decreasing rates of PT. Molten sulfur stocks containing Sn, Cu, Bi sulfides and Pb-barite exposed on the inner crater slope were presumed to be extinct molten sulfur pools at the former lake bottom. This was strongly supported by the barite precipitation temperature estimated through the consideration of the temperature dependence of Pb-chlorocomplex formation.     

Multicomponent ionic diffusion in porewaters: Coulombic effects revisited

The diffusion of an ion in porewaters cannot occur independently of the other ions in solution as a result of Coulombic coupling, as well as from other effects not considered here. Unfortunately, a longstanding disagreement exists about the correct form and meaning of the equations that describe Coulombic coupling in porewaters, i.e., Ben-Yaakov [Am. J. Sci. 281 (1981) 974] vs. Lasaga [Am. J. Sci. 281 (1981) 981]. This paper re-examines this controversy by reformulating the problem starting from fundamental concepts of mass and charge conservation. We show that these antagonistic formulations are both valid and, in fact, equivalent, when the different interpretations of charge balance are resolved. Most of the disagreements between Ben-Yaakov and Lasaga are then shown to result from differing methods of solution, not fundamental disparities in their models. We note, however, that the explanation for the concept of “stationary” gradients of nonreacting ions as given Ben-Yaakov is inaccurate, and such gradients do lead to diffusive fluxes that are counterbalanced by electrochemical migrational fluxes to produce no net flux (excluding advective flux). We further find that the bicarbonate diffusive flux will not balance the diffusional charge flux of sulfate during its reduction if advection is present. This latter imbalance generates compensating fluxes in the other nonreacting ions. We have applied our theory to a simplified case of sulfate reduction in a marine sediment. The results show that nonreacting ions do diffuse and that with normally expected values of porewater advection, the ratio of the bicarbonate to the sulfate flux can be far different than the ideal value of −2.     

The Okavango giant mafic dyke swarm (NE Botswana): its structural significance within the Karoo Large Igneous Province

The structural organization of a giant mafic dyke swarm, the Okavango complex, in the northern Karoo Large Igneous Province (LIP) of NE Botswana is detailed. This N110°E-oriented dyke swarm extends for 1500 km with a maximum width of 100 km through Archaean basement terranes and Permo-Jurassic sedimentary sequences. The cornerstone of the study is the quantitative analysis of N>170 (exposed) and N>420 (detected by ground magnetics) dykes evidenced on a ca. 80-km-long section lying in crystalline host-rocks, at high-angle to the densest zone of the swarm (Shashe area). Individual dykes are generally sub-vertical and parallel to the entire swarm. Statistical analysis of width data indicates anomalous dyke frequency (few data <5.0 m) and mean dyke thickness (high value of 17 m) with respect to values classically obtained from other giant swarms. Variations of mean dyke thicknesses from 17 (N110°E swarm) to 27 m (adjoining and coeval N70°E giant swarm) are assigned to the conditions hosting fracture networks dilated as either shear or pure extensional structures, respectively, in response to an inferred NNW–SSE extension. Both fracture patterns are regarded as inherited brittle basement fabrics associated with a previous (Proterozoic) dyking event. The Okavango N110°E dyke swarm is thus a polyphase intrusive system in which total dilation caused by Karoo dykes (estimated frequency of 87%) is 12.2% (6315 m of cumulative dyke width) throughout the 52-km-long projected Shashe section. Assuming that Karoo mafic dyke swarms in NE Botswana follow inherited Proterozoic fractures, as similarly applied for most of the nearly synchronous giant dyke complexes converging towards the Nuanetsi area, leads us to consider that the resulting triple junction-like dyke/fracture pattern is not a definitive proof for a deep mantle plume in the Karoo LIP.     

Quantum-thermodynamic treatment of intrinsic anharmonicity; Wallace’s theorem revisited

Wallace (in Thermodynamics of crystals, 1972) developed a theorem, rooted in rigid lattice dynamics, which incorporates intrinsic anharmonic effects in solids. The practical application of this theorem in mineral physics is computationally involved and this is the main reason for the theorem not getting the attention it deserves. Because intrinsic anharmonicity is an important issue at the extreme conditions in planetary mantles, we derived a method which removes the computational obstacles in applying this theorem. We extended the theorem to incorporate details of the phonon spectrum and tested our algorithm on forsterite (Mg₂SiO₄). Using a least squares inversion technique applied to all available experimental data, we show that it results in an accurate representation of thermodynamic properties and sound wave velocities of Mg₂SiO₄ in its complete pressure–temperature stability range. We also show that the accuracy of our results is not significantly affected by the use of a different equation of state.     

Effect of Fe on garnet–melt trace element partitioning: experiments in FCMAS and quantification of crystal-chemical controls in natural systems

Garnet–melt trace element partitioning experiments were performed in the system FeO–CaO–MgO–Al₂O₃–SiO₂ (FCMAS) at 3 GPa and 1540°C, aimed specifically at studying the effect of garnet Fe²⁺ content on partition coefficients (D^Grt/Melt). D^Grt/Melt, measured by SIMS, for trivalent elements entering the garnet X-site show a small but significant dependence on garnet almandine content. This dependence is rationalised using the lattice strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal–melt partition coefficients from elastic moduli. Nature 372, 452–454], which describes partitioning of an element i with radius r_i and valency Z in terms of three parameters: the effective radius of the site r₀(Z), the strain-free partition coefficient D₀(Z) for a cation with radius r₀(Z), and the apparent compressibility of the garnet X-site given by its Young's modulus E_X(Z). Combination of these results with data in Fe-free systems [Van Westrenen, W., Blundy, J.D., Wood, B.J., 1999. Crystal-chemical controls on trace element partitioning between garnet and anhydrous silicate melt. Am. Mineral. 84, 838–847] and crystal structure data for spessartine, andradite, and uvarovite, leads to the following equations for r₀(3+) and E_X(3+) as a function of garnet composition (X) and pressure (P):

r₀(3+) [Å]=0.930X_Py+0.993X_Gr+0.916X_Alm+0.946X_Spes+1.05(X_And+X_Uv)−0.005(P [GPa]−3.0)(±0.005 Å)

E_X(3+) [GPa]=3.5×10¹²(1.38+r₀(3+) [Å])^−26.7(±30 GPa)

Accuracy of these equations is shown by application to the existing garnet–melt partitioning database, covering a wide range of P and T conditions (1.8 GPa<P<5.0 GPa; 975°C<T<1640°C). D^Grt/Melt for all 3+ elements entering the X-site (REE, Sc and Y) are predicted to within 10–40% at given P, T, and X, when D^Grt/Melt for just one of these elements is known. In the absence of such knowledge, relative element fractionation (e.g. D_Sm^Grt/Melt/D_Nd^Grt/Melt) can be predicted. As an example, we predict that during partial melting of garnet peridotite, group A eclogite, and garnet pyroxenite, r₀(3+) for garnets ranges from 0.939±0.005 to 0.953±0.009 Å. These values are consistently smaller than the ionic radius of the heaviest REE, Lu. The above equations quantify the crystal-chemical controls on garnet–melt partitioning for the REE, Y and Sc. As such, they represent a major advance en route to predicting D^Grt/Melt for these elements as a function of P, T and X.     

The effect of Ca-Tschermaks component on trace element partitioning between clinopyroxene and silicate melt

We have studied the influence of Ca-Tschermaks (Calcium Tschermaks or CaTs) content of clinopyroxene on the partitioning of trace elements between this phase and silicate melt at fixed temperature and pressure. Ion probe analyses of experiments carried out in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂, at 0.1 MPa and 1218°C, produced crystal-melt partition coefficients (D) of 36 trace elements (Li, Cl, Sc, Ti, V, Cr, Fe, Co, Ge, Sr, Y, Zr, Nb, Mo, Ru, Rh, In, Sn, Sb, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta and W), for clinopyroxene compositions between 10 and 32 mol% CaTs. Partition coefficients for 2+ to 5+ cations show, for each charge, a near parabolic dependence of log D on ionic radius of the substituting cation, for partitioning into both the M1 and M2 sites of clinopyroxene. Fitting the results to the elastic strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal-melt partition coefficients from elastic moduli. Nature 372, 452–454] we obtain results for the strain-free partition coefficients of theoretical cations (D₀), with site radius r₀, and for the site's Young's Modulus (E).

In agreement with earlier data our results show that increasing ^ivAl concentration in cpx is matched by increasing D, E_M1, E_M2 and D₀ for tri-, tetra- and pentavalent cations. The degree of fractionation between chemically similar elements (i.e. Ta/Nb, Zr/Hf) also increases. In contrast, D values for mono-, di- and hexavalent cations decrease with increasing ^ivAl in the cpx. The large suite of trace elements used has allowed us to study the effects of cation charge on D₀, r₀ and E. We have found that D₀ and r₀ decrease with increasing cation charge, e.g. r₀=0.66 Å for 4+ cations and 0.59 Å for 5+ cations substituting into M1. Values of E_M1 and E_M2 increase with cation charge as well as with increasing ^ivAl content. The increase in E_M2 is linear and close to the trend set by Hazen and Finger [Hazen, R.M., Finger, L.W., 1979. Bulk modulus-volume relationship for cation–anion polyhedra. J. Geophys. Res. 84 (10) 6723–6728] for oxides. E_M1 values are much higher and do not fit the trend predicted by the Hazen and Finger relationship.     

Micro- and nano-textural evidence of Ti(–Ca–Fe) mobility during fluid–rock interactions in carbonaceous lawsonite-bearing rocks from New Zealand

Understanding the mobility of chemical elements during fluid–rock interactions is critical to assess the geochemical evolution of a rock undergoing burial and metamorphism and, more generally, to constrain the geochemical budget of the subduction factory. In particular, determining the behavior and mobility of Ti in aqueous fluids constitutes a great challenge that is still under scrutiny. Here, we study plant fossils preserved in blueschist metasedimentary rocks from the Marybank Formation (New Zealand). Using scanning and transmission electron microscopies (SEM and TEM), we show that the carbonaceous material (CM) composing the fossils contains abundant nano-inclusions of Ti- and Fe-oxides. These nanocrystals are mainly anatase, rutile, and Fe–Ti oxides. The mineral composition observed within the fossils is significantly different from that detected in the surrounding rock matrix. We propose that Ti and Fe might have been mobilized by the alteration of a detrital Ti–Fe-rich protolith during an early diagenetic event under acidic and reducing conditions. Aqueous fluids rich in organic ligands released by the degradation of organic matter may have been involved. Moreover, using mass balance and petrological observations, we show that the contrasted mineralogy between the rock matrix and the fossil CM might be the consequence of the chemical isolation of fossil CM during the prograde path of the rock. Such an isolation results from the early formation of quartz and Fe-rich phyllosilicate layers enclosing the fossil as characterized by SEM and TEM investigations. Overall, this study shows that investigating minerals associated with CM down to the nanometer scale in metamorphic rocks can provide a precious record of early prograde geochemical conditions.