Oxygen isotope heterogeneity of the mantle beneath the Canary Islands: a discussion of the paper of Gurenko et al.

Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) report laser-assisted fluorination (LF) and secondary ionization mass spectrometry (SIMS) ¹⁸O/¹⁶O datasets for olivine grains from the Canary Islands of Gran Canaria, Tenerife, La Gomera, La Palma and El Hierro. As with prior studies of oxygen isotopes in Canary Island lavas (e.g. Thirlwall et al. Chem Geol 135:233–262, 1997; Day et al. Geology 37:555–558, 2009, Geochim Cosmochim Acta 74:6565–6589, 2010), these authors find variations in δ¹⁸O_ol (~4.6–6.0 ‰) beyond that measured for mantle peridotite olivine (Mattey et al. Earth Planet Sci Lett 128:231–241, 1994) and interpret this variation to reflect contributions from pyroxenite-peridotite mantle sources. Furthermore, Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) speculate that δ¹⁸O_ol values for La Palma olivine grains measured by LF (Day et al. Geology 37:555–558, 2009, Geochim Cosmochim Acta 74:6565–6589, 2010) may be biased to low values due to the presence of altered silicate, possibly serpentine. The range in δ¹⁸O_ol values for Canary Island lavas are of importance for constraining their origin. Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) took a subset (39 SIMS analyses from 13 grains from a single El Hierro lava; EH4) of a more extensive dataset (321 SIMS analyses from 110 grains from 16 Canary Island lavas) to suggest that δ¹⁸O_ol is weakly correlated (R ² = 0.291) with the parameter used by Gurenko et al. (Earth Planet Sci Lett 277:514–524, 2009) to describe the estimated weight fraction of pyroxenite-derived melt (Xpx). With this relationship, end-member δ¹⁸O values for HIMU-peridotite (δ¹⁸O = 5.3 ± 0.3 ‰) and depleted pyroxenite (δ¹⁸O = 5.9 ± 0.3 ‰) were defined. Although the model proposed by Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) implicates similar pyroxenite-peridotite mantle sources to those proposed by Day et al. (Geology 37:555–558, 2009, Geochim Cosmochim Acta 74:6565–6589, 2010) and Day and Hilton (Earth Planet Sci Lett 305:226–234, 2011), there are significant differences in the predicted δ¹⁸O values of end member components in the two models. In particular, Day et al. (Geochim Cosmochim Acta 74:6565–6589, 2010) proposed a mantle source for La Palma lavas with low-δ¹⁸O (<5 ‰), rather than higher-δ¹⁸O (c.f. the HIMU-peridotite composition of Gurenko et al. in Contrib Mineral Petrol 162:349–363, 2011). Here we question the approach of using weakly correlated variations in δ¹⁸O_ol and the Xpx parameter to define mantle source oxygen isotope compositions, and provide examples of why this approach appears flawed. We also provide reasons why the LF datasets previously published for Canary Island lavas remain robust and discuss why LF and SIMS data may provide complementary information on oxygen isotope variations in ocean island basalts (OIB), despite unresolved small-scale uncertainties associated with both techniques.     

From peat bog to lignite seam: a new method to calculate the consolidation coefficient of lignite seams,Wielkopolska region in central Poland

Only one quantitative method is described in the literature to estimate the consolidation coefficient of lignite seams. More studies concentrate on the compaction of plant tissues, which explains why the obtained results are overestimated. Moreover, most studies do commonly not determine the consolidation of the whole peat bog but usually those of some of its elements only. We propose a new approach, which conceptually is fairly close to the Hager et al. (Fortschr Geol Rheinld Westf 29:319–352, 1981) method. Our method has been tested on an example of the first Middle-Polish Lignite Seam and the second Lusatian Lignite Seam from some Wielkopolska lignite deposits in central Poland. The consolidation coefficients, C _c, range between 2.34 and 2.56 for the second Lusatian Lignite Seam and between 1.80 and 2.14 for the first Middle-Polish Lignite Seam.     

Ca–Mg diffusion in diopside: tracer and chemical inter-diffusion coefficients

