Large-scale DNA sequencing in taxonomy and conservation : a case study with the mayfly family Heptageniidae (Ephemeroptera) in the Alps and Madagascar

Les larves aquatiques d'éphémères (Ephemeroptera) colonisent toutes les eaux douces du monde et sont couramment utilisées comme bio-indicateurs de la qualité de l'eau. Le genre Rhithrogena (Heptageniidae) est le deuxième plus diversifié chez les éphémères, et plusieurs espèces européennes ont une distribution restreinte dans des environnements alpins sensibles. Les espèces de Rhithrogena ont été classées en "groupes d'espèces" faciles à identifier. Cependant, malgré leur importance écologique et en terme de conservation, beaucoup d'espèces présentent des différences morphologiques ambiguës, suggérant que lataxonomie actuelle ne refléterait pas correctement leur diversité évolutive. De plus, aucune information sur leurs relations, leur origine, le taux de spéciation ou les mécanismes ayant provoqué leur remarquable diversification dans les Alpes n'est disponible. Nous avons d'abord examiné le statut spécifique d'environ 50% des espèces européennes de Rhithrogena en utilisant un large échantillonnage de populations alpines incluant 22 localités typiques, ainsi qu'une analyse basée sur le modèle général mixte de Yule et de coalescence (GMYC) appliqué à un gène mitochondrial standard (coxl) et à un gène nucléaire développé spécifiquement pour cette étude. Nous avons observé un regroupement significatif des séquences coxl en 31 espèces potentielles, et nos résultats ont fortement suggéré la présence d'espèces cryptiques et de fractionnements taxonomiques excessifs chez les Rhithrogena. Nos analyses phylogénétiques ont démontré la monophylie de quatre des six groupes d'espèces reconnus présents dans notre échantillonnage. La taxonomie ADN développée dans cette étude pose les bases d'une future révision de ce genre important mais cryptique en Europe. Puis nous avons mené une étude phylogénétique multi-gènes entre les espèces européennes de Rhithrogena. Les données provenant de trois gènes nucléaires et de deux gènes mitochondriaux ont été largement concordantes, et les relations entre les espèces bien résolues au sein de la plupart des groupes d'espèces dans une analyse combinant tous les gènes. En l'absence de points de calibration extérieurs tels que des fossiles, nous avons appliqué à nos données mitochondriales une horloge moléculaire standard pour les insectes, suggérant une origine des Rhithrogena alpins à la limite Oligocène / Miocène. Nos résultats ont montré le rôle prépondérant qu'ont joué les glaciations du quaternaire dans leur diversification, favorisant la spéciation d'au moins la moitié des espèces actuelle dans les Alpes. La biodiversité et le taux d'endémisme à Madagascar, notamment au niveau de la faune des eaux douces, sont parmi les plus extraordinaires et les plus menacés au monde. On pense que beaucoup d'espèces d'éphémères sont restreintes à un seul bassin versant (microendémisme) dans les zones forestières, ce qui les rendrait particulièrement sensibles à la réduction et à la dégradation de leur habitat. Mis à part deux espèces décrites, Afronurus matitensis et Compsoneuria josettae, les Heptageniidae sont pratiquement inconnus à Madagascar. Les deux genres ont une distribution discontinue en Afrique, à Madagascar et en Asie du Sud-Est, et leur taxonomie complexe est régulièrement révisée. L'approche standard pour comprendre leur diversité, leur endémisme et leur origine requerrait un échantillonnage étendu sur plusieurs continents et des années de travaux taxonomiques. Pour accélérer le processus, nous avons utilisé des collections de musées ainsi que des individus fraîchement collectés, et appliqué une approche combinant taxonomie ADN et phylogénie. L'analyses GMYC du gène coxl a délimité 14 espèces potentielles à Madagascar, dont 70% vraisemblablement microendémiques. Une analyse phylogénique incluant des espèces africaines et asiatiques portant sur deux gènes mitochondriaux et quatre gènes nucléaires a montré que les Heptageniidae malgaches sont monophylétiques et groupe frère des Compsoneuria africains. L'existence de cette lignée unique, ainsi qu'un taux élevé de microendémisme, mettent en évidence leur importance en terme de conservation. Nos résultats soulignent également le rôle important que peuvent jouer les collections de musées dans les études moléculaires et en conservation. - Aquatic nymphs of mayflies (Ephemeroptera) colonize all types of freshwaters throughout the world and are extensively used as bio-indicators of water quality. Rhithrogena (Heptageniidae) is the second most species-rich genus of mayflies, and several European species have restricted distributions in sensitive Alpine environments and therefore are of conservation interest. The European Rhithrogena species are arranged into "species groups" that are easily identifiable. However, despite their ecological and conservation importance, ambiguous morphological differences among many species suggest that the current taxonomy may not accurately reflect their evolutionary diversity. Moreover, no information about their relationships, origin, timing of speciation and mechanisms promoting their successful diversification in the Alps is available. We first examined the species status of ca. 50% of European Rhithrogena diversity using a widespread sampling scheme of Alpine species that included 22 type localities, general mixed Yule- coalescent (GMYC) model analysis of one standard mitochondrial (coxl) and one newly developed nuclear marker. We observed significant clustering of coxl into 31 GMYC species, and our results strongly suggest the presence of both cryptic diversity and taxonomic oversplitting in Rhithrogena. Phylogenetic analyses recovered four of the six recognized species groups in our samples as monophyletic. The DNA taxonomy developed here lays the groundwork for a future revision of this important but cryptic genus in Europe. Then we conducted a species-level, multiple-gene phylogenetic study of European Rhithrogena. Data from three nuclear and two mitochondrial loci were broadly congruent, and species-level relationships were well resolved within most species groups in a combined analysis. In the absence of external calibration points like fossils, we applied a standard insect molecular clock hypothesis to our mitochondrial data, suggesting an origin of Alpine Rhithrogena in the Oligocene / Miocene boundary. Our results highlighted the preponderant role that quaternary glaciations played in their diversification, promoting speciation of at least half of the current diversity in the Alps. Madagascar's biodiversity and endemism are among the most extraordinary and endangered in the world. This includes the island's freshwater biodiversity, although detailed knowledge of the diversity, endemism, and biogeographic origin of freshwater invertebrates is lacking. Many mayfly species are thought to be restricted to single river basins (microendemic species) in forested areas, making them particularly sensitive to habitat reduction and degradation. The Heptageniidae are practically unknown in Madagascar except for two described species, Afronurus matitensis and Compsoneuria josettae. Both genera have a disjunct distribution in Africa, Madagascar and Southeast Asia, and a complex taxonomic status still in flux. The standard approach to understanding their diversity, endemism, and origin would require extensive field sampling on several continents and years of taxonomic work. Here we circumvent this using museum collections and freshly collected individuals in a combined approach of DNA taxonomy and phylogeny. The cox/-based GMYC analysis revealed 14 putative species on Madagascar, 70% of which potentially microendemics. A phylogenetic analysis that included African and Asian species and data from two mitochondrial and four nuclear loci indicated the Malagasy Heptageniidae are monophyletic and sister to African Compsoneuria. The observed monophyly and high microendemism highlight their conservation importance. Our results also underline the important role that museum collections can play in molecular studies, especially in critically endangered biodiversity hotspots like Madagascar.

Conservation of the Eurasian Lynx (Lynx lynx) in a fragmented landscape - habitat models, dispersal and potential distribution

