Ocean-floor hydrothermal metamorphism in the Limousin ophiolites (western French Massif Central): evidence of a rare preserved Variscan oceanic marker

The Limousin ophiolite is located at the suture zone between two major thrust sheets in the western French Massif Central. This ophiolitic section comprises mantle‐harzburgite, mantle‐dunite, wehrlites, troctolites and layered gabbros. It has recorded a static metamorphic event transforming the gabbros into undeformed amphibolites and the magmatic ultramafites into serpentinites and/or pargasite‐bearing chloritites. With various thermobarometric methods, it is possible to show that the different varieties of amphibole have registered low‐P (c. 0.2 GPa) conditions with temperature ranging from high‐T, late‐magmatic conditions to greenschist–zeolite metamorphic facies. The abundance of undeformed metamorphic rocks (which is typical of the lower oceanic crust), the occurrence of Ca–Al (–Mg) metasomatism illustrated by the growth of Ca–Al silicates in veins or replacing the primary magmatic minerals, the P–T conditions of the metamorphism and the numerous similarities with oceanic crustal rocks from Ocean Drilling Program and worldwide ophiolites are the main arguments for an ocean‐floor hydrothermal metamorphism in the vicinity of a palaeo‐ridge. Among the West‐European Variscan ophiolites, the Limousin ophiolites constitute an extremely rare occurrence that has not been involved in any HP (subduction‐related) or MP (orogenic) metamorphism as observed in other ophiolite occurrences (i.e. France, Spain and Germany).     

P-T-t-deformation paths recorded by kinzigites during diapirism in the western Variscan belt (Golfe du Morbihan, southern Brittany, France)

Abstract The high-grade metamorphic rocks of southern Brittany underwent a complex tectonic evolution under various P-T conditions (high-P, high-T), related to stacking of nappes during Palaeozoic continentcontinent collision. The east to west thrusting observed in the whole belt is strongly perturbed by vertical movements attributed to the ascent of anatectic granites in the high-T area. The field reconstruction of subvertical, closed elliptical structures in gneisses and migmatites, associated with the subhorizontal, doubly radial pattern of stretching lineation in the mica schists, suggests the existence of an elliptical diapiric body buried at depth beneath the present erosion level. Deformation is associated with a complex P-T evolution partly recorded in aluminous gneisses (kinzigites, e.g. morbihanites). A chronology of successive episodes of mineral growth at different compositions is established by detailed studies of the mineral-microstructure relationships in X-Z sections, using the deformation-partitioning concept (low- and high-strain zones). Several thermometric and barometric calibrations are applied to mineral pairs either in contact or not in contact but in equivalent microstructiiral positions with respect to the deformation history. This methodology provides a continuous microstructural control of P-T variations through time and leads to three P-T-t-d paths constructed from numerous successive P-T estimations. Path 1 is a clockwise retrograde path preserved in low-strain zones, which records general exhumation movements after crustal thickening. Paths 2 and 3 are clockwise prograde/retrograde paths from high-strain zones; they are interpreted and discussed in the light of models of crustal anatexis and upward movement of magma (diapirism). Deformation and P-T effects induced by diapirism can be distinguished from the general deformation-metamorphic history of a belt, and would seem to be produced during a late stage of its history. The present microstructural-petrological approach to defining successive mineral equilibria in relation to progressive deformation steps provides a far more accurate evaluation of the metamorphic evolution than is possible by ‘standard’thermobarometry.     

Crustal structure at the easternmost termination of the Variscan belt based on CELEBRATION 2000 and ALP 2002 data

The eastern margin of the Variscan belt in Europe comprises plate boundaries between continental blocks and terranes formed during different tectonic events. The crustal structure of that complicated area was studied using the data of the international refraction experiments CELEBRATION 2000 and ALP 2002. The seismic data were acquired along SW–NE oriented refraction and wide-angle reflection profiles CEL10 and ALP04 starting in the Eastern Alps, passing through the Moravo-Silesian zone of the Bohemian Massif and the Fore-Sudetic Monocline, and terminating in the TESZ in Poland. The data were interpreted by seismic tomographic inversion and by 2-D trial-and-error forward modelling of the P waves. Velocity models determine different types of the crust–mantle transition, reflecting variable crustal thickness and delimiting contacts of tectonic units in depth. In the Alpine area, few km thick LVZ with the Vp of 5.1 km s^− 1 dipping to the SW and outcropping at the surface represents the Molasse and Helvetic Flysch sediments overthrust by the Northern Calcareous Alps with higher velocities. In the Bohemian Massif, lower velocities in the range of 5.0–5.6 km s^− 1 down to a depth of 5 km might represent the SE termination of the Elbe Fault Zone. The Fore-Sudetic Monocline and the TESZ are covered by sediments with the velocities in the range of 3.6–5.5 km s^− 1 to the maximum depth of 15 km beneath the Mid-Polish Trough. The Moho in the Eastern Alps is dipping to the SW reaching the depth of 43–45 km. The lower crust at the eastern margin of the Bohemian Massif is characterized by elevated velocities and high Vp gradient, which seems to be a characteristic feature of the Moravo-Silesian. Slightly different properties in the Moravian and Silesian units might be attributed to varying distances of the profile from the Moldanubian Thrust front as well as a different type of contact of the Brunia with the Moldanubian and its northern root sector. The Moho beneath the Fore-Sudetic Monocline is the most pronounced and is interpreted as the first-order discontinuity at a depth of 30 km.     