We have experimentally determined the tracer diffusion coefficients (D*) of ⁴⁴Ca and ²⁶Mg in a natural diopside (~Di₉₆) as function of crystallographic direction and temperature in the range of 950–1,150 °C at 1 bar and f(O₂) corresponding to those of the WI buffer. The experimental data parallel to the a*, b, and c crystallographic directions show significant diffusion anisotropy in the a–c and b–c planes, with the fastest diffusion being parallel to the c axis. With the exception of logD*(²⁶Mg) parallel to the a* axis, the experimental data conform to the empirical diffusion “compensation relation”, converging to logD ~ −19.3 m²/s and T ~ 1,155 °C. Our data do not show any change of diffusion mechanism within the temperature range of the experiments. Assuming that D* varies roughly linearly as a function of angle with respect to the c axis in the a–c plane, at least within a limited domain of ~20° from the c-axis, our data do not suggest any significant difference between D*(//c) and D*(⊥(001)), the latter being the diffusion data required to model compositional zoning in the (001) augite exsolution lamellae in natural clinopyroxenes. Since the thermodynamic mixing property of Ca and Mg is highly nonideal, calculation of chemical diffusion coefficient of Ca and Mg must take into account the effect of thermodynamic factor (TF) on diffusion coefficient. We calculate the dependence of the TF and the chemical interdiffusion coefficient, D(Ca–Mg), on composition in the diopside–clinoenstatite mixture, using the available data on mixing property in this binary system. Our D*(Ca) values parallel to the c axis are about 1–1.5 log units larger than those Dimanov et al. (1996). Incorporating the effect of TF, the D(Ca–Mg) values calculated from our data at 1,100–1,200 °C is ~0.6–0.7 log unit greater than the experimental quasibinary D((Ca–Mg + Fe)) data of Fujino et al. (1990) at 1 bar, and ~0.6 log unit smaller than that of Brady and McCallister (1983) at 25 kb, 1,150 °C, if our data are normalized to 25 kb using activation volume (~4 and ~6 cm³/mol for Mg and Ca diffusion, respectively) calculated from theoretical considerations.     

Relation Between Normalized Axial Stress and Failure Strain in Heterogeneous Carbonate Rocks Exhibiting Large Axial Strains

A semi-analytical equation for the modeling of stress–strain relationship for heterogeneous carbonate rocks exhibiting large axial strains (ε_af > 1%) is formulated. The equation is derived by modifying the stress–strain model based on Haldane’s distribution proposed by Palchik (2006) for carbonate rocks exhibiting ε _af ≤ 1%. The developed exponential model is used to relate normalized axial stress (σ _a/σ _c) over the whole pre-failure strain range to current axial strain (ε _a) and failure strain (ε _af). For carbonate rocks exhibiting ε _af > 1%, the value of pre-calculated parameter δ involved in the stress–strain model is not constant, but dependent on the failure strain value (ε _af). The normalized stress–strain model can be used to calculate the failure strain in terms of uniaxial compressive strength and stress–strain measurement at one point only. The advantages of the failure strain model and ways of its use in engineering practice are discussed.     

On the Use of Skempton’s Compression Index Equation

In this study, an alternate approach to establish the e-log p relationships for clayey soils within a vertical pressure range of 10–1,000 kPa is discussed. Skempton’s compression index equation correlating the liquid limit, w _L, and the compression index, C _c, and the reported equation correlating the void ratio at liquid limit, e _L, and the void ratio at a vertical pressure of 100 kPa, e ₁₀₀, by Burland (1990), were used to establish the e-log p relationships for several reconstituted normally consolidated clayey soils. Consolidation test results of 13 clayey soils covering a sufficiently wide range of liquid limit were selected from the literature. Also, consolidation tests were carried out on two highly expansive soils in this study. A comparison of the experimental consolidation test results with the calculated e-log p relationships in the current study indicated that in general, the agreements between the calculated relationships and the experimental results are good. The agreements were found to be slightly better for soils with liquid limits less than about 70%. A comparison of the calculated e-log p relationships in the current study with that determined following methods suggested by Nagaraj and Srinivasa Murthy (1983) and Burland (1990) showed that all the three methods yielded very similar results for soils with liquid limit less than 70%. For soils with liquid limits greater than 70%, the difference between the e-log p relationships calculated in this study and that following Burland (1990)’s method was insignificant, whereas Nagaraj and Srinivasa Murthy (1983)’s method slightly over-predicted the void ratios at larger vertical stresses.     

What’s Hiding Behind Ontogenetic δ<Superscript>13</Superscript>C Variations in Mollusk Shells? New Insights from the Great Scallop (<Emphasis Type="Italic">Pecten maximus</Emphasis>)

Mollusk shells contain geochemical information about environmental conditions that prevailed at the time of formation. We investigated ontogenetic and seasonal variations of δ¹³C in calcitic shells of Pecten maximus. Ontogenetic variations of δ¹³C_shell in three large specimens collected in Norway, France, and Spain exhibited a similar linear decrease with increasing shell height. We removed this linear drift (detrending). These three residual time series displayed variations that could be linked to environmental fluctuations. To check it, we reanalyzed the isotopic datasets of Lorrain et al. (Journal of Experimental Marine Biology and Ecology 275:47–61, 2002, Geochimica et Cosmochimica Acta 68:3509–3519, 2004), who worked on three scallops harvested in 2000 in the bay of Brest (France), a well-monitored ecosystem. Lowest values of δ¹³C_{shell detrended} were recorded in all shells in late spring–early summer, most likely reflecting corresponding variations in food availability. Our results indicate that ontogenetic and seasonal variations of δ¹³C_shell cannot be used as a proxy for past δ¹³C_DIC variations but should be considered as promising tools for ecophysiological studies.     

Hydraulic Conductivity and Coefficient of Consolidation of Two Sulphide Clays in Sweden

The variation and anisotropy in hydraulic conductivity and the coefficient of consolidation was investigated for two Swedish sulphide clays. A series of constant rate of strain oedometer tests was performed on samples trimmed in the vertical and horizontal direction. A methodology to evaluate the horizontal coefficients of consolidation c _h via the horizontal hydraulic conductivity k _h and the vertical compression modulus M _v is proposed. Laboratory evaluations of c _h are also compared with determinations of c _h from in situ piezometer measurements in vertically drained sulphide clay. Furthermore, the validity of the empirical correlation between hydraulic conductivity change index C _k and initial void ratio e ₀, C _k  = 0.5e ₀ (Tavenas et al. in Can Geotech J 20(4):645–660, 1983b), was investigated for the sulphide clays. The results from the investigation show large ranges in measured hydraulic conductivities and coefficients of consolidation. However, the results indicate that the correlation C _k  = 0.5e ₀ is valid. The anisotropy in hydraulic conductivity and the coefficient of consolidation of the sulphide clays tested seems to be small. For design purposes, multiple tests for assessment of hydraulic conductivity and the coefficient of consolidation should be made, and a partial factor of safety, depending on the requisite level of safety and the spatial variability of the parameters, should be introduced. For design purposes in this type of clay, k _h  = k _v and c _h  = c _v are suggested.     

Elasticity and equation of state of Li<Subscript>2</Subscript>B<Subscript>4</Subscript>O<Subscript>7</Subscript>

A single crystal X-ray diffraction study on lithium tetraborate Li₂B₄O₇ (diomignite, space group I4₁ cd) has been performed under pressure up to 8.3 GPa. No phase transitions were found in the pressure range investigated, and hence the pressure evolution of the unit-cell volume of the I4₁ cd structure has been described using a third-order Birch–Murnaghan equation of state (BM-EoS) with the following parameters: V ₀  = 923.21(6) ų, K ₀  = 45.6(6) GPa, and K′ = 7.3(3). A linearized BM-EoS was fitted to the axial compressibilities resulting in the following parameters a ₀  = 9.4747(3) Å, K _0a  = 73.3(9) GPa, K′ _a  = 5.1(3) and c ₀  = 10.2838(4) Å, K _0c  = 24.6(3) GPa, K′ _c  = 7.5(2) for the a and c axes, respectively. The elastic anisotropy of Li₂B₄O₇ is very large with the zero-pressure compressibility ratio β _0c /β _0a  = 3.0(1). The large elastic anisotropy is consistent with the crystal structure: A three-dimensional arrangement of relatively rigid tetraborate groups [B₄O₇]²⁻ forms channels occupied by lithium along the polar c–axis, and hence compression along the c axis requires the shrinkage of the lithium channels, whereas compression in the a direction depends mainly on the contraction of the most rigid [B₄O₇]²⁻ units. Finally, the isothermal bulk modulus obtained in this work is in general agreement with that derived from ultrasonic (Adachi et al. in Proceedings-IEEE Ultrasonic Symposium, 228–232, 1985; Shorrocks et al. in Proceedings-IEEE Ultrasonic Symposium, 337–340, 1981) and Brillouin scattering measurements (Takagi et al. in Ferroelectrics, 137:337–342, 1992).     

Reply to “Oxygen isotope heterogeneity of the mantle beneath the Canary Islands: a discussion of the paper of Gurenko et al.”

The comment by Day et al. (Contrib Mineral Petrol, 2012) (1) discusses the validity of the previously obtained oxygen isotope data for El Hierro and La Palma (Canary Island) olivines, (2) questions the approach by Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) of using weakly correlated variations of δ¹⁸O_olivine values with X _px (proportion of pyroxenite-derived melt in the parental magma), and (3) provides reasons why oxygen isotope data by secondary ion mass spectrometry (SIMS) “offer sensitive means for detecting melt-crust interactions.” We respond these comments and report a new set of oxygen isotope measurements performed by SIMS and single-grain laser fluorination methods. These measurements confirm our previous data and conclusions and demonstrate the ability of the SIMS technique to analyze O isotopes in terrestrial samples with 2-sigma uncertainty better than ±0.25 ‰.     

