La détection de cavités en milieu urbain est importante pour prévenir différentes causes d'accidents liés à des possibles effondrements. Les cavités sont aussi des cibles d'intérêts pour les archéologues, car les cavités oubliées sont de potentielles sources de matériel révélant des usages passés. Ces cavités sont de tailles différentes, d'origine anthropique ou non, en milieu extérieur ou sous des bâtiments. Leur taille, ainsi que les propriétés physiques du milieu extérieur dans lequel elles se situent, permettent l'utilisation de différentes méthodes géophysiques. Nous nous sommes concentrés sur l'utilisation de deux méthodes géophysiques, le radar de sol et la tomographie par mesures de résistivité électrique, pour localiser et déterminer les cavités métriques à sub-métrique dans le proche sous-sol (6 premiers mètres). Les mesures de radar de sol sont sensibles aux variations de permittivité diélectrique entre la cavité et le milieu extérieur. Nous montrons par des modélisations numériques un effet sur l'amplitude de la réfléchie en fonction du déport qui permet de discriminer entre une cavité vide et une cavité pleine d'eau. Nous étudions aussi l'amplitude de la réfléchie à incidence normale sur le toit d'une cavité à section carrée en fonction de sa profondeur et de sa taille. Nous mettons en évidence une relation logarithmique profondeur versus taille de cavité pour laquelle l'amplitude de la réfléchie est maximum pour les fréquences de prospection typiques du radar de sol. Par ailleurs nous confirmons qu'alors que les mesures radar permettent de déterminer avec précision les dimensions d'une anomalie dans un sous-sol homogène et peu conducteur, les mesures de résistivité électrique permettent elles de déterminer des zones de hautes résistivités à l'emplacement des cavités. Nous couplons ces deux méthodes géophysiques dans deux études de cas, en utilisant la profondeur des interfaces déterminées sur des radargrammes pour contraindre les modèles de résistivité inversés par l'ajout d'information a priori. / The detection of cavities in urban areas is important to prevent different causes of accidents related to possible collapse. The cavities are also interesting targets to archaeologists because forgotten cavities are potential sources of material revealing past uses. These cavities are of different sizes, of anthropogenic origin or not, in an outdoor setting or under buildings. Their size and the physical properties of the external environment in which they are located, allow the use of different geophysical methods. We focused on the use of two of them, Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT), to locate and determine cavities in the near subsurface (first 6 meters). GPR data are sensitive to variations in dielectric permittivity between the cavity and the external environment. We show by numerical modelling an effect on the amplitude of the reflected signal depending on the offset which could enable discrimination between an empty cavity and a cavity filled with water. We also study the amplitude of the reflected wave at normal incidence on the roof of a cavity of square cross section in terms of its depth and size. We show a logarithmic relationship between the cavity size and its depth at which the amplitude of the reflection is maximum for frequencies of typical exploration with GPR. Furthermore, we confirm that while GPR data determine accurately the size of an anomaly in homogeneous low conductive medium, ERT helps to determine areas with high resistivity at the location of cavities. We combine these two geophysical methods in two case studies, using the depth of interfaces detected on radargrams as a priori information to constrain the inversion of electrical resistivity models.

