The trapped fluid phase in upper mantle xenoliths from Victoria,Australia: implications for mantle metasomatism

Mantle-derived xenoliths of spinel lherzolite, spinel pyroxenite, garnet pyroxenite and wehrlite from Bullenmerri and Gnotuk maars, southwestern Victoria, Australia contain up to 3 vol.% of fluids trapped at high pressures. The fluid-filled cavities range in size from fluid inclusions (1–100 m) up to vugs 11/2 cm across, lined with euhedral high-pressure phases. The larger cavities form an integral part of the mosaic microstructure. Microthermometry and Raman laser microprobe analysis show that the fluids are dominantly CO₂. Small isolated inclusions may have densities 1.19 g/cm³, but most inclusions show microstructural evidence of partial decrepitation during eruption, and these have lower fluid densities. Mass-spectrometric analysis of gases released by crushing or heating shows the presence of He, N₂, Ar, H₂S, CO_s and SO₂ in small quantities; these may explain the small freezing-point depressions observed in some inclusions. Petrographic, SEM and microprobe studies show that the trapped fluids have reacted with the cavity walls (in clinopyroxene grains) to produce secondary amphiboles and carbonates. The trapped CO₂ thus represents only a small residual proportion of an original volatile phase, which has undergone at least two stages of modification — first by equilibration with spinel lherzolite to form amphibole (±mica±apatite), then by limited reaction with the walls of the fluid inclusions. The inferred original fluid was a CO₂-H₂O mixture, with significant contents of (at least) Cl and sulfur species. Generation of this fluid phase in the garnet-peridotite stability field, followed by its migration to the spinel peridotite stability field, would provide an efficient mechanism for metasomatic enrichment of the upper mantle in LIL elements. This migration could involve either a volatile flux or transport in small volumes of silicate melt that crystallize in the spinel peridotite field. These observations suggest that some portions of the subcontinental upper mantle contain large reservoirs of free fluid CO₂, which may be liberated during episodes of rifting or magmatism, to induce granulite-facies metamorphism of the lower crust.     

Petrography of low rank coals with reference to liquefaction potential

Ten sub-bituminous coals were analyzed petrographically and chemically as part of a program by the Alberta Research Council to investigate the compositional properties of Alberta coals and relationship to liquefaction potential. In addition, four lignites, two from Saskatchewan and one each from Texas and North Dakota, were studied for comparative purposes. Liquefaction experiments were carried out on three of the Alberta coals and on density fractions of the two U.S. lignites. Petrographically the ten sub-bituminous coals could be divided into two populations; one with high huminite and liptinite contents (>90%) and one with substantially lower contents (70–75%). The two U.S. lignites showed huminite plus liptinite contents of less than 70% while the Saskatchewan coals were somewhat higher. Reflectance measurements on the sub-bituminous coals showed a range of 0.33 to 0.51 with a rough correlation to the ASTM rank designations of these coals. Reflectivities determined on the lignites were generally lower. Float-sink fractionation of the two U.S. lignites and two sub-bituminous coals showed enrichment of huminite and liptinite in the lighter fractions and inertinite and mineral matter in the heavier splits. Examination of the relationship between liquefaction yields and petrographic composition showed a rough positive correlation with huminite-liptinite content. Examination of liquefaction residues suggested total conversion of the liptinite, extensive conversion of huminite and possible partial conversion and reactivity of the semifusinite maceral of the inertinite group.     

Notes on late Cenomanian and Turonianammonites from Touraine,Western France

Recent collecting in the type area of the Turonian stage and adjacent parts of Touraine in western France provides new geographic and stratigraphic records as well as additional systematic information on species of the ammonite genera Puzosia, Watinoceras, Metoicoceras, Mammites, Prohauericeras, Collignoniceras, Lecointriceras, Subprionocyclus and Sciponoceras.     

Sulphur saturation in the lower and critical zones of the Eastern Bushveld Complex