Prediction of the thermodynamic properties of metal–arsenate and metal–arsenite aqueous complexes to high temperatures and pressures and some geological consequences

The standard thermodynamic properties at 25°C, 1 bar (ΔG _f^o, ΔH _f^o, S ^o, C _P^o, V ^o, ω) and the coefficients of the revised Helgeson–Kirkham–Flowers equations of state were evaluated for several aqueous complexes formed by dissolved metals and either arsenate or arsenite ions. The guidelines of Shock and Helgeson (Geochim Cosmochim Acta 52:2009–2036, 1988) and Sverjensky et al. (Geochim Cosmochim Acta 61:1359–1412, 1997) were followed and corroborated with alternative approaches, whenever possible. The SUPCRT92 computer code was used to generate the log K of the destruction reactions of these metal–arsenate and metal–arsenite aqueous complexes at pressures and temperatures required by the EQ3/6 software package, version 7.2b. Apart from the AlAsO₄^o and FeAsO₄^o complexes, our log K at 25°C, 1 bar are in fair agreement with those of Whiting (MS Thesis, Colorado School of Mines, Golden, CO, 1992). Moreover, the equilibrium constants evaluated in this study are in good to fair agreement with those determined experimentally for the Ca–dihydroarsenate and Ca–hydroarsenate complexes at 40°C (Mironov et al., Russ J Inorg Chem 40:1690, 1995) and for Fe(III)–hydroarsenate complex at 25°C (Raposo et al., J Sol Chem 35:79–94, 2006), whereas the disagreement with the log K measured for the Ca–arsenate complex at 40°C (Mironov et al., Russ J Inorg Chem 40:1690, 1995) might be due to uncertainties in this measured value. The implications of aqueous complexing between dissolved metals and arsenate/arsenite ions were investigated for seawater, high-temperature geothermal liquids and acid mine drainage and aqueous solutions deriving from mixing of acid mine waters and surface waters. Electronic Supplementary Material The online version of this article () contains supplementary material, which is available to authorized users.     