Cette thèse porte sur la dynamique des grandes échelles des écoulements géophysiques turbulents, en particulier sur leur organisation en écoulements parallèles orientés dans la direction est-ouest (jets zonaux). Ces structures ont la particularité d'évoluer sur des périodes beaucoup plus longues que la turbulence qui les entoure. D'autre part, on observe dans certains cas, sur ces échelles de temps longues, des transitions brutales entre différentes configurations des jets zonaux (multistabilité). L'approche proposée dans cette thèse consiste à moyenner l'effet des degrés de liberté turbulents rapides de manière à obtenir une description effective des grandes échelles spatiales de l'écoulement, en utilisant les outils de moyennisation stochastique et la théorie des grandes déviations. Ces outils permettent d'étudier à la fois les attracteurs, les fluctuations typiques et les fluctuations extrêmes de la dynamique des jets. Cela permet d'aller au-delà des approches antérieures, qui ne décrivent que le comportement moyen des jets.Le premier résultat est une équation effective pour la dynamique lente des jets, la validité de cette équation est étudiée d'un point de vue théorique, et les conséquences physiques sont discutées. De manière à décrire la statistique des évènements rares tels que les transitions brutales entre différentes configurations des jets, des outils issus de la théorie des grandes déviations sont employés. Des méthodes originales sont développées pour mettre en œuvre cette théorie, ces méthodes peuvent par exemple être appliquées à des situations de multistabilité. / This thesis deals with the dynamics of geophysical turbulent flows at large scales, more particularly their organization into east-west parallel flows (zonal jets). These structures have the particularity to evolve much slower than the surrounding turbulence. Besides, over long time scales, abrupt transitions between different configurations of zonal jets are observed in some cases (multistability). Our approach consists in averaging the effect of fast turbulent degrees of freedom in order to obtain an effective description of the large scales of the flow, using stochastic averaging and the theory of large deviations. These tools provide theattractors, the typical fluctuations and the large fluctuations of jet dynamics. This allows to go beyond previous studies, which only describe the average jet dynamics. Our first result is an effective equation for the slow dynamics of jets, the validityof this equation is studied from a theoretical point of view, and the physical consequences are discussed. In order to describe the statistics of rare events such as abrupt transitions between different jet configurations, tools from large deviation theory are employed. Original methods are developped in order to implement this theory, those methods can be applied for instance in situations of multistability.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Devido a sua grande extensão latitudinal, a costa brasileira é influenciada por diferentes regimes climáticos e oceanográficos. Adicionalmente, a distribuição da população brasileira é caracterizada por uma alta concentração nas capitais litorâneas. Todos esses fatores levam à construção de inúmeras estruturas de engenharia que podem de alguma forma impactar o transporte de sedimento e consequentemente o balanço sedimentar de algumas praias. Uma das formas mais eficientes de recuperar esse balanço sedimentar é a alimentação artificial de praias, com sedimentos provenientes da plataforma continental com características semelhantes. Os métodos geofísicos acústicos permitem mapear de forma eficiente o fundo e o subfundo marinho para a busca de áreas fontes para a lavra de forma eficiente. O objetivo deste estudo é a identificação de padrões geoacústicos a partir de analises quantitativas e qualitativas, para a caracterização de áreas fontes de areias siliciclásticas compatíveis com sedimentos de praias em erosão na cidade do Rio de Janeiro. O mapeamento foi realizado na plataforma continental interna do Rio de Janeiro, em área adjacente à praia de Itaipuaçu (Maricá, RJ). As análises evidenciaram depósitos, localizados entre 19 e 30 metros, de areias com granulometria, textura e selecionamento em condições ideias para a recomposição de praias do Rio de Janeiro. / The Brazilian coast is afected by different oceanographic characteristics, mainly due to its large latitudinal extention. Additionally, the distribution of the Brazilian population is characterized by a high concentration in the coastal capital. These facts generate a large number of coastal engineering structures that may impact the sediment transport and consequently the sediment balance in some beaches. All these factors impact some sand beaches and the best choice for its mitigation and recovering is the beach nourishment - with sediments from the continental shelf with similar characteristics. However to find a good place for dredging in the shelf, we need to use geophysical methods that allow to map the bottom and sub bottom environments more efficiently. The objective of this study is the identification of geoacoustics patterns, following qualitative and quantitative analyzes, for characterizing source areas of siliciclastic sands in ideal conditions for recovering eroded beaches in Rio de Janeiro. The mapping was done in the shelf in front of Itaipuaçu beach (Maricá, RJ). The quantitative results showed well sorted median sand placers seaward of the closure depth, located between 19 and 30 meters, ideal for Rios beaches recovering.

