Méthodes sismiques actives adaptées à l'étude de milieux hétérogènes : application aux glissements de terrain

Résumé Les glissements de terrain représentent un des principaux risques naturels dans les régions montagneuses. En Suisse, chaque année les glissements de terrains causent des dégâts qui affectent les infrastructures et ont des coûts financiers importants. Une bonne compréhension des mécanismes des glissements peut permettre d'atténuer leur impact. Celle-ci passe notamment par la connaissance de la structure interne du glissement, la détermination de son volume et de son ou ses plans de glissement. Dans un glissement de terrain, la désorganisation et la présence de fractures dans le matériel déplacé engendre un changement des paramètres physiques et en particulier une diminution des vitesses de propagation des ondes sismiques ainsi que de la densité du matériel. Les méthodes sismiques sont de ce fait bien adaptées à l'étude des glissements de terrain. Parmi les méthodes sismiques, l'analyse de la dispersion des ondes de surface est une méthode simple à mettre en oeuvre. Elle présente l'avantage d'estimer les variations des vitesses de cisaillement avec la profondeur sans avoir spécifiquement recours à l'utilisation d'une source d'onde S et de géophones horizontaux. Sa mise en oeuvre en trois étapes implique la mesure de la dispersion des ondes de surface sur des réseaux étendus, la détermination des courbes de dispersion pour finir par l'inversion de ces courbes. Les modèles de vitesse obtenus à partir de cette procédure ne sont valides que lorsque les milieux explorés ne présentent pas de variations latérales. En pratique cette hypothèse est rarement vérifiée, notamment pour un glissement de terrain dans lequel les couches remaniées sont susceptibles de présenter de fortes hétérogénéités latérales. Pour évaluer la possibilité de déterminer des courbes de dispersion à partir de réseaux de faible extension des mesures testes ont été effectuées sur un site (Arnex, VD) équipé d'un forage. Un profil sismique de 190 m de long a été implanté dans une vallée creusée dans du calcaire et remplie par des dépôts glacio-lacustres d'une trentaine de mètres d'épaisseur. Les données acquises le long de ce profil ont confirmé que la présence de variations latérales sous le réseau de géophones affecte l'allure des courbes de dispersion jusqu'à parfois empêcher leur détermination. Pour utiliser l'analyse de la dispersion des ondes de surface sur des sites présentant des variations latérales, notre approche consiste à déterminer les courbes de dispersions pour une série de réseaux de faible extension, à inverser chacune des courbes et à interpoler les différents modèles de vitesse obtenus. Le choix de la position ainsi que de l'extension des différents réseaux de géophones est important. Il tient compte de la localisation des hétérogénéités détectées à partir de l'analyse de sismique réfraction, mais également d'anomalies d'amplitudes observées sur des cartes qui représentent dans le domaine position de tir - position du récepteur, l'amplitude mesurée pour différentes fréquences. La procédure proposée par Lin et Lin (2007) s'est avérée être une méthode efficace permettant de déterminer des courbes de dispersion à partir de réseaux de faible extension. Elle consiste à construire à partir d'un réseau de géophones et de plusieurs positions de tir un enregistrement temps-déports qui tient compte d'une large gamme de distances source-récepteur. Au moment d'assembler les différentes données une correction de phase est appliquée pour tenir compte des hétérogénéités situées entre les différents points de tir. Pour évaluer cette correction nous suggérons de calculer pour deux tir successif la densité spectrale croisée des traces de même offset: Sur le site d'Arnex, 22 courbes de dispersions ont été déterminées pour de réseaux de géophones de 10 m d'extension. Nous avons également profité du forage pour acquérir un profil de sismique verticale en ondes S. Le modèle de vitesse S déduit de l'interprétation du profil de sismique verticale est utilisé comme information à priori lors l'inversion des différentes courbes de dispersion. Finalement, le modèle en deux dimension qui a été établi grâce à l'analyse de la dispersion des ondes de surface met en évidence une structure tabulaire à trois couches dont les limites coïncident bien avec les limites lithologiques observées dans le forage. Dans celui-ci des argiles limoneuses associées à une vitesse de propagation des ondes S de l'ordre de 175 m/s surmontent vers 9 m de profondeur des dépôts de moraine argilo-sableuse caractérisés par des vitesses de propagation des ondes S de l'ordre de 300 m/s jusqu'à 14 m de profondeur et supérieur ou égal à 400 m/s entre 14 et 20 m de profondeur. Le glissement de la Grande Combe (Ballaigues, VD) se produit à l'intérieur du remplissage quaternaire d'une combe creusée dans des calcaires Portlandien. Comme dans le cas du site d'Arnex les dépôts quaternaires correspondent à des dépôts glacio-lacustres. Dans la partie supérieure la surface de glissement a été localisée à une vingtaine de mètres de profondeur au niveau de l'interface qui sépare des dépôts de moraine jurassienne et des dépôts glacio-lacustres. Au pied du glissement 14 courbes de dispersions ont été déterminées sur des réseaux de 10 m d'extension le long d'un profil de 144 m. Les courbes obtenues sont discontinues et définies pour un domaine de fréquence de 7 à 35 Hz. Grâce à l'utilisation de distances source-récepteur entre 8 et 72 m, 2 à 4 modes de propagation ont été identifiés pour chacune des courbes. Lors de l'inversion des courbes de dispersion la prise en compte des différents modes de propagation a permis d'étendre la profondeur d'investigation jusqu'à une vingtaine de mètres de profondeur. Le modèle en deux dimensions permet de distinguer 4 couches (Vs1 < 175 m/s, 175 m/s < Vs2 < 225 m/s, 225 m/s < Vs3 < 400 m/s et Vs4 >.400 m/s) qui présentent des variations d'épaisseur. Des profils de sismiques réflexion en ondes S acquis avec une source construite dans le cadre de ce travail, complètent et corroborent le modèle établi à partir de l'analyse de la dispersion des ondes de surface. Un réflecteur localisé entre 5 et 10 m de profondeur et associé à une vitesse de sommation de 180 m/s souligne notamment la géométrie de l'interface qui sépare la deuxième de la troisième couche du modèle établi à partir de l'analyse de la dispersion des ondes de surface. Abstract Landslides are one of the main natural hazards in mountainous regions. In Switzerland, landslides cause damages every year that impact infrastructures and have important financial costs. In depth understanding of sliding mechanisms may help limiting their impact. In particular, this can be achieved through a better knowledge of the internal structure of the landslide, the determination of its volume and its sliding surface or surfaces In a landslide, the disorganization and the presence of fractures in the displaced material generate a change of the physical parameters and in particular a decrease of the seismic velocities and of the material density. Therefoe, seismic methods are well adapted to the study of landslides. Among seismic methods, surface-wave dispersion analysis is a easy to implement. Through it, shearwave velocity variations with depth can be estimated without having to resort to an S-wave source and to horizontal geophones. Its 3-step implementation implies measurement of surface-wave dispersion with long arrays, determination of the dispersion curves and finally inversion of these curves. Velocity models obtained through this approach are only valid when the investigated medium does not include lateral variations. In practice, this assumption is seldom correct, in particular for landslides in which reshaped layers likely include strong lateral heterogeneities. To assess the possibility of determining dispersion curves from short array lengths we carried out tests measurements on a site (Arnex, VD) that includes a borehole. A 190 m long seismic profile was acquired in a valley carved into limestone and filled with 30 m of glacio-lacustrine sediments. The data acquired along this profile confirmed that the presence of lateral variations under the geophone array influences the dispersion-curve shape so much that it sometimes preventes the dispersion curves determination. Our approach to use the analysis of surface-wave dispersion on sites that include lateral variations consists in obtaining dispersion curves for a series of short length arrays; inverting each so obtained curve and interpolating the different obtained velocity model. The choice of the location as well as the geophone array length is important. It takes into account the location of the heterogeneities that are revealed by the seismic refraction interpretation of the data but also, the location of signal amplitude anomalies observed on maps that represent, for a given frequency, the measured amplitude in the shot position - receiver position domain. The procedure proposed by Lin and Lin (2007) turned out to be an efficient one to determine dispersion curves using short extension arrays. It consists in building a time-offset from an array of geophones with a wide offset range by gathering seismograms acquired with different source-to-receiver offsets. When assembling the different data, a phase correction is applied in order to reduce static phase error induced by lateral variation. To evaluate this correction, we suggest to calculate, for two successive shots, the cross power spectral density of common offset traces. On the Arnex site, 22 curves were determined with 10m in length geophone-arrays. We also took advantage of the borehole to acquire a S-wave vertical seismic profile. The S-wave velocity depth model derived from the vertical seismic profile interpretation is used as prior information in the inversion of the dispersion-curves. Finally a 2D velocity model was established from the analysis of the different dispersion curves. It reveals a 3-layer structure in good agreement with the observed lithologies in the borehole. In it a clay layer with a shear-wave of 175 m/s shear-wave velocity overlies a clayey-sandy till layer at 9 m depth that is characterized down to 14 m by a 300 m/s S-wave velocity; these deposits have a S-wave velocity of 400 m/s between depths of 14 to 20 m. The La Grand Combe landslide (Ballaigues, VD) occurs inside the Quaternary filling of a valley carved into Portlandien limestone. As at the Arnex site, the Quaternary deposits correspond to glaciolacustrine sediments. In the upper part of the landslide, the sliding surface is located at a depth of about 20 m that coincides with the discontinuity between Jurassian till and glacio-lacustrine deposits. At the toe of the landslide, we defined 14 dispersion curves along a 144 m long profile using 10 m long geophone arrays. The obtained curves are discontinuous and defined within a frequency range of 7 to 35 Hz. The use of a wide range of offsets (from 8 to 72 m) enabled us to determine 2 to 4 mode of propagation for each dispersion curve. Taking these higher modes into consideration for dispersion curve inversion allowed us to reach an investigation depth of about 20 m. A four layer 2D model was derived (Vs1< 175 m/s, 175 m/s <Vs2< 225 m/s, 225 m/s < Vs3 < 400 m/s, Vs4> 400 m/s) with variable layer thicknesses. S-wave seismic reflection profiles acquired with a source built as part of this work complete and the velocity model revealed by surface-wave analysis. In particular, reflector at a depth of 5 to 10 m associated with a 180 m/s stacking velocity image the geometry of the discontinuity between the second and third layer of the model derived from the surface-wave dispersion analysis.