On the Ratios between Elastic Modulus and Uniaxial Compressive Strength of Heterogeneous Carbonate Rocks

The ratios M _R = E/σ _c for 11 heterogeneous carbonate (dolomites, limestones and chalks) rock formations collected from different regions of Israel were examined. Sixty-eight uniaxial compressive tests were conducted on weak-to-strong (5 MPa < σ _c < 100 MPa) and very strong (σ _c > 100 MPa) rock samples exhibiting wide ranges of elastic modulus (E = 6100–82300 MPa), uniaxial compressive strength (σ _c = 14–273.9 MPa), Poisson's ratio (ν = 0.13–0.49), and dry bulk density (ρ = 1.7–2.7 g/cm³). The observed range of M _R = 60.9–1011.4 and mean value of M _R = 380.5 are compared with the results obtained by Deere (Rock mechanics in engineering practice, Wiley, London, pp 1–20, 1968) for limestones and dolomites, and the statistical analysis of M _R distribution is performed. Mutual relations between E, σ _c, ρ, M _R for all studied rocks, and separately for concrete rock formations are revealed. Linear multiple correlations between E on the one hand and σ _c and ρ on the other for Nekorot and Bina limestone and Aminadav dolomite are obtained. It is established that the elastic modulus and M _R in very strong carbonate samples are more correlated with ρ−σ _c combination and ε _a max, respectively, than in weak to strong samples. The relation between M _R and maximum axial strain (ε _a max) for all studied rock samples (weak-to-strong and very strong) is discussed.     