Evaluation of approaches to calculate debris-flow parameters for hazard assessment

Many different runout prediction methods can be applied to estimate the mobility of future debris flows during hazard assessment. The present article reviews the empirical, analytical, simple flow routing and numerical techniques. All these techniques were applied to back-calculate a debris flow, which occurred in 1982 at La Guingueta catchment, in the Eastern Pyrenees. A sensitivity analysis of input parameters was carried out, while special attention was paid to the influence of rheological parameters. We used the Voellmy fluid rheology for our analytical and numerical modelling, since this flow resistance law coincided best with field observations. The simulation results indicated that the “basal” friction coefficients rather affect the runout distance, while the “turbulence” terms mainly influence flow velocity. A comparison of the velocity computed on the fan showed that the analytical model calculated values similar to the numerical ones. The values of our rheological parameters calibrated at La Guingueta agree with data back-calculated for other debris flows. Empirical relationships represent another method to estimate total runout distance. The results confirmed that they contain an important uncertainty and they are strictly valid only for the conditions, which were the basis for their development. With regards to the simple flow routing algorithm, this methods could satisfactorily simulate the total area affected by the 1982 debris flow, but it was not able to directly calculate total runout distance and velocity. Finally, a suggestion on how different runout prediction methods can be applied to generate debris-flow hazard maps is presented. Taking into account the definition of hazard and intensity, the best choice would be to divide the resulting hazard maps into two types: “final hazard maps” and “preliminary hazard maps”. Only the use of numerical models provided final hazard maps, because they could incorporate different event magnitudes and they supplied output-values for intensity calculation. In contrast, empirical relationships and flow routing algorithms, or a combination of both, could be applied to create preliminary hazard maps. The present study only focussed on runout prediction methods. Other necessary tasks to complete the hazard assessment can be looked up in the “Guidelines for landslide susceptibility, hazard and risk zoning” included in this Special Issue.     

Variscan basic dykes in the Pelagonian (Northern Greece and south FYROM): Geodynamic significance based on petrological, geochemical and geochronological studies

Dolerite dykes intruding Variscan plutonites were studied in terms of mineralogy, petrology, geochemistry and geochronology. The main mineral constituents were studied and the sequence of crystallization has been derived. The geochemical characteristic indicate mantle origin of the dolerites and magma sources different from the hosting granitoids. From SHRIMP analyses of five spots on four different zircon crystals, resulted a 292.0±4.1 Ma age that is interpreted as the time of crystallization of the dolerite. The hosting granitoids are probably the result of mixing between two possible end-members: enriched mantle and acid metaigneous or lower crustal metasediments.

The Variscan age of the dolerites, in combination with the geochemical characteristics, indicated that the enriched mantle basaltic material should be the source of the dolerite veins. These mantle-derived basaltic melts may represent the underplated material, which probably provided the necessary thermal input to the dehydration melting in the lower crust. The dolerites should have intruded the newly formed batholiths before or at the first stages of their uplift, recording the last events of the Variscan subduction.     

The Meso-Cenozoic thermo-tectonic evolution of the Eastern Pyrenees: an <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar fission track and (U–Th)/He thermochronological study of the Canigou and Mont-Louis massifs

Vertical displacements on the SW–NE Têt fault (Eastern Pyrenees Axial Zone, France), which separates the Variscan Canigou-Carança and Mont-Louis massifs, were constrained using a thermochronologic multi-method approach. ⁴⁰Ar/³⁹Ar data from the granitic Mont-Louis massif record its Variscan cooling history and reveal no ages younger than Early Cretaceous, while the Canigou-Carança gneiss massif records systematically younger ⁴⁰Ar/³⁹Ar ages. These younger ⁴⁰Ar/³⁹Ar ages in the Canigou-Carança gneiss massif are the result of partial to total rejuvenation of argon isotopic systems related to a thermal flow coeval with the Cretaceous HT-BP metamorphism in the North Pyrenean Zone. Only the deepest rocks from the Canigou-Carança suffered this extensive Mid-Cretaceous thermal overprint probably due to differential burial around 4 km at that time. The post Mid-Cretaceous vertical displacements along the Têt fault are recorded by “low” temperature thermochronology using K-feldspar ⁴⁰Ar/³⁹Ar, zircon and apatite fission track and (U–Th)/He datings. The Mont-Louis granite samples experienced a long period of protracted cooling reflecting a lack of thermo-tectonic activity in this area from Late Palaeozoic to Early Cenozoic, followed by cooling from 55–60 Ma to Late Eocene at a mean rate of 15–20°C/Ma in the final stage. This cooling stage corresponds to Têt fault reactivation with a reversed component, promoting exhumation of the Mont-Louis roof zone contemporaneously with the south-vergent Pyrenean thrusting. In the Canigou-Carança massif, the main cooling event occurred from 32 to 18 Ma at a maximum rate of 30°C/Ma during Early Oligocene followed by a more moderate rate of 3°C/Ma from Late Oligocene to Early Burdigalian, coeval with the normal reactivation of the Têt fault in brittle conditions that accommodated the final exhumation of the massif during the opening of the Gulf of Lion.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