Investigación geofísica de los paleovalles pleistocenos del NE de Cataluña: el ejemplo de la Riera de Tossa (Costa Brava meridional)

The geometry of Riera de Tossa (Costa Brava, Catalonia) Pleistocene paleovalley has been studied by two geophysical prospecting methods: vertical electric sounding (VES-DC) and seismic refraction. The results show the existence of a sudden slope change of the thalweg that permits to differentiate two sections where erosion processes operated in unequal intensity. Based on geological and geomorphological available data, it can be inferred that this slope rupture is caused by the litological control imposed by the alineation of resistant materials, which limited the headward erosion extend of the creek during the Upper Pleistocene

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

Um estudo das condições de percolação e estabilidade em barragens de terra mediante métodos geofísicos: caso do dique de Sant Llorenç de Montgai - Espanha

The monitoring of Earth dam makes use of visual inspection and instrumentation to identify and characterize the deterioration that compromises the security of earth dams and associated structures. The visual inspection is subjective and can lead to misinterpretation or omission of important information and, some problems are detected too late. The instrumentation are efficient but certain technical or operational issues can cause restrictions. Thereby, visual inspections and instrumentation can lead to a lack of information. Geophysics offers consolidated, low-cost methods that are non-invasive, non-destructive and low cost. They have a strong potential and can be used assisting instrumentation. In the case that a visual inspection and strumentation does not provide all the necessary information, geophysical methods would provide more complete and relevant information. In order to test these theories, geophysical acquisitions were performed using Georadar (GPR), Electric resistivity, Seismic refraction, and Refraction Microtremor (ReMi) on the dike of the dam in Sant Llorenç de Montgai, located in the province of Lleida, 145 km from Barcelona, Catalonia. The results confirmed that the geophysical methods used each responded satisfactorily to the conditions of the earth dike, the anomalies present and the geological features found, such as alluvium and carbonate and evaporite rocks. It has also been confirmed that these methods, when used in an integrated manner, are able to reduce the ambiguities in individual interpretations. They facilitate improved imaging of the interior dikes and of major geological features, thus inspecting the massif and its foundation. Consequently, the results obtained in this study demonstrated that these geophysical methods are sufficiently effective for inspecting earth dams and they are an important tool in the instrumentation and visual inspection of the security of the dams

