Étude des sédiments quaternaires, de la molasse et sa tectonique, dans le Grand Lac (Léman) à partir de données sismiques 2D et 3D

RESUME L'Institut de Géophysique de l'Université de Lausanne a développé au cours de ces dernières années un système d'acquisition de sismique réflexion multitrace à haute résolution 2D et 3D. L'objectif de cette thèse était de poursuivre ce développement tout améliorant les connaissances de la géologie sous le lac Léman, en étudiant en particulier la configuration des grands accidents sous-lacustres dans la Molasse (Tertiaire) qui forme l'essentiel du substratum des formations quaternaires. En configuration 2D, notre système permet d'acquérir des profils sismiques avec une distance inter-CDP de 1,25 m. La couverture varie entre 6 et 18 selon le nombre de traces et la distance inter-tir. Le canon à air (15/15 eu. in.), offre une résolution verticale de 1,25 ni et une pénétration maximale de 300 m sous le fond de l'eau. Nous avons acquis au total plus de 400 km de sections 2D dans le Grand Lac et le Haut Lac entre octobre 2000 et juillet 2004. Une campagne de sismique 3D a fourni des données au large d'Evian sur une surface de 442,5 m sur 1450 m, soit 0,64 km2. La navigation ainsi que le positionnement des hydrophones et de la source ont été réalisés avec des GPS différentiels. Nous avons utilisé un traitement sismique conventionnel, sans appliquer d'AGC et en utilisant une migration post-stack. L'interprétation du substratum antéquaternaire est basée sur l'identification des sismofaciès, sur leurs relations avec les unités géologiques adjacentes au lac, ainsi que sur quelques données de forages. Nous obtenons ainsi une carte des unités géologiques dans le Grand Lac. Nous précisons la position du chevauchement subalpin entre la ville de Lausanne, sur la rive nord, et le bassin de Sciez, sur la rive sud. Dans la Molasse de Plateau, nous avons identifié les décrochements de Pontarlier et de St. Cergue ainsi que plusieurs failles non reconnues jusqu'ici. Nous avons cartographié les accidents qui affectent la Molasse subalpine ainsi que le plan de chevauchement du flysch sur la Molasse près de la rive sud du lac. Une nouvelle carte tectonique de la région lémanique a ainsi pu être dressée. L'analyse du substratum ne montre pas de failles suggérant une origine tectonique de la cuvette lémanique. Par contre, nous suggérons que la forme du creusement glaciaire, donc de la forme du lac Léman, a été influencée par la présence de failles dans le substratum antéquaternaire. L'analyse des sédiments quaternaires nous a permis de tracer des cartes des différentes interfaces ou unités qui les composent. La carte du toit du substratum antéquaternaire montre la présence de chenaux d'origine glaciaire dont la profondeur maximale atteint la cote -200 ni. Leur pente est dirigée vers le nord-est, à l'inverse du sens d'écoulement actuel des eaux. Nous expliquons cette observation par l'existence de circulations sous-glaciaires d'eau artésienne. Les sédiments glaciaires dont l'épaisseur maximale atteint 150 ni au centre du lac ont enregistré les différentes récurrences glaciaires. Dans la zone d'Evian, nous mettons en évidence la présence de lentilles de sédiments glaciolacustres perchées sur le flanc de la cuvette lémanique. Nous avons corrélé ces unités avec des données de forage et concluons qu'il s'agit du complexe inférieur de la pile sédimentaire d'Evian. Celui-ci, âgé de plus de 30 000 ans, serait un dépôt de Kame associé à un lac périglaciaire. La sismique réflexion 3D permet de préciser l'orientation de l'alimentation en matériel détritique de l'unité. La finesse des images obtenues nous permet également d'établir quels types d'érosion ont affecté certaines unités. Les sédiments lacustres, dont l'épaisseur maximale imagée atteint plus de 225 m et sans doute 400 ni sous le delta du Rhône, indiquent plusieurs mécanismes de dépôts. A la base, une mégaturbidite, épaisse d'une trentaine de mètres en moyenne, s'étend entre l'embouchure de la Dranse et le delta du Rhône. Au-dessus, la décantation des particules en suspension d'origine biologique et détritique fournit l'essentiel des sédiments. Dans la partie orientale du lac, les apports détritiques du Rhône forment un delta qui prograde vers l'ouest en s'imbriquant avec les sédiments déposés par décantation. La structure superficielle du delta a brutalement évolué, probablement à la suite de l'évènement catastrophique du Tauredunum (563 A.D.). Sa trace probable se marque par la présence d'une surface érosive que nous avons cartographiée. Le delta a ensuite changé de géométrie, avec notamment un déplacement des chenaux sous-lacustres. Sur l'ensemble de nos sections sismiques, nous n'observons aucune faille dans les sédiments quaternaires qui attesterait d'une tectonique postglaciaire du substratum. ABSTRACT During the last few years the institute of Geophysics of the University of Lausanne cleveloped a 2D and 3D high-resolution multichannel seismic reflection acquisition system. The objective of the present work was to carry on this development white improving our knowledge of the geology under Lake Geneva, in particular by studying the configuration of the large accidents affecting the Tertiary Molasse that makes up the basement of most Quaternary deposits. In its 2D configuration, our system makes it possible to acquire seismic profiles with a CDP interval of 1.25 m. The fold varies from 6 to 18 depending on the number of traces and the shooting interval. Our air gun (15/15 cu. in.) provides a vertical resolution of 1.25 m and a maximum penetration depth of approximately 300 m under water bottom. We acquired more than 400 km of 2D sections in the Grand Lac and the Haut Lac between October 2000 and July 2004. A 3D seismic survey off the city of Evian provided data on a surface of 442.5 m x 1450 m (0.64 km2). Ship's navigation as well as hydrophone- and source positioning were carried out with differential GPS. The seismic data were processed following a conventional sequence without .applying AGC and using post-stack migration. The interpretation of the pre-Quaternary substratum is based on sismofacies, on their relationships with terrestrial geological units and on some borehole data. We thus obtained a map of the geological units in the Grand Lac. We defined the location of the subalpine thrust from Lausanne, on the north shore, to the Sciez Basin, on the south shore. Within the Molasse de Plateau, we identified the already know Pontarlier and St Cergue transforms Fault as well as faults. We mapped faults that affect subalpine Molasse as well as the thrust fault plane between alpine flysch and Molasse near the lake's south shore. A new tectonic map of the Lake Geneva region could thus be drawn up. The substratum does not show faults indicating a tectonic origin for the Lake Geneva Basin. However, we suggest that the orientation of glacial erosion, and thus the shape of Lake Geneva, vas influenced by the presence of faults in the pre-Quaternary basement. The analysis of Quaternary sediments enabled us to draw up maps of various discontinuities or internal units. The top pre-Quaternary basement map shows channels of glacial origin, the deepest of them reaching an altitude of 200 m a.s.l. The channel's slopes are directed to the North-East, in opposite direction of the present water flow. We explain this observation by the presence of artesian subglacial water circulation. Glacial sediments, the maximum thickness of which reaches 150 m in the central part of the lake, record several glacial recurrences. In the Evian area, we found lenses of glacio-lacustrine sediments set high up on the flank of the Lake Geneva Bassin. We correlated these units with on-land borehole data and concluded that they represent the lower complex of the Evian sedimentary pile. The lower complex is aider than 30 000 years, and it could be a Kame deposit associated with a periglacial lake. Our 3D seismic reflexion survey enables us to specify the supply direction of detrital material in this unit. With detailed seismic images we established how some units were affected by different erosion types. The lacustrine sediments we imaged in Lake Geneva are thicker than 225 m and 400 m or more Linder the Rhone Delta. They indicate several depositional mechanisms. Their base is a major turbidite, thirty meters thick on average, that spreads between the Dranse mouth and the Rhone delta. Above this unit, settling of suspended biological and detrital particles provides most of the sediments. In the eastern part of the lake, detrital contribution from the Rhone builds a delta that progrades to the west and imbricates with the settling sediments. The shallow structure of the Rhone delta abruptly evolved, probably after the catastrophic Tauredunum event (563 A.D.). It probably coincides with an erosive surface that we mapped. As a result, the delta geometry changed, in particular associated with a displacement of water bottom channels. In all our seismic sections, we do not observe fault in the Quaternary sediments that would attest postglacial tectonic activity in the basement.