The age of the Kagenfels granite (northern Vosges) and its bearing on the intrusion scheme of late Variscan granitoids

In the Saxothuringian part of the Vosges (France), a first series of Variscan plutonic rocks (diorites to granites) has been intruded by several younger granites. Rocks of both the older generations have been cross-cut by the late orogenic Kagenfels granite. The averages of the hitherto published mineral ages of the earlier rock generations are 331 and 334 Ma, respectively, whereas Rb-Sr and K-Ar dates around 290 Ma have been reported for the Kagenfels granite. Because of the unlikely large age hiatus, a redetermination of the intrusion age of the Kagenfels granite formation appeared to be irrevocable. The newly obtained mineral ages on the Kagenfels granite (K-Ar and ⁴⁰Ar/³⁹Ar biotite ages as well as single zircon radiogenic ²⁰⁷Pb/²⁰⁶Pb data: 331 ± 5 Ma) are about 40 Ma older than the previous results. They are interpreted as giving the time of emplacement of the Kagenfels granite during the latest Visan. The mineral ages of the earlier plutonic rocks in this part of the Variscan Orogeny in all probability are not significantly different from their ages of intrusion. Therefore the age concordance of all three granitoid generations constrains a rather narrow time interval of orogenic magmatism close to the Lower-Upper Carboniferous boundary.     

HP–LT Variscan metamorphism in the Cubito-Moura schists (Ossa-Morena Zone, southern Iberia)

Multi-equilibrium thermobarometry shows that low-grade metapelites (Cubito-Moura schists) from the Ossa–Morena Zone underwent HP–LT metamorphism from 340–370 °C at 1.0–0.9 GPa to 400–450 °C at 0.8–0.7 GPa. These HP–LT equilibriums were reached by parageneses including white K mica, chlorite and chloritoid, which define the earliest schistosity (S₁) in these rocks. The main foliation in the schists is a crenulation cleavage (S₂), which developed during decompression from 0.8–0.7 to 0.4–0.3 GPa at increasing temperatures from 400–450 °C to 440–465 °C. Fe³⁺ in chlorite decreased greatly during prograde metamorphism from molar fractions of 0.4 determined in syn-S₁ chlorites down to 0.1 in syn-S₂ chlorites. These new data add to previous findings of eclogites in the Moura schists indicating that a pile of allochtonous rocks situated next to the Beja-Acebuches oceanic amphibolites underwent HP–LT metamorphism during the Variscan orogeny. To cite this article: G. Booth-Rea et al., C. R. Geoscience 338 (2006).     

Aftershocks series monitoring of the September 18, 2004 M = 4.6 earthquake at the western Pyrenees: A case of reservoir-triggered seismicity?

On September 18, 2004, a 4.6 mbLg earthquake was widely felt in the region around Pamplona, at the western Pyrenees. Preliminary locations reported an epicenter less than 20 km ESE of Pamplona and close to the Itoiz reservoir, which started impounding in January 2004. The area apparently lacks of significant seismic activity in recent times. After the main shock, which was preceded by series of foreshocks reaching magnitudes of 3.3 mbLg, a dense temporal network of 13 seismic stations was deployed there to monitor the aftershocks series and to constrain the hypocentral pattern. Aftershock determinations obtained with a double-difference algorithm define a narrow epicentral zone of less than 10 km², ESE–WNW oriented. The events are mainly concentrated between 3 and 9 km depth. Focal solutions were computed for the main event and 12 aftershocks including the highest secondary one of 3.8 mbLg. They show mainly normal faulting with some strike-slip component and one of the nodal planes oriented NW–SE and dipping to the NE. Cross-correlation techniques applied to detect and associate events with similar waveforms, provided up to 33 families relating the 67% of the 326 relocated aftershocks. Families show event clusters grouped by periods and migrating from NW to SE. Interestingly, the narrow epicentral zone inferred here is located less than 4 km away from the 111-m high Itoiz dam. These hypocentral results, and the correlation observed between fluctuations of the reservoir water level and the seismic activity, favour the explanation of this foreshock–aftershock series as a rapid response case of reservoir-triggered seismicity, burst by the first impoundment of the Itoiz reservoir. The region is folded and affected by shallow dipping thrusts, and the Itoiz reservoir is located on the hangingwall of a low angle southward verging thrust, which might be a case sensible to water level fluctuations. However, continued seismic monitoring in the coming years is mandatory in this area to infer more reliable seismotectonic and hazard assessments.