Over a decade ago, it was discovered that the earth’s natural seismic wavefields, propagating as seismic noise, can be processed using correlation methods to produce surface waves, similar to those generated by earthquakes. This discovery represents a paradigm shift in seismology and has led to several tomographic studies of earth structures, at different scales and resolutions, in previously difficult-to-study areas around the world. This PhD dissertation presents research results on multi-scale and multi-purpose applications of ambient seismic noise wavefields under three topics: (1) Imaging of sedimentary basins and sub-basin structures in eastern and western Canada using ambient seismic noise, (2) Combining measurements from ambient seismic noise with earthquake datasets for imaging crustal and mantle structures, and (3) Temporal variation in cultural seismic noise and noise correlation functions (NCFs) during the COVID-19 lockdown in Canada. The first topic involved imaging the sedimentary basins in eastern and western Canada using shear wave velocities derived from ambient noise group velocities. The results show that the basins are characterized by varying depths, with maximums along the studied cross-sections in excess of 10 km, in eastern and western Canada. Characteristics of accreted terranes in eastern and western Canada are also revealed in the results. A seismically distinct basement is imaged in eastern Canada and is interpreted to be a vestige of the western African crust trapped beneath eastern Canada at the opening of the Atlantic Ocean. In western Canada, the 3D variation of the Moho and sedimentary basin depths is imaged. The thickest sediments in eastern Canada are found beneath the Queen Charlotte, Williston and the Alberta Deep basins, while the Moho is the deepest beneath the Williston basin and parts of Alberta basin and northern British Columbia. For the second topic, I worked on improving the seismological methodology to construct broadband (period from 2 to 220 s) dispersion curves by combining the dispersion measurements derived from ambient seismic noise with those from earthquakes. The broadband dispersion curves allow for imaging earth structures spanning the shallow crust to the upper mantle. For the third topic, I used ambient seismic data from the earlier stages of the COVID-19 pandemic to study the temporal variation of seismic power spectra and the potential impacts of COVID-19 lockdown on ambient NCFs in four cities in eastern and western Canada. The results show mean seismic power drops of 24% and 17% during the lockdown in eastern Canada, near Montreal and Ottawa respectively and reductions of 27% and 17% near Victoria and Sidney respectively. NCF signal quality within the secondary microseism band reached maximum before the lockdown, minimum during lockdown and at intermediate levels during the gradual reopening phase for the western Canada station pair. / Graduate

In this study, ambient noise wavefield was used for the first time to image spatial and temporal upper crustal seismic structures in NW-Bohemia/Vogtland region. The data come from 111 stations and were collected from continuous recordings of the permanent station networks of Germany and Czech Academy of Sciences as well as temporary stations of the BOHEMA and PASSEQ experiments. Rayleigh and Love waves travelling between each station-pair are extracted by cross-correlating long time series of ambient noise data recorded at the stations. Group velocity dispersion curves are obtained by time-frequency analysis of cross-correlation functions between 0.1 and 1 Hz, and are tomographically inverted to provide 2-D group velocity maps. At shorter periods Rayleigh wave group velocity maps are in good agreement with surface geology where low velocity anomalies appear along Mariánské Lázně Fault and Eger rift. A low velocity zone is observed at the northern edge of Mariánské Lázně Fault which shifts slightly to the south with increasing period and correlates well with the main focal zone of the earthquake swarms at 5 s period. We invert the 2-D group velocity maps into a 3-D shear wave velocity model. In this step Love waves were excluded from further analysis because of their high level of misfit to modelled dispersion curves. Horizontal and vertical sections through the model reveal a clear low velocity zone above the Nový Kostel seismic focal zone which narrows towards the top of the seismic activity and ends above the shallowest hypocenters at 7 km depth. We investigate temporal variation of seismic velocity within and around the Nový Kostel associated with 2008 and 2011 earthquake swarms by employing Passive Image Interferometry method using 7 continuous seismograms recorded by the WEBNET network. The results reveals stable seismic velocities without a clear post seismic velocity change during earthquake swarms in the Nový Kostel area.