Análise da compacidade do solo na área urbana de Bauru/ Estado de São Paulo utilizando Ensaios de Penetração Padronizado e sísmica de refração

This study was done to determine compactness of soil strata in the urban area of Bauru, central São Paulo State using standard penetration test (SPT) and seismic refraction methods. Frequency distribution within boreholes was measured every meter to a depth of 30 meters. The geometric mean of standard penetration test values was computed every meter. A biplot graph showing the relationship of SPT with depth indicates soil stratification. In addition, five sections of seismic refraction values, which also indicate the compactness of soil strata, were obtained. A comparison of the two methods shows significant correlation between the results obtained from each.

Inversão de dados de sísmica de refração profunda a partir da curva tempo-distância

O trabalho em pauta tem como objetivo o modelamento da crosta, através da inversão de dados de refração sísmica profunda, segundo camadas planas horizontais lateralmente homogêneas, sobre um semi-espaço. O modelo direto é dado pela expressão analítica da curva tempo-distância como uma função que depende da distância fonte-estação e do vetor de parâmetros velocidades e espessuras de cada camada, calculado segundo as trajetórias do raio sísmico, regidas pela Lei de Snell. O cálculo dos tempos de chegada por este procedimento, exige a utilização de um modelo cujas velocidades sejam crescentes com a profundidade, de modo que a ocorrência das camadas de baixa velocidade (CBV) é contornada pela reparametrização do modelo, levando-se em conta o fato de que o topo da CBV funciona apenas como um refletor do raio sísmico, e não como refrator. A metodologia de inversão utilizada tem em vista não só a determinação das soluções possíveis, mas também a realização de uma análise sobre as causas responsáveis pela ambiguidade do problema. A região de pesquisa das prováveis soluções é vinculada segundo limites superiores e inferiores para cada parâmetro procurado, e pelo estabelecimento de limites superiores para os valores de distâncias críticas, calculadas a partir do vetor de parâmetros. O processo de inversão é feito utilizando-se uma técnica de otimização do ajuste de curvas através da busca direta no espaço dos parâmetros, denominado COMPLEX. Esta técnica apresenta a vantagem de poder ser utilizada com qualquer função objeto, e ser bastante prática na obtenção de múltiplas soluções do problema. Devido a curva tempo-distância corresponder ao caso de uma multi-função, o algoritmo foi adaptado de modo a minimizar simultaneamente várias funções objetos, com vínculos nos parâmetros. A inversão é feita de modo a se obter um conjunto de soluções representativas do universo existente. Por sua vez, a análise da ambiguidade é realizada pela análise fatorial modo-Q, através da qual é possível se caracterizar as propriedades comuns existentes no elenco das soluções analisadas. Os testes com dados sintéticos e reais foram feitos tendo como aproximação inicial ao processo de inversão, os valores de velocidades e espessuras calculados diretamente da interpretação visual do sismograma. Para a realização dos primeiros, utilizou-se sismogramas calculados pelo método da refletividade, segundo diferentes modelos. Por sua vez, os testes com dados reais foram realizados utilizando-se dados extraídos de um dos sismogramas coletados pelo projeto Lithospheric Seismic Profile in Britain (LISPB), na região norte da Grã-Bretanha. Em todos os testes foi verificado que a geometria do modelo possui um maior peso na ambiguidade do problema, enquanto os parâmetros físicos apresentam apenas suaves variações, no conjunto das soluções obtidas.

