A new method to estimate the infilling of alluvial sediment of glacial valleys using a sloping local base level

A new method is used to estimate the volumes of sediments of glacial valleys. This method is based on the concept of sloping local base level and requires only a digital terrain model and the limits of the alluvial valleys as input data. The bedrock surface of the glacial valley is estimated by a progressive excavation of the digital elevation model (DEM) of the filled valley area. This is performed using an iterative routine that replaces the altitude of a point of the DEM by the mean value of its neighbors minus a fixed value. The result is a curved surface, quadratic in 2D. The bedrock surface of the Rhone Valley in Switzerland was estimated by this method using the free digital terrain model Shuttle Radar Topography Mission (SRTM) (~92 m resolution). The results obtained are in good agreement with the previous estimations based on seismic profiles and gravimetric modeling, with the exceptions of some particular locations. The results from the present method and those from the seismic interpretation are slightly different from the results of the gravimetric data. This discrepancy may result from the presence of large buried landslides in the bottom of the Rhone Valley.

Preliminary report of the first season work of the International Pluridisciplinary Archaeological Expedition to Bactria 2007

A prior long-term and complex evaluation of the already available data on the geophysical prospecting during the first season work carried out at 2006, at the archaeological site of Tchinguiz Tepe of Termez, took place to decide the strategy to follow during the campaign of 2007. This previous evaluation of the information, on one hand, leaded to the decision to increase the geophysical prospecting at Tchinguiz Tepe, on the other hand, to decide the exact location of areas where the archaeological interventions.would carry out. The main objective at the beginning of this new season was to crosscheck the reliabilityof the measurements and, at the same time, to establish the unknown up to the present archaeologicaland chronological sequence of Tchinguiz Tepe. Meanwhile, the geophysical prospecting also wasextended to the outskirts of the city were the localisation of an unknown up to now Buddhist Monasterywas possible.

Preliminary report of the first season work of the International Pluridisciplinary Archaeological Expedition to Bactria 2006

The International Pluridisciplinary Archaeological Expedition in Bactria (IPAEB) was created in 2006.The name underlines the international character of the team (which includes Uzbeks, Spanish, French, British and Greek members), the presence of specialists from various fields apart from archaeology and the fame of Bactria1.

Preliminary report of the first season work of the International Pluridisciplinary Archaeological Expedition to Bactria 2008

This year the development of our project can be divided into two main clearly different parts, on one hand the laboratory work, where the sampled ceramic individuals has been prepared and analyzed and the elaboration of the data obtained during the excavation of 2007 has been finished and, on the other hand, the field work developed at the archaeological site during this specific year (2008).In the mark of the analytical work a significant number of ceramic individuals (144) from the different stratigraphical units from various areas of the excavation of Termez and Tchinguiz Tepe sampled duringthe new and previous field works has been archaeometrically characterized. This specific material includedindividuals dated into the Hellenistic and Sassanian period, which has been confirmed by C14dating upon organic samples.At the same time, in the mark of the field work of 2008 the archaeological record, already started tobe under study during the excavation of 2007, has been completed and two new archaeological recordshave been registered on of which is located in the area of Tchinguiz Tepe. For the archeological studythe information of the previous geophysical prospecting has been indisputably taken into considerationand the same methodology has been applied to crosscheck the latter archaeological results.

Preliminary Report of the work of the International Pluridisciplinary Archaeological Expedition to Bactria (2009)

Institute of Archaeology & Institute of Fine Arts. Academy of Science of the Republic of Uzbekistan- Universitat de Barcelona- Ministerio de Cultura (Gobierno de España)- Ministerio de Ciencia e Innovación (Gobierno de España)

Geochemical exploration in the Montseny Mountains (NE Spain)

A five year program of systematic multi-element geochemical exploration of the Catalonian Coastal Ranges has been initiated by the Geological Survey of Autonomic Government of Catalonia (Generalitat de Catalunya) and the Department of Geological and Geophysical Exploration (University of Barcelona). This paper reports the first stage results of this regional survey, covering an area of 530 km2 in the Montseny Mountains, NE of Barcelona (Spain). Stream sediments for metals and stream waters for fluoride were chosen because of the regional characteristics. Four target areas for future tactic survey were recognized after the prospect. The most important is a 40 km* zone in the Canoves-Vilamajor area, with high base metal values accompanied by Cd, Ni, Co, As and Sb anomalies. Keywords: Catalanides. Geochemical exploration. Stream sediments. Base metal anomalies. Principal Component Analysis.