Seismic records contain information that allows geoscientists to make inferences about the structure and properties of the Earth’s interior. Traditionally, seismic imaging and tomography methods require wavefields to be generated and recorded by identifiable sources and receivers, and use these directly-recorded signals to create models of the Earth’s subsurface. However, in recent years the method of seismic interferometry has revolutionised earthquake seismology by allowing unrecorded signals between pairs of receivers, pairs of sources, and source-receiver pairs to be constructed as Green’s functions using either cross-correlation, convolution or deconvolution of wavefields. In all of these formulations, seismic energy is recorded and emitted by surrounding boundaries of receivers and sources, which need not be active and impulsive but may even constitute continuous, naturally-occurring seismic ambient noise. In the first part of this thesis, I provide a comprehensive overview of seismic interferometry, its background theory, and examples of its application. I then test the theory and evaluate the effects of approximations that are commonly made when the interferometric formulae are applied to real datasets. Since errors resulting from some approximations can be subtle, these tests must be performed using almost error-free synthetic data produced with an exact waveform modelling method. To make such tests challenging the method and associated code must be applicable to multiply-scattering media. I developed such a modelling code specifically for interferometric tests and applications. Since virtually no errors are introduced into the results from modelling, any difference between the true and interferometric waveforms can safely be attributed to specific origins in interferometric theory. I show that this is not possible when using other, previously available methods: for example, the errors introduced into waveforms synthesised by finite-difference methods due to the modelling method itself, are larger than the errors incurred due to some (still significant) interferometric approximations; hence that modelling method can not be used to test these commonly-applied approximations. I then discuss the ability of interferometry to redatum seismic energy in both space and time, allowing virtual seismograms to be constructed at new locations where receivers may not have been present at the time of occurrence of the associated seismic source. I present the first successful application of this method to real datasets at multiple length scales. Although the results are restricted to limited bandwidths, this study demonstrates that the technique is a powerful tool in seismologists’ arsenal, paving the way for a new type of ‘retrospective’ seismology where sensors may be installed at any desired location at any time, and recordings of seismic events occurring at any other time can be constructed retrospectively – even long after their energy has dissipated. Within crustal seismology, a very common application of seismic interferometry is ambient-noise tomography (ANT). ANT is an Earth imaging method which makes use of inter-station Green’s functions constructed from cross-correlation of seismic ambient noise records. It is particularly useful in seismically quiescent areas where traditional tomography methods that rely on local earthquake sources would fail to produce interpretable results due to the lack of available data. Once constructed, interferometric Green’s functions can be analysed using standard waveform analysis techniques, and inverted for subsurface structure using more or less traditional imaging methods. In the second part of this thesis, I discuss the development and implementation of a fully non-linear inversion method which I use to perform Love-wave ANT across the British Isles. Full non-linearity is achieved by allowing both raypaths and model parametrisation to vary freely during inversion in Bayesian, Markov chain Monte Carlo tomography, the first time that this has been attempted. Since the inversion produces not only one, but a large ensemble of models, all of which fit the data to within the noise level, statistical moments of different order such as the mean or average model, or the standard deviation of seismic velocity structures across the ensemble, may be calculated: while the ensemble average map provides a smooth representation of the velocity field, a measure of model uncertainty can be obtained from the standard deviation map. In a number of real-data and synthetic examples, I show that the combination of variable raypaths and model parametrisation is key to the emergence of previously-unobserved, loop-like uncertainty topologies in the standard deviation maps. These uncertainty loops surround low- or high-velocity anomalies. They indicate that, while the velocity of each anomaly may be fairly well reconstructed, its exact location and size tend to remain uncertain; loops parametrise this location uncertainty, and hence constitute a fully non-linearised, Bayesian measure of spatial resolution. The uncertainty in anomaly location is shown to be due mainly to the location of the raypaths that were used to constrain the anomaly also only being known approximately. The emergence of loops is therefore related to the variation in raypaths with velocity structure, and hence to 2nd and higher order wave-physics. Thus, loops can only be observed using non-linear inversion methods such as the one described herein, explaining why these topologies have never been observed previously. I then present the results of fully non-linearised Love-wave group-velocity tomography of the British Isles in different frequency bands. At all of the analysed periods, the group-velocity maps show a good correlation with known geology of the region, and also robustly detect novel features. The shear-velocity structure with depth across the Irish Sea sedimentary basin is then investigated by inverting the Love-wave group-velocity maps, again fully non-linearly using Markov chain Monte Carlo inversion, showing an approximate depth to basement of 5 km. Finally, I discuss the advantages and current limitations of the fully non-linear tomography method implemented in this project, and provide guidelines and suggestions for its improvement.