Afin de cartographier la structure souterraine des volcans et détecter des ressources géothermiques de haute température, on utilise souvent l’imagerie de résistivité électrique. La résistivité électrique des volcans est affectée par plusieurs facteurs: volume et salinité de l’eau interstitielle, abondance de minéraux conducteurs, température de la roche et présence de magma. Ce travail de thèse tente de contraindre l'interprétation des structures de résistivité électrique autour des volcans actifs, afin de développer des outils innovants pour l'exploration des ressources géothermiques. La contribution des minéraux conducteurs est au cœur de la thèse: conducteurs ioniques solides (minéraux argileux, en particulier la smectite) ou semi-conducteurs électroniques (pyrite, oxydes de fer), mais l’influence de la porosité, de la salinité, de la température et de la présence de magma est aussi étudiée. La thèse utilise le volcan Krafla comme terrain d’étude pour affiner les interprétations des structures de résistivité électriques, du fait de la disponibilité de carottes, de données, de bibliographie et d’infrastructure. La smectite et la pyrite sont formées par altération hydrothermale des roches volcaniques et témoignent ainsi des convections hydrothermales. Les oxydes de fer en revanche sont plutôt formés lors de la cristallisation du magma et sont dissous lors des circulations hydrothermales. La contribution de la smectite à la conductivité électrique de roches volcaniques, saturées en eau à différentes salinités, est d'abord étudiée en laboratoire (à température ambiante) par spectroscopie d’impédance électrique « résistivité complexe ». Des variations non linéaires de la conductivité électrique à 1 kHz avec la salinité sont observées et discutées. La conduction interfoliaire est suggérée comme un mécanisme important par lequel la smectite conduit le courant électrique. L'influence de la pyrite et des oxydes de fer sur les effets de polarisation provoquée est ensuite analysée en utilisant l'angle de phase de l'impédance, qui dépend de la fréquence. Un angle de phase maximal supérieur à 20 mrad est attribué à la pyrite si la roche est conductrice et aux oxydes de fer si la roche est résistive. L'angle de phase maximal augmente d'environ 22 mrad pour chaque pourcent de pyrite ou d'oxyde de fer. Ces résultats de laboratoire en domaine fréquentiel sont appliqués à l’interprétation de tomographies de résistivité complexe sur le terrain en domaine temporel. Smectite, pyrite et oxydes de fer ont pu être identifiés jusqu'à 200 m de profondeur. La température in-situ, plus élevées qu’en laboratoire, semble augmenter la conductivité de la smectite. De manière générale, la tomographie de résistivité complexe est recommandée comme méthode complémentaire aux sondages électromagnétiques pour l'exploration géothermique. / Electromagnetic soundings are widely used to image the underground structure of volcanoes and look for hightemperature geothermal resources. The electrical resistivity of volcanoes is affected by several characteristics of rocks: volume and salinity of pore fluid, abundance of conductive minerals, rock temperature and presence of magma. This thesis aims at improving the interpretation of electrical resistivity structures around active volcanoes, in order to develop innovative tools for the assessment of geothermal resources. I focus on conductive minerals, which can either be solid ionic conductors (clay minerals, in particular smectite) or electronic semi-conductors (pyrite and iron-oxides), but I also investigate the effects of porosity, salinity, temperature and presence of magma. I use Krafla volcano as a laboratory area, where extensive literature, borehole data, core samples, surface soundings and infrastructures are available. Smectite and pyrite are formed upon hydrothermal alteration of volcanic rocks and thus witness hydrothermal convection. On the other hand, iron-oxides are mostly formed during the primary crystallization of magma and dissolved by hydrothermal fluids. The contribution of smectite to the electrical conductivity of volcanic rocks saturated with pore water at different salinity is first investigated in the laboratory (room temperature) by electrical impedance spectroscopy “complex resistivity”. Non-linear variations of the conductivity at 1 kHz with salinity are observed and discussed. Interfoliar conduction is suggested as an important mechanism by which smectite conducts electrical current. The influence of pyrite and iron-oxides on induced polarization effects is then analyzed, using the frequency-dependent phase-angle of the impedance. A maximum phase-angle higher than 20 mrad is attributed to pyrite if the rock is conductive and to ironoxides if the rock is resistive. The maximum phase-angle increases by about 22 mrad for each additional per cent of pyrite or iron-oxide. These laboratory frequency-domain findings are partly upscaled to interpret field time-domain complex resistivity tomography at Krafla: smectite, pyrite and iron-oxides can be identified down to 200 m. The in-situ temperature, higher than in laboratory conditions, appears to significantly increase the conductivity associated to smectite. In general, time-domain complex resistivity measurements are recommended as a complementary method to electromagnetic soundings for geothermal exploration.

