Laboratory Investigations on Extremely High Suction Measurements for Fine-grained Soils

A Dewpoint PotentiaMeter (WP4) was used to measure suction of two fine-grained soils: a locally available silty soil and commercially available white clay, rapidly. Using these results, efforts were made to check the suitability and efficiency of various fitting functions, for defining the soil–water characteristic curve, SWCC, for high suction ranges (0–80 MPa). In addition to this, a knowledge-based database SoilVision 3.34 was used to estimate the SWCC using Pedo-transfer functions, PTFs. The study brings out that the Fredlund et al. [1997, Proceedings of the 3rd Symposium on Unsaturated Soil, Rio de Janeiro, Brazil, pp. 13–23] PTF yields the best estimate of SWCC for fine-grained soils. The influence of the soil type and dry unit weight, on suction and the SWCC fitting parameters, have also been studied.     

Experimental Validation of an Method for In Situ Estimation of Hydraulic Conductivity of Water Treatment Granular Materials

Domestic wastewater-treatment system is today a widely held technique. In order to evaluate the efficiency of these systems, it is necessary to determinate the in situ hydraulic conductivity of the water-treatment granular materials constituting them. The in situ measurement of the hydraulic conductivity of soils has proven to be imprecise, take a long time. Empirical equations based on the physical properties of soil have been proposed to overcome these difficulties. In this context, the originality of this paper is to propose an in situ method to obtain reliable input parameters for the predictive equations of Chapuis (Can Geotech J 41(5): 787–795, 2004) and Carrier (J Geotech Geoenviron Eng 129(11): 1054–1056, 2003) by using both a penetrometer and a geoendoscope. This method is described in the first part of this paper. Then, the validation of the method based on laboratory tests performed under controlled conditions for three kinds of soils is presented. Hydraulic conductivity obtained when applying the method is compared to that measured by a Darcy permeameter. The difference between these two hydraulic conductivities is less than 25%. Finally, the precision of the results is discussed.     

Residual and Fully Softened Strength Evaluation of Soils using Artificial Neural Networks

A backpropagation artificial neural network (ANN) model is developed to predict the secant friction angle of residual and fully softened soils, using data reported by Stark et al. (J Geotech Geoenviron Eng ASCE 131:575–588, 2005). In the ANN model, index properties such as liquid limit, plastic limit, activity, clay fraction and effective normal stress are used as input variables while secant residual friction angle is used as output variable. The model is verified using data that were not used for model training and testing. The results also indicate that the secant residual friction angle of cohesive soils can be predicted quite accurately using liquid limit, clay fraction and effective normal stress as input variables with R ² = 0.93. The sensitivity analysis results indicate that plastic limit and activity have no appreciable effect on ANN predicted secant friction angles. The secant friction angle predictions of the ANN model were also compared with those of Stark’s et al. (2005) curves and the empirical formulas suggested for the same data sets by Wright (Evaluation of soil shear strengths for slope and retaining wall stability with emphasis on high plasticity clays, 2005). The comparison shows that the ANN model predictions are very close to those suggested by the Stark et al. (2005) curves but much better than the prediction of Wright’s (2005) empirical equations. The results also show that ANN is an alternative powerful tool to predict the secant friction angle of soils.     

Reply to Oczlon’s (2006) comment

The comment of Martin Oczlon contains some significant contributions to the topics discussed in the article of Huckriede et al. (Int J Earth Sci 93:414–431, 2004). Contrary to Oczlon’s comment, the central results of Huckriede et al. (Int J Earth Sci 93:414–431, 2004) are clearly different from the tectonic model of Oczlon (Geol Rundsch 83:20–31, 1994). Additionally, there is no reason for a new interpretation of the ⁴⁰K/⁴⁰Ar muscovite cooling-ages from allochthonous units.     

Crossover of cation partitioning in olivines: a combination of ab initio and Monte Carlo study

We report studies based on a combination of ab initio electronic structure and Monte Carlo (MC) technique on the problem of cation partitioning among inequivalent octahedral sites, M1 and M2 in mixed olivines containing Mg²⁺ and Fe²⁺ ions. Our MC scheme uses interactions derived out of ab initio, density functional calculations carried out on measured crystal structure data. Our results show that there is no reversal of the preference of Fe for M1 over M2 as a function of temperature. Our findings do not agree with the experimental findings of Redfern et al. (Phys Chem Miner 27:630–637, 2000), but are in agreement with those of Heinemann et al. (Eur J Mineral 18:673–689, 2006) and Morozov et al. (Eur J Mineral 17:495–500, 2005).     

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.     