Este trabalho consiste no estudo da estrutura sedimentar da Bacia do Paraná da região de Bebedouro utilizando função do receptor, ondas de superfície e sísmica de refração. Funções do Receptor (FR) mostram a resposta da estrutura geológica abaixo de uma estação sismográfica. Para diminuir a não unicidade na inversão do traço de uma FR, foram usadas curvas de dispersão de ondas de superfície como vínculos adicionais. Ondas de superfície foram também obtidas com a correlação cruzada de ruído sísmico ambiental. Esse método passivo permite a obtenção de dispersão em períodos intermediários entre os dados de telessismos e os dados de sísmica rasa auxiliando na determinação estruturas sedimentares. Foi utilizada a inversão conjunta de FR de alta e baixa frequência, dispersão de fase e grupo de períodos menores que 2 s e dispersão de ondas Love e Rayleigh continentais (períodos entre 10 e 100 s). Na inversão conjunta, foram usados vários modelos iniciais diferentes para garantir melhor cobertura do espaço de soluções possíveis. Os resultados mostram uma camada de basalto com uma espessura entre 200 e 400 m, sob uma camada superficial de arenito de 50 a 100m. O embasamento da bacia está torno de 2.5 3.0 km. A espessura e razão Vp/Vs crustal da região foi estimada em 40.0 (1.0) km e 1.78 (0.02). Foi possível identificar zonas de baixa velocidade dentro do pacote de basalto, possivelmente relacionadas a zonas de fraturas ou camadas de basalto alterado. Estas camadas de baixa velocidade estão na parte sul, próximas à área de maior sismicidade induzida e também relativamente mais perto dos poços profundos de maior vazão. Isto corrobora com a hipótese de que a sismicidade local está associada à exploração de poços e estrutura de falhas pré -existentes. / A study of the seismic structure of the Paraná Basin in Bebedouro, SP, was carried out with Receiver Functions (RF). Both low frequency (with ~0.5 Hz low pass filter) and high-frequency (~10Hz) RFs were jointly inverted with surface wave dispersion curves. Surface-wave data included: long-period group velocities of Rayleigh and Love waves from continental-scale tomography in the period range 10-100 s (useful to control crustal scale structure and Moho depth), intermediate period group velocities near 1s period obtained from cross-correlation of ambient noise (useful to control sedimentary structure), phase velocities near 0.2-1.0 s obtained from local earthquakes, and phase velocities of Rayleigh waves from shallow seismic refraction data near ~0.1 s period (useful to control shallow layers). Joint inversion was obtained with several different initial models to better cover all possible solutions. The results indicate that the basalt layer is 200 to 400m thick, beneath sandstones of about 50 100m thick. The Basin basement was detected between 2.5 3.0 km depth, consistent with expected values from regional bore-hole data in the Paraná Basin. Beneath several stations a low-velocity zone was detected in the middle of the basalt pack which is interpreted as a zone of fractured or altered basalt. This anomalous low-velocity layer was detected near the most seismically active zone in the Andes district. It is also relatively close to the deep wells with larger outflow. The existence of this low-velocity zone in the middle of the basalt layer is consistent with the model proposed for the water-well induced seismicity.