The Gumeshevo skarn-porphyry copper deposits in the Central Urals,Russia: Evolution of the ore-magmatic system as deduced from isotope geochemistry (Sr,Nd, C,O, H)

The Early Devonian Gumeshevo deposit is one of the largest ore objects pertaining to the dioritic model of the porphyry copper system paragenetically related to the low-K quartz diorite island-arc complex. The (⁸⁷Sr/⁸⁶Sr)_t and (ɛ_Nd)_t of quartz diorite calculated for t = 390 Ma are 0.7038–0.7045 and 5.0–5.1, respectively, testifying to a large contribution of the mantle component to the composition of this rock. The contents of typomorphic trace elements (ppm) are as follows: 30–48 REE sum, 5–10 Rb, 9–15 Y, and 1–2 Nb. The REE pattern is devoid of Eu anomaly. Endoskarn of low-temperature and highly oxidized amphibole-epidote-garnet facies is surrounded by the outer epidosite zone. Widespread retrograde metasomatism is expressed in replacement of exoskarn and marble with silicate (chlorite, talc, tremolite)-magnetite-quartz-carbonate mineral assemblage. The ⁸⁷Sr/⁸⁶Sr ratios of epidote in endoskarn and carbonate in retrograde metasomatic rocks (0.7054–0.7058 and 0.7053–0.7065, respectively) are intermediate between the Sr isotope ratios of quartz dioritic rocks and marble (⁸⁷Sr/⁸⁶Sr = 0.70784 ± 2). Isotopic parameters of the fluid equilibrated with silicates of skarn and retrograde metasomatic rocks replacing exoskarn at 400°C are δ¹⁸O = +7.4 to +8.5‰ and δD = −49 to −61‰ (relative to SMOW). The δ¹³C and δ¹⁸O of carbonates in retrograde metasomatic rocks after marble are −5.3 to +0.6 (relative to PDB) and +13.0 to +20.2% (relative to SMOW), respectively. Sulfidation completes metasomatism, nonuniformly superimposed on all metasomatic rocks and marbles with formation of orebodies, including massive sulfide ore. The δ³⁴S of sulfides is 0 to 2‰ (relative to CDT);⁸⁷Sr/⁸⁶Sr of calcite from the late calcite-pyrite assemblage replacing marble is 0.704134 ± 6. The δ¹³C and ⁸⁷Sr/⁸⁶Sr of postore veined carbonates correlate positively (r = 0.98; n = 6). The regression line extends to the marble field. Its opposite end corresponds to magmatic (in terms of Bowman, 1998b) calcite with minimal δ¹³C, δ¹⁸O, and ⁸⁷Sr/⁸⁶Sr values (−6.9 ‰, +6.7‰, and 0.70378 ± 4, respectively). The aforementioned isotopic data show that magmatic fluid was supplied during all stages of mineral formation and interacted with marble and other rocks, changing its Sr, C, and O isotopic compositions. This confirms the earlier established redistribution of major elements and REE in the process of metasomatism. A contribution of meteoric and metamorphic water is often established in quartz from postore veins.     

Habitat Use by Coastal Birds Inferred from Stable Carbon and Nitrogen Isotopes

Conservationists need to know the degree of habitat fidelity for species of conservation concern. Stable Isotope Analysis in R quantified the contribution of terrestrial vs. saltmarsh primary production sources to terrestrial passerine food webs from four habitats of Sapelo Island, Georgia (USA), saltmarsh, maritime scrub–shrub, maritime broadleaf (oak), and maritime narrowleaf (pine) forests, using δ ¹³C and δ ¹⁵N. Models suggested Northern Parula (Parula americana) in oak forests, White-eyed Vireos (Vireo griseus) in shrub, and Brown-headed Nuthatches (Sitta pusilla) in pine forests derived most of their food from habitats they occupied (53–100%). Saltmarsh provided 47–94% of Painted Bunting (Passerina ciris) food sources, supporting previous findings by Springborn and Meyers (2005). Thus, Painted Bunting conservation in the Southeastern USA should focus on Springborn and Meyers’ suggestion of maritime scrub–shrub habitat and forests with <75% canopy, >50% ground cover, and patches of shrubs that are within 700 m of saltmarsh.     