Closed-Form Solutions for a Circular Tunnel in Elastic-Brittle-Plastic Ground with the Original and Generalized Hoek–Brown Failure Criteria

The paper presents a closed-form solution for the convergence curve of a circular tunnel in an elasto-brittle-plastic rock mass with both the Hoek–Brown and generalized Hoek–Brown failure criteria, and a linear flow rule, i.e., the ratio between the minor and major plastic strain increments is constant. The improvement over the original solution of Brown et al. (J Geotech Eng ASCE 109(1):15–39, 1983) consists of taking into account the elastic strain variation in the plastic annulus, which was assumed to be fixed in the original solution by Brown et al. The improvement over Carranza-Torres’ solution (Int J Rock Mech Min Sci 41(Suppl 1):629–639, 2004) consists of providing a closed-form solution, rather than resorting to numerical integration of an ordinary differential equation. The presented solution, by rigorously following the theory of plasticity, takes into account that the elastic strain components change with radial and circumferential stress changes within the plastic annulus. For the original Hoek–Brown failure criterion, disregarding the elastic strain change leads to underestimate the convergence by up to 55%. For a rock mass failing according to the generalized Hoek–Brown failure criterion, using the original failure criterion leads to a high probability (97%) of underestimating the convergence by up to 100%. As a consequence, the onset or degree of squeezing may be underestimated, and the loading on the support/reinforcement calculated with the convergence/confinement method may be largely underestimated.     

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.     

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).     

Focal parameters of seismic sources during the 1981 and 1983 eruption at Mt. Etna volcano (Sicily,Italy)

In this study 50 seismic events, preceding and accompanying the eruptions occurring in 1981 and 1983, have been considered. Seismic moments, fault radii, stress drops and seismic energies have been calculated using Brune’s model (J Geophys Res 75:4997–5009, 1970; J Geophys Res 76:5002, 1971); site, anelastic attenuation along the propagation path, geometrical spreading and interaction with the free surface effects are taken into account. For each event we have also estimated the equivalent Wood–Anderson magnitude (MWAeq) (Scherbaum and Stoll in Bull Seism Soc Am 73:1321–1343, 1983); relations among all these source parameters have been determined. Furthermore, the hypothesis of self-similarity (Aki in J Geophys Res 72:1217–1231, 1967) is not verified for events with seismic moments <10¹² N-m: in fact the relationship between log-stress drop and log-moment is linear up to a moment of 10¹² N-m (events of 1981 eruption), while for higher moments (events of 1983 eruption) the slope of the regression line is not significantly different from zero. We suppose that such a behaviour is related to a heterogeneous medium with barriers on the faults. Finally, the main conclusion is that eruptions of 1981 and 1983 differ from one another both in eruptive and seismic aspects; analysis of seismic energies indicates an increase in Mt. Etna’s activity, confirmed by studies performed on the following lateral eruption of 1991–1993 (Patanè et al. in Bull Volcanol 47:941–952, 1995), occurring on the same structural trend.     

Numerical investigation of undrained cavity expansion in fine-grained soils

Based on the spherical cavity expansion (SCE) problem, Cudmani and Osinov (Can Geotech J 38:622–638, 2001), Osinov and Cudmani (Int J Numer Anal Method Geomech 25:473–495, 2001) developed a semi-empirical method of interpretation of CPT for coarse-grained soils (sand, gravel) using a hypoplastic constitutive model. Using a material-independent shape factor, the cone penetration resistance was related to the limit pressure required to expand a spherical cavity. The shape factor was observed to be a function of the soil state only, in particular the pressure-dependent relative density. This paper presents an analogous interpretation technique for CPT in fine-grained soils using the shape factor concept, Cavity Expansion approach, and a hypoplastic constitutive model. Relations for the shape factor and the limit pressure have been proposed based on the parameters affecting these quantities. A validation of the proposed interpretation technique with experimental results has also been performed.     

Paraelasticity

A 3D hysteretic, fully reversible constitutive model with rate-independent damping and with variable secant stiffness is proposed. A reversible dilatancy–contractancy effect is an optional feature as described in the companion paper (Niemunis et al. in Acta Geotech, 2011). The present paper describes the basic aspects of the model only. A strain path reversal is defined, and a treatment of the past reversals using a stack structure is proposed.     