The use of geophysical testing methods has considerable potential to be a cost effective and accurate technique to assess near-surface soil conditions. Multi channel analysis of surface waves (MASW) test is a geophysical non-intrusive test that uses the dispersive characteristic of Rayleigh waves to estimate low strain shear modulus and damping coefficient of near-surface soil. Also, this technique is used to detect underground voids. Recently, MASW technique has gained more attention, partly because of its ease of use and partly because of the significant improvements in data acquisition systems. The theories of MASW test consider the effect of horizontal soil layering, though the effect of lateral inhomogeneities (i. e. cavities and voids), inclined layering and inverse layering (i. e. a layered system in which the top layers are stiffer than the bottom ones) are not addressed properly in these theories. <br /><br /> The objective of this dissertation is to investigate the effect of lateral inhomogeneities on the propagation of Rayleigh waves in an elastic half-space excited by a transient loading. The results can be applied to locate underground cavities using MASW test and to improve the MASW analysis techniques. In lieu of theoretical solutions, two and three dimensional numerical models are constructed to simulate the MASW test. To assure the quality of the obtained data, numerical models are calibrated with Lamb solution. Voids with different sizes and embedment depths are inserted in the medium. Responses along the surface as well as inside the medium are recorded and analyzed in time, frequency, spatial and frequency-wave number domains. Different material types and sources are used to generalize the results. Afterwards, the combined effect of void and layered systems on the surface responses are studied. To verify the results, experimental field and laboratory data are presented and the trends are compared to the numerical results. <br /><br /> It is found that the void starts to vibrate in response to the Rayleigh wave excitation. Due to the vibration of the void energy partitioning occurs. Part of the incident energy is reflected in the form of Rayleigh wave. Another part is converted to body waves, and spread into the medium. The transferred part of the energy is attenuated and has smaller amplitudes. Finally, a part of energy is trapped in the void region and bounces back and forth between the void boundaries, until it damps. The trapped energy is associated to higher modes of Rayleigh waves and excited Lamb waves. The effect of trapped energy is seen as a region in the vicinity of the void with concentrated energy, in frequency domain. The extents of this region depends on the void size, and the frequency content of the incident energy. Thus, in some cases it is possible to correspond the size of the model to the extents of the region with energy concentration. <br /><br /> A new technique is proposed to determine the location of a void, and estimate its embedment depth. The technique is called Attenuation Analysis of Rayleigh Waves (AARW), and is based on the observed damping effect of the void on the surface responses. For verification, the results are compared to experimental field and laboratory data. The observations are in good agreement with the observed numerical results. Further, the AARW technique showed to be a promising tool for void detection.