Métodos geofísicos aplicados a prospecção de água subterrânea no município de Ponta de Pedras - Ilha de Marajó - Pará

Este trabalho foi realizado numa área de 25 km² em Ponta de Pedras, Ilha de Marajó. Tem por objetivo o dimensionamento geométrico e seleção de zonas favoráveis à água subterrânea. Foram utilizados dois métodos geofísicos: sísmica de refração e eletroresistividade. O método sísmico foi empregado em caráter experimental, visando verificar a viabilidade de sua aplicação na área. Os resultados da interpretação confirmaram a suspeita inicial, da existência de horizontes geológicos pouco profundos, não possíveis de serem distinguidos pelo método sísmico. Com o método de eletroresistividade foram feitas sondagens verticais e perfis horizontais. Os resultados da interpretação indicaram áreas promissoras na parte central e numa faixa ao norte da área prospectada. A espessura do aquífero raso varia entre 0 e 50 m.

Geophysical Analysis of Seasonal Montandon Gravel Ridge Water Table Fluctuation and Moisture Gradient Variation due to Storm Events

The purpose of this research project is to continue exploring the Montandon Long-Term Hydrologic Research Site(LTHR) by using multiple geophysical methods to obtain more accurate and precise information regarding subsurface hydrologic properties of a local gravel ridge,which are important to both the health of surrounding ecosystems and local agriculture. Through using non-invasive geophysical methods such as seismic refraction, Direct Current resistivity and ground penetrating radar (GPR) instead of invasive methods such as boreholedrilling which displace sediment and may alter water flow, data collection is less likely to bias the data itself. In addition to imaging the gravel ridge subsurface, another important researchpurpose is to observe how both water table elevation and the moisture gradient (moisture content of the unsaturated zone) change over a seasonal time period and directly after storm events. The combination of three types of data collection allows the strengths of each method combine together and provide a relatively strongly supported conclusions compared to previous research. Precipitation and geophysical data suggest that an overall increase in precipitation during the summer months causes a sharp decrease in subsurface resistivity within the unsaturated zone. GPR velocity data indicate significant immediate increase in moisture content within the shallow vadose zone (< 1m), suggesting that rain water was infiltrating into the shallow subsurface. Furthermore, the combination of resistivity and GPR results suggest that the decreased resistivity within the shallow layers is due to increased ion content within groundwater. This is unexpected as rainwater is assumed to have a DC resistivity value of 3.33*105 ohm-m. These results may suggest that ions within the sediment must beincorporated into the infiltrating water.