Heat capacity,entropy and phase equilibria of stishovite

The low-temperature heat capacity (C _P) of stishovite (SiO₂) synthesized with a multi-anvil device was measured over the range of 5–303 K using the heat capacity option of a physical properties measurement system (PPMS) and around ambient temperature using a differential scanning calorimeter (DSC). The entropy of stishovite at standard temperature and pressure calculated from DSC-corrected PPMS data is 24.94 J mol⁻¹ K⁻¹, which is considerably smaller (by 2.86 J mol⁻¹ K⁻¹) than that determined from adiabatic calorimetry (Holm et al. in Geochimica et Cosmochimica Acta 31:2289–2307, 1967) and about 4% larger than the recently reported value (Akaogi et al. in Am Mineral 96:1325–1330, 2011). The coesite–stishovite phase transition boundary calculated using the newly determined entropy value of stishovite agrees reasonably well with the previous experimental results by Zhang et al. (Phys Chem Miner 23:1–10, 1996). The calculated phase boundary of kyanite decomposition reaction is most comparable with the experimental study by Irifune et al. (Earth Planet Sci Lett 77:245–256, 1995) at low temperatures around 1,400 K, and the calculated slope in this temperature range is mostly consistent with that determined by in situ X-ray diffraction experiments (Ono et al. in Am Mineral 92:1624–1629, 2007).     

In search of the mixed derivative ?2 M /? P ? T ( M = G , K ): joint analysis of ultrasonic data for polycrystalline pyrope from gas- and solid-medium apparatus

Elastic wave velocities for dense (99.8% of theoretical density) isotropic polycrystalline specimens of synthetic pyrope (Mg₃Al₂Si₃O₁₂) were measured to 1,000 K at 300 MPa by the phase comparison method of ultrasonic interferometry in an internally heated gas-medium apparatus. The temperature derivatives of the elastic moduli [(∂Ks/∂T) _P = −19.3(4); (∂G/∂T) _P = −10.4(2) MPa K⁻¹] measured in this study are consistent with previous acoustic measurements on both synthetic polycrystalline pyrope in a DIA-type cubic anvil apparatus (Gwanmesia et al. in Phys Earth Planet Inter 155:179–190, 2006) and on a natural single crystal by the rectangular parallelepiped resonance (RPR; Suzuki and Anderson in J Phys Earth 31:125–138, 1983) method but |(∂Ks/∂T) _P | is significantly larger than from a Brillouin spectroscopy study of single-crystal pyrope (Sinogeikin and Bass in Phys Earth Planet Inter 203:549–555, 2002). Alternative approaches to the retrieval of mixed derivatives of the elastic moduli from joint analysis of data from this study and from the solid-medium data of Gwanmesia et al. in Phys Earth Planet Inter 155:179–190 (2006) yield ∂² G/∂P∂T = [0.07(12), 0.20(14)] × 10⁻³ K⁻¹ and ∂² K _S /∂P∂T = [−0.20(24), 0.22(26)] × 10⁻³ K⁻¹, both of order 10⁻⁴ K⁻¹ and not significantly different from zero. More robust inference of the mixed derivatives will require solid-medium acoustic measurements of precision significantly better than 1%.     

Coupling of oceanic and continental crust during Eocene eclogite-facies metamorphism: evidence from the Monte Rosa nappe,western Alps