A suite of ultramafic and mafic rocks from the lower, critical and lower portion of the main zones of the Bushveld Complex has been analysed for Th, Cs, Zr, Ni, Cr and Au by INAA and XRF spectrometry. The incompatible elements Th, Cs, and Zr correlate positively, and show a gradual upward increase in abundance. Assuming constant average proportion of intercumulus material, this upward increase implies that the zones of the Complex studied represent crystallization of a single magma type some 3600 m thick. Pyroxenites dominate the lower portion of the section studied and their Ni content shows an initial rapid decrease from 850 ppm in the lowermost rocks, to around 500 ppm, with considerable scatter. This distribution is most likely to have resulted from bottom crystallization with superimposed convective overturn near the transient floor of the chamber. Gold abundances are generally higher in chromitites, and correlate positively with Ni, indicating the presence of significant amounts of cumulus immiscible sulphide. In the silicate rocks, Au does not correlate with any of the analysed elements, and it is concluded that Au was trapped in small quantities of immiscible sulphide which precipitated continuously during crystallization. There is an upward increase in the amount of cumulus immiscible sulphide, indicating a progressive increase in sulphur solubility in the magma.     

Overstepping the garnet isograd: a comparison of QuiG barometry and thermodynamic modeling

The consequences of overstepping the garnet isograd reaction have been investigated by comparing the composition of garnet formed at overstepped P–T conditions (the overstep or “OS” model) with the P–T conditions that would be inferred by assuming garnet nucleated in equilibrium with the matrix assemblage at the isograd (the equilibrium or “EQ” model). The garnet nucleus composition formed at overstepped conditions is calculated as the composition that produces the maximum decrease in Gibbs free energy from the equilibrated, garnet-absent, matrix assemblage for the bulk composition under study. Isopleths were then calculated for this garnet nucleus composition assuming equilibrium with the matrix assemblage (the EQ model). Comparison of the actual P–T conditions of nucleation (the OS model) with those inferred from the EQ model reveals considerable discrepancy between the two. In general, the inferred garnet nucleation P–T conditions (the EQ model) are at a lower temperature and higher or lower pressure (depending on the coexisting calcic phase(s)) than the actual (OS model) nucleation conditions. Moreover, the degree of discrepancy increases with the degree of overstepping. Independent estimates of the pressure of nucleation of garnet were made using the Raman shift of quartz inclusions in garnet (quartz-in-garnet or QuiG barometry). To test the validity of this method, an experimental synthesis of garnet containing quartz inclusions was made at 800 °C, 20 kbar, and the measured Raman shift reproduced the synthesis conditions to within 120 bars. Raman band shifts from three natural samples were then used to calculate an isochore along which garnet was presumed to have nucleated. Model calculations were made at several temperatures along this isochore (the OS model), and these P–T conditions were compared to those computed assuming equilibrium nucleation (the EQ model) to estimate the degree of overstepping displayed by these samples. A sample from the garnet isograd in eastern Vermont is consistent with overstepping of around 10 degrees and 0.6 kbar (affinities of around 2 kJ/mole garnet). A sample from the staurolite–kyanite zone in the same terrane requires overstepping of around 50 °C and 2–5 kbar (affinities of around 10–18 kJ/mole garnet). A similar amount of overstepping was inferred for a blueschist sample from Sifnos, Greece. These results indicate that overstepping of garnet nucleation reactions may be common and pronounced in regionally metamorphosed terranes, and that the P–T conditions and paths inferred from garnet zoning studies may be egregiously in error.     

Temporal evolution and compositional signatures of two supervolcanic systems recorded in zircons from Mangakino volcanic centre,New Zealand

