Cumulating processes at the crust-mantle transition zone inferred from Permian mafic-ultramafic xenoliths (Puy Beaunit,France)

Ultramafic and mafic xenoliths of magmatic origin, sampled in the Beaunit vent (northern French Massif Central), derive from the Permian (257 Ma) Beaunit layered complex (BLC) that was emplaced at the crust-mantle transition zone (∼1 GPa). These plutonic xenoliths are linked to a single fractional crystallisation process in four steps: peridotitic cumulates; websteritic cumulates; Al-rich mafic cumulates (plagioclase, pyroxenes, garnet, amphibole and spinel) and finally low-Al mafic cumulates. This sequence of cumulates can be related to the compositional evolution of hydrous Mg basaltic magma that evolved to high-Al basalt and finally to andesitic basalt. Sr and Nd isotopic compositions confirm the co-genetic character of the various magmatic xenoliths and argue for an enriched upper mantle source comparable to present mantle wedges above subduction zones. LILE, LREE and Pb enrichment are a common feature of all xenoliths and argue for an enriched sub-alkaline transitional parental magma. The existence of a Permian magma chamber at 30 km depth suggests that the low-velocity zone observed locally beneath the Moho probably does not represent an anomalous mantle but rather a sequence of mafic/ultramafic cumulates with densities close to those of mantle rocks.     

Volcanic edifice stability during cryptodome intrusion

Limit equilibrium analyses were applied to the 1980 Mount St. Helens and 1956 Bezymianny failures in order to examine the influence on stability of structural deformation produced by cryptodome emplacement. Weakening structures associated with the cryptodome include outward-dipping normal faults bounding a summit graben and a flat shear zone at the base of the bulged flank generated by lateral push of the magma. Together with the head of the magmatic body itself, these structures serve directly to localize failure along a critical surface with low stability deep within the interior of the edifice. This critical surface, with the safety coefficient reduced by 25-30%, is then very sensitive to stability condition variation, in particular to the pore-pressure ratio (r_u) and seismicity coefficient (n). For r_u=0.3, or n=0.2, the deep surface suffers catastrophic failure, removing a large volume of the edifice flank. In the case of Mount St. Helens, failure occurred within a material with angle of friction ~40°, cohesion in the range 10⁵-10⁶ Pa, and probably significant water pore pressure. On 18 May 1980, detachment of slide block I occurred along a newly formed rupture surface passing through the crest of the bulge. Although sliding of block I may have been helped by the basal shear zone, significant pore pressure and a triggering earthquake were required (r_u=0.3 and n=0.2). Detachment of the second block was guided by the summit normal fault, the front of the cryptodome, and the basal shear zone. This occurred along a deep critical surface, which was on the verge of failure even before the 18 May 1980 earthquake. The stability of equivalent surfaces at Bezymianny Volcano appears significantly higher. Thus, although magma had already reached the surface, weaker materials, or higher pore pressure and/or seismic conditions were probably required to reach the rupture threshold. From our analysis, we find that deep-seated sector collapses formed by removing the edifice summit cannot generally result from a single slide. Cryptodome-induced deformation does, however, provide a deep potential slip surface. As previously thought, it may assist deep-seated sector collapse because it favors multiple retrogressive slides. This leads to explosive depressurization of the magmatic and hydrothermal systems, which undermines the edifice summit and produces secondary collapses and explosive blasts.     

Constraints in 3D gravity inversion

A three-dimensional (3D) inversion program is developed to interpret gravity data using a selection of constraints. This selection includes minimum distance, flatness, smoothness and compactness constraints, which can be combined using a Lagrangian formulation. A multigrid technique is also implemented to resolve separately large and short gravity wavelengths. The subsurface in the survey area is divided into rectangular prismatic blocks and the problem is solved by calculating the model parameters, i.e. the densities of each block. Weights are given to each block depending on depth, a priori information on density and the density range allowed for the region under investigation. The present computer code is tested on modelled data for a dipping dike and multiple bodies. Results combining different constraints and a weight depending on depth are shown for the dipping dike. The advantages and behaviour of each method are compared in the 3D reconstruction. Recovery of geometry (depth, size) and density distribution of the original model is dependent on the set of constraints used. From experimentation, the best combination of constraints for multiple bodies seems to be flatness and a minimum volume for the multiple bodies. The inversion method is tested on real gravity data from the Rouyn-Noranda (Quebec) mining camp. The 3D inversion model for the first 10 km is in agreement with the known major lithological contacts at the surface; it enables the determination of the geometry of plutons and intrusive rocks at depth.     

Strain partitioning in a pluton during emplacement in transpressional regime: the example of the Néouvielle granite (Pyrenees)