The use of geophysical testing methods has considerable potential to be a cost effective and accurate technique to assess near-surface soil conditions. Multi channel analysis of surface waves (MASW) test is a geophysical non-intrusive test that uses the dispersive characteristic of Rayleigh waves to estimate low strain shear modulus and damping coefficient of near-surface soil. Also, this technique is used to detect underground voids. Recently, MASW technique has gained more attention, partly because of its ease of use and partly because of the significant improvements in data acquisition systems. The theories of MASW test consider the effect of horizontal soil layering, though the effect of lateral inhomogeneities (i. e. cavities and voids), inclined layering and inverse layering (i. e. a layered system in which the top layers are stiffer than the bottom ones) are not addressed properly in these theories. <br /><br /> The objective of this dissertation is to investigate the effect of lateral inhomogeneities on the propagation of Rayleigh waves in an elastic half-space excited by a transient loading. The results can be applied to locate underground cavities using MASW test and to improve the MASW analysis techniques. In lieu of theoretical solutions, two and three dimensional numerical models are constructed to simulate the MASW test. To assure the quality of the obtained data, numerical models are calibrated with Lamb solution. Voids with different sizes and embedment depths are inserted in the medium. Responses along the surface as well as inside the medium are recorded and analyzed in time, frequency, spatial and frequency-wave number domains. Different material types and sources are used to generalize the results. Afterwards, the combined effect of void and layered systems on the surface responses are studied. To verify the results, experimental field and laboratory data are presented and the trends are compared to the numerical results. <br /><br /> It is found that the void starts to vibrate in response to the Rayleigh wave excitation. Due to the vibration of the void energy partitioning occurs. Part of the incident energy is reflected in the form of Rayleigh wave. Another part is converted to body waves, and spread into the medium. The transferred part of the energy is attenuated and has smaller amplitudes. Finally, a part of energy is trapped in the void region and bounces back and forth between the void boundaries, until it damps. The trapped energy is associated to higher modes of Rayleigh waves and excited Lamb waves. The effect of trapped energy is seen as a region in the vicinity of the void with concentrated energy, in frequency domain. The extents of this region depends on the void size, and the frequency content of the incident energy. Thus, in some cases it is possible to correspond the size of the model to the extents of the region with energy concentration. <br /><br /> A new technique is proposed to determine the location of a void, and estimate its embedment depth. The technique is called Attenuation Analysis of Rayleigh Waves (AARW), and is based on the observed damping effect of the void on the surface responses. For verification, the results are compared to experimental field and laboratory data. The observations are in good agreement with the observed numerical results. Further, the AARW technique showed to be a promising tool for void detection.

The purpose of this study is to assess the engineering geological and geotechnical characteristics and to perform seismic hazard studies of the Upper Pliocene to Quaternary (Plio-Quaternary) deposits located towards the north of Ankara through surface wave testing methods. Based on a general engineering geological and seismic site characterization studies, site classification systems are assigned in seismic hazard assessments. The objective of the research is to determine the regional and local seismic soil conditions (i.e., shear wave velocities, soil predominant periods and soil amplification factors) and to characterize the soil profile of the sites in this region by the help of surface geophysical methods. These studies have been supported by engineering geological and geotechnical field studies carried out prior to and during this study. By integrating these studies, local soil conditions and dynamic soil characteristics for the study area have been assessed by detailed soil characterization in the region. As a result, seismic hazard assessments have been performed for &Ccedil / ubuk and its close vicinity with the aid of Geographical Information Systems (GIS) through establishing seismic characterization and local soil conditions of the area.

<p>Determination of the subsurface geology is very important for the rock quarry industry. This is primarily done by drilling and mapping of the surface geology. However in Sweden the bedrock is often completely covered by Quaternary sediments making the prediction of subsurface geology quite difficult. Incorrect prediction of the rock-mass quality can lead to economic problems for the quarry. By performing geophysical measurements a more complete understanding of the subsurface geology can be determined. This study shows that by doing 2D-parallel data sampling a 3D inversion of the dataset is possible, which greatly enhances the visualization of the subsurface. Furthermore the electrical resistivity technique together with the induced polarization method proved to be very efficient in detecting fracture frequency, identification of major fracture zones, and variations in rock-mass quality all of which can affect the aggregate quality. With this technique not only the rock-mass quality is determined but also the thickness of the overburden. Implementation of geophysics can be a valuable tool for the quarry industry, resulting in substantial economic benefits.  </p>