Paleohydrological Reconstruction from Late Holocene Records in Interdune Lakes (N'Guigmi, Northern Bank of the Lake Chad, Niger)

An old erg covers the northern part of the Lake Chad basin. This dune landform allowed the formation of many inter- dune ponds of various sizes. Still present in certain zones where the groundwater level is high (e.g. Kanem, southern Manga), these ponds formed in the past a vast network of lacustrine microsystems, as shown by the nature and the dis- tribution of their deposits. In the Manga, these interdune deposits represent the main sedimentary records of the Holo- cene environmental succession. Their paleobiological (pollens, diatoms, ostracods) and geochemical (δ18O, δ13C, Sr/ Ca) contents are often the basis for paleoenvironmental reconstruction. On the other hand, their sedimentological char- acters are rarely exploited. This study of palustro-lacustrine deposits of the Holocene N'Guigmi lake (northern bank of the Lake Chad; Niger) is based on the relationships between the sedimentological features and the climato-hydrological fluctuations. The mineralogical parameters (e.g. calcium carbonate content, clay mineralogy) and the nature of autoch- thonous mineralization (i.e. amorphous silica, clays, calcium carbonates) can be interpreted using a straightforward hy- dro-sedimentary model. Established to explain the geochemical dynamics of Lake Chad, this model is based on a bio- geochemical cycle of the main elements (i.e. silicium, calcium) directly controlled by the local hydrological balance (i.e. rainfall/evaporation ratio). All these results show that a detailed study of sedimentological features can provide impor- tant paleohydrological informations about the regional aridification since ca 6500 14C BP.

An evaluation of the inorganic and organic geochemistry of the San Vicente Mississippi Valley-type zinc-lead district, central Peru: Implications for ore fluid composition, mixing processes, and sulfate reduction

Mississippi Tialley-type zinc-lead deposits and ore occurrences in the San Vicente belt are hosted in dolostones of the eastern Upper Triassic to Lower Jurassic Pucara basin, central Peru. Combined inorganic and organic geochemical data from 22 sites, including the main San Vicente deposit, minor ore occurrences, and barren localities, provide better understanding of fluid pathways and composition, ore precipitation mechanisms, Eh-pH changes during mineralization, and relationships between organic matter and ore formation. Ore-stage dark replacement dolomite and white sparry dolomite are Fe and rare earth element (REE) depleted, and Mn enriched, compared to the host dolomite. In the main deposit, they display significant negative Ce and probably Eu anomalies. Mixing of an incoming hot, slightly oxidizing, acidic brine (H2CO3 being the dominant dissolved carbon species), probably poor in REE and Fe, with local intraformational, alkaline, reducing waters explains the overall carbon and oxygen isotope variation and the distributions of REE and other trace elements in the different hydrothermal carbonate generations. The incoming ore fluid flowed through major aquifers, probably basal basin detrital units, with limited interaction with the carbonate host rocks. The hydrothermal carbonates show a strong regional chemical homogeneity, indicating access of the ore fluids by interconnected channelways near the ore occurrences. Negative Ce anomalies in the main deposit, that are absent at the district scale, indicate local ore-fluid chemical differences. Oxidation of both migrated and indigenous hydrocarbons by the incoming fluid provided the local reducing conditions necessary for sulfate reduction to H2S, pyrobitumen precipitation, and reduction of Eu3+ to Eu2+. Fe-Mn covariations, combined with the REE contents of the hydrothermal carbonates, are consistent with the mineralizing system shifting from reducing/rock-dominated to oxidizing/fluid-dominated conditions following ore deposition. Sulfate and sulfide sulfur isotopes support sulfide origin from evaporite-derived sulfate by thermochemical organic reduction; further evidence includes the presence of C-13-depleted calcite cements (similar to-12 parts per thousand delta(13)C) as sulfate pseudomorphs, elemental sulfur, altered organic matter in the host dolomite, and isotopically heavier, late, solid bitumen. Significant alteration of the indigenous and extrinsic hydrocarbons, with absent bacterial membrane biomarkers (hopanes) is observed. The light delta(34)S of sulfides from small mines and occurrences compared to the main deposit reflect a local contribution of isotopically light sulfur, evidence of local differences in the ore-fluid chemistry.

Olivine growth-rates in a tholeiitic basalt - an experimental study of melt inclusions in plagioclase

Basaltic glass inclusions trapped in plagioclase phenocrysts (An84) are remnant of their parent magmatic liquid. They can be used as natural reactors for the experimental investigation of olivine growth rate as a function of temperature. The growth of one olivine nucleus can be observed at constant temperature. Supercooling from 15-degrees to 150-degrees-C have been investigated. Growth habits vary from equant to feather in qualitative agreement with previous studies. Growth rates vary from < 10(-10) m s-1 to 6.10(-7) m s-1; they vary with the direction, the growth process (planar or dendritic) and the degree of supercooling. Chemical analysis of crystal overgrowth and the remaining liquid composition enables a mass-balance calculation which confirms the rates determined optically. The small number of results obtained so far does not permit to draw inferences on the growth mechanisms.

The birth of the Rheic Ocean - Early Palaeozoic subsidence patterns and subsequent tectonic plate scenarios

New plate-tectonic reconstructions of the Gondwana margin suggest that the location of Gondwana-derived terranes should not only be guided by the models, but should also consider the possible detrital input from some Asian blocks (Hunia), supposed to have been located along the Cambrian Gondwana margin, and accreted in the Silurian to the North-Chinese block. Consequently, the Gondwana margin has to be subdivided into a more western domain, where the future Avalonian blocks will be separated from Gondwana by the opening Rheic Ocean, whereas in its eastern continuation, hosting the future basement areas of Central Europe, different periods of crustal extension should be distinguished. Instead of applying a rather cylindrical model, it is supposed that crustal extension follows a much more complex pattern, where local back-arcs or intra-continental rifts are involved. Guided by the age data of magmatic rocks and the pattern of subsidence curves, the following extensional events can be distinguished: During the early to middle Cambrian, a back-arc setting guided the evolution at the Gondwana margin. Contemporaneous intra-continental rift basins developed at other places related to a general post-PanAfrican extensional phase affecting Africa Upper Cambrian formation of oceanic crust is manifested in the Chamrousse area, and may have lateral cryptic relics preserved in other places. This is regarded as the oceanisation of some marginal basins in a context of back-arc rifting. These basins were closed in a mid-Ordovician tectonic phase, related to the subduction of buoyant material (mid-ocean ridge?) Since the Early Ordovician, a new phase of extension is observed, accompanied by a large-scale volcanic activity, erosion of the rift shoulders generated detritus (Armorican Quartzite) and the rift basins collected detrital zircons from a wide hinterland. This phase heralded the opening of Palaeotethys, but it failed due to the Silurian collision (Eo-Variscan phase) of an intra-oceanic arc with the Gondwana margin. During this time period, at the eastern wing of the Gondwana margin begins the drift of the future Hunia microcontinents, through the opening of an eastern prolongation of the already existing Rheic Ocean. The passive margin of the remaining Gondwana was composed of the Galatian superterranes, constituents of the future Variscan basement areas. Remaining under the influence of crustal extension, they will start their drift to Laurussia since the earliest Devonian during the opening of the Palaeotethys Ocean. (C) 2008 Elsevier B.V. All rights reserved.

Magmatic fluids in the Breccia-hosted epithermal Au-Ag deposit of Rosia Montana, Romania

The breccia-hosted epithermal Au-Ag deposit of Rosia Montana is located 7 kin northeast of Abrud, in the northern part of the South Apuseni Mountains, Romania. Estimated total reserves of 214.91 million metric toils (Mt) of ore at 1.46 g/t An and 6.9 g/t Ag (10.1 Moz of An and 47.6 Moz of Ag) make Rosia Montana one of the largest gold deposits in Europe. At this location, Miocene calc-alkaline magmatic and hydrothermal activity was associated with local extensional tectonics within a strike-slip regime related to the indentation of the Adriatic microplate into the European plate during the Carpathian orogenesis. The host rocks of the magmatic complex consist of pre-Mesozoic metamorphosed continental crust covered by Cretaceous turbiditic sediment (flysch). Magmatic activity at Rosia Montana and its surroundings occurred in several pulses and lasted about 7 m.y, Rosia Montana is a breccia-hosted epithermal system related to strong phreatomagmatic activity due to the shallow emplacement of the Montana dacite. The Montana dacite intruded Miocene volcaniclastic material (volcaniclastic breccias) and crops out at Cetate and Carnic Hills. Current mining is focused primarily on the Cetate open pit, which was mapped in detail, leading to the recognition of three distinct breccia bodies: the dacite breccia with a dominantly hydrothermal matrix, the gray polymict breccia with a greater proportion of sand-sized matrix support, and the black polymict breccia, which reached to the surface, contains carbonized tree trunks and has a dominantly barren elastic matrix. The hydrothermal alteration is pervasive. Adularia alteration with a phyllic overprint is ubiquitous; silicification and argillic alteration occur locally. Mineralization consists of quartz, adularia, carbonates (commonly Mn-rich), pyrite, Fe-poor sphalerite, galena, chalcopyrite, tetrahedrite, and native gold and occurs as disseminations, as well as in veins and filling vugs within the Montana dacite and the different breccias. The age of mineralization (12.85 +/- 0.07 Ma) was determined by Ar-40- Ar-39 dating on hydrothermal adularia crystals from vugs in the dacite breccia in the Cetate open pit. Microthermometric measurements of fluid inclusions in quartz phenocrysts from the Montana dacite revealed two fluid types that are absent from the hydrothermal breccia and must have been trapped at depth prior to dacite dome emplacement: brine inclusions (32-55 -wt % NaCl equiv, homogenizing at T-h > 460 degrees C) and intermediate density fluids (4.9-15.6 wt % NaCl equiv, T, between 345 degrees-430 degrees C). Secondary aqueous fluid inclusion assemblages in the phenocrysts have salinities of 0.2 to 2.2 wt percent NaCl equiv and T-h of 200 degrees to 280 degrees C. Fluid inclusion assemblages in hydrothermal quartz from breccias and veins have salinities of 0.2 to 3.4 wt percent NaCl equiv and T-h, from 200 degrees to 270 degrees C. The oxygen isotope composition of several zones of an ore-related epithermal quartz crystal indicate a very constant delta O-18 of 4.5 to 5.0 per mil for the mineralizing fluid, despite significant salinity and temperature variation over time. Following microthermometry, selected fluid inclusion assemblages were analyzed by laser ablation-inductively coupled-plasma mass spectrometry (LA-ICMS). Despite systematic differences in salinity between phenocryst-hosted fluids trapped at depth and fluids from quartz in the epithermal breccias, all fluids have overlapping major and trace cation ratios, including identical Na/K/Rb/Sr/Cs/Ba. Consistent with the constant near-magmatic oxygen isotope composition of the hydrothermal fluids, these data strongly indicate a common magmatic component of these chemically conservative solutes in all fluids. Cu, Pb, Zn, and Mn show variations in concentration relative to the relatively non-reactive alkalis, reflecting the precipitation of sulfide minerals together with An in the epithermal breccia, and possibly of Cu in an inferred subjacent porphyry environment. The magmatic-hydrothermal processes responsible for epithermal Au-Ag mineralization at Rosia Montana are, however, not directly related to the formation of the spatially associated porphyry Cu-Au deposit of Rosia Poieni, which occurred lout 3 m.y. later.

Silver-base metal epithermal vein and listwaenite types of deposit Crnac, Rogozna Mts., Kosovo. Part I: Ore mineral geochemistry and sulfur isotope study

Pb-Zn-Ag vein and listwaenite types of mineralization in Crnac deposit, Western Vardar zone, were deposited within several stages: (i) the pre-ore stage comprises pyrite, arsenopyrite, pyrrhotite, quartz, kaolinite and is followed by magnetite-pyrite; (ii) the syn-ore stage is composed of galena, sphalerite, tetrahedrite and stefanite; and (iii) the post-ore stage is composed of carbonates, pyrite, arsenopyrite and minor galena. The vein type mineralization is hosted by Jurassic amphibolites and veins terminate within overlying serpentinites. Mineralized listwaenites are developed along the serpentinite-amphibolite interface. The reserves are estimated to 1.7 Mt of ore containing in average 7.6% lead, 2.9% zinc, and 102 g/t silver. Sulfides from the pre- and syn-mineralization assemblage of the vein- and listwaenite-types of mineralization from the Crnac Pb-Zn-Ag deposit have been analyzed using microprobe, crush-leachates and sulfur isotopes. The pre-ore assemblage precipitated under high sulfur fugacities (f(S(2)) = 10(-8)-10(-6) bar) from temperatures ranging between 350 degrees C and 380 degrees C. Most likely water-rock reactions, boiling and/or increase of pH caused an increase of delta(34)S of pyrite toward upper levels within the deposit. The decomposition of pre-ore pyrrhotite to a pyrite-magnetite mixture occurred at a fugacity of sulfur from f(S(2)) = 8.7 x 10(-10) to 9.6 x 10(-9) bar and fugacity of oxygen from f(O(2)) = 2.4 x 10(-30) to 3.1 x 10(-28) bars, indicating a contribution of an oxidizing fluid, i.e. meteoric water during pre-ore stages of hydrothermal activity. The crystallization temperatures obtained by the sphalerite-galena isotope geothermometer range from 230 to 310 degrees C. The delta(34)S values of pre- and syn-ore sulfides (pyrite, galena, sphalerite, delta(34)S = 0.3-5.9 parts per thousand) point to magmatic sulfur. Values of delta(34)S of galena and sphalerite are decreasing upwards due to precipitation of early formed sulfide minerals. Post-ore assemblage precipitated at temperature below 190 degrees C. Based on data presented above, we assume two fluid sources: (i) a magmatic source, supported by sulfur isotopic compositions within pre- and syn-ore minerals and a high mol% of fluorine found within pre- and syn-ore leachates, and (ii) a meteoric source, deduced by coincident pyrite-magnetite intergrowth, sulfur isotopic trends within syn-ore minerals and decrease of crystallization temperatures from the pre-ore stage (380-350 degrees C), towards the syn-ore (310-215 degrees C) and post-ore stages (<190 degrees C). Post-ore fluids are Na-Ca-Mg-K-Li chlorine rich and were modified via water-rock reactions. Simple mineral assemblage and sphalerite composition range from 1.5 to 10.1 mol% of FeS catalog Crnac to a group of intermediate sulfidation epithermal deposit. (C) 2011 Elsevier B.V. All rights reserved.

Mid- to Late Holocene shoreline reconstruction and human occupation in Ancient Eretria (South Central Euboea, Greece)

Few studies have aimed to reconstruct landscape change in the area of Eretria (South Central Euboea, Greece) during the last 6000 years. The aim of this paper is to partially fill in this gap by examining the interaction be- tween Mid- to Late Holocene shoreline evolution and human occupation, which is documented in the area from the Late Neolithic to the Late Roman period (with discontinuities). Evidence of shoreline displacements is derived from the study of five boreholes (maximum depth of 5.25 m below the surface) drilled in the lowlands of Eretria. Based on sedimentological analyses and micro/macrofaunal identifications, different facies have been identified in the cores and which reveal typical features of deltaic progradation with marine, lagoonal, fluvio- deltaic and fluvial environments. In addition, a chronostratigraphy has been obtained based on 20 AMS 14C radio- carbon dates performed on samples of plant remains and marine/lagoonal shells found in situ. The main sequences of landscape reconstruction in the plain of Eretria can be summarized as follows: a marine environ- ment predominated from ca. 4000 to 3200 cal. BC and a gradual transition to shallow marine conditions is ob- served ca. 3200-3000 cal. BC due to the general context of deltaic progradation west of the ancient city. Subsequently, from ca. 3000 to 2000 cal. BC, a lagoon occupied the area in the vicinity of the Temple of Apollo and the settlement's development was restricted to several fluvio-deltaic levees, thus severely limiting human activities in the plain. From ca. 2000 to 800 cal. BC, a phase of shallow marine presence prevailed and constrained settlement on higher ground, forcing abandonment of the major part of the plain. Finally, since the eighth century BC, the sea has regressed southward and created the modern landscape.

Genesis of the Jurassic carbonate-hosted Pb-Zn deposits of Jebel Ressas (North-eastern Tunisia): Evidence from mineralogy, petrography and trace metal contents and isotope (O, C, S, Pb) geochemistry

The Jebel Ressas Pb-Zn deposits in North-Eastern Tunisia occur mainly as open-space fillings (lodes, tectonic breccia cements) in bioclastic limestones of the Upper Jurassic Ressas Formation and along the contact of this formation with Triassic rocks. The galena-sphalerite association and their alteration products (cerussite, hemimorphite, hydrozincite) are set within a calcite gangue. The Triassic rocks exhibit enrichments in trace metals, namely Pb, Co and Cd enrichment in clays and Pb, Zn, Cd, Co and Cr enrichment in carbonates, suggesting that the Triassic rocks have interacted with the ore-bearing fluids associated with the Jebel Ressas Pb-Zn deposits. The delta(18)O content of calcite associated with the Pb-Zn mineralization suggests that it is likely to have precipitated from a fluid that was in equilibrium with the Triassic dolostones. The delta(34)S values in galenas from the Pb-Zn deposits range from -1.5 to +11.4%, with an average of 5.9% and standard deviation of 3.9%. These data imply mixing of thermochemically-reduced heavy sulfur carried in geothermal- and fault-stress-driven deep-seated source fluid with bacterially-reduced light sulfur carried in topography-driven meteoric fluid. Lead isotope ratios in galenas from the Pb-Zn deposits are homogenous and indicate a single upper crustal source of base-metals for these deposits. Synthesis of the geochemical data with geological data suggests that the base-metal mineralization at Jebel Ressas was formed during the Serravallian-Tortonian (or Middle-Late Miocene) Alpine compressional tectonics.

Développement de l'interprétation quantitative de la sismique réflexion lacustre : application de l'AVO aux sédiments quaternaires dans le lac Léman

Résumé Des développements antérieurs, au sein de l'Institut de Géophysique de Lausanne, ont permis de développer des techniques d'acquisition sismique et de réaliser l'interprétation des données sismique 2D et 3D pour étudier la géologie de la région et notamment les différentes séquences sédimentaires du Lac Léman. Pour permettre un interprétation quantitative de la sismique en déterminant des paramètres physiques des sédiments la méthode AVO (Amplitude Versus Offset) a été appliquée. Deux campagnes sismiques lacustres, 2D et 3D, ont été acquises afin de tester la méthode AVO dans le Grand Lac sur les deltas des rivières. La géométrie d'acquisition a été repensée afin de pouvoir enregistrer les données à grands déports. Les flûtes sismiques, mises bout à bout, ont permis d'atteindre des angles d'incidence d'environ 40˚ . Des récepteurs GPS spécialement développés à cet effet, et disposés le long de la flûte, ont permis, après post-traitement des données, de déterminer la position de la flûte avec précision (± 0.5 m). L'étalonnage de nos hydrophones, réalisé dans une chambre anéchoïque, a permis de connaître leur réponse en amplitude en fonction de la fréquence. Une variation maximale de 10 dB a été mis en évidence entre les capteurs des flûtes et le signal de référence. Un traitement sismique dont l'amplitude a été conservée a été appliqué sur les données du lac. L'utilisation de l'algorithme en surface en consistante a permis de corriger les variations d'amplitude des tirs du canon à air. Les sections interceptes et gradients obtenues sur les deltas de l'Aubonne et de la Dranse ont permis de produire des cross-plots. Cette représentation permet de classer les anomalies d'amplitude en fonction du type de sédiments et de leur contenu potentiel en gaz. L'un des attributs qui peut être extrait des données 3D, est l'amplitude de la réflectivité d'une interface sismique. Ceci ajoute une composante quantitative à l'interprétation géologique d'une interface. Le fond d'eau sur le delta de l'Aubonne présente des anomalies en amplitude qui caractérisent les chenaux. L'inversion de l'équation de Zoeppritz par l'algorithme de Levenberg-Marquardt a été programmée afin d'extraire les paramètres physiques des sédiments sur ce delta. Une étude statistique des résultats de l'inversion permet de simuler la variation de l'amplitude en fonction du déport. On a obtenu un modèle dont la première couche est l'eau et dont la seconde est une couche pour laquelle V P = 1461 m∕s, ρ = 1182 kg∕m3 et V S = 383 m∕s. Abstract A system to record very high resolution (VHR) seismic data on lakes in 2D and 3D was developed at the Institute of Geophysics, University of Lausanne. Several seismic surveys carried out on Lake Geneva helped us to better understand the geology of the area and to identify sedimentary sequences. However, more sophisticated analysis of the data such as the AVO (Amplitude Versus Offset) method provides means of deciphering the detailed structure of the complex Quaternary sedimentary fill of the Lake Geneva trough. To study the physical parameters we applied the AVO method at some selected places of sediments. These areas are the Aubonne and Dranse River deltas where the configurations of the strata are relatively smooth and the discontinuities between them easy to pick. A specific layout was developed to acquire large incidence angle. 2D and 3D seismic data were acquired with streamers, deployed end to end, providing incidence angle up to 40˚ . One or more GPS antennas attached to the streamer enabled us to calculate individual hydrophone positions with an accuracy of 50 cm after post-processing of the navigation data. To ensure that our system provides correct amplitude information, our streamer sensors were calibrated in an anechoic chamber using a loudspeaker as a source. Amplitude variations between the each hydrophone were of the order of 10 dB. An amplitude correction for each hydrophone was computed and applied before processing. Amplitude preserving processing was then carried out. Intercept vs. gradient cross-plots enable us to determine that both geological discontinuities (lacustrine sediments/moraine and moraine/molasse) have well defined trends. A 3D volume collected on the Aubonne river delta was processed in order ro obtain AVO attributes. Quantitative interpretation using amplitude maps were produced and amplitude maps revealed high reflectivity in channels. Inversion of the water bottom of the Zoeppritz equation using the Levenberg-Marquadt algorithm was carried out to estimate V P , V S and ρ of sediments immediately under the lake bottom. Real-data inversion gave, under the water layer, a mud layer with V P = 1461 m∕s, ρ = 1182 kg∕m3 et V S = 383 m∕s.

Robust Nonlinear Mapping of Soil Contamination Using Support Regression

Spatial data analysis mapping and visualization is of great importance in various fields: environment, pollution, natural hazards and risks, epidemiology, spatial econometrics, etc. A basic task of spatial mapping is to make predictions based on some empirical data (measurements). A number of state-of-the-art methods can be used for the task: deterministic interpolations, methods of geostatistics: the family of kriging estimators (Deutsch and Journel, 1997), machine learning algorithms such as artificial neural networks (ANN) of different architectures, hybrid ANN-geostatistics models (Kanevski and Maignan, 2004; Kanevski et al., 1996), etc. All the methods mentioned above can be used for solving the problem of spatial data mapping. Environmental empirical data are always contaminated/corrupted by noise, and often with noise of unknown nature. That's one of the reasons why deterministic models can be inconsistent, since they treat the measurements as values of some unknown function that should be interpolated. Kriging estimators treat the measurements as the realization of some spatial randomn process. To obtain the estimation with kriging one has to model the spatial structure of the data: spatial correlation function or (semi-)variogram. This task can be complicated if there is not sufficient number of measurements and variogram is sensitive to outliers and extremes. ANN is a powerful tool, but it also suffers from the number of reasons. of a special type ? multiplayer perceptrons ? are often used as a detrending tool in hybrid (ANN+geostatistics) models (Kanevski and Maignank, 2004). Therefore, development and adaptation of the method that would be nonlinear and robust to noise in measurements, would deal with the small empirical datasets and which has solid mathematical background is of great importance. The present paper deals with such model, based on Statistical Learning Theory (SLT) - Support Vector Regression. SLT is a general mathematical framework devoted to the problem of estimation of the dependencies from empirical data (Hastie et al, 2004; Vapnik, 1998). SLT models for classification - Support Vector Machines - have shown good results on different machine learning tasks. The results of SVM classification of spatial data are also promising (Kanevski et al, 2002). The properties of SVM for regression - Support Vector Regression (SVR) are less studied. First results of the application of SVR for spatial mapping of physical quantities were obtained by the authorsin for mapping of medium porosity (Kanevski et al, 1999), and for mapping of radioactively contaminated territories (Kanevski and Canu, 2000). The present paper is devoted to further understanding of the properties of SVR model for spatial data analysis and mapping. Detailed description of the SVR theory can be found in (Cristianini and Shawe-Taylor, 2000; Smola, 1996) and basic equations for the nonlinear modeling are given in section 2. Section 3 discusses the application of SVR for spatial data mapping on the real case study - soil pollution by Cs137 radionuclide. Section 4 discusses the properties of the modelapplied to noised data or data with outliers.

The thermal circuits of the Argentera Massif (western Alps, Italy): An example of low-enthalpy geothermal resources controlled by Neogene alpine tectonics

This work presents geochemistry and structural geology data concerning the low enthalpy geothermal circuits of the Argentera crystalline Massif in northwestern Italian Alps. I n this area some thermal springs (50-60 degreesC), located in the small Bagni di Vinadio village, discharge mixtures made up of a Na-Cl end-member and a Na-SO4 component. The latter is also discharged by the thermal springs of Terme di Valdieri located some kilometres apart within the same tectonic complex. Both end-members share the same meteoric origin and the same reservoir temperature, which is close to 150 degreesC. Explanations are thus required to understand how they reach the surface and how waters of the same origin and circulating in similar rocks can attain such different compositions. Sodium-sulphate waters discharged at both sites, likely represent the common interaction product of meteoric waters with the widespread granitic-migmatitic rocks of the Argentera Massif, whereas Na-CI waters originate through leaching of mineralised cataclastic rocks, which are rich in phyllosilicatic minerals and fluid inclusions, both acting as Cl- sources. Due to the relatively low inferred geothermal gradient of the region, -25C/km, meteoric waters have to descend to depths of 5.5-6 km to attain temperatures of similar to 150 degreesC. These relevant depths can be reached by descending meteoric waters, due to the recent extensional stress field, which allows the development of geothermal circulations at greater depths than in other sectors of the Alps by favouring a greater fractures aperture. The ascent of the thermal waters rakes place along brittle shear zones. In both sites, the thermal waters emerge at the bottoms of the valleys, close to either the lateral termination of a brittle shear zone at Terme di Valdieri, or a step-over between two en-echelon brittle shear zones at Bagni di Vinadio. These observations attest to a strong control operated on the location of outlet regions by both brittle tectonics and the minima in hydraulic potential inside the fractured massif.