High precision U–Pb geochronology of rutile from quartz–carbonate–white mica–rutile veins that are hosted within eclogite and schist of the Monte Rosa nappe, western Alps, Italy, indicate that the Monte Rosa nappe was at eclogite-facies metamorphic conditions at 42.6 ± 0.6 Ma. The sample area [Indren glacier, Furgg zone; Dal Piaz (2001) Geology of the Monte Rosa massif: historical review and personal comments. SMPM] consists of eclogite boudins that are exposed inside a south-plunging overturned synform within micaceous schist. Associated with the eclogite and schist are quartz–carbonate–white mica–rutile veins that formed in tension cracks in the eclogite and along the contact between eclogite and surrounding schist. Intrusion of the veins at about 42.6 Ma occurred at eclogite-facies metamorphic conditions (480–570°C, >1.3–1.4 GPa) based on textural relations, oxygen isotope thermometry, and geothermobarometry. The timing of eclogite-facies metamorphism in the Monte Rosa nappe determined in this study is identical to that of the Gran Paradiso nappe [Meffan-Main et al. (2004) J Metamorphic Geol 22:261–281], confirming that these two units have shared the same Alpine metamorphic history. Furthermore, the Gran Paradiso and Monte Rosa nappes underwent eclogite-facies metamorphism within the same time interval as the structurally overlying Zermatt-Saas ophiolite [∼50–40 Ma; e.g., Amato et al. (1999) Earth Planet Sci Lett 171:425–438; Mayer et al. (1999) Eur Union Geosci 10:809 (abstract); Lapen et al. (2003) Earth Planet Sci Lett 215:57–72]. The nearly identical P–T–t histories of the Gran Paradiso, Monte Rosa, and Zermatt-Saas units suggest that these units shared a common Alpine tectonic and metamorphic history. The close spatial and temporal associations between high pressure (HP) ophiolite and continental crust during Alpine orogeny indicates that the HP internal basement nappes in the western Alps may have played a key role in exhumation and preservation of the ophiolitic rocks through buoyancy-driven uplift. Coupling of oceanic and continental crust may therefore be critical in preventing permanent loss of oceanic crust to the mantle.     

Zr diffusion in titanite

Chemical diffusion of Zr under anhydrous, pO₂-buffered conditions has been measured in natural titanite. The source of diffusant was either zircon powder or a ZrO₂–Al₂O₃–titanite mixture. Experiments were run in sealed silica glass capsules with solid buffers (to buffer at NNO or QFM). Rutherford Backscattering Spectrometry (RBS) was used to measure diffusion profiles. The following Arrhenius parameters were obtained for Zr diffusion parallel to c over the temperature range 753–1,100°C under NNO-buffered conditions: D _Zr = 5.33 × 10⁻⁷ exp(−325 ± 30 kJ mol⁻¹/RT) m² s⁻¹ Diffusivities are similar for experiments buffered at QFM. These data suggest that titanite should be moderately retentive of Zr chemical signatures, with diffusivities slower than those for O and Pb in titanite, but faster than those for Sr and the REE. When applied in evaluation of the relative robustness of the recently developed Zr-in-titanite geothermometer (Hayden and Watson, Abstract, 16th V.M. Goldschmidt Conference 2006), these findings suggest that Zr concentrations in titanite will be less likely to be affected by later thermal disturbance than the geothermometer based on Zr concentrations in rutile (Zack et al. in Contrib Mineral Petrol 148:471–488, 2004; Watson et al. in Contrib Mineral. Petrol, 2006), but much less resistant to diffusional alteration subsequent to crystallization than the Ti-in-Zircon geothermometer (Watson and Harrison in Science 308:841–844, 2005).     

A fluid inclusion and stable isotope study of 200 Ma of fluid evolution in the Galway Granite, Connemara, Ireland