Mangakino, the oldest rhyolitic caldera centre delineated in the Taupo Volcanic Zone of New Zealand, generated two very large (super-sized) ignimbrite eruptions, the 1.21 ± 0.04 Ma >500 km³ Ongatiti and ~1.0 Ma ~1,200 km³ Kidnappers events, the latter of which was followed after a short period of erosion by the ~200 km³ Rocky Hill eruption. We present U/Pb ages and trace-element analyses on zircons from pumice clasts from these three eruptions by Secondary Ion Mass Spectrometry (SIMS) using SHRIMP-RG instruments to illustrate the evolution of the respective magmatic systems. U–Pb age spectra from the Ongatiti imply growth of the magmatic system over ~250 kyr, with a peak of crystallisation around 1.32 Ma, ~100 kyr prior to eruption. The zircons are inferred to have then remained stable in a mush with little crystallisation and/or dissolution before later rejuvenation of the system at the lead-in to eruption. The paired Kidnappers and Rocky Hill eruptions have U–Pb zircon ages and geochemical signatures that suggest they were products of a common system grown over ~200 kyr. The Kidnappers and Rocky Hill samples show similar weakly bimodal age spectra, with peaks at 1.1 and 1.0 Ma, suggesting that an inherited antecrystic population was augmented by crystals grown at ages within uncertainty of the eruption age. In the Kidnappers, this younger age peak is dominantly seen in needle-shaped low U grains with aspect ratios of up to 18. In all three deposits, zircon cores show larger ranges and higher absolute concentrations of trace elements than zircon rims, consistent with zircon crystallisation from evolving melts undergoing crystal fractionation involving plagioclase and amphibole. Abundances and ratios of many trace elements frequently show variations between different sectors within single grains, even where there is no visible sector zoning in cathodoluminescence (CL) imaging. Substitution mechanisms, as reflected in the molar (Sc + Y + REE³⁺)/P ratio, differ in the same growth zone between the sides (along a-axis and b-axis: values approaching 1.0) and tips (c-axis: values between 1.5 and 5.0) of single crystals. These observations have implications for the use of zircons for tracking magmatic processes, particularly in techniques where CL zonation within crystals is not assessed and small analytical spot sizes cannot be achieved. These observations also limit applicability of the widely used Ti-in-zircon thermometer. The age spectra for the Ongatiti and Kidnappers/Rocky Hill samples indicate that both magmatic systems were newly built in the time-breaks after respective previous large eruptions from Mangakino. Trace element variations defining three-component mixing suggest that zircons, sourced from multiple melts, contributed to the population in each system.     

Early Neoproterozoic granulite facies metamorphism of mafic dykes from the Vestfold Block,east Antarctica

Proterozoic mafic dykes from the southwestern Vestfold Block experienced heterogeneous granulite facies metamorphism, characterized by spotted or fractured garnet‐bearing aggregates in garnet‐absent groundmass. The garnet‐absent groundmass typically preserves an ophitic texture composed of lathy plagioclase, intergranular clinopyroxene and Fe–Ti oxides. Garnet‐bearing domains consist mainly of a metamorphic assemblage of garnet, clinopyroxene, orthopyroxene, hornblende, biotite, plagioclase, K‐feldspar, quartz and Fe–Ti oxides. Chemical compositions and textural relationships suggest that these metamorphic minerals reached local equilibrium in the centre of the garnet‐bearing domains. Pseudosection calculations in the model system NCFMASHTO (Na₂O–CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃) yield P–T estimates of 820–870 °C and 8.4–9.7 kbar. Ion microprobe U–Pb zircon dating reveals that the NW‐ and N‐trending mafic dykes were emplaced at 1764 ± 25 and 1232 ± 12 Ma, respectively, whereas their metamorphic ages cluster between 957 ± 7 and 938 ± 9 Ma. The identification of granulite facies mineral inclusions in metamorphic zircon domains is also consistent with early Neoproterozoic metamorphism. Therefore, the southwestern margin of the Vestfold Block is inferred to have been buried to depths of ~30–35 km beneath the Rayner orogen during the late stage of the late Mesoproterozoic/early Neoproterozoic collision between the Indian craton and east Antarctica (i.e. the Lambert Terrane or the Ruker craton including the Lambert Terrane). The lack of penetrative deformation and intensive fluid–rock interaction in the rigid Vestfold Block prevented the nucleation and growth of garnet and resulted in the heterogeneous granulite facies metamorphism of the mafic dykes.     

Three-Dimensional Numerical Investigations of the Failure Mechanism of a Rock Disc with a Central or Eccentric Hole