Velocities and elastic moduli of basalts from Pacific DSDP sites

Compressional (Vp) and shear (Vs) wave velocities have been measured to 10 kb in 32 cores of basalt from 14 Pacific sites of the Deep Sea Drilling Project. Both Vp and V s show wide ranges (3.70 to 6.38 km/sec for Vp and 1.77 to 3.40 km/sec for V s at 0.5 kb) which are linearly related to density and sea floor age, confirming earlier findings by Christensen and Salisbury of decreasing velocity with progressive submarine weathering based on studies of basalts from five sites in the Atlantic. Combined Pacific and Atlantic data give rates of decreasing velocity of -1.89 and -1.35 km/sec per 100 my for Vp and Vs respectively. New analyses of oceanic seismic refraction data indicate a decrease in layer 2 velocities with age similar to that observed in the laboratory, suggesting that weathering penetrates to several hundred meters in many regions and is largely responsible for the extreme range and variability of layer 2 refraction velocities.

Compressional and shear wave velocities and elastic constants at DSDP Leg 59 Holes

Prior to the Deep Sea Drilling Project the composition of the oceanic crust could only be inferred from seismic-refraction and gravity data and the recovery of a wide variety of dredged rocks. Through the success of the Deep Sea Drilling Project, it is now clear that the top of oceanic Layer 2 usually consists of basalt. Several laboratory studies (e.g., Fox et al., 1972; Christensen and Shaw, 1970; Hyndman and Drury, 1976) have demonstrated that the seismic velocities of oceanic basalt are similar to velocities reported from refraction studies of Layer 2 and that the variability in Layer 2 velocities has many causes, the most important being fracturing and sea-floor alteration produced by the interaction of basalt and sea water (Christensen and Salisbury, 1973). To date, most reported measurements of velocities in oceanic basalts are from samples obtained from the main ocean basins. With the exception of an earlier study of velocities and related elastic properties of a suite of rocks from DSDP Sites 292, 293, 294, and 296 located in the Philippine Sea (Christensen et al., 1975; Fountain et al., 1975), elastic properties have not been determined for oceanic rocks from marginal basins. In this chapter compressional- and shear-wave velocities and elastic constants are reported at elevated confining pressures for basalt and volcanic breccias from Holes 447A, 448, and 448A.