Abstract: The expansion of a recovering population - whether re-introduced or spontaneously returning - is shaped by (i) biological (intrinsic) factors such as the land tenure system or dispersal, (ii) the distribution and availability of resources (e.g. prey), (iii) habitat and landscape features, and (iv) human attitudes and activities. In order to develop efficient conservation and recovery strategies, we need to understand all these factors and to predict the potential distribution and explore ways to reach it. An increased number of lynx in the north-western Swiss Alps in the nineties lead to a new controversy about the return of this cat. When the large carnivores were given legal protection in many European countries, most organizations and individuals promoting their protection did not foresee the consequences. Management plans describing how to handle conflicts with large predators are needed to find a balance between "overabundance" and extinction. Wildlife and conservation biologists need to evaluate the various threats confronting populations so that adequate management decisions can be taken. I developed a GIS probability model for the lynx, based on habitat information and radio-telemetry data from the Swiss Jura Mountains, in order to predict the potential distribution of the lynx in this mountain range, which is presently only partly occupied by lynx. Three of the 18 variables tested for each square kilometre describing land use, vegetation, and topography, qualified to predict the probability of lynx presence. The resulting map was evaluated with data from dispersing subadult lynx. Young lynx that were not able to establish home ranges in what was identified as good lynx habitat did not survive their first year of independence, whereas the only one that died in good lynx habitat was illegally killed. Radio-telemetry fixes are often used as input data to calibrate habitat models. Radio-telemetry is the only way to gather accurate and unbiased data on habitat use of elusive larger terrestrial mammals. However, it is time consuming and expensive, and can therefore only be applied in limited areas. Habitat models extrapolated over large areas can in turn be problematic, as habitat characteristics and availability may change from one area to the other. I analysed the predictive power of Ecological Niche Factor Analysis (ENFA) in Switzerland with the lynx as focal species. According to my results, the optimal sampling strategy to predict species distribution in an Alpine area lacking available data would be to pool presence cells from contrasted regions (Jura Mountains, Alps), whereas in regions with a low ecological variance (Jura Mountains), only local presence cells should be used for the calibration of the model. Dispersal influences the dynamics and persistence of populations, the distribution and abundance of species, and gives the communities and ecosystems their characteristic texture in space and time. Between 1988 and 2001, the spatio-temporal behaviour of subadult Eurasian lynx in two re-introduced populations in Switzerland was studied, based on 39 juvenile lynx of which 24 were radio-tagged to understand the factors influencing dispersal. Subadults become independent from their mothers at the age of 8-11 months. No sex bias neither in the dispersal rate nor in the distance moved was detected. Lynx are conservative dispersers, compared to bear and wolf, and settled within or close to known lynx occurrences. Dispersal distances reached in the high lynx density population - shorter than those reported in other Eurasian lynx studies - are limited by habitat restriction hindering connections with neighbouring metapopulations. I postulated that high lynx density would lead to an expansion of the population and validated my predictions with data from the north-western Swiss Alps where about 1995 a strong increase in lynx abundance took place. The general hypothesis that high population density will foster the expansion of the population was not confirmed. This has consequences for the re-introduction and recovery of carnivores in a fragmented landscape. To establish a strong source population in one place might not be an optimal strategy. Rather, population nuclei should be founded in several neighbouring patches. Exchange between established neighbouring subpopulations will later on take place, as adult lynx show a higher propensity to cross barriers than subadults. To estimate the potential population size of the lynx in the Jura Mountains and to assess possible corridors between this population and adjacent areas, I adapted a habitat probability model for lynx distribution in the Jura Mountains with new environmental data and extrapolated it over the entire mountain range. The model predicts a breeding population ranging from 74-101 individuals and from 51-79 individuals when continuous habitat patches < 50 km2 are disregarded. The Jura Mountains could once be part of a metapopulation, as potential corridors exist to the adjoining areas (Alps, Vosges Mountains, and Black Forest). Monitoring of the population size, spatial expansion, and the genetic surveillance in the Jura Mountains must be continued, as the status of the population is still critical. ENFA was used to predict the potential distribution of lynx in the Alps. The resulting model divided the Alps into 37 suitable habitat patches ranging from 50 to 18,711 km2, covering a total area of about 93,600 km2. When using the range of lynx densities found in field studies in Switzerland, the Alps could host a population of 961 to 1,827 residents. The results of the cost-distance analysis revealed that all patches were within the reach of dispersing lynx, as the connection costs were in the range of dispersal cost of radio-tagged subadult lynx moving through unfavorable habitat. Thus, the whole Alps could once be considered as a metapopulation. But experience suggests that only few disperser will cross unsuitable areas and barriers. This low migration rate may seldom allow the spontaneous foundation of new populations in unsettled areas. As an alternative to natural dispersal, artificial transfer of individuals across the barriers should be considered. Wildlife biologists can play a crucial role in developing adaptive management experiments to help managers learning by trial. The case of the lynx in Switzerland is a good example of a fruitful cooperation between wildlife biologists, managers, decision makers and politician in an adaptive management process. This cooperation resulted in a Lynx Management Plan which was implemented in 2000 and updated in 2004 to give the cantons directives on how to handle lynx-related problems. This plan was put into practice e.g. in regard to translocation of lynx into unsettled areas. Résumé: L'expansion d'une population en phase de recolonisation, qu'elle soit issue de réintroductions ou d'un retour naturel dépend 1) de facteurs biologiques tels que le système social et le mode de dispersion, 2) de la distribution et la disponibilité des ressources (proies), 3) de l'habitat et des éléments du paysage, 4) de l'acceptation de l'espèce par la population locale et des activités humaines. Afin de pouvoir développer des stratégies efficaces de conservation et de favoriser la recolonisation, chacun de ces facteurs doit être pris en compte. En plus, la distribution potentielle de l'espèce doit pouvoir être déterminée et enfin, toutes les possibilités pour atteindre les objectifs, examinées. La phase de haute densité que la population de lynx a connue dans les années nonante dans le nord-ouest des Alpes suisses a donné lieu à une controverse assez vive. La protection du lynx dans de nombreux pays européens, promue par différentes organisations, a entraîné des conséquences inattendues; ces dernières montrent que tout plan de gestion doit impérativement indiquer des pistes quant à la manière de gérer les conflits, tout en trouvant un équilibre entre l'extinction et la surabondance de l'espèce. Les biologistes de la conservation et de la faune sauvage doivent pour cela évaluer les différents risques encourus par les populations de lynx, afin de pouvoir rapidement prendre les meilleuresmdécisions de gestion. Un modèle d'habitat pour le lynx, basé sur des caractéristiques de l'habitat et des données radio télémétriques collectées dans la chaîne du Jura, a été élaboré afin de prédire la distribution potentielle dans cette région, qui n'est que partiellement occupée par l'espèce. Trois des 18 variables testées, décrivant pour chaque kilomètre carré l'utilisation du sol, la végétation ainsi que la topographie, ont été retenues pour déterminer la probabilité de présence du lynx. La carte qui en résulte a été comparée aux données télémétriques de lynx subadultes en phase de dispersion. Les jeunes qui n'ont pas pu établir leur domaine vital dans l'habitat favorable prédit par le modèle n'ont pas survécu leur première année d'indépendance alors que le seul individu qui est mort dans l'habitat favorable a été braconné. Les données radio-télémétriques sont souvent utilisées pour l'étalonnage de modèles d'habitat. C'est un des seuls moyens à disposition qui permette de récolter des données non biaisées et précises sur l'occupation de l'habitat par des mammifères terrestres aux moeurs discrètes. Mais ces méthodes de- mandent un important investissement en moyens financiers et en temps et peuvent, de ce fait, n'être appliquées qu'à des zones limitées. Les modèles d'habitat sont ainsi souvent extrapolés à de grandes surfaces malgré le risque d'imprécision, qui résulte des variations des caractéristiques et de la disponibilité de l'habitat d'une zone à l'autre. Le pouvoir de prédiction de l'Analyse Ecologique de la Niche (AEN) dans les zones où les données de présence n'ont pas été prises en compte dans le calibrage du modèle a été analysée dans le cas du lynx en Suisse. D'après les résultats obtenus, la meilleure mé- thode pour prédire la distribution du lynx dans une zone alpine dépourvue d'indices de présence est de combiner des données provenant de régions contrastées (Alpes, Jura). Par contre, seules les données sur la présence locale de l'espèce doivent être utilisées pour les zones présentant une faible variance écologique tel que le Jura. La dispersion influence la dynamique et la stabilité des populations, la distribution et l'abondance des espèces et détermine les caractéristiques spatiales et temporelles des communautés vivantes et des écosystèmes. Entre 1988 et 2001, le comportement spatio-temporel de lynx eurasiens subadultes de deux populations réintroduites en Suisse a été étudié, basé sur le suivi de 39 individus juvéniles dont 24 étaient munis d'un collier émetteur, afin de déterminer les facteurs qui influencent la dispersion. Les subadultes se sont séparés de leur mère à l'âge de 8 à 11 mois. Le sexe n'a pas eu d'influence sur le nombre d'individus ayant dispersés et la distance parcourue au cours de la dispersion. Comparé à l'ours et au loup, le lynx reste très modéré dans ses mouvements de dispersion. Tous les individus ayant dispersés se sont établis à proximité ou dans des zones déjà occupées par des lynx. Les distances parcourues lors de la dispersion ont été plus courtes pour la population en phase de haute densité que celles relevées par les autres études de dispersion du lynx eurasien. Les zones d'habitat peu favorables et les barrières qui interrompent la connectivité entre les populations sont les principales entraves aux déplacements, lors de la dispersion. Dans un premier temps, nous avons fait l'hypothèse que les phases de haute densité favorisaient l'expansion des populations. Mais cette hypothèse a été infirmée par les résultats issus du suivi des lynx réalisé dans le nord-ouest des Alpes, où la population connaissait une phase de haute densité depuis 1995. Ce constat est important pour la conservation d'une population de carnivores dans un habitat fragmenté. Ainsi, instaurer une forte population source à un seul endroit n'est pas forcément la stratégie la plus judicieuse. Il est préférable d'établir des noyaux de populations dans des régions voisines où l'habitat est favorable. Des échanges entre des populations avoisinantes pourront avoir lieu par la suite car les lynx adultes sont plus enclins à franchir les barrières qui entravent leurs déplacements que les individus subadultes. Afin d'estimer la taille de la population de lynx dans le Jura et de déterminer les corridors potentiels entre cette région et les zones avoisinantes, un modèle d'habitat a été utilisé, basé sur un nouveau jeu de variables environnementales et extrapolé à l'ensemble du Jura. Le modèle prédit une population reproductrice de 74 à 101 individus et de 51 à 79 individus lorsque les surfaces d'habitat d'un seul tenant de moins de 50 km2 sont soustraites. Comme des corridors potentiels existent effectivement entre le Jura et les régions avoisinantes (Alpes, Vosges, et Forêt Noire), le Jura pourrait faire partie à l'avenir d'une métapopulation, lorsque les zones avoisinantes seront colonisées par l'espèce. La surveillance de la taille de la population, de son expansion spatiale et de sa structure génétique doit être maintenue car le statut de cette population est encore critique. L'AEN a également été utilisée pour prédire l'habitat favorable du lynx dans les Alpes. Le modèle qui en résulte divise les Alpes en 37 sous-unités d'habitat favorable dont la surface varie de 50 à 18'711 km2, pour une superficie totale de 93'600 km2. En utilisant le spectre des densités observées dans les études radio-télémétriques effectuées en Suisse, les Alpes pourraient accueillir une population de lynx résidents variant de 961 à 1'827 individus. Les résultats des analyses de connectivité montrent que les sous-unités d'habitat favorable se situent à des distances telles que le coût de la dispersion pour l'espèce est admissible. L'ensemble des Alpes pourrait donc un jour former une métapopulation. Mais l'expérience montre que très peu d'individus traverseront des habitats peu favorables et des barrières au cours de leur dispersion. Ce faible taux de migration rendra difficile toute nouvelle implantation de populations dans des zones inoccupées. Une solution alternative existe cependant : transférer artificiellement des individus d'une zone à l'autre. Les biologistes spécialistes de la faune sauvage peuvent jouer un rôle important et complémentaire pour les gestionnaires de la faune, en les aidant à mener des expériences de gestion par essai. Le cas du lynx en Suisse est un bel exemple d'une collaboration fructueuse entre biologistes de la faune sauvage, gestionnaires, organes décisionnaires et politiciens. Cette coopération a permis l'élaboration du Concept Lynx Suisse qui est entré en vigueur en 2000 et remis à jour en 2004. Ce plan donne des directives aux cantons pour appréhender la problématique du lynx. Il y a déjà eu des applications concrètes sur le terrain, notamment par des translocations d'individus dans des zones encore inoccupées.

The Paleozoic magmatic history of the Maggia and Sambuco nappes, Lower Penninic, Central Lepontine Alps