The sinking of dense shelf waters down the continental slope (or “cascading”) contributes to oceanic water mass formation and carbon cycling. Cascading is therefore of significant importance for the global overturning circulation and thus climate. The occurrence of cascades is highly intermittent in space and time and observations of the process itself (rather than its outcomes) are scarce. Global climate models do not typically resolve cascading owing to numerical challenges concerning turbulence, mixing and faithful representation of bottom boundary layer dynamics. This work was motivated by the need to improve the representation of cascading in numerical ocean circulation models. Typical 3-D hydrostatic ocean circulation models are employed in a series of numerical experiments to investigate the process of dense water cascading in both idealised and realistic model setups. Cascading on steep bottom topography is modelled using POLCOMS, a 3-D ocean circulation model using a terrain-following s-coordinate system. The model setup is based on a laboratory experiment of a continuous dense water flow from a central source on a conical slope in a rotating tank. The descent of the dense flow as characterised by the length of the plume as a function of time is studied for a range of parameters, such as density difference, speed of rotation, flow rate and (in the model) diffusivity and viscosity. Very good agreement between the model and the laboratory results is shown in dimensional and non-dimensional variables. It is confirmed that a hydrostatic model is capable of reproducing the essential physics of cascading on a very steep slope if the model correctly resolves velocity veering in the bottom boundary layer. Experiments changing the height of the bottom Ekman layer (by changing viscosity) and modifying the plume from a 2-layer system to a stratified regime (by enhancing diapycnal diffusion) confirm previous theories, demonstrate their limitations and offer new insights into the dynamics of cascading outside of the controlled laboratory conditions. In further numerical experiments, the idealised geometry of the conical slope is retained but up-scaled to oceanic dimensions. The NEMO-SHELF model is used to study the fate of a dense water plume of similar properties to the overflow of brine-enriched shelf waters from the Storfjorden in Svalbard. The overflow plume, resulting from sea ice formation in the Storfjorden polynya, cascades into the ambient stratification resembling the predominant water masses of Fram Strait. At intermediate depths between 200-500m the plume encounters a layer of warm, saline AtlanticWater. In some years the plume ‘pierces’ the Atlantic Layer and sinks into the deep Fram Strait while in other years it remains ‘arrested’ at Atlantic Layer depths. It has been unclear what parameters control whether the plume pierces the Atlantic Layer or not. In a series of experiments we vary the salinity ‘S’ and the flow rate ‘Q’ of the simulated Storfjorden overflow to investigate both strong and weak cascading conditions. Results show that the cascading regime (piercing, arrested or ‘shaving’ - an intermediate case) can be predicted from the initial values of S and Q. In those model experiments where the initial density of the overflow water is considerably greater than of the deepest ambient water mass we find that a cascade with high initial S does not necessarily reach the bottom if Q is low. Conversely, cascades with an initial density just slightly higher than the deepest ambient layer may flow to the bottom if the flow rate Q is high. A functional relationship between S/Q and the final depth level of plume waters is explained by the flux of potential energy (arising from the introduction of dense water at shallow depth) which, in our idealised setting, represents the only energy source for downslope descent and mixing. Lastly, the influence of tides on the propagation of a dense water plume is investigated using a regional NEMO-SHELF model with realistic bathymetry, atmospheric forcing, open boundary conditions and tides. The model has 3 km horizontal resolution and 50 vertical levels in the sh-coordinate system which is specially designed to resolve bottom boundary layer processes. Tidal effects are isolated by comparing results from model runs with and without tides. A hotspot of tidally-induced horizontal diffusion leading to the lateral dispersion of the plume is identified at the southernmost headland of Spitsbergen which is in close proximity to the plume path. As a result the lighter fractions in the diluted upper layer of the plume are drawn into the shallow coastal current that carries Storfjorden water onto the Western Svalbard Shelf, while the dense bottom layer continues to sink down the slope. This bifurcation of the plume into a diluted shelf branch and a dense downslope branch is enhanced by tidally-induced shear dispersion at the headland. Tidal effects at the headland are shown to cause a net reduction in the downslope flux of Storfjorden water into deep Fram Strait. This finding contrasts previous results from observations of a dense plume on a different shelf without abrupt topography. The dispersive mechanism which is induced by the tides is identified as a mechanism by which tides may cause a relative reduction in downslope transport, thus adding to existing understanding of tidal effects on dense water overflows.

No description available.