Hydroxyl in mantle olivine xenocrysts from the Udachnaya kimberlite pipe

The incorporation of hydrogen in mantle olivine xenocrysts from the Udachnaya kimberlite pipe was investigated by Fourier-transform infrared spectroscopy and secondary ion mass spectrometry (SIMS). IR spectra were collected in the OH stretching region on oriented single crystals using a conventional IR source at ambient conditions and in situ at temperatures down to −180°C as well as with IR synchrotron radiation. The IR spectra of the samples are complex containing more than 20 strongly polarized OH bands in the range 3,730–3,330 cm⁻¹. Bands at high energies (3,730–3,670 cm⁻¹) were assigned to inclusions of serpentine, talc and the 10 Å phase. All other bands are believed to be intrinsic to olivine. The corresponding point defects are (a) associated with vacant Si sites (3,607 cm⁻¹ E || a, 3,597 E || a, 3,571 cm⁻¹ E || c, 3,567 E || c, and 3,556 E || b), and (b) with vacant M1 sites (most of the bands polarized parallel to a). From the pleochroic behavior and position of the OH bands associated with the vacant M1 sites, we propose two types of hydrogen—one bonded to O1 and another to O2, so that both OH vectors are strongly aligned parallel to a. The O2–H groups may be responsible for the OH bands at higher wavenumbers than those for the O1–H groups. The multiplicity of the corresponding OH bands in the spectra can be explained by different chemical environments and by slightly different distortions of the M1 sites in these high-pressure olivines. Four samples were investigated by SIMS. The calculated integral molar absorption coefficient using the IR and SIMS results of 37,500±5,000 L mol H₂O cm⁻² is within the uncertainties slightly higher than the value determined by Bell et al. (J Geophys Res 108(B2):2105–2113, 2003) (28,450±1,830 L mol H₂O cm⁻²). The reason for the difference is the different distributions of the absorption intensity of the spectra of both studies (mean wavenumber 3,548 vs. 3,570 cm⁻¹). Olivine samples with a mean wavenumber of about 3,548 cm⁻¹ should be quantified with the absorption coefficient as determined in this study; those containing more bands at higher wavenumber (mean wavenumber 3,570 cm⁻¹) should be quantified using the value determined by Bell et al. (J Geophys Res 108(B2):2105–2113, 2003).

Insights into the high-pressure behavior of kaolinite from infrared spectroscopy and quantum-mechanical calculations

The high-pressure behavior of Keokuk kaolinite has been studied to 9.5 GPa by infrared spectroscopy using synchrotron radiation. The kaolinite-I → kaolinite-II and kaolinite-II → kaolinite-III transformations have clear spectroscopic expression, with discontinuities coinciding with the transformation pressures bracketed by X-ray diffraction (Welch and Crichton in Am Mineral 95:651–654, 2010). The experimental spectra have been interpreted from band assignments derived from density functional theory for the structures of kaolinite-II and kaolinite-III, using as starting models the ab initio structures reported by Mercier and Le Page (Acta Crystallogr A B64:131–143, 2008, Mater Sci Technol 25:437–442, 2009) and unit-cell parameters from Welch and Crichton (Am Mineral 95:651–654, 2010). The relaxed theoretical structures are very similar to those reported by Mercier and Le Page (Acta Crystallogr A B64:131–143, 2008, Mater Sci Technol 25:437–442, 2009) in their theoretical investigation of kaolinite polytypes at high pressure. The vibrational spectra calculated from the quantum-mechanical analysis allow band assignments of the IR spectra to be made and provide insights into the behavior of different OH environments in the two high-pressure polytypes. The single perpendicular-interlayer OH group of kaolinite-III has a distinctive spectroscopic signature that is diagnostic of this polytype (ν = 3,595 cm⁻¹ at 9.5 GPa) and is sensitive to the compression/expansion of the interlayer space. This OH group also has a distinctive signature in the calculated spectra. The spectra collected on decompression are those of kaolinite-III and persist largely unchanged to 4.6 GPa, except for a continuous blue shift of the 3,595 cm⁻¹ band to 3,613 cm⁻¹. Finally, kaolinite-I is recovered at 0.6 GPa, confirming the kaolinite-III → kaolinite-I transformation previously observed by X-ray diffraction, and the irreversibility of the kaolinite-II → kaolinite-III transformation. The ambient spectra collected at the start and finish of the experiment are those of kaolinite-I, and start/finish band frequencies agree to within 6 cm⁻¹.     