Résumé : Les corps magmatiques sont des indicateurs essentiels dans toute reconstitution paléogéographique et/ou géodynamique d'un cycle orogénique, en particulier en contexte polycyclique, où la plupart des autres indices ont été oblitérés. Ils sont aisément datables et leurs caractéristiques géochimiques permettent de contraindre leur contexte tectonique de mise en place. Cette approche a été appliquée aux socles pré-mésozoïques des nappes penniques inférieures de Sambuco et de la Maggia, dans les Alpes centrales lepontines. Plusieurs événements magmatiques ont été identifiés dans le socle de Sambuco et datés par la méthode U-Pb sur zircon couplée à la technique LA-ICPMS. La suite calco-alcaline mafique rubanée de Scheggia est datée du Cambrien inférieur à 540-530 Ma ; le métagranite alumineux oeillé de Sasso Nero a un âge de 480-470 Ma, tout comme bien d'autres «older orthogneisses» des socles alpins. Il contient des zircons hérités d'âge panafricain à 630-610 Ma, indicateur d'une affiliation gondwanienne de ces terrains. Le pluton calco-alcalin du Matorello est daté à environ 300-310 Ma, et les filons lamprophyriques qu'il abrite à 300 Ma. La granodiorite de Cocco et le leucogranite de Ruscada, tous deux intrudés dans le socle de la nappe adjacente de la Maggia, ont des âges similaires à celui du Matorello. Ceci ajouté aux similitudes magmatiques observées entre Cocco et Matorello suggère une proximité paléogéographique des deux nappes au Permien-Carbonifère. Or ces dernières sont actuellement considérées appartenir à deux domaines paléogéographiques mésozoïques distincts : helvétique pour Sambuco et briançonnais pour Maggia, séparés par un bassin océanique. Si tel fut le cas, aucun mouvement décrochant ne doit avoir décalé les marges continentales de l'océan, retrouvées en parfaite coïncidence lors de sa fermeture. Le Matorello est un pluton recristallisé en faciès amphibolite et plissé par cinq phases successives de déformation non-coaxiales, qui ont conduit à son renversement complet, attesté par des indicateurs de paléogravité. Il préserve de spectaculaires phénomènes de coexistence liquide de magmas (essaims d'enclaves et Bills composites). Ce pluton était originellement tabulaire, construit par l'accumulation de multiples injections de magma en feuillets d'épaisseur métrique à décamétrique. Suivant le rythme de mise en place, les injections successives ont rapidement cristallisé avec des contours nets et bien définis (Bills composites) ou se sont mélangées avec les précédentes pour former une couche non consolidée de plusieurs dizaines de mètres d'épaisseur (granodiorite principale). Les injections individuelles sont délimitées par de subtils contrastes en granulométrie, proportions modales ou ségrégation de minéraux (schlieren), ou par des phénomènes d'érosion le long des surfaces de contact. Deux couches métriques à contour sinueux consistent en une accumulation compacte d'enclaves mafiques arrondies dans une matrice granodioritique fine. Le granoclassement des enclaves, la présence de figures de charge et de phénomènes érosifs en base de couche, ainsi que des schlieren de biotite entrecroisés évoquent l'injection de coulées de magma chargé d'enclaves et de faible viscosité en régime hydrodynamique turbulent dans un encaissant granodioritique encore largement liquide. La nature hybride des roches implique une chambre magmatique sous-jacente, en cours de différenciation et périodiquement réalimentée. Les magmas sont des liquides mafiques dérivés du manteau et des liquides anatectiques d'origine crustale, comme l'indique la gamme mesurée des rapports isotopiques initiaux du Sr (0.704 à 0.709) et des valeurs epsilon Nd (-2.1 à -4.7). Ces données montrent également que la contribution crustale est dominante, en accord avec les isotopes du plomb. Les phénomènes d'hybridation ont vraisemblablement eu lieu en base de croûte et dans la chambre magmatique sous-jacente au laccolite du Matorello. Les indicateurs de paléogravité du Matorello contribuent accessoirement à la compréhension de l'architecture actuelle de la nappe de Sambuco. Des plis isoclinaux à surface axiale verticale peuvent être mis en évidence par le contact entre les faciès dioritique et granodioritique. L'antiforme dont le Matorello forme le coeur est un synclinal, ce qui le positionne dans le Flanc inverse du grand pli couché que forme la nappe de Sambuco. Par ailleurs, des blocs de gneiss retrouvés dans le wildflysch sommital de la couverture de la nappe d'Antigorio ont été affiliés dans cette étude au pluton du Matorello. Ceci implique que le front de la nappe de Sambuco chevauchait déjà la partie est du bassin d'Antigorio au moment de sa fermeture. Par conséquent, ce n'est qu'en position externe que la nappe du Lebendun chevauche directement la nappe d'Antigorio. Abstract Magmatic bodies are important markers in paleo-geographic or geodynamic reconstructions of orogenic cycles, even more so in the case of polycyclic events where many of the other markers have been overwritten or destroyed. Plutons are relatively easy to date and their geochemical properties help constrain the tectonic context in which they were emplaced. This study focuses on the pre-mesozoic basement in the Sambuco and Maggia lower Penninic nappes located in the central Lepontine domain of the Alps. A number of magmatic events have been identified in the Sambuco basement. These events were dated using LA-ICPMS U/Pb on zircon grains. The mafic calc-alkaline banded Scheggia suite is dated as lower Cambrian, 540-530 Ma. The Al-rich Sasso-Nero lenticular gneiss is 480-470 Ma old (similarly to many older orfhogneisses of the Alpine basement) and contains 630-610 Ma old pan-African inherited zircons that illustrate the Gondwanian origin of these terranes.The calc-alkaline Matorello pluton is dated as 310-300 Ma whereas the lamprophyric bodies it contains are of 300 Ma. The Cocco granodiorite and the Ruscada leucogranite both intrude the basement of the adjacent Maggia nappe and are of similar ages to the Matorello. The ages as well as the geochemical similarities between the Cocco, Rucada and Matorello plutons suggest their paleo-geographic proximity at the Permian-Carboniferous boundary. However, these nappes are currently considered as belonging to two different Mesozoic paleo-geographic domains. Indeed, the Sambuco is considered as Helvetic whereas the Maggia is said to be Briançonnais, both separated by an oceanic basin. If this is the case, then it is essential that nostrike-slip movement has misaligned both continental margins since these coincide perfectly now that the oceanic domain closed. The Matorello pluton was originally a tabular intrusion, built up by the accumulation of multiple, several meter-thick, subhorizontal sheet-like injections of magma. Depending on their emplacement rate, the successive magma injections either solidified rapidly with sharp and rather well-defined boundaries (like the composite sills) or mingled with previous injections generating a thick molten layer up to several tens to hundred meters thick, like in the main granodioritic facies. These coalesced injections are hardly distinguishable, however subtle contrasts in granulometry, mineral modal proportions or mineral sorting (cross-bedded biotite-rich schlieren), as well as erosional features and/or crystal entrapment along contact surfaces allow to distinguish between the different injections. Two exceptional meter-thick layers display sinuous boundaries with the host granodiorite and consist of a densely packed accumulation of mafic enclaves in a granodioritic matrix. Gravitational sorting of the enclaves with load cast features at the base of the layers and sinuous biotite schlieren point to injection of low viscosity turbulent composite magma flows in the still largely molten granodiorite host. The hybrid nature of these rocks implies the existence of á periodically replenished and differentiated underlying magma chamber. Magmas are mafic liquids derived from the mantle and anatectic liquids of crustal origin, as shown by the (87Sr/86Sr), and epsilon Nd values (0.704-0.709 and -2.1 to -4.7 respectively. These data show that the crustal contribution is important, as confirmed by the Pb isotopes. The hybridisation processes seem to have occurred in the lower crust in magma chambers underlying the Matorello laccolith. The paleo-gravity markers in the Matorello help understand the architecture of the Sambuco nappe. Isoclinal folds with a vertical axial plane can be seen at the contact between dioritic and granodioritic facies. The antiform structure of which the Matorello is the heart is in fact a syncline. This places it in the inverse flanc of the large recumbent fold that constitutes the Sambuco nappe. The gneiss blocs found in the summital wildflysh cover of the Antigorio nappe have been linked to the Matorello pluton. This means that the front of the Sambuco nappe already overlapped the Antigorio basin when it closed. This implies that the Lebendun nappe can only overlap the Antigorio nappe in it's external position. Résumé grand public La chaîne alpine est la conséquence de la collision tertiaire entre deux masses continentales, l'Europe au nord et la péninsule apulienne africaine au sud, originellement séparées par l'océan mésozoïque téthysien. Cette collision a fermé un espace large de plusieurs centaines de km avec pour résultat l'écaillage de la croûte terrestre en unités tectoniques de dimensions variables, qui se sont empilées, imbriquées, éventuellement replissées en nappes de géométrie complexe. Cet amoncellement de 40 km d'épaisseur a vu sa température et sa pression lithostatique internes augmenter jusqu'à des valeurs de l'ordre de 680 °C et 6000 bars, induisant une recristallisation métamorphique des roches. L'un des objectifs de la géologie alpine est de reconstituer la géographie de la région aux temps mésozoïques de l'océan téthysien, en d'autres termes, de replacer chacune des unités tectoniques identifiées au sein de l'empilement alpin dans sa position originelle. Le défi est de taille et peut être comparé à celui de la reconstitution d'un vaste puzzle, dont certaines pièces seraient endommagées au niveau de leur contour ou leurs couleurs (métamorphisme), dissimulées par d'autres (enfouissement), voire tombées de la table de jeu (subduction, échappement latéral). Plusieurs approches ont été mises en oeuvre au cours du siècle écoulé. On citera en particulier la stratigraphie, la tectonique et le paléomagnétisme. Dans ce travail, nous avons essentiellement utilisé des techniques de datation isotopique absolue des roches (U/Pb sur zircon) qui, sur la base des connaissances acquises par l'ensemble des autres disciplines géologiques, nous ont permis de mieux contraindre ta paléogéographie mésozoïque du domaine «pennique inférieur » des Alpes centrales lépontines. Et au-delà? Nous savons tous que la disposition des continents à la surface de la Terre évolue constamment. Il est donc tentant d'essayer de remonter plus loin encore dans le temps et de reconstituer la physionomie de la marge sud européenne, tout au moins certains éléments de son histoire, au cours de l'ère paléozoïque. Les traces de ces événements très anciens sont naturellement ténues et dans ce contexte, les techniques de datation mentionnées ci-dessus deviennent les outils les plus performants. Ainsi, des datations u/Pb sur zircon nous ont permis de recenser plusieurs intrusions magmatiques, attribuées à quatre événements orogéniques anté-alpins. Des âges néoprotérozoïques (630-610 millions d'années ou Ma), cambrien inférieur (540-530 Ma), ordovicien inférieur (480-470 Ma) et carbonifère supérieur-permien inférieur (310-285 Ma) ont été obtenus dans le socle de la nappe de Sambuco. Des âges similaires à 300 Ma ont été obtenus dans la nappe voisine de la Maggia, qui permettent de relier ces deux unités. Aujourd'hui côte à côte, ces deux nappes devaient également se trouver proches l'une de l'autre il y a 300 Ma, lors de l'extension post-varisque. Les structures magmatiques spectaculaires préservées dans le pluton du Matorello (300 Ma) contraignent la géométrie actuelle de la nappe de Sambuco dans laquelle l'intrusion s'est mise en place. La forme originelle du pluton, aujourd'hui retourné et replissé plusieurs fois, s'avère être tabulaire, faite d'intrusions de faible épaisseur (1-300 m) s'étalant en forme de disque (30m à 2 km de diamètre). Les injections successives de magma se sont accumulées sous un toit dioritique précoce; elles sont issues, par le refais de fractures, d'une chambre magmatique plus profonde, périodiquement réalimentée par des magmas calco-alcalins d'origine mantellique contaminés parla croûte continentale profonde (εNd = -2.1 à -4.7). Des accumulations d'enclaves magmatiques arrondies et granoclassées dans des paléo-chenaux à fond érosif témoignent de conditions de mise en place hydrodynamiques à haute énergie. Ces enclaves sont emmenées de la chambre magmatique sous-jacente à la faveur d'épisodes de fracturation hydraulique liés à l'injection de magmas matelliques chauds dans des liquides différenciés riches en eau. Cette hypothèse est étayée par l'existence de filons composites. Une paléohorizontale a pu être déduite au sein du pluton, indiquant que cette partie de la nappe de Sambuco est verticalisée et isoclinalement replissée par la déformation alpine. Finalement, des blocs érodés du socle Sambuco ont été retrouvés dans le wildflysch sommital de la couverture sédimentaire mésozoïque de la nappe d'Antigorio sous-jacente. Ceci suggère que les blocs ont été fournis parle front de la nappe de Sambuco en train de chevaucher sur la nappe d'Antigorio au moment de la fermeture du bassin sédimentaire de cette dernière.

The Mont Collon mafic complex (Austroalpine Dent Blanche nappe) : permian evolution of the Western European mantle

Résumé: Le complexe du Mont Collon (nappe de la Dent Blanche, Austroalpin) est l'un des exemples les mieux préservés du magmatisme mafique permien des Alpes occidentales. Il est composé d'affleurements discontinus et d'une stratification magmatique en son centre (Dents de Bertol) et est composé à 95% de roches mafiques cumulatives (gabbros à olivine et/ou cpx, anorthositiques, troctolites, wehrlites et wehrlites à plagioclase) et localement de quelques gabbros pegmatitiques. Ces faciès sont recoupés par de nombreux filons acides (aphtes, pegmatites quartziques, microgranodiorites et filons anorthositiques) et mafiques tardifs (dikes mélanocrates riches en Fe et Ti). Les calculs thermométriques (équilibre olivine-augite) montrent des températures de 1070-1120 ± 6°C, tandis que le thermomètre amphibole-plagioclase indique une température de 740 ± 40°C à 0.5 GPa pour les amphiboles magmatiques tardives. La geobarométrie sur pyroxène donne des pressions moyennes de 0.3-0.6 GPa, indiquant un emplacement dans la croûte moyenne. De plus, les températures obtenues sur des amphiboles coronitiques indiquent des températures de l'ordre de 700 ± 40°C confirmant que les réactions coronitiques apparaissent dans des conditions subsolidus. Les âges concordants U/Pb sur zircons de 284.2 ± 0.6 et 282.9 ± 0.6 Ma obtenus sur un gabbro pegmatitique et une pegmatitique quartzique, sont interprétés comme des âges de cristallisation. Les datations 40Ar/39Ar sur amphiboles des filons mélanocrates donnent un âge plateau de 260.2 ± 0.7 Ma, qui est probablement très proche de l'âge de cristallisation. Ainsi, cet age 40Ar/39Ar indique un second évènement magmatique au sein du complexe. Les compositions des roches totales en éléments majeurs et traces montrent peu de variations, ainsi que le Mg# (75-80). Les éléments traces enregistrent le caractère cumulatif des roches (anomalie positive en Eu) et révèlent des anomalies négatives systématiques en Nb, Ta, Zr, Hf et Ti dans les faciès basiques. Le manque de corrélation entre éléments majeurs et traces est caractéristique d'un processus de cristallisation in situ impliquant une quantité variable de liquide interstitiel (L) entre les phases cumulus. Les distributions des éléments traces dans les minéraux sont homogènes, indiquant une rééquilibration .subsolidus entre cristaux et liquide interstitiel. Un modèle quantitatif basé sur les équations de cristallisation in situ de Langmuir reproduisent correctement les concentrations en terres rares légères des minéraux cumulatifs montrant la présence de 0 à 35% de liquide interstitiel L pour des degrés de différenciation F de 0 à 45%, par rapport au faciès les moins évolués du complexe. En outre, les valeurs de L sont bien corrélées avec les proportions modales d'amphibole interstitielle et les concentrations en éléments incompatibles des roches (Zr, Nb). Le liquide parental calculé des cumulats du Mont Collon est caractérisé par un enrichissement relatif en terres rares légères et Th, un appauvrissement en terres rares lourdes typique d'une affinité transitionnelle (T-MORB) et une forte anomalie négative en Nb-Ta. Les roches cumulatives montrent des compositions isotopiques en Nd-Sr proches de la terre globale silicatée (BSE), soit 0.6<εNdi<+3.2, 0.7045<87Sr/86Sri<0.7056. Les rapports initiaux en Pb indiquent une source dans le manteau enrichi subcontinental lithosphérique, préalablement contaminé par des sédiments océaniques. Les dikes mélanocrates Fe-Ti sont représentatifs de liquides et ont des spectres de terres rares enrichis, une anomalie positive en Nb-Ta et des εNdi de +7, des 87Sr/86Sri de 0.703 et des rapports initiaux en Pb, similaires à ceux des basaltes d'île océanique, indiquant une source asthénosphérique modérément appauvrie. Ainsi, la fusion partielle du manteau lithosphérique subcontinental est induite par l'amincissement post-orogénique et la remontée de l'asthénosphère. Les filons mélanocrates proviennent, après délamination du manteau lithosphérique, de la fusion de l'asthénosphère. Abstract The early Permian Mont Collon mafic complex (Dent Blanche nappe, Austroalpine nappe system) is one of the best preserved examples of the Permian mafic magmatism in the Western Alps. It is composed of discontinuous exposures and a well-preserved magmatic layering (the Dents de Bertol cliff) crops out in the center part of the complex. It mainly consists of cumulative mafic rocks, which represent 95 vol-% of the mafic complex (ol- and cpx-bearing gabbros and rare anorthositic layers, troctolites, wehrlites and plagioclase-wehrlites) and locally pegmatitic gabbros. All these facies are crosscut by widespread acidic (aplites, quartz-rich pegmatites, microgranodiorites) and late mafic Fe-Ti melanocratic dikes. Olivine-augite thermometric calculations yield a range of 1070-1120 ± 6°C, while amphibole-plagioclase thermometer yields a temperature of 740 ± 40°C at 0.5 GPa. Pyroxene geobarometry points to a pressure of 0.3-0.6 GPa, indicating a middle crustal level of emplacement. Moreover, temperature calculations on the Mont Conon coronitic amphiboles indicate temperatures of 700 ± 40°C, close to those calculated for magmatic amphiboles. These temperatures confirm that coronitic reactions occurred at subsolidus conditions. ID-TIMS U/Pb zircon ages of 284.2 ± 0.6 and 282.9 ± 0.6 Ma obtained on a pegmatitic gabbro and a quartz-pegmatitic dike, respectively, were interpreted as the crystallization ages of these rocks. 40Ar/39Ar dating on amphiboles from Fe-Ti melanocratic dikes yields a plateau age of 260.2 ± 0.7 Ma, which is probably very close to the crystallization age. Consequently, this 40Ar/P39Ar age indicates a second magmatic event. Whole-rock major- and trace-element compositions show little variation across the whole intrusion and Mg-number stays within a narrow range (75-80). Trace-element concentrations record the cumulative nature of the rocks (e.g. positive Eu anomaly) and reveal systematic Nb, Ta, Zr, Hf and Ti negative anomalies for all basic facies. The lack of correlation between major and trace elements is characteristic of an in situ crystallization process involving variable amounts of interstitial liquid (L) trapped between the cumulus mineral phases. LA-ICPMS measurements show that trace-element distributions in minerals are homogeneous, pointing to subsolidus re-equilibration between crystals and interstitial melts. A quantitative modeling based on Langmuir's in situ crystallization equation successfully reproduced the Rare Earth Element (REE) concentrations in cumulitic minerals. The calculated amounts of interstitial liquid L vary between 0 and 35% for degrees of differentiation F of 0 to 45%, relative to the least evolved facies of the intrusion. Furthermore, L values are well correlated with the modal proportions of interstitial amphibole and whole-rock incompatible trace-element concentrations (e.g. Zr, Nb) of the tested samples. The calculated parental melt of the Mont Collon cumulates is characterized by a relative enrichment in Light REE and Th, a depletion in Heavy REE, typical of a transitional affinity (T-MORB), and strong negative Nb-Ta anomaly. Cumulative rocks display Nd-Sr isotopic compositions close to the BSE (-0.6 < εNdi < +3.2, 0.7045 < 87Sr/86Sri < 0.7056). Initial Pb ratios point to an origin from the melting of an enriched subcontinental lithospheric mantle source, previously contaminated at the source by oceanic sediments. The contrasted alkaline Fe-Ti melanocratic dikes are representative of liquids. They display enriched fractionated REE patterns, a positive Nb-Ta anomaly and εNdi of +7, 87Sr/86Sri of 0.703 and initial Pb ratios, all reminiscent of Ocean Island Basalt-type rocks, pointing to a moderately

Événements et déformations tardi-métamorphiques dans le segment Ossola-Ticino (Val Vigezzo-Centovalli, Italie-Suisse)

RESUME: Une zone tectonique large et complexe, connue sous le nom de ligne des Centovalli, traverse le secteur des Alpes Centrales compris entre Domodossola et Locarno. Cette région, formée par le Val Vigezzo et la vallée des Centovalli, constitue la terminaison méridionale du dôme Lepontin et représente une portion de la zone des racines des nappes alpines. Elle fait partie d’une grande et complexe zone de cisaillement, en partie associée à des phénomènes hydrothermaux d’âge alpin (<20 Ma), qui comprend le système tectonique Insubrien et celui du Simplon. Le Val Vigezzo et les Centovalli constituent un vrai carrefour entre les principaux accidents tectoniques des Alpes ainsi qu'une zone de juxtaposition du socle Sudalpin avec la zone des racines de l’Austroalpin et du Pennique. Les phases de déformation et les structures géologiques qui peuvent être étudiées s'étalent sur une période comprise entre environ 35 Ma et l'actuel. L’étude détaillée de terrain a mis en évidence la présence de nombreuses roches et structures de déformation de type ductile et cassant tels que des mylonites, des cataclasites, des pseudotachylites, des kakirites, des failles minéralisées, des gouges de faille et des plis. Sur le terrain on a pu distinguer au moins quatre générations de plis liés aux différentes phases de déformation. Le nombre et la complexité de ces structures indiquent une histoire très compliquée, selon plusieurs étapes distinctes, parfois liées, voire même superposées. Une partie de ces structures de déformation affectent aussi les dépôts sédimentaires d’âge quaternaire, notamment des limons et des sables lacustres. Ces sédiments constituent les restes d'un bassin lacustre attribué à l'époque interglaciaire Riss/Würm (éemien, 67.000-120.000 ans) et ils affleurent dans la partie centrale de la zone étudiée, à l'Est de la plaine de Santa Maria Maggiore. Ces sédiments montrent en leur sein toute une série de structures de déformation tels que des plans de faille inverses, des structures conjuguées de raccourcissement et des véritables plis. Ces failles et ces plis représenteraient les évidences de surface d’une déformation probablement active en époque quaternaire. Une autre formation rocheuse a retenu tout notre attention; il s'agit d'un corps de brèches péridotitiques monogéniques qui affleure en discontinuité le long du versant méridional et le long du fond de la vallée Vigezzo sur environ 20 km. Ces brèches se posent indifféremment sur le socle (unités Finero, Orselina) ou sur les sédiments lacustres. Elles sont traversées par des plans de failles qui développent des véritables stries de faille et des gouges de faille; l’orientation de ces plans est la même que celle affectant les failles à gouges du socle. La genèse de cette brèche est liée à l'altération et au modelage glacier (rock-glaciers) d'une brèche tectonique originelle qui borde la partie externe du Corps de Finero. Les structures de déformation de cette brèche, pareillement à celles des sédiments lacustres, ont été considérées comme les évidences de surface d'une tectonique quaternaire active dans la région. La dernière phase de déformation cassante qui affecte cette région peut donc être considérée comme active en époque quaternaire. Une vue d’ensemble de la région étudiée nous permet de reconnaître à l’échelle régionale une zone de cisaillement complexe orientée E-W, parallèlement à l’axe de la vallée Centovalli-Val Vigezzo. Les données de terrain, indiquent que cette zone de cisaillement débute sous conditions ductiles et évolue en plusieurs étapes jusqu’à des conditions de failles cassantes de surface. La reconstruction de l'évolution géodynamique de la région a permis de définir trois étapes distinctes qui marquent le passage, de ce secteur de socle cristallin, de conditions P-T profondes à des conditions de surface. Dans ce contexte, on a reconnu trois phases principales de déformation à l’échelle régionale qui caractérisent ces trois étapes. La phase la plus ancienne est constituée par des mylonites en faciès amphibolite, associées à des mouvements de cisaillement dextre, qui sont ensuite remplacés par des mylonites en faciès schistes verts et des plis rétrovergentes liés au rétrocharriage des nappes alpines. Une deuxième étape est identifiée par le développement d’une phase hydrothermale liée à un système de failles extensives et décrochantes dextres à direction principale E-W, NE-SW et NW-SE. Leur caractérisation minéralogique a permis la mise en évidence des phases cristallines de néoformation liées à cet événement constituées par : K-feldspath (microcline), chlorites (Fe+Mg), épidotes, prehnite, zéolites (laumontite), sphène, calcite. Dans ce contexte, pour obtenir une meilleure caractérisation de cet événement hydrothermal on a utilisé des géothermomètres sur chlorites, sensible aussi à la pression et a la a(H2O), qui ont donné des valeurs descendantes comprises entre 450-200°C. Les derniers mouvements sont mis en évidence par le développement d’une série de plans majeurs de failles à gouge, qui forment une structure en sigmoïdes d’épaisseur kilométrique reconnaissable à l’échelle de la vallée et caractérisée par des mouvements transpressifs avec une composante décrochante dextre toujours importante. Cette phase de déformation forme un système conjugué de failles avec direction moyenne E-W qui coupent la zone des racines des nappes alpines, la zone du Canavese et le corps ultramafique de Finero. Ce système se déroule de manière subparallèle à l'axe de la vallée le long de plusieurs dizaines de kilomètres. Une analyse complète et détaillée des gouges de faille par XRD a montré que la fraction argileuse (<2 µm) de ces gouges contient une partie de néoformation très importante constituée par, des illites, des chlorites et des interstratifiés de type illite/smectite ou chlorite/smectite. Des datations avec méthode K-Ar sur ces illites ont donné des valeurs comprises entre 12 et 4 Ma qui représentent l'âge de cette dernière déformation cassante. L'application de la méthode de la cristallinité de l'illite (C.I.) a permis d'évaluer les conditions thermiques qui caractérisent le déroulement de cette dernière phase tectonique qui se produit sous conditions de température caractéristiques de l'anchizone et de la diagenèse. L'ensemble des structures de déformation qu'on vient de décrire s'insère parfaitement dans le contexte de convergence oblique entre la plaque adriatique et celle européenne qui à produit l'orogène alpin. On peut considérer les structures tectoniques du Val Vigezzo-Centovalli comme l'expression d'une zone majeure de cisaillement "Simplo-Insubrienne". L'empilement structural et les structures tectoniques affleurantes dans la région sont le résultat de l'interaction entre un régime tectonique transpressif et un régime transtensif. Ces deux champs de tension sont antagonistes entre eux mais sont reliés, de toute façon, à une seule phase décrochante dextre principale, due à une convergence oblique entre deux plaques. À l'échelle de l'évolution géodynamique on peut distinguer différentes étapes au cours desquelles les structures de ces deux régimes tectoniques interagissent en manière différente. En accord avec les données géophysiques et les reconstructions paléodynamiques prises dans la littérature on considère que la ligne Rhône-Simplon-Centovalli représente l'évidence de surface de la suture majeure profonde entre la plaque Adriatique et celle Européenne. Les vitesses de soulèvement qui ont été calculées dans cette étude pour cette région des Alpes donnent une valeur moyenne de 0.8 mm/a qui est tout à fait comparable avec les données proposées par la littérature sur cette zone. La zone Val Vigezzo-Centovalli peut être donc considérée comme un carrefour géologique où se croisent différentes phases tectoniques qui représentent les évidences de surface d'une suture profonde majeure entre deux plaques dans un contexte de collision continentale. ABSTRACT: A wide and complex tectonic zone known as Centovalli line, crosses the Central Alps sector between Domodossola and Locarno. This area, formed by the Vigezzo Valley and Centovalli valley, constitutes the southernmost termination of the Lepontin dome and represents a portion of the alpine nappes root zone. It belongs to a large and complex shear-zone, partly associated with hydrothermal phenomena of alpine age (<20 My), which includes the Insubric Line and the Simplon fault zone. Vigezzo Valley and Centovalli constitute a real crossroads between the mains alpines tectonics lines as well as a zone of juxtaposition of the Southalpine basement with the Austroalpin and Pennique root zone. The deformation phases and the geological structures that can be studied between approximately 35 My and the present. The detailed field study showed the presence of many brittle and ductile deformation structures and fault rocks such as mylonites, cataclasites, pseudotachylites, kakirites, mineralized faults, fault gouges and folds. In the field we could distinguish at least four folds generations related to the various deformation phases. The number and the complexity of these structures indicate a very complicated history, comprising several different stages, that sometimes are related and even superimposed. Part of these deformation structures affect also the sedimentary deposits of quaternary age, in particular the silts and sands lake deposit. These sediments constitute the remainders of a lake basin ascribed to the interglacial Riss/Würm (Eemien, 67.000-120.000 years) and outcroping in the central part of the studied area, in the Eastern part of Santa Maria Maggiore plain. These sediments show a whole series of deformation structures such as inverse fault planes, combined shortening structures and true folds. These faults and folds would represent the surface evidence of a probably active tectonic deformation in quaternary time. Another rock formation attracted all our attention. It is a body of monogenic peridotite breccia which outcrops in discontinuity along the southernmost slope and the bottom of the Vigezzo valley on approximately 20 km. This breccia lies indifferently on the basement (Finero and Orselina units) or on the lake sediments. They are crossed by fault planes which developed slikenside and fault gouges whose orientation is the same of the faults gouges in the alpine basement. This breccia results from the weathering and the surface modelling of an original tectonic breccia which borders the external part of Finero peridotite body. This breccia deformation structures, like those of the lake sediments, were regarded as the surface interaction of active quaternary tectonics in the area. So the last brittle deformation phases which affects this area seems to be actives in quaternary time. Theoverall picture of the studied area on a regional scale enables us to point out a complex shear-zone directed E-W, parallel to the axis of the Centovalli and Vigezzo Valley. The field analysis indicates that this shear-zone began under ductile conditions and evolved in several stages to brittle faulting under surface conditions. The analysis of the geodynamic evolution of the area allows to define three different stages which mark the transition of this alpine basement root zone, from deep P-T conditions to P-T surface conditions. In this context on regional scale three principal deformation phases, which characterize these three stages can be distinguished. The oldest phase consisted of the amphibolitie facies mylonites, associated to dextral strikeslip movements. They are then replaced by green-schists facies mylonites and backfolds related to the backthrusting of the alpines nappes. A second episode is caracterized by the development of an hydrothermal phase bound to an extensive fault and dextral strike-slip fault system, with E-W, NW-SE and SE-NW principal directionsThe principal neoformed mineral phases related to this event are: K-feldspar (microcline), chlorites (Fe+Mg), epidotes prehnite, zéolites (laumontite), sphene and calcite. In this context, to obtain a better characterization of this hydrothermal event, we have used an chlorite geothermometer, sensitive also to the pressure and has the a(H2O), which gave downward values ranging between 450-200°C. The last movements are caracterized by the development of important gouge fault plans, which form a sigmoid structure of kilometric thickness which is recognizable at the valley scale, and is characterized by transpressive movements always with a significant dextral strike-slip component. This deformation phase forms a combined faults system with an average E-W direction, which cuts trough the alpine root zone, the Canavese zone and the Finero ultramafic body. This fault system takes place subparallel to the axis of the valley over several tens of kilometers. A complete and detailed XRD analysis of the gouges fault showed that the clay fraction (<2µm) contains a very significant neo-formation of illite, chlorites and mixed layered clays such as illite/smectite or chlorite/smectite. The K-Ar datings of the illite fraction <2µm gave values ranging between 12 and 4 My and the illite fraction <0.2µm gave more recents values until to 2,4-0 My.This values represent the age of this last brittle deformation. The application of the illite crystallinity method (C.I.) allowed evaluating the thermal conditions which characterize this tectonic phase that occured under temperature conditions of the anchizone and diagenesis. The whole set of deformation structures which we just described, perfectly fit the context of oblique convergence between the Adriatic and the European plate that produced the alpine orogen. We can regard the Vigezzo valley and Centovalli tectonic structures as the expression of a major "Simplo-Insubric" shear-zone. Structural stacking and tectonic structures that outcrop in the studied area, are the result of the interaction between a transpressive and a transtensve tectonic phases. These two tension fields are antagonistic but they are also connected, in any event, with only one principal dextral strike-slip movement, caused by an oblique convergence between two plates. On the geodynamic evolution scale we can distinguish various stages during which these two tectonic structures fields interact in various ways. In agreement with the geophysical data and the paleodynamic recostructions taken in the literature we considers that the Rhone-Simplon-Centovalli line are the surface feature of the major collision between the Adriatique and the European plate at depth. The uplift speeds we calculated in this study for this Alpine area give an average value of 0.8 mm/a, which is in good agreement with the data suggested by the literature on this zone. TheVigezzo Valley and Centovalli zone can therefore be regarded as a geological crossroad where various tectonic phases are superimposed. They represent the evidences of a major and deeper suture between two plates in a continental collision context.

Linking landscape- and population ecology for large population management modelling: the case of Ibex (capra ibex) in Switzerland

La biologie de la conservation est communément associée à la protection de petites populations menacées d?extinction. Pourtant, il peut également être nécessaire de soumettre à gestion des populations surabondantes ou susceptibles d?une trop grande expansion, dans le but de prévenir les effets néfastes de la surpopulation. Du fait des différences tant quantitatives que qualitatives entre protection des petites populations et contrôle des grandes, il est nécessaire de disposer de modèles et de méthodes distinctes. L?objectif de ce travail a été de développer des modèles prédictifs de la dynamique des grandes populations, ainsi que des logiciels permettant de calculer les paramètres de ces modèles et de tester des scénarios de gestion. Le cas du Bouquetin des Alpes (Capra ibex ibex) - en forte expansion en Suisse depuis sa réintroduction au début du XXème siècle - servit d?exemple. Cette tâche fut accomplie en trois étapes : En premier lieu, un modèle de dynamique locale, spécifique au Bouquetin, fut développé : le modèle sous-jacent - structuré en classes d?âge et de sexe - est basé sur une matrice de Leslie à laquelle ont été ajoutées la densité-dépendance, la stochasticité environnementale et la chasse de régulation. Ce modèle fut implémenté dans un logiciel d?aide à la gestion - nommé SIM-Ibex - permettant la maintenance de données de recensements, l?estimation automatisée des paramètres, ainsi que l?ajustement et la simulation de stratégies de régulation. Mais la dynamique d?une population est influencée non seulement par des facteurs démographiques, mais aussi par la dispersion et la colonisation de nouveaux espaces. Il est donc nécessaire de pouvoir modéliser tant la qualité de l?habitat que les obstacles à la dispersion. Une collection de logiciels - nommée Biomapper - fut donc développée. Son module central est basé sur l?Analyse Factorielle de la Niche Ecologique (ENFA) dont le principe est de calculer des facteurs de marginalité et de spécialisation de la niche écologique à partir de prédicteurs environnementaux et de données d?observation de l?espèce. Tous les modules de Biomapper sont liés aux Systèmes d?Information Géographiques (SIG) ; ils couvrent toutes les opérations d?importation des données, préparation des prédicteurs, ENFA et calcul de la carte de qualité d?habitat, validation et traitement des résultats ; un module permet également de cartographier les barrières et les corridors de dispersion. Le domaine d?application de l?ENFA fut exploré par le biais d?une distribution d?espèce virtuelle. La comparaison à une méthode couramment utilisée pour construire des cartes de qualité d?habitat, le Modèle Linéaire Généralisé (GLM), montra qu?elle était particulièrement adaptée pour les espèces cryptiques ou en cours d?expansion. Les informations sur la démographie et le paysage furent finalement fusionnées en un modèle global. Une approche basée sur un automate cellulaire fut choisie, tant pour satisfaire aux contraintes du réalisme de la modélisation du paysage qu?à celles imposées par les grandes populations : la zone d?étude est modélisée par un pavage de cellules hexagonales, chacune caractérisée par des propriétés - une capacité de soutien et six taux d?imperméabilité quantifiant les échanges entre cellules adjacentes - et une variable, la densité de la population. Cette dernière varie en fonction de la reproduction et de la survie locale, ainsi que de la dispersion, sous l?influence de la densité-dépendance et de la stochasticité. Un logiciel - nommé HexaSpace - fut développé pour accomplir deux fonctions : 1° Calibrer l?automate sur la base de modèles de dynamique (par ex. calculés par SIM-Ibex) et d?une carte de qualité d?habitat (par ex. calculée par Biomapper). 2° Faire tourner des simulations. Il permet d?étudier l?expansion d?une espèce envahisseuse dans un paysage complexe composé de zones de qualité diverses et comportant des obstacles à la dispersion. Ce modèle fut appliqué à l?histoire de la réintroduction du Bouquetin dans les Alpes bernoises (Suisse). SIM-Ibex est actuellement utilisé par les gestionnaires de la faune et par les inspecteurs du gouvernement pour préparer et contrôler les plans de tir. Biomapper a été appliqué à plusieurs espèces (tant végétales qu?animales) à travers le Monde. De même, même si HexaSpace fut initialement conçu pour des espèces animales terrestres, il pourrait aisément être étndu à la propagation de plantes ou à la dispersion d?animaux volants. Ces logiciels étant conçus pour, à partir de données brutes, construire un modèle réaliste complexe, et du fait qu?ils sont dotés d?une interface d?utilisation intuitive, ils sont susceptibles de nombreuses applications en biologie de la conservation. En outre, ces approches peuvent également s?appliquer à des questions théoriques dans les domaines de l?écologie des populations et du paysage.<br/><br/>Conservation biology is commonly associated to small and endangered population protection. Nevertheless, large or potentially large populations may also need human management to prevent negative effects of overpopulation. As there are both qualitative and quantitative differences between small population protection and large population controlling, distinct methods and models are needed. The aim of this work was to develop theoretical models to predict large population dynamics, as well as computer tools to assess the parameters of these models and to test management scenarios. The alpine Ibex (Capra ibex ibex) - which experienced a spectacular increase since its reintroduction in Switzerland at the beginning of the 20th century - was used as paradigm species. This task was achieved in three steps: A local population dynamics model was first developed specifically for Ibex: the underlying age- and sex-structured model is based on a Leslie matrix approach with addition of density-dependence, environmental stochasticity and culling. This model was implemented into a management-support software - named SIM-Ibex - allowing census data maintenance, parameter automated assessment and culling strategies tuning and simulating. However population dynamics is driven not only by demographic factors, but also by dispersal and colonisation of new areas. Habitat suitability and obstacles modelling had therefore to be addressed. Thus, a software package - named Biomapper - was developed. Its central module is based on the Ecological Niche Factor Analysis (ENFA) whose principle is to compute niche marginality and specialisation factors from a set of environmental predictors and species presence data. All Biomapper modules are linked to Geographic Information Systems (GIS); they cover all operations of data importation, predictor preparation, ENFA and habitat suitability map computation, results validation and further processing; a module also allows mapping of dispersal barriers and corridors. ENFA application domain was then explored by means of a simulated species distribution. It was compared to a common habitat suitability assessing method, the Generalised Linear Model (GLM), and was proven better suited for spreading or cryptic species. Demography and landscape informations were finally merged into a global model. To cope with landscape realism and technical constraints of large population modelling, a cellular automaton approach was chosen: the study area is modelled by a lattice of hexagonal cells, each one characterised by a few fixed properties - a carrying capacity and six impermeability rates quantifying exchanges between adjacent cells - and one variable, population density. The later varies according to local reproduction/survival and dispersal dynamics, modified by density-dependence and stochasticity. A software - named HexaSpace - was developed, which achieves two functions: 1° Calibrating the automaton on the base of local population dynamics models (e.g., computed by SIM-Ibex) and a habitat suitability map (e.g. computed by Biomapper). 2° Running simulations. It allows studying the spreading of an invading species across a complex landscape made of variously suitable areas and dispersal barriers. This model was applied to the history of Ibex reintroduction in Bernese Alps (Switzerland). SIM-Ibex is now used by governmental wildlife managers to prepare and verify culling plans. Biomapper has been applied to several species (both plants and animals) all around the World. In the same way, whilst HexaSpace was originally designed for terrestrial animal species, it could be easily extended to model plant propagation or flying animals dispersal. As these softwares were designed to proceed from low-level data to build a complex realistic model and as they benefit from an intuitive user-interface, they may have many conservation applications. Moreover, theoretical questions in the fields of population and landscape ecology might also be addressed by these approaches.

Phylogéographie des principales vipères européennes (Vipera ammodytes, V. aspis et V. berus), structuration génétique et multipaternité chez Vipera berus

Résumé: Les vipères du genre Vipera sont des serpents venimeux distribués dans la totalité du Paléarctique. Malgré cette répartition considérable, elles sont extrêmement menacées, leur déclin étant principalement dû à la destruction et à la fragmentation de leur habitat ainsi qu'à la persécution humaine. Afin d'apporter de nouveaux éléments dans le contexte de la protection de ce groupe de reptiles, nous avons utilisé durant ce travail de thèse différents marqueurs moléculaires pour étudier la structuration génétique à petite et à large échelle chez trois espèces appartenant au genre Vipera. La première étude, une phylogéographie moléculaire de la vipère ammodytes (Vipera ammodytes), a montré dans l'ensemble de l'aire de répartition une forte structuration génétique provenant d'isolements antérieures au Pléistocène. La présence d'un nombre important de clades dans le centre des Balkans suggère que cette région a fourni de nombreux refuges isolés durant les glaciations. Ces dernières ont également eu un impact considérable sur la diversité génétique au sein de la majorité des clades, suite à d'importants goulots d'étranglement durant le Pléistocène. L'étude de la phylogéographie de la vipère aspic (Vipera aspis) a montré une différenciation génétique entre les populations présentes de chaque côté des Alpes, mais également une forte structuration interne avec la mise en évidence d'un refuge en France. Cette étude est la première à établir clairement l'utilisation d'un refuge français pour un vertébré terrestre. La troisième partie de cette thèse a étudié la phylogéographie de la vipère péliade (Vipera berus), espèce cible de ce travail. En plus de la mise en évidence d'un groupe génétique inattendu (localisé dans le nord de l'Italie, le sud de l'Autriche, le nord de la Slovénie et l'extrême sud-est de la Suisse), la variabilité génétique au sein du groupe nordique (comprenant les animaux de l'entier de l'aire de répartition de l'espèce à l'exception des individus du groupe italien et les animaux provenant des Balkans) est suffisamment importante pour conclure à l'utilisation de refuges glaciaires nordiques durant les dernières glaciations, en complément des refuges habituellement décrits pour la majorité des espèces animales (soit les péninsules ibérique, italienne et balakanique). Ces résultats nous ont conduit à effectuer une étude morphologique (quatrième partie) comparant les vipères péliades du "clade italien" et du "clade nordique" décrits ci-dessus. Seules de petites différences morphologiques ont pu être mises en évidence, malgré une séparation de ces groupes estimée à plus d'un million d'années. Une étude à plus petite échelle, centrée sur le Massif jurassien et certaines populations alpines et françaises, a été entreprise afin d'estimer leur diversité génétique et d'évaluer la structuration génétique entre les populations à l'aide de marqueurs microsatellites (cinquième partie). Une importante structuration a été observée entre les populations distantes de plus de 3 kilomètres, la structuration entre les populations plus proches étant plus limitée. De plus, une diversité génétique plus faible dans les populations jurassiennes et alpines comparativement aux populations du massif central et de la côte atlantique a été constatée, probablement due à une perte de diversité génétique lors de la recolonisation post-glaciaire. La sixième étude s'est intéressée au succès reproducteur des mâles de vipères péliades en conditions naturelles. Une corrélation entre la taille des mâles et leur succès reproducteur a été relevée, les individus de plus grande taille ayant un succès reproducteur plus élevé. Le taux de multipaternité a aussi été investigué, démontrant que la proportion de pontes issues de plusieurs pères est élevée (69%) malgré la faible densité de vipères observée sur le site étudié. Finalement, aucun lien entre le nombre de pères au sein d'une ponte et la mortalité des jeunes à la naissance n'a pu être mis en évidence, contrastant avec des travaux précédents. En conclusion, l'observation de la structuration très marquée chez les vipères péliades devrait permettre d'affiner les méthodes de protection de l'espèce dans le massif jurassien. A plus large échelle, l'importante structuration génétique observée chez les vipères ammodytes, aspic et péliade résultant de l'utilisation de nombreux refuges glaciaires, complémentaires aux refuges habituellement utilisés par les espèces animales, démontre l'intérêt de l'analyse phylogéographique des reptiles pour la compréhension des phénomènes de colonisation et d' extinction des populations durant la fin du Tertiaire et le Quaternaire. La mise en évidence chez les différentes espèces de vipères étudiées de nombreux groupes génétiques distincts (ESUs) devrait conduire à des modifications de la taxonomie ainsi qu'au statut de protection de ces espèces. Abstract: The vipers of the genus Vipera are venomous snakes widespread throughout the Palaearctic regions. Despite a large distribution area, several species are extremely threatened, especially due to the destruction and fragmentation of their habitats, as well as by human persecution. In order to increase the knowledge on these species and to improve their protection, several molecular markers have been used to investigate the genetic structure on small and large scales, within three species of the genus Vipera. The first study, a molecular phylogeography of the nose-horned viper (Vipera ammodytes), showed a considerable structuring throughout the distribution area, due to isolation into refugia before the Pleistocene. A high number of clades in the centre of the Balkans suggests that this region harboured numerous isolated glacial refugia during the last glaciation. Moreover, low genetic diversity within several clades implies that most populations of nose-horned vipers have suffered bottlenecks during the Pleistocene. The study of the phylogeography of the asp viper (Vipera aspis) showed genetic differentiation between populations on each side of the Alps, as well as considerable internal genetic structure, suggesting the use of a glacial refugium in France. This study is the first to establish firmly the occurrence of a French refugia for a terrestrial vertebrate. The third part of this work involved a phylogeographic study of the adder (Vipera berus), the target species of this thesis. Three clades were revealed: a Balkan clade (corresponding to the subspecies V. b. sachalinensis), an unexpected Italian clade (limited to northern Italy, southern Austria, northern Slovenia and southeasternmost corner of Switzerland) and a Northern clade clade (including adders of the whole distribution area excepted animals from the Balkan and the Italian clades). The genetic variability within the Northern clade is sufficiently high to conclude that a northern glacial refugia during the last glaciation, in addition to those refugia already described for the main species (Iberian, Italian and Balkan peninsula). These results motivated a morphological study (part four) comparing the adders from the Italian and the Northern clades describe above. Only small morphological differences have been found, despite the split between these two clades have taken place more than 1 million years ago. A study on a local scale, focused on the Jura Mountains, on a few populations in the Alps and France was, performed to estimate the genetic diversity and the genetic structure between populations using microsatellite markers (part five). Considerable structure was observed between populations separated by more than 3 kilometres, whereas the structure between closer populations is less marked. Moreover, lower genetic diversity in the populations from Jura Mountains and Alps was noticed compared to populations from Massif Central of Atlantic coast. Such loss of genetic variation probably followed post-glacial recolonisation. The sixth study focused on the reproductive success of male adders in the wild. A positive correlation between body length and reproductive success was observed. Multiple paternity was also observed in most of clutches (69%) despite the low density of adders in the study area. Finally, no relationship was found between the number of fathers in a clutch and the survival of offspring at birth, contradicting previous studies. To conclude, the observation of a significant genetic structure in Vipera berus will enable recommendations to be made to improve protection of this species in the Jura Mountain. On a larger scale, the considerable genetic structure found within Vipera ammdoytes, V. aspis and V. berus, resulting from isolation in additional glacial refugia to those already described for other species, demonstrates the relevance of phylogeographic studies of reptiles to better understand the colonisation and disappearance during the last Tertiary and the Quaternary. The observation of several groups of evolutionary significant units (ESUs) within the three studied species might lead to a revision of the taxonomy, as well as their conservation status.

Three-dimensional hydrodynamic numerical modelling of fold nappe formation, basement-cover deformation and slab detachment with applications to the helvetic nappe system (W Switzerland)

Many three-dimensional (3-D) structures in rock, which formed during the deformation of the Earth's crust and lithosphere, are controlled by a difference in mechanical strength between rock units and are often the result of a geometrical instability. Such structures are, for example, folds, pinch-and-swell structures (due to necking) or cuspate-lobate structures (mullions). These struc-tures occur from the centimeter to the kilometer scale and the related deformation processes con-trol the formation of, for example, fold-and-thrust belts and extensional sedimentary basins or the deformation of the basement-cover interface. The 2-D deformation processes causing these structures are relatively well studied, however, several processes during large-strain 3-D defor-mation are still incompletely understood. One of these 3-D processes is the lateral propagation of these structures, such as fold and cusp propagation in a direction orthogonal to the shortening direction or neck propagation in direction orthogonal to the extension direction. Especially, we are interested in fold nappes which are recumbent folds with amplitudes usually exceeding 10 km and they have been presumably formed by ductile shearing. They often exhibit a constant sense of shearing and a non-linear increase of shear strain towards their overturned limb. The fold axes of the Morcles fold nappe in western Switzerland plunges to the ENE whereas the fold axes in the more eastern Doldenhorn nappe plunges to the WSW. These opposite plunge direc-tions characterize the Rawil depression (Wildstrubel depression). The Morcles nappe is mainly the result of layer parallel contraction and shearing. During the compression the massive lime-stones were more competent than the surrounding marls and shales, which led to the buckling characteristics of the Morcles nappe, especially in the north-dipping normal limb. The Dolden-horn nappe exhibits only a minor overturned fold limb. There are still no 3-D numerical studies which investigate the fundamental dynamics of the formation of the large-scale 3-D structure including the Morcles and Doldenhorn nappes and the related Rawil depression. We study the 3-D evolution of geometrical instabilities and fold nappe formation with numerical simulations based on the finite element method (FEM). Simulating geometrical instabilities caused by sharp variations of mechanical strength between rock units requires a numerical algorithm that can accurately resolve material interfaces for large differences in material properties (e.g. between limestone and shale) and for large deformations. Therefore, our FE algorithm combines a nu-merical contour-line technique and a deformable Lagrangian mesh with re-meshing. With this combined method it is possible to accurately follow the initial material contours with the FE mesh and to accurately resolve the geometrical instabilities. The algorithm can simulate 3-D de-formation for a visco-elastic rheology. The viscous rheology is described by a power-law flow law. The code is used to study the 3-D fold nappe formation, the lateral propagation of folding and also the lateral propagation of cusps due to initial half graben geometry. Thereby, the small initial geometrical perturbations for folding and necking are exactly followed by the FE mesh, whereas the initial large perturbation describing a half graben is defined by a contour line inter-secting the finite elements. Further, the 3-D algorithm is applied to 3-D viscous nacking during slab detachment. The results from various simulations are compared with 2-D resulats and a 1-D analytical solution. -- On retrouve beaucoup de structures en 3 dimensions (3-D) dans les roches qui ont pour origines une déformation de la lithosphère terrestre. Ces structures sont par exemple des plis, des boudins (pinch-and-swell) ou des mullions (cuspate-lobate) et sont présentés de l'échelle centimétrique à kilométrique. Mécaniquement, ces structures peuvent être expliquées par une différence de résistance entre les différentes unités de roches et sont généralement le fruit d'une instabilité géométrique. Ces différences mécaniques entre les unités contrôlent non seulement les types de structures rencontrées, mais également le type de déformation (thick skin, thin skin) et le style tectonique (bassin d'avant pays, chaîne d'avant pays). Les processus de la déformation en deux dimensions (2-D) formant ces structures sont relativement bien compris. Cependant, lorsque l'on ajoute la troisiéme dimension, plusieurs processus ne sont pas complètement compris lors de la déformation à large échelle. L'un de ces processus est la propagation latérale des structures, par exemple la propagation de plis ou de mullions dans la direction perpendiculaire à l'axe de com-pression, ou la propagation des zones d'amincissement des boudins perpendiculairement à la direction d'extension. Nous sommes particulièrement intéressés les nappes de plis qui sont des nappes de charriage en forme de plis couché d'une amplitude plurikilométrique et étant formées par cisaillement ductile. La plupart du temps, elles exposent un sens de cisaillement constant et une augmentation non linéaire de la déformation vers la base du flanc inverse. Un exemple connu de nappes de plis est le domaine Helvétique dans les Alpes de l'ouest. Une de ces nap-pes est la Nappe de Morcles dont l'axe de pli plonge E-NE tandis que de l'autre côté de la dépression du Rawil (ou dépression du Wildstrubel), la nappe du Doldenhorn (équivalent de la nappe de Morcles) possède un axe de pli plongeant O-SO. La forme particulière de ces nappes est due à l'alternance de couches calcaires mécaniquement résistantes et de couches mécanique-ment faibles constituées de schistes et de marnes. Ces différences mécaniques dans les couches permettent d'expliquer les plissements internes à la nappe, particulièrement dans le flanc inver-se de la nappe de Morcles. Il faut également noter que le développement du flanc inverse des nappes n'est pas le même des deux côtés de la dépression de Rawil. Ainsi la nappe de Morcles possède un important flanc inverse alors que la nappe du Doldenhorn en est presque dépour-vue. A l'heure actuelle, aucune étude numérique en 3-D n'a été menée afin de comprendre la dynamique fondamentale de la formation des nappes de Morcles et du Doldenhorn ainsi que la formation de la dépression de Rawil. Ce travail propose la première analyse de l'évolution 3-D des instabilités géométriques et de la formation des nappes de plis en utilisant des simulations numériques. Notre modèle est basé sur la méthode des éléments finis (FEM) qui permet de ré-soudre avec précision les interfaces entre deux matériaux ayant des propriétés mécaniques très différentes (par exemple entre les couches calcaires et les couches marneuses). De plus nous utilisons un maillage lagrangien déformable avec une fonction de re-meshing (production d'un nouveau maillage). Grâce à cette méthode combinée il nous est possible de suivre avec précisi-on les interfaces matérielles et de résoudre avec précision les instabilités géométriques lors de la déformation de matériaux visco-élastiques décrit par une rhéologie non linéaire (n>1). Nous uti-lisons cet algorithme afin de comprendre la formation des nappes de plis, la propagation latérale du plissement ainsi que la propagation latérale des structures de type mullions causé par une va-riation latérale de la géométrie (p.ex graben). De plus l'algorithme est utilisé pour comprendre la dynamique 3-D de l'amincissement visqueux et de la rupture de la plaque descendante en zone de subduction. Les résultats obtenus sont comparés à des modèles 2-D et à la solution analytique 1-D. -- Viele drei dimensionale (3-D) Strukturen, die in Gesteinen vorkommen und durch die Verfor-mung der Erdkruste und Litosphäre entstanden sind werden von den unterschiedlichen mechani-schen Eigenschaften der Gesteinseinheiten kontrolliert und sind häufig das Resulat von geome-trischen Istabilitäten. Zu diesen strukturen zählen zum Beispiel Falten, Pich-and-swell Struktu-ren oder sogenannte Cusbate-Lobate Strukturen (auch Mullions). Diese Strukturen kommen in verschiedenen Grössenordungen vor und können Masse von einigen Zentimeter bis zu einigen Kilometer aufweisen. Die mit der Entstehung dieser Strukturen verbundenen Prozesse kontrol-lieren die Entstehung von Gerbirgen und Sediment-Becken sowie die Verformung des Kontaktes zwischen Grundgebirge und Stedimenten. Die zwei dimensionalen (2-D) Verformungs-Prozesse die zu den genannten Strukturen führen sind bereits sehr gut untersucht. Einige Prozesse wäh-rend starker 3-D Verformung sind hingegen noch unvollständig verstanden. Einer dieser 3-D Prozesse ist die seitliche Fortpflanzung der beschriebenen Strukturen, so wie die seitliche Fort-pflanzung von Falten und Cusbate-Lobate Strukturen senkrecht zur Verkürzungsrichtung und die seitliche Fortpflanzung von Pinch-and-Swell Strukturen othogonal zur Streckungsrichtung. Insbesondere interessieren wir uns für Faltendecken, liegende Falten mit Amplituden von mehr als 10 km. Faltendecken entstehen vermutlich durch duktile Verscherung. Sie zeigen oft einen konstanten Scherungssinn und eine nicht-lineare zunahme der Scherverformung am überkipp-ten Schenkel. Die Faltenachsen der Morcles Decke in der Westschweiz fallen Richtung ONO während die Faltenachsen der östicher gelegenen Doldenhorn Decke gegen WSW einfallen. Diese entgegengesetzten Einfallrichtungen charakterisieren die Rawil Depression (Wildstrubel Depression). Die Morcles Decke ist überwiegend das Resultat von Verkürzung und Scherung parallel zu den Sedimentlagen. Während der Verkürzung verhielt sich der massive Kalkstein kompetenter als der Umliegende Mergel und Schiefer, was zur Verfaltetung Morcles Decke führ-te, vorallem in gegen Norden eifallenden überkippten Schenkel. Die Doldenhorn Decke weist dagegen einen viel kleineren überkippten Schenkel und eine stärkere Lokalisierung der Verfor-mung auf. Bis heute gibt es keine 3-D numerischen Studien, die die fundamentale Dynamik der Entstehung von grossen stark verformten 3-D Strukturen wie den Morcles und Doldenhorn Decken sowie der damit verbudenen Rawil Depression untersuchen. Wir betrachten die 3-D Ent-wicklung von geometrischen Instabilitäten sowie die Entstehung fon Faltendecken mit Hilfe von numerischen Simulationen basiert auf der Finite Elemente Methode (FEM). Die Simulation von geometrischen Instabilitäten, die aufgrund von Änderungen der Materialeigenschaften zwischen verschiedenen Gesteinseinheiten entstehen, erfortert einen numerischen Algorithmus, der in der Lage ist die Materialgrenzen mit starkem Kontrast der Materialeigenschaften (zum Beispiel zwi-schen Kalksteineinheiten und Mergel) für starke Verfomung genau aufzulösen. Um dem gerecht zu werden kombiniert unser FE Algorithmus eine numerische Contour-Linien-Technik und ein deformierbares Lagranges Netz mit Re-meshing. Mit dieser kombinierten Methode ist es mög-lich den anfänglichen Materialgrenzen mit dem FE Netz genau zu folgen und die geometrischen Instabilitäten genügend aufzulösen. Der Algorithmus ist in der Lage visko-elastische 3-D Ver-formung zu rechnen, wobei die viskose Rheologie mit Hilfe eines power-law Fliessgesetzes beschrieben wird. Mit dem numerischen Algorithmus untersuchen wir die Entstehung von 3-D Faltendecken, die seitliche Fortpflanzung der Faltung sowie der Cusbate-Lobate Strukturen die sich durch die Verkürzung eines mit Sediment gefüllten Halbgraben bilden. Dabei werden die anfänglichen geometrischen Instabilitäten der Faltung exakt mit dem FE Netz aufgelöst wäh-rend die Materialgranzen des Halbgrabens die Finiten Elemente durchschneidet. Desweiteren wird der 3-D Algorithmus auf die Einschnürung während der 3-D viskosen Plattenablösung und Subduktion angewandt. Die 3-D Resultate werden mit 2-D Ergebnissen und einer 1-D analyti-schen Lösung verglichen.

Accessory phases as a tool to investigate magmatic underplating in polymetamorphic metasediments : an example from a fossil lower crust : Ivrea Verbano zne, NW Italy

The Ivrea and the Strona-Ceneri zones, NW italy and S Switzerland, offer the possibility to study the continental crust of the Southern Alps. Because of its high metamorphic degree and the abundant Permo- Carboniferous mafic intrusions, the Ivrea Zone is classically interpreted an exposed section trough the Permian lower crust. The present work is focused here on metasedimentary slices (septa) intercalated within Permian gabbro (mafic complex). In particular I studied the evolution of accessory phases such as rutile and zircon and the chemistry of the metasediments. The septa build an irregular and discontinuous band that cut obliquely the mafic complex from its deepest part (N) to its roof (S). The chemistry of the metasediments evolves along the band and the chemical evolution can be compared with that observed in the country-rock surrounding the mafic intrusion to the NE and overprinted by a main regional metamorphic event. This suggests that the degree of chemical depletion of the septa was mainly established during the same regional metamorphic event. Moreover it suggests that incorporation of the septa within the gabbro did not modify their original stratigraphie distribution within the crust. It implies that the mafic complex has been emplaced following a dynamic substantially different from the classic model of « gabbro glacier » (Quick et al., 1992; Quick et al., 1994). It is more likely that it has been emplaced by repeated injections of sills at different depths during a protracted period of time. Zircon trace elements and U-Pb ages suggest that regional metamorphism occurred 330-320Ma, the first sills in the deepest part of the Mafic Complex are injected at ~300Ma, the mafic magmas reached higher levels in the crust at 285Ma and the magmatic activity continued locally until 275Ma. The ages of detrital cores in zircons fix the maximal sedimentation age at ~370Ma, this age corresponds therefore with the maximal age of the incorporation of the Ivrea zone within the lower crust. I propose that the Ivrea zone has been accreted to the lower crust during the Hercynian orogeny sensu lato. The analysis of detrital ages suggests that the source terrains for the Ivrea zone and those for the Strona-Ceneri zone have a completely different Palaeozoic history. The systematic analysis of rutile in partially molten metasediments of the Ivrea zone reveals the occurrence of two generations. The two generations are characterized by a different chemistry and textural distribution. A first generation is formed during pro-grade metamorphism in the restitic counterpart. The second generation is formed in the melts during cooling at the same time that part of the first generation re-equilibrate. Re-equilibration of the first generation seems to be spatially controlled by the presence of fluids. Locally the second generation forms overgrowths on the first generation. Considered the different diffusivity of U and Pb in rutile, U heterogeneities have important implication for U-Pb dating of rutile. ID-TIMS and LA-ICPMS dating coupled with a careful textural investigation (SEM) suggest that rutile grains are characterized by multiple path along which Pb diffusion can occur: volume diffusion is an important process, but intragrain and subgrain boundaries provide additional high diffusivity pathways for Pb escape and reduce drastically the effective diffusion length. -- La zone d'Ivrea et la zone de Strona-Ceneri, en Italie nord-occidentale et Suisse méridionale, offrent la possibilité d'étudier la croûte continentale des Alpes du Sud. En raison du haut degré métamorphique et l'abondance d'intrusions mafiques d'âge Permo-Carbonifère [complexe mafique), la zone d'Ivrea est interprétée classiquement comme de la croûte inférieure permienne. Ce travail ce concentre sur des bandes metasédimentaires (septa) incorporées dans les magmas mafiques lors de l'intrusion. Les septa forment une bande irrégulière qui coupe obliquement le complexe mafique du bas (N) vers le haut (S). La chimie des septa évolue du bas vers le haut et l'évolution chimique se rapproche de l'évolution observé dans la roche encaissante l'intrusion affecté par un événement métamorphique régionale. Cette relation suggère que le degré d'appauvrissement chimique des septa a été établit principalement lors de l'événement métamorphique régional. De plus l'incorporation dans les gabbros n'a pas perturbée la distribution stratigraphique originelle des septa. Ces deux observations impliquent que le métamorphisme dans la roche encaissante précède la mise en place du gabbro et que cette dernière ne se fait pas selon le modèle classique (« gabbro glacier » de Quick et al., 1992, 1994), mais se fait plutôt par injections répétées de sills a différentes profondeurs. Les âges U-Pb et les éléments traces des zircons suggèrent que le métamorphisme régionale a eu lieu 330-320Ma, alors que les premiers sills dans la partie profonde du Mafic Complex s'injectent à ~300Ma, le magmatisme mafique atteigne des niveaux supérieurs à 285Ma et continue localement jusqu'à 270Ma. Les âges des coeurs détritiques des zircons permettent de fixer l'âge maximale de sédimentation à ~370Ma ce qui correspond donc à l'âge maximale de l'incorporation de la zone d'Ivrea dans la croûte inférieur. L'analyse systématique des rutiles, nous a permit de montrer l'existence de plusieurs générations qui ont une répartition texturale et une chimie différente. Une génération se forme lors de l'événement UHT dans les restites, une autre génération se forme dans les liquides lors du refroidissement, au même temps qu'une partie de la première génération se rééquilibre au niveau du Zr. Localement la deuxième génération peut former des surcroissances autour de la première génération. Dans ces cas, des fortes différences en uranium entre les deux générations ont des importantes implications pour la datation U-Pb sur rutile. Classiquement les ratios Pb/U dans le rutile sont interprétés comme indiquant l'âges du refroidissement du minéral sous une température à la quelle la diffusion du Pb dans le minéral n'est plus détectable et la diffusion à plus hautes températures est assumée se faire par «volume diffusion» dans le grain (Mezger et al., 1989). Par des datations ID-TIMS (sur grain entier) et LA-ICPMS (in-situ) et une analyse texturale (MEB) approfondie nous montrons que cette supposition est trop simpliste et que le rutile est repartie en sous-domaines. Chacun de ces domaines a ça propre longueur ou chemin de diffusion spécifique. Nous proposons donc une nouvelle approche plus cohérente pour l'interprétation des âges U-Pb sur rutile.