Este trabalho consiste no estudo da estrutura sedimentar da Bacia do Paraná da região de Bebedouro utilizando função do receptor, ondas de superfície e sísmica de refração. Funções do Receptor (FR) mostram a resposta da estrutura geológica abaixo de uma estação sismográfica. Para diminuir a não unicidade na inversão do traço de uma FR, foram usadas curvas de dispersão de ondas de superfície como vínculos adicionais. Ondas de superfície foram também obtidas com a correlação cruzada de ruído sísmico ambiental. Esse método passivo permite a obtenção de dispersão em períodos intermediários entre os dados de telessismos e os dados de sísmica rasa auxiliando na determinação estruturas sedimentares. Foi utilizada a inversão conjunta de FR de alta e baixa frequência, dispersão de fase e grupo de períodos menores que 2 s e dispersão de ondas Love e Rayleigh continentais (períodos entre 10 e 100 s). Na inversão conjunta, foram usados vários modelos iniciais diferentes para garantir melhor cobertura do espaço de soluções possíveis. Os resultados mostram uma camada de basalto com uma espessura entre 200 e 400 m, sob uma camada superficial de arenito de 50 a 100m. O embasamento da bacia está torno de 2.5 3.0 km. A espessura e razão Vp/Vs crustal da região foi estimada em 40.0 (1.0) km e 1.78 (0.02). Foi possível identificar zonas de baixa velocidade dentro do pacote de basalto, possivelmente relacionadas a zonas de fraturas ou camadas de basalto alterado. Estas camadas de baixa velocidade estão na parte sul, próximas à área de maior sismicidade induzida e também relativamente mais perto dos poços profundos de maior vazão. Isto corrobora com a hipótese de que a sismicidade local está associada à exploração de poços e estrutura de falhas pré -existentes. / A study of the seismic structure of the Paraná Basin in Bebedouro, SP, was carried out with Receiver Functions (RF). Both low frequency (with ~0.5 Hz low pass filter) and high-frequency (~10Hz) RFs were jointly inverted with surface wave dispersion curves. Surface-wave data included: long-period group velocities of Rayleigh and Love waves from continental-scale tomography in the period range 10-100 s (useful to control crustal scale structure and Moho depth), intermediate period group velocities near 1s period obtained from cross-correlation of ambient noise (useful to control sedimentary structure), phase velocities near 0.2-1.0 s obtained from local earthquakes, and phase velocities of Rayleigh waves from shallow seismic refraction data near ~0.1 s period (useful to control shallow layers). Joint inversion was obtained with several different initial models to better cover all possible solutions. The results indicate that the basalt layer is 200 to 400m thick, beneath sandstones of about 50 100m thick. The Basin basement was detected between 2.5 3.0 km depth, consistent with expected values from regional bore-hole data in the Paraná Basin. Beneath several stations a low-velocity zone was detected in the middle of the basalt pack which is interpreted as a zone of fractured or altered basalt. This anomalous low-velocity layer was detected near the most seismically active zone in the Andes district. It is also relatively close to the deep wells with larger outflow. The existence of this low-velocity zone in the middle of the basalt layer is consistent with the model proposed for the water-well induced seismicity.

Near surface measurements of seismic velocity and Q are useful in a number of situations, for example for when carrying out re-datuming and migration for depth images, or when analysing ground conditions for building. This thesis concentrates on the estimation of surface wave group and phase velocity as well as Q structure through the use of cross correlations of ambient noise recordings. Linearised tomography estimates are made for the British Isles, the Permain Basin of Texas and New Mexico, the Ekofisk Life of Field Seismic (LoFS) array and the Aquistore CO2 storage site. The results correspond well with the known geological structure and/or structure observed in velocity maps by other researchers. For the Ekofisk array a non-linearised non-linear method was also applied and the results estimated by these two methods for the Ekofisk LoFS array are compared. By non-linearised non-linear it is meant that the inversion method is fully non-linear and no linearisation has taken place in the method, this term will be used throughout this thesis for all methods which fall into this category. The tomography results from the two methods had similarities in their general structure but differences in the finer details, and so suggest that the substantial increase in time required for the non-linearised non-linear method is not justified. Linearised tomographic inversion of the Aquistore array was used to determine the potential of using ambient noise tomography for monitoring of CO2 storage sites. It was found that the repeatability of the tomography at the Aquistore site was not good enough to allow ambient noise tomography to be used for monitoring; however, it may be possible to apply the method at other sites. A Q and phase velocity inversion of the Ekofisk array is also presented, with results mostly showing excellent correlation with known geological features. It is shown that the higher frequencies are more sensitive to the effects of sea floor subsidence at the site, while lower frequencies are more sensitive to the effects of faulting. A final near surface method called ambient noise gradiometry was applied, this concentrates on estimating locations of sources of seismic energy within receiver arrays. Ambient noise gradiometry is applied to synthetic and real data for this purpose. It was found that using ambient noise gradiometry allows internal sources of energy to be identified but they produce a bias in the phase velocity tomography result. Two methods of reducing this bias are presented, both of which also provide an estimate of the source term for different sections of time of the recording.