This microstructural and anisotropy of magnetic susceptibility study of the internal structures of the Hercynian Néouvielle granite pluton (100 km²) provides new data indicating that the pluton was emplaced during the main Hercynian tectonic event recognized in the Pyrenees. It also provides new data about the later Alpine deformation localized along narrow mylonitic bands. These bands acted as reverse faults and have not rotated the Hercynian structures that define the main part of the pluton. The pluton is composed of two structural domains: the northern half of the pluton displays a beak shape in map view, with subhorizontal E-W trending lineations of magmatic origin; the southern half is semi-circular and displays rather steeply northward plunging lineations corresponding to magmatic and high temperature (HT) solid-state microstructures. These features are associated with magma deformation during emplacement. Magma deformation corresponds, in the northern half of the pluton, to an E-W strike-slip deformation recognized in the enveloping pelitic metasediments of Carboniferous age and, in the southern half of the pluton, to southward overthrusting recognized in the enveloping quartzites of Devonian age. Juxtaposition in a single granite body of transcurrent and compressive domains is viewed as a strain partitioning in the magma. This strain partitioning is linked to both the transpressive character of the main regional deformation event and the rheological contrast between the pelitic country rocks and quartzose country rocks.     

Empirical relationships between diversity of invertebrate communities and altitude in rivers: Application to biomonitoring

The relationships between the diversity of invertebrate communities and the altitude of sampling sites were analysed in 438 benthic samples, collected between 1982 and 1991, in 56 rivers of western Switzerland. Diversity, estimated from total number of taxa (genus or family) and from number of taxa intolerant of pollution, was positively correlated with increasing altitude. In contrast, density of human population and the level of organic pollution were negatively correlated with increasing altitude. Therefore, the upstream increase of invertebrate diversity was attributed to the decrease of human population which is the main source of organic pollution. In this study, altitude was used, instead of organic pollution, to predict diversity. Empirical relationships between diversity and altitude were applied to surveys of water quality to describe the general altitudinal pattern characteristic for each region and to single out anomalous sites and rivers. In addition, changes in the altitudinal patterns of diversity can be used to monitor the recovery of rivers from pollution.     

The use of transplanted cultured tropical oysters (Saccostrea commercialis) to monitor Cd levels in North Queensland coastal waters (Australia)

Bivalves are commonly used to detect metal pollution in the marine environment. Commercially cultured Milky oysters (Saccostrea commercialis) were transplanted in various sites along the North Queensland coast and analyzed for two metals of potentially anthropogenic origin (Cd, Zn). To provide additional information, naturally occurring Black Lip oysters (Saccostrea echinata) were also collected at the transplantation sites. The study demonstrated that the oysters species transplanted are good bioindicators of these metal concentrations in tropical waters, sensitive to variations in the environment at concentrations which are much smaller than pollution signals commonly reported for temperate waters. Three transplant experiments were carried out from May 1999 to February 2000. Milky oysters transplanted to offshore areas (Orpheus Is., Kelso Reef) accumulated Cd in the soft parts whereas oysters sampled from cages placed in Ross Creek and the Herbert River estuaries showed a decrease in Cd concentration, which resulted from an increase in dry weight. Dry weight appeared to be an important covariant affecting Cd concentration in the oysters whereas it does not unambiguously affect Zn concentrations. For the duration of the experiments, oysters sampled from the Magnetic Is. reference site showed effectively constant Cd concentrations and total Cd contents which indicates that any seasonal cycle affecting metal concentration is weak. It was found that Cd accumulation in oysters increased as ambient dissolved Cd concentration decreased, from which it was concluded that for these oysters, the predominant source of Cd was from the particulate phase rather than the dissolved phase.     

Mid-Holocene to present-day evolution of the Indian monsoon in transient global simulations

The low-frequency evolution of Indian rainfall mean-state and associated interannual-to-decadal variability is discussed for the last 6000 years from a multi-configuration ensemble of fully coupled global transient simulations. This period is marked by a shift of Indian Summer Monsoon Rainfall (ISMR) distribution towards drier conditions, including extremes, and a contraction of the rainy season. The drying is larger in simulations with higher horizontal resolution of the atmosphere and revised land surface hydrology. Vegetation–climate interactions and the way runoff is routed to ocean modulate the timing of the monsoon onset but have negligible effects on the evolution of seasonal rainfall amounts in our modeling framework in which carbon cycling is always active. This drying trend is accompanied by changes in ISMR interannual-to-decadal variability decreasing over north and south India but increasing over central India (20°–25° N). The ISMR interannual-to-decadal variability is decomposed into six physically consistent regimes using a clustering technique to further characterize its changes and associated teleconnections. From 6 to 3.8 kyr bp, the century-to-century modulations in the frequency of occurrence associated to the regimes are asynchronous between the simulations. Orbitally-driven trends can only be detected for two regimes over the whole 6–0 kyr bp period. These two regimes reflect increased influence of ENSO on both ISMR and Indian Ocean Dipole as the inter-hemispheric energy gradient weakens. Severe long-term droughts are also shown to be a combination of long-term drying and internally generated low-frequency modulations of the interannual-to-decadal variability.