The quantitative analysis of OH in vesuvianite: a polarized FTIR and SIMS study

A well-characterized suite of vesuvianite samples from the volcanic ejecta (skarn or syenites) from Latium (Italy) was studied by single-crystal, polarized radiation, Fourier-transform infrared (FTIR) spectroscopy and secondary-ion mass-spectrometry (SIMS). OH-stretching FTIR spectra consist of a rather well-defined triplet of broad bands at higher-frequency (3,700–3,300 cm^–1) and a very broad composite absorption below 3,300 cm^–1. Measurements with E//c or Ec show that all bands are strongly polarized with maximum absorption for E//c. They are in agreement with previous band assignments (Groat et al. Can Mineral 33:609, 1995) to the two O(11)–H(1) and O(10)–H(2) groups in the structure. Pleochroic measurements with changing direction of the E vector of the incident radiation show that the orientation of the O(11)–H(1) dipole is OHc~35°, in excellent agreement with the neutron data of Lager et al. (Can Mineral 37:763, 1999). A SIMS-based calibration curve at ~10% rel. accuracy has been worked out and used as reference for the quantitative analysis of H₂O in vesuvianite by FTIR. Based on previous SIMS results for silicate minerals (Ottolini and Hawthorne in J Anal At Spectrom 16:1266, 2001; Ottolini et al. in Am Mineral 87:1477, 2002) the SiO₂ and FeO content of the matrix were assumed as the major factors to be considered at a first approximation in the selection of the standards for H. The lack of vesuvianite standards for quantitative SIMS analysis of H₂O has been here overcome by selecting low-silica elbaite crystals (Ottolini et al. in Am Mineral 87:1477, 2002). The resulting integrated molar absorption FTIR coefficient for vesuvianite is _i=100.000±2.000 l mol^–1 cm^–2. SIMS data for Li, B, F, Sr, Y, Be, Ba REE, U and Th are also provided in the paper.     

Upper-bound multi-rigid-block solutions for bearing capacity of two-layered soils

The limit analysis method has been widely used in the stability analysis of geotechnical problems including the bearing capacity of foundations. Two main approaches have been followed in the limit analysis to improve the calculation of the bearing capacity of foundations. One approach is to combine limit analysis with the finite element method and linear/nonlinear programming. The other is to use a multi-rigid-block mechanism to obtain an upper-bound solution. In this paper, the multi-rigid-block upper-bound method with a modified failure mechanism from that of Florkiewicz [Florkiewicz A. Upper bound to bearing capacity of layered soils. Can Geotech J 1989;26(4):730–6.] was employed to calculate the bearing capacity of foundations. Attention was paid particularly to the bearing capacity of strip footings over a two-layered soil. In order to verify the effectiveness of the modified mechanism, comparisons were made with other well-known solutions. The results showed improvements over the best available multi-rigid-block upper-bound solutions given by Michalowski and Shi [Michalowski RL, Shi L. Bearing capacity of footings over two-layer foundation soils. J Geotech Eng ASCE 1995;121(5):421–8.], and fair consistence with the results from the finite element limit analysis in Shiau et al. [Shiau JS, Lyamin AV, Sloan SW. Bearing capacity of a sand layer on clay by finite element limit analysis. Can Geotech J 2003;40(5):900–15.].     

Synchrotron IR study of hydrous ringwoodite (γ-Mg2SiO4) up to 30 GPa

High-pressure synchrotron infrared (IR) absorption spectra were collected between 650 and 4,000 cm⁻¹ at ambient temperature for hydrous Mg-ringwoodite (γ-Mg₂SiO₄) up to 30 GPa. The main feature in the OH⁻ stretching region is an extremely broad band centred at 3,150 cm⁻¹. The hydrogen bond is strong for most protons and the most probable site for protonation is the tetrahedral edge. With increasing pressure, this band shifts downward while decreasing its integrated intensity until disappearance at a pressure of 25 GPa. Only one band at 2,450 cm⁻¹ and an absorption plateau persist with a maximum wavenumber of 3,800 cm⁻¹. This behaviour is reversible upon pressure release. We interpret this as a second-order phase transition occurring in hydrated Mg-ringwoodite at high pressure (beyond ∼ 25 GPa). This result is compatible with the observation by Kleppe et al. (Phys Chem Miner 29:473–476, 2002a) who suggested the presence of Si–O–Si linkages and/or partial increase in the coordination of Si. Beyond the phase transition, the protons are delocalized and their environment on the ringwoodite structure is probably quite different from that at low pressure. Data obtained in situ at high pressures and temperatures are needed to better understand the effect of protonation on the structure and to better constrain this phase transition.     