Surface-generalized ambient noise in a shallow ocean waveguide with a sediment layer possessing a specific class of density and sound speed distributions capable of describing a realistic seabed environment is considered in this analysis. This class of non-uniform sediment layer has the density and sound speed distributions varying with respect to depth as a genearlized-exponential and an inverse-square function, respectively. The study invokes a formulation developed by Kuperman and Ingenito for surface noise generation, in conjunction with the analytical solutions for the Helmholtz equation corresponding to the sediment layer, to arrive at an analytical expression convenient for numerical implementation. The intensity and spatial correlation of the noise sound field are analyzed with respect to the variations of the system parameters, including frequency, sediment layer thickness, sound speed gradient, with emphasis on the effects of sediment properties on the ambient noise field. The results have demonstrated that the intensity of the noise field is relatively sensitive to the variations of the paramters, while the spatial correlation is not, suggesting that the energy distribution, rather than the spatial structure, of the noise field is susceptible to the environmental variation.

A rapid growth of installed wind power capacity is expectedin the next few years. However, the siting of wind turbines ona large scale raises concerns about their environmental impact,notably with respect to noise. To this end, variable speed windturbines offer a promising solution for applications in denselypopulated areas like the European countries, as this designwould enable an efficient utilisation of the masking effect dueto ambient noise. In rural and recreational areas where windturbines are sited, the ambient noise originates from theaction of wind on the vegetation and about the listener's ear(pseudo-noise). It shows a wind speed dependence similar tothat of the noise from a variable speed wind turbine and cantherefore mask the latter for a wide range of conditions.However, a problem inherent to the design of these machines istheir proclivity to pure tone generation, because of theenhanced difficulty of avoiding structural resonances in themechanical parts. Pure tones are deemed highly annoying and areseverely regulated by most noise policies. In relation to thisproblem, the vibration transmission of structure-borne sound tothe tower of the turbine is investigated, in particular whenthe tower is stiffened at its upper end. Furthermore, sincenoise annoyance due to wind turbine is mostly a masking issue,the wind-related sources of ambient noise are studied and theirmasking potentials assessed. With this aim, prediction modelsfor wind-induced vegetation noise and pseudo-noise have beendeveloped. Finally, closely related to the effect of masking,is the difficulty, regularly encountered by local authoritiesand wind farm developers, to measure noise immission from windturbines. A new measurement technique has thus been developedin the course of this work. Through improving thesignal-to-noise ratio between wind turbine noise and ambientnoise, the new technique yields more accurate measurementresults. Keywords: Masking, vibration transmission, diffraction,ambient noise, pseudo-noise, cylindrical shell, perturbationmethods, structural mobility, acoustic outdoor measurement. / QC 20100616