The diametrical compression of a circular disc (Brazilian test) or cylinder with a small eccentric hole is a simple but important test to determine the tensile strength of rocks. This paper studies the failure mechanism of circular disc with an eccentric hole by a 3D numerical model (RFPA3D). A feature of the code RFPA3D is that it can numerically simulate the evolution of cracks in three-dimensional space, as well as the heterogeneity of the rock mass. First, numerically simulated Brazilian tests are compared with experimental results. Special attention is given to the effect of the thickness to radius ratio on the failure modes and the peak stress of specimens. The effects of the compressive strength to tensile strength ratio (C/T), the loading arc angle (2α), and the homogeneity index (m) are also studied in the numerical simulations. Secondly, the failure process of a rock disc with a central hole is studied. The effects of the ratio of the internal hole radius (r) to the radius of the rock disc (R) on the failure mode and the peak stress are investigated. Thirdly, the influence of the vertical and horizontal eccentricity of an internal hole on the initiation and propagation of cracks inside a specimen are simulated. The effect of the radius of the eccentric hole and the homogeneity index (m) are also investigated.     

Laboratory Estimation of Rock Joint Stiffness and Frictional Parameters

Numerical modeling of complex rock engineering problems involves the use of various input parameters which control usefulness of the output results. Hence, it is of utmost importance to select the right range of input physical and mechanical parameters based on laboratory or field estimation, and engineering judgment. Joint normal and shear stiffnesses are two popular input parameters to describe discontinuities in rock, which do not have specific guidelines for their estimation in literature. This study attempts to provide simple methods to estimate joint normal and shear stiffnesses in the laboratory using the uniaxial compression and small-scale direct shear tests. Samples have been prepared using rocks procured from different depths, geographical locations and formations. The study uses a mixture of relatively smooth natural joints and saw-cut joints in the various rock samples tested. The results indicate acceptable levels of uncertainty in the calculation of the stiffness parameters and provide a database of good first estimates and empirical relations which can be used for calculating values for joint stiffnesses when laboratory estimation is not possible. Joint basic friction angles have also been estimated as by-products in the small scale direct shear tests.     

Effect of antecedent-hydrological conditions on rainfall triggering of debris flows in ash-fall pyroclastic mantled slopes of Campania (southern Italy)

Mountainous areas surrounding the Campanian Plain and the Somma-Vesuvius volcano (southern Italy) are among the most risky areas of Italy due to the repeated occurrence of rainfall-induced debris flows along ash-fall pyroclastic soil-mantled slopes. In this geomorphological framework, rainfall patterns, hydrological processes taking place within multi-layered ash-fall pyroclastic deposits and soil antecedent moisture status are the principal factors to be taken into account to assess triggering rainfall conditions and the related hazard. This paper presents the outcomes of an experimental study based on integrated analyses consisting of the reconstruction of physical models of landslides, in situ hydrological monitoring, and hydrological and slope stability modeling, carried out on four representative source areas of debris flows that occurred in May 1998 in the Sarno Mountain Range. The hydrological monitoring was carried out during 2011 using nests of tensiometers and Watermark pressure head sensors and also through a rainfall and air temperature recording station. Time series of measured pressure head were used to calibrate a hydrological numerical model of the pyroclastic soil mantle for 2011, which was re-run for a 12-year period beginning in 2000, given the availability of rainfall and air temperature monitoring data. Such an approach allowed us to reconstruct the regime of pressure head at a daily time scale for a long period, which is representative of about 11 hydrologic years with different meteorological conditions. Based on this simulated time series, average winter and summer hydrological conditions were chosen to carry out hydrological and stability modeling of sample slopes and to identify Intensity-Duration rainfall thresholds by a deterministic approach. Among principal results, the opposing winter and summer antecedent pressure head (soil moisture) conditions were found to exert a significant control on intensity and duration of rainfall triggering events. Going from winter to summer conditions requires a strong increase of intensity and/or duration to induce landslides. The results identify an approach to account for different hazard conditions related to seasonality of hydrological processes inside the ash-fall pyroclastic soil mantle. Moreover, they highlight another important factor of uncertainty that potentially affects rainfall thresholds triggering shallow landslides reconstructed by empirical approaches.