Fluid inclusions in granite quartz and three generations of veins indicate that three fluids have affected the Caledonian Galway Granite. These fluids were examined by petrography, microthermometry, chlorite thermometry, fluid chemistry and stable isotope studies. The earliest fluid was a H₂O-CO₂-NaCl fluid of moderate salinity (4–10 wt% NaCl eq.) that deposited late-magmatic molybdenite mineralised quartz veins (V₁) and formed the earliest secondary inclusions in granite quartz. This fluid is more abundant in the west of the batholith, corresponding to a decrease in emplacement depth. Within veins, and to the east, this fluid was trapped homogeneously, but in granite quartz in the west it unmixed at 305–390 °C and 0.7–1.8 kbar. Homogeneous quartz δ¹⁸O across the batholith (9.5 ± 0.4‰n = 12) suggests V₁ precipitation at high temperatures (perhaps 600 °C) and pressures (1–3 kbar) from magmatic fluids. Microthermometric data for V₁ indicate lower temperatures, suggesting inclusion volumes re-equilibrated during cooling. The second fluid was a H₂O-NaCl-KCl, low-moderate salinity (0–10 wt% NaCl eq.), moderate temperature (270–340 °C), high δD (−18 ± 2‰), low δ¹⁸O (0.5–2.0‰) fluid of meteoric origin. This fluid penetrated the batholith via quartz veins (V₂) which infill faults active during post-consolidation uplift of the batholith. It forms the most common inclusion type in granite quartz throughout the batholith and is responsible for widespread retrograde alteration involving chloritization of biotite and hornblende, sericitization and saussuritization of plagioclase, and reddening of K-feldspar. The salinity was generated by fluid-rock interactions within the granite. Within granite quartz this fluid was trapped at 0.5–2.3 kbar, having become overpressured. This fluid probably infiltrated the Granite in a meteoric-convection system during cooling after intrusion, but a later age cannot be ruled out. The final fluid to enter the Granite and its host rocks was a H₂O-NaCl-CaCl₂-KCl fluid with variable salinity (8–28 wt% NaCl eq.), temperature (125–205 °C), δD (−17 to −45‰), δ¹⁸O (−3 to + 1.2‰), δ¹³C_CO2 (−19 to 0‰) and δ³⁴S_sulphate (13–23‰) that deposited veins containing quartz, fluorite, calcite, barite, galena, chalcopyrite sphalerite and pyrite (V₃). Correlations of salinity, temperature, δD and δ¹⁸O are interpreted as the result of mixing of two fluid end-members, one a high-δD (−17 to −8‰), moderate-δ¹⁸O (1.2–2.5‰), high-δ¹³C_CO2 (> −4‰), low-δ³⁴S_sulphate (13‰), high-temperature (205–230 °C), moderate-salinity (8–12 wt% NaCl eq.) fluid, the other a low-δD (−61 to −45‰), low-δ¹⁸O (−5.4 to −3‰), low-δ¹³C (<−10‰), high-δ³⁴S_sulphate (20–23‰) low-temperature (80–125 °C), high-salinity (21–28 wt% NaCl eq.) fluid. Geochronological evidence suggests V₃ veins are late Triassic; the high-δD end-member is interpreted as a contemporaneous surface fluid, probably mixed meteoric water and evaporated seawater and/or dissolved evaporites, whereas the low-δD end-member is interpreted as a basinal brine derived from the adjacent Carboniferous sequence. This study demonstrates that the Galway Granite was a locus for repeated fluid events for a variety of reasons; from expulsion of magmatic fluids during the final stages of crystallisation, through a meteoric convection system, probably driven by waning magmatic heat, to much later mineralisation, concentrated in its vicinity due to thermal, tectonic and compositional properties of granite batholiths which encourage mineralisation long after magmatic heat has abated. Received: 3 April 1996 / Accepted: 5 May 1997     

Suspended Matter,Chl-a,CDOM, Grain Sizes,and Optical Properties in the Arctic Fjord-Type Estuary,Kangerlussuaq, West Greenland During Summer

Optical constituents as suspended particulate matter (SPM), chlorophyll (Chl-a), colored dissolved organic matter (CDOM), and grain sizes were obtained on a transect in the arctic fjord-type estuary Kangerlussuaq (66°) in August 2007 along with optical properties. These comprised diffuse attenuation coefficient of downwelling PAR (K _d(PAR)), upwelling PAR (K _u(PAR)), particle beam attenuation coefficient (c _p), and irradiance reflectance R(−0, PAR). PAR is white light between 400 and 700 nm. The estuary receives melt water from the Greenland Inland Ice and stations covered a transect from the very high turbid melt water outlet to clear marine waters. Results showed a strong spatial variation with high values as for suspended matter concentrations, CDOM, diffuse attenuation coefficient K _d(PAR), particle beam attenuation coefficients (c _p), and reflectance R(−0, PAR) at the melt water outlet. Values of optical constituents and properties decreased with distance from the melt water outlet to a more or less constant level in central and outer part of the estuary. There was a strong correlation between inorganic suspended matter (SPMI) and diffuse attenuation coefficient K _d(PAR) (r ² = 0.92) and also for particle beam attenuation coefficient (c _p; r ² = 0.93). The obtained SPMI specific attenuation—K _d^*(PAR) = 0.13 m² g⁻¹ SPMI—and the SPMI specific particle beam attenuation—c _p^* = 0.72 m² g⁻¹—coefficients were about two times higher than average literature values. Irradiance reflectance R(−0, PAR) was comparatively high (0.09−0.20) and showed a high (r ² = 0.80) correlation with K _u(PAR). Scattering dominated relative to absorption—b(PAR)/a(PAR) = 12.3. Results strongly indicated that the high values in the optical properties were related to the very fine particle sizes (mean = 2–6 μm) of the suspended sediment. Data and results are discussed and compared to similar studies from both temperate and tropical estuaries.