Basanite–nephelinite suite from early Kilauea: carbonated melts of phlogopite–garnet peridotite at Hawaii’s leading magmatic edge

A basanite–nephelinite glass suite from early submarine Kilauea defines a continuous compositional array marked by increasing concentrations of incompatible components with decreasing SiO₂, MgO, and Al₂O₃. Like peripheral and post-shield strongly alkalic Hawaiian localities (Clague et al. in J Volcanol Geotherm Res 151:279–307, 2006; Dixon et al. in J Pet 38:911–939, 1997), the early Kilauea basanite–nephelinite glasses are interpreted as olivine fractionation products from primary magnesian alkalic liquids. For early Kilauea, these were saturated with a garnet–phlogopite–sulfide peridotite assemblage, with elevated dissolved CO₂ contents responsible for the liquids’ distinctly low-SiO₂ concentrations. Reconstructed primitive liquids for early Kilauea and other Hawaiian strongly alkalic localities are similar to experimental 3 GPa low-degree melts of moderately carbonated garnet lherzolite, and estimated parent magma temperatures of 1,350–1,400°C (olivine–liquid geothermometry) match the ambient upper mantle geotherm shortly beneath the base of the lithosphere. The ~3 GPa source regions were too hot for stable crystalline carbonate and may have consisted of ambient upper mantle peridotite containing interstitial carbonate–silicate or carbonatitic liquid, possibly (Dixon et al. in Geochem Geophys Geosyst 9(9):Q09005, 2008), although not necessarily, from the Hawaiian mantle plume. Carbonate-enriched domains were particularly susceptible to further melting upon modest decompression during upward lithospheric flexure beneath the advancing Hawaiian Arch, or by conductive heating or upward drag by the Hawaiian mantle plume. The early Kilauea basanite–nephelinite suite has a HIMU-influenced isotopic character unlike other Hawaiian magmas (Shimizu et al. in EOS Tran Amer Geophys Union 82(47): abstr V12B-0962, 2001; Shimizu et al. in Geochim Cosmochim Acta 66(15A):710, 2002) but consistent with oceanic carbonatite involvement (Hoernle et al. in Contrib Mineral Petrol 142:520–542, 2002). It may represent the melting products of a fertile domain in the ambient upper mantle impinged upon and perturbed by the sustained plume source that feeds later shield-stage magmatism.     

Scaling and multivariate analysis of medium to large landslide events: Haida Gwaii, British Columbia

Gimbarzevsky (1988) collected an exceptional landsliding inventory for Haida Gwaii, British Columbia that included over 8,000 landsliding vectors covering an area of approximately 10,000 km². This database was never published in the referred literature, despite its regional significance. It was collected prior to widespread application of GIS technologies in landsliding studies, limiting the analyses undertaken at the time. Gimbarzevsky identified landslides using 1:50,000 aerial photographs, and transferred the information to NTS map sheets. In our study, we digitized the landslide vectors from these original map sheets and connected each landslide to a digital elevation model. Lengths of landslide vectors are compared to the landsliding inventory for Haida Gwaii analyzed in Rood (1984), Martin Y et al. BC Can J Earth Sci 39:289–305 (2002); the latter inventory is based on larger-scale aerial photographs (~1:12,000). Rood’s database contains a more complete record of smaller landslides, while the inventory of Gimbarzevsky provides improved statistical representation of less frequent, medium to large landslides. It is suggested that combined landslide delineation at different scales could provide a more complete landslide record. Discriminant analysis was undertaken to assess which of nine predictor variables, chosen on the basis of mechanical theory, best predict failed versus unfailed locations. Seven of the nine variables were found to be statistically significant in discriminating amongst failed and unfailed locations. Results show that 81.7% of original grouped cases were correctly classified.     

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.