Velocity-density systematics for basalts at DSPD Leg 54 Holes

Precise velocity and density measurements at atmospheric and elevated pressures have been obtained on basalts drilled in 13 holes during Leg 54. The vp-sigma data show systematics which are controlled by the degree of crystallinity (or conversely, glassiness), microstructure, and original mineralogy and chemistry. Extensively fractionated basalts with marked iron enrichment produce anomalously low measured velocities at densities above 2.90 g/cm**3. Also, the effective in situ pressure acting on Leg 54 basalts is less than hydrostatic, and perhaps close to zero. At zero effective pressure, the measured velocities average 2.5 km/s higher than East Pacific Rise upper crustal velocities determined by seismic refraction. This implies that the in situ velocities are undoubtedly a result of the highly fragmented nature of East Pacific Rise crust.

Basalt density, basement age, and intrusive/extrusive relations from DSDP Legs 2 through 7, 9, and 14

Densities of layer 2 basalt recovered during the Deep Sea Drilling Project have been found to decrease steadily with age, a finding ascribed to progressive submarine weathering in the context of sea-floor spreading. The least-squares solution for 52 density measurements gives a rate of decrease in density of (Delta p)/(Delta t) = -0.0046 g per ccm m.y. = -16 percent per 100 m.y., which is in excellent agreement with earlier estimates based on observed chemical depletion rates of dredged oceanic basalt. Weathering of sea-floor basalt, should it penetrate to any considerable depth in layer 2, will decrease layer 2 seismic refraction velocities, act as a source of geothermal heat, and substantially influence the chemistry of sea water and the overlying column of sediment.

Magnetic and gravity measurements on North-South profiles across the Grosse Meteor Bank, North Atlantic