Prospección geoquímica de flúor en las aguas superficiales del macizo Montseny (Barcelona y Girona)

En el Macizo Montseny-Guilleries se conocen antiguas explotaciones de fluorita y fluorita-baritina (Osor, Tagamanent, Espinelves, Sant Marcal, etc). Es por el1o que se ha utilizado el método de prospección geoquímica de flúor en aguas superficiales para detectar posibles zonas de interés hasta ahora desconocidas.

É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.

Lake dwellers occupation gap in Lake Geneva (France-Switzerland) possibly explained by an earthquake-mass movement-tsunami event during Early Bronze Age

High-resolution seismic and sediment core data from the 'Grand Lac' basin of Lake Geneva reveal traces of repeated slope instabilities with one main slide-evolved mass-flow (minimum volume 0.13 km3) that originated from the northern lateral slope of the lake near the city of Lausanne. Radiocarbon dating of organic remains sampled from the top of the main deposit gives an age interval of 1865-1608 BC. This date coincides with the age interval for a mass movement event described in the 'Petit Lac' basin of Lake Geneva (1872-1622 BC). Because multiple mass movements took place at the same time in different parts of the lake, we consider the most likely trigger mechanism to be a strong earthquake (Mw 6) that occurred in the period between 1872 and 1608 BC. Based on numerical simulations, we show the major deposit near Lausanne would have generated a tsunami with local wave heights of up to 6 m. The combined effects of the earthquake and the following tsunami provide a possible explanation for a gap in lake dwellers occupation along the shores of Lake Geneva revealed by dendrochronological dating of two palafitte archaeological sites.

Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Le Briançonnais, terrain exotique dans les Alpes ?

New reconstructions of the Western Alps from late Early Jurassic till early Tertiary are proposed. These reconstructions use deep lithospheric data gathered through recent seismic surveys and tomographic studies carried out in the Alps. The present day position, under the Po plain, of the southern limit of the European plate (fig. 1), allows to define the former geometry of the Brianconnais peninsula. The Brianconnais domain is regarded as an exotic terrane formerly belonging to the European margin until Late Jurassic, then transported eastward during the drift of Iberia (fig. 5). Therefore, on a present day Western Alps cross section, a duplication of the European continental margin can be recognized (fig. 10). Stratigraphic and sedimentological data along a zone linking the Pyrenean fracture zone to the Brianconnais, can be related to a rifting event starting in Oxfordian time. This event is responsible for the Late Jurassic till mid-Cretaceous drift of Iberia opening, first the northern Atlantic, then the Gulf of Biscay. Simultaneously, the drift of the Brianconnais will open the Valais ocean and close the Piemontese ocean. The resulting oblique collision zone between the Brianconnais and the Apulian margin generates HP/LT metamorphism starting in Early Cretaceous. The eastward drift of the Brianconnais peninsula will eventually bring it in front of a more northerly segment of the former European margin. The thrusting of the Brianconnais unto that margin takes place in early Tertiary (fig. 9), following the subduction of the Valais ocean. The present nappe pile results not only from continent/continent frontal collision, but also from important lateral displacement of terranes, the most important one being the Brianconnais. The dilemma of `'en echelon'' oceanic domains in the Alps is an outcome of these translations. A solution is found when considering the opening of a Cretaceous Valais ocean across the European margin, running out eastward into the Piemontese ocean, where the drift is taken up along a former transform fault and compensated by subduction under the Apulian margin (fig. 8). In the Western Alps we are then dealing with two oceans, the Piemontese and the Valaisan and a duplicated European margin. In the Eastern Alps the single Piemontese ocean is cut by newly created oceanic crust. All these elements will be incorporated into the Penninic structural domain which does not represent a former unique paleogeographic area, it is a composite accretionary domain squeezed between Europe and Apulia.

Induced-feeding bioassays for detection of Bacillus thuringiensis insecticidal activity against Epilachna paenulata (Coleoptera)

The objective of this work was to establish and test the induced-feeding bioassay in order to detect Bacillus thuringiensis insecticidal activity against Epilachna paenulata (Coleoptera: Coccinellidae). Larvae were induced to swallow high concentrations of spore-crystal suspensions of seven exotic and 30 Argentine B. thuringiensis strains. The great majority of strains showed no toxicity to E. paenulata larvae, and observed mortality was lower than 30%. Induced-feeding bioassay is feasible, and should be used for prospecting strains that produce right combinations of Cry proteins needed to an efficient pest control.

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.