Le développement de l’énergie géothermique a conduit à l’exploitation de ressources établies dans des contextes géologiques et géodynamiques très variés. L’exploration géophysique de ces réservoirs complexes nécessite l’utilisation de plusieurs méthodes d’imagerie complémentaire. Ce travail de thèse porte sur l’exploration d’une ressource géothermique située en contexte de socle fracturé dans le Massif Central français par magnétotellurique, tomographie de bruit ambiant et gravimétrie.La magnétotellurique est une méthode d’imagerie 3D résolvante qui est sensible à la présence d’eau et aux argiles d’altération hydrothermale mais limitée par sa couverture spatiale. La tomographie de bruit sismique présente une bonne résolution verticale mais ne résout pas les variations horizontales de vitesse. Cette méthode est sensible aux variations des propriétés mécaniques des roches et donc aux milieux fracturés. Enfin la gravimétrie apporte une contrainte sur les variations lithologiques et possède une bonne résolution latérale mais une faible résolution verticale.Nous présentons une méthode d’inversion conjointe des données sismiques et gravimétriques sous contrainte d’un modèle de résistivité obtenu par inversion magnétotellurique indépendante. L’inversion conjointe nécessite de définir des couplages entre modèles. Par absence de connaissance a priori de relations pétrophysiques, nous avons couplé les modèles de densité, de résistivité et de vitesse avec une loi qui contraint les paramètres à être corrélés en moyenne. Cette stratégie vise à faire ressortir des relations caractéristiques des objets géologiques de la ressource géothermique.Cette méthodologie d’inversion conjointe a été testée sur des modèles synthétiques. L’application aux données réelles acquises dans le Massif Central a permis de définir une zone en profondeur de forte corrélation interprétée comme la transition ductile fragile. La partie intermédiaire des modèles, plus homogène, permet de distinguer différentes unités géologiques séparées par une zone de faille. Enfin la partie superficielle se distingue par une forte hétérogénéité des paramètres résultants probablement de processus d’altération de surface. / The development of geothermal energy has led to the exploitation of resources established in varied geological and geodynamic contexts. Geophysical exploration of these complex reservoirs requires the use of several complementary imaging methods. This PhD thesis focuses on the exploration of a geothermal resource located within the fractured basement in the French Massif Central using magnetotelluric, ambient noise tomography and gravimetry. Magnetotelluric is a 3D imaging method with a good resolution power that is sensitive to the presence of water and hydrothermal weathering clays but is limited by its spatial coverage. Seismic noise tomography has a good vertical resolution but does not resolve well horizontal velocity variations. This method is sensitive to variations of the mechanical properties of rocks and thus to fractured media. Finally gravimetry brings constraint on the lithological variations and has a good lateral resolution but lacks vertical resolution.We present a method of joint inversion of seismic and gravimetric data under the constraint of a resistivity model obtained by independent magnetotelluric inversion. Joint inversion requires defining model couplings. By lack of prior knowledge of petrophysical relationships, we have coupled the density, resistivity and velocity models with a law that constraints the parameters to be correlated on average.This strategy aims to bring out the characteristic relationships of the geological objects of the geothermal resource. This joint inversion methodology has been tested on synthetic models. The application to the real data acquired in the Massif Central has made it possible to define a deep zone of high correlation interpreted as the fragile ductile transition. The intermediate part of the models, more homogeneous, allows to distinguish different geological units separated by a fault zone. Finally the superficial part is distinguished by strong heterogeneity of the parameters resulting probably from surface alteration process.

We have performed 3D ambient noise tomography of the Bohemian Massif. We invert adopted inter-station dispersion curves of both Love and Rayleigh waves in periods 4-20 s, which were extracted from ambient noise cross-correlations, using a two-step approach. In the first step, the inter-station dispersion curves are localized for each period into the so-called dispersion maps. To account for finite-frequency effects, gradient method employing Fréchet kernels is used. Assuming membrane wave approximation of the surface wave propagation at each period, the kernels were calculated using the adjoint method. To reduce the effect of data noise, the kernels were regularized by Gaussian smoothing. The proper level of regularization is assessed on synthetic tests. In the second step, the phase-velocity dispersion maps are inverted into a 3D S-wave velocity model using the Bayesian approach. The posterior probability density function describing the solution is sampled by more than one million models obtained by Monte-Carlo approach (parallel tempering). The calculated variance of the model shows that the well resolved part corresponds to the upper crust (i.e., upper 20 km). The mean velocity model contains mainly large scale structures that show good correlation with the main geologic domains of the Bohemian...