During the first section of the "Meteor" cruise No. 2 a profile was run from the Azores to the south across the flanks of the Mid-Atlantic Ridge with a chain of seamounts. The profile extended between the Cruiser (living) and the Hyeres seamounts, which, according to our soundings, form a connected massif, and across the centre of the Grosse Meteor Bank (30°N, 28.5 °W). These seamounts rise from a depth of more than 4000 m up till close to the surface of the sea forming there a large almost flat plateau. In the case of the Grosse Meteor Bank, this plateau has a N-S extension of approx. 30 nautical miles and an E-W extension of approx. 20 nautical miles and reaches a height of 275 m in water depth. The gravity measurements yielded a density of the topographic masses of 2.6 g/cm**3 for the Grosse Meteor Bank. Magnitude and shape of the measured free-air anomaly are very well shown in a model computation with this density. The theoretical gravity effects of the seismically detected swell of cristalline rock and of the Moho depression (mountain root) are not indicated by the observational data. It can, therefore, be assumed that the latter two neutralize each other. It seems, accordingly, that there is no local isostatic compensation of the topographic masses. Hence, the density of 2.6 g/cm**3 obtained would be about the true density of rock. In connection with the mean velocity of P waves (Aric et al., 1968) obtained by seismic refraction methods it must be concluded that the material of the 1200-4000 m thick surface layer of the Grosse Meteor Bank consists of consolidated sediments. This finding is supported by the total intensity of the Earth's magnetic field over the Grosse Meteor Bank. On the assumption of a homogeneous magnetization in the direction of the present Earth's field, the computed anomaly of the massif deviates considerably from the measured anomaly while the magnetic field of the seismically detected crystalline body is capable of interpreting the observed data. Deviating from the prevailing interpretation of the seamounts' plateau as a volcanic cone with submarine abrasion, the Grosse Meteor Bank and the seamounts in the vicinity are assumed to be of continental origin. The questions whether these seamounts submerged later on or whether the sealevel has risen subsequently are, therefore, largely nonexistent.

Permafrost mapping on glaciers in the Ötztal, Austria

The occurrence of permafrost in the region of the Hochebenkar rock glaciers has been mapped in detail. For this purpose basal temperatures of the winter snow cover were measured at over 250 sites (BTS-method), II refraction-seismic profiles were taken in frozen and unfrozen unconsolidated sediments and 12 springs were investigated in terms of their summer temperature 128 W. Haeberli und G. PatzeIt variations. The combination of seismic refraction and the BTS-method allows rapid and reliable mapping of alpine permafrost, and at the same time enables differences in active-Iayer thickness to be established. The resuIts of the observations confirm the ideas developed in the Swiss Alps in recent years about the relationships between permafrost and rock glaciers: rock glaciers are creep phenomena of discontinuous alpine permafrost.

The TOMO-ETNA experiment: an imaging active campaign at Mt. Etna volcano. Context, main objectives, working-plans and involved research projects

The TOMO-ETNA experiment was devised to image of the crust underlying the volcanic edifice and, possibly, its plumbing system by using passive and active refraction/reflection seismic methods. This experiment included activities both on-land and offshore with the main objective of obtaining a new high-resolution seismic tomography to improve the knowledge of the crustal structures existing beneath the Etna volcano and northeast Sicily up to Aeolian Islands. The TOMO ETNA experiment was divided in two phases. The first phase started on June 15, 2014 and finalized on July 24, 2014, with the withdrawal of two removable seismic networks (a Short Period Network and a Broadband network composed by 80 and 20 stations respectively) deployed at Etna volcano and surrounding areas. During this first phase the oceanographic research vessel “Sarmiento de Gamboa” and the hydro-oceanographic vessel “Galatea” performed the offshore activities, which includes the deployment of ocean bottom seismometers (OBS), air-gun shooting for Wide Angle Seismic refraction (WAS), Multi-Channel Seismic (MCS) reflection surveys, magnetic surveys and ROV (Remotely Operated Vehicle) dives. This phase finished with the recovery of the short period seismic network. In the second phase the Broadband seismic network remained operative until October 28, 2014, and the R/V “Aegaeo” performed additional MCS surveys during November 19-27, 2014. Overall, the information deriving from TOMO-ETNA experiment could provide the answer to many uncertainties that have arisen while exploiting the large amount of data provided by the cutting-edge monitoring systems of Etna volcano and seismogenic area of eastern Sicily.