<p>The resolution of near-surface seismic reflection data is often limited by attenuation and scattering in the shallow subsurface which reduces the high frequencies in the data. Compensating for attenuation and scattering, as well as removing the propagating source wavelet in a time-variant manner can be used to improve the resolution. Here we investigate continuous non-stationary iterative time-domain deconvolution (CNS-ITD), where the seismic wavelet is allowed to vary along the seismic trace. The propagating seismic wavelet is then a combination of the source wavelet and the effects of attenuation and scattering effects, and can be estimated in a data-driven manner by performing a Gabor decomposition of the data. For each Gabor window, the autocorrelation is estimated and windowed about zero lag to estimate the propagating wavelet. Using the matrix-vector equations, the estimated propagating wavelets are assigned to the related columns of a seismic wavelet matrix, and these are then interpolated to the time location where the maximum of the envelope of the trace occurs within the iterative time-domain deconvolution. Advantages of using this data-driven, time-varying approach include not requiring prior knowledge of the attenuation and scattering structure and allowing for the sparse estimation of the reflectivity within the iterative deconvolution. We first apply CNS-ITD to synthetic data with a time-varying attenuation, where the method successfully identified the reflectors and increased the resolution of the data. We then applied CNS-ITD to two observed shallow seismic reflection datasets where improved resolution was obtained. </p>

Using a specially constructed Gas Hydrate Resonant Column (GHRC), the University of Southampton explored different methods of hydrate synthesis and measured the properties of the resulting sediments, such as shear wave velocity (Vs), compressional wave velocity (Vp) and their respective attenuation measurements (Qs -1 and Qp -1). Two approaches were considered. The first utilises an excess gas technique, where known water volume in the pore space dictates the quantity of hydrate. The second approach uses a known quantity of methane gas within the water saturated pore space to constrain the volume of hydrate. Results from the two techniques show that hydrates formed in excess gas environments cause stiffening of the sediment structure at low concentrations (3%), whereas, even at high concentrations of hydrate (40%) in excess water environments, only moderate increase in stiffness was observed. Additionally, attenuation results show a peak in damping at approximately 5% hydrate in excess gas tests, whereas in excess water tests, damping continues to increase with increasing hydrate content in the pore space. By considering the results from the two approaches, it becomes apparent that formation method has an influence on the properties of the hydrate bearing sand, and must therefore influence the morphology of the hydrate in the pore space.

Conventional AVO analysis has been used for the past 4 decades to aid in locating oil and gas reservoirs for extraction. It is, however, not possible to use it to acquire information on the porosity of the rock, the fluid saturation or other important petrophysical parameters. In this thesis, I study the effects of attenuation on seismic waveform signatures, due to wave induced fluid-flow. In the first part of the thesis, 2 models were used to synthetically model the attenuation caused by the wave induced fluid-flow: White's model and the double-porosity dual-permeability (DPDP) model. The parameters used for modeling were both synthetic and acquired from real well data of a known oil field off the coast of Norway. White's model was found to model relatively high attenuation (5%) for intermediately consolidated gas reservoirs while oil saturated intermediately consolidated reservoirs showed such low attenuation (0.3%) that it is easy to say that for the real-world situation it would not be detected. The DPDP model seemed to be able to better describe the attenuation and gave attenuations up to 10% for an intermediately consolidated oil reservoir, but due to lack of parameters from well data it was not sufficiently able to model the real-world situation. The synthetic data, however, show interesting characteristics and it is therefore recommended that more and detailed well parameters be acquired if the research should continue. For the second part, Svenska Petroleum Exploration AB and Det Norske Oljeselskap ASA provided stacked seismic data that were spectrally analyzed for hints of attenuation variation with frequency (using Fourier Transform and Complex Spectral Decomposition). Twelve locations, on the stacked seismic cube, were analyzed; six oil saturated; and six (assumed) water saturated. At each location, a main trace was selected along with the two nearby traces on each side of it (five in total). The Complex Spectral Decomposition method seemed unable to correctly break down the stacked section's signal, which is why Fourier Transform was used for further analysis. The frequency analysis showed a peak at ~30 Hz for both oil and water saturated reservoirs which seems like a characteristic frequency of the source, but that was unfortunately not confirmed and not enough time was available to test the assumption. The Fourier transform seems to show some difference between oil and water saturated traces, but that could well be because of lithological differences and not the pore fluid. It is therefore recommended, if research is to be continued, that 4D seismic data is used to analyze the same location with respect to time. It is also recommended that pre-stack or shot data be used as information is lost in stacked data. / Konventionell AVO-analys har använts under fyra deceniär som ett hjälpmedel för att finna olje- och gasreserver, men tekniken kan även användas för att erhålla information om bergets porositet, vätskemättnaden och andra viktiga petrofysiska parametrar. I denna avhandlingen har jag studerat hur våginducerat vätskeflöde påverkar dämpningen av den seismiska vågformssignaturen. I den första delen av avhandlingen användes två metoder för att syntetisk modellera dämpning orsakad av våginducerat vätskeflöde: "White's modell" och "double-porosity dual-permeability (DPDP) modellen". Både syntetiska parametrar och verkliga parametrar från borrhålsdata från ett känt norskt oljefält användes vid modelleringen. White's modell visade sig modellera relativt kraftig dämpning (5%) för medelstarkt konsoliderade gasreservoarer medan för oljereservoarer med motsvaranda konsolidering dämpningen var så låg (0.3%) att det är uppenbart att i en verklig situation skulle dämpningen inte vara mätbar. DPDP modelleringen verkar vara bättre på att beskriva dämpningen och gav dämpningar upp till 10% för en medelstarkt konsoliderad oljereservoar. Brist på parametrar från borrhålsdata gjorde att det inte var möjligt att på ett tillfredställande sätt modellera en verklig situation.Dock visade syntetisk data intressant karaktäristik och det rekommenderas därför att mer och detaljerade borrhålsparametrar mäts om ytterligare forskning om detta ska genomföras. För den andra delen av avhandlingen har Svenska Petroleum Exploration AB och Det Norske Oljeselskap ASA bidragit med stackad seismisk data som även var spectralanalyserad för indikationer på frekvensberoende dämpningsvariationer (utfört med fouriertransform och komplex spectraldekomposition). Tolv områden på den stackade kuben analyserades; sex oljemättade och sex som antogs vara vattenmättade. I varje område valdes en huvudtracé och de två närmaste tracéerna på vardera sida (totalt fem tracéer). Metoden med komplex spectraldekomposition klarade inte att analysera signalen från den stackade sektionen, varför fouriertransform användes för vidare analys. Frekvensanalysen gav en topp vid ~30 Hz för både olje- och vattenmättade reservoarer vilket tycks vara en karaktäristisk frekvens för källan. Detta kunde tyvärr inte bekräftas och tiden räckte inte till för att testa antagandet. Fouriertransformen tycks visa en viss skillnad mellan olje- och vattenmättade tracéer, men det kan också bero på skillnad i litologin snarare än porvätskan. Där för rekommenderas vid fortsättning på denna forskning att 4D seismisk data används för att analysera samma område men med data från olika tidpunkter. Det rekommenderas även att ostackad eller råa skott-data används eftersom väsentlig information kan försvinna när data stackas. / <p>Advisor present: Dr. Chris Juhlin</p><p>Examiner: Dr. Milovan Urosevic</p><p>Opponent: Álvaro Polín Tornero</p>

Biot's theory when applied to homogeneous media (involving the macroscopic flow mechanism) cannot explain the high level of attenuation observed in natural porous media over the seismic frequency range. However, several successful mesocopic inhomogeneity models have been developed to account for P wave attenuation. In this thesis I further develop the approaches to tackle S wave velocity and attenuation, to simulate transient wave propagation in poroelastic media, and to construct new models for determining the effective parameters of porous media containing mesoscopic inhomogeneities. As an important application of the double-porosity dual-permeability (DPDP) model, I have reformulated the effective Biot model using the total-field variables. This gives rise to new and more general governing equations than the previous approach based on the host phase field variables (which become a special case of the more general treatment). The analytical transient solution and dispersion characteristics for the double-porosity model and also for a poro-viscoacoustic model are derived over the entire frequency range for a homogeneous medium. The comparison between the results of the two models shows that dissipation by local mesoscopic flow of the double porosity model is very hard to fit by a single Zener element over a broad band. I chose the relaxation function to approximate the dispersion behaviour of the double porosity model just around the source centre frequency. It is shown that for most water-filled sandstones having a double porosity structure, wave propagation can be well described by the poro-viscoaoustic model with a single Zener element in the seismic frequency range. The transient solution for heterogeneous double porosity media is obtained by a numerical pseudospectral time splitting technique. This method is extended to 2.5-D poro-viscoelastic media to capture both P and S wave behaviour. I also demonstrate that if the frequency is below several Hz, then a single Kelvin-Voigt element gives an even better result than a single Zener element. I propose a two-phase permeability spherical inclusion model and obtain the dispersion curves of phase velocity and dissipation factor for the composite. I then determine the effective dynamic permeability of porous media with mesoscopic heterogeneities over the whole frequency range. This result is used to check the validity of other measures of effective dynamic permeability, deduced from the effective hydraulic permeability by replacing the permeability of the components with their dynamic values as determined from the Johnson, Koplik and Dashen (JKD) model. I also investigate the scattering of plane transverse waves by a spherical porous inclusion embedded in an infinite poroelastic medium. The vector displacement wave equations of Biot’s theory are solved as an infinite series of vector spherical harmonics for the case of a plane S-wave incidence. Then, the non-self-consistent theory is used to derive the dispersion characteristics of shear wave velocity and attenuation for a porous rock having mesoscopic spherical inclusions which are designed to represent either the patchy saturation model or the double porosity model with dilute concentrations of identical inclusions. / http://proxy.library.adelaide.edu.au/login?url= http://library.adelaide.edu.au/cgi-bin/Pwebrecon.cgi?BBID=1457632 / Thesis (Ph.D.) -- University of Adelaide, School of Chemistry and Physics, 2009

Au cours de leur propagation, les ondes sismiques sont atténuées par deux phénomènes : l'absorption causée par les propriétés anélastiques des matériaux, d'une part, et la diffusion ou " scattering " causée par la présence d'hétérogénéités de petite échelle dans le milieu d'autre part. L'objectif de cette thèse est de cartographier les propriétés de diffusion et d'absorption des ondes sismiques dans deux contextes géophysiques extrêmes présentant des échelles spatiales très différentes. La première partie du manuscrit est consacrée à la stratification d'hétérogénéité dans la Lune. À l'aide d'un modèle original de diffusion en géométrie sphérique, nous avons inversé les mesures de temps d'arrivée du maximum d'énergie et de décroissance de la coda sismique réalisée sur les données des missions Apollo. Nos inversions mettent en évidence un très fort contraste des propriétés de scattering entre le mégarégolithe très atténuant et le manteau lunaire profond transparent. L'atténuation est très largement dominée par le scattering et suggère la présence de fracturation jusqu'à environ 100 km de profondeur, affectant ainsi le manteau. Une nouvelle méthode d'estimation de la profondeur des séismes superficiels fondée sur les signaux diffus a été développée et permet de confirmer l'existence de failles actives autour de 50 km de profondeur. La deuxième partie de la thèse est consacrée à la structure d'atténuation de Taïwan, une région qui présente des structures géologiques très variées dans un contexte tectonique de double subduction. On utilise la MLTWA (Multiple Lapse Time Window Analysis) -une méthode fondée sur le rapport entre énergie cohérente et incohérente du signal sismique- pour imager les variations latérales d'atténuation. Dans un premier temps nous avons travaillé dans l'hypothèse classique de diffusion isotrope dans un demi-espace. Nos résultats mettent en évidence un niveau d'atténuation globale très élevé ainsi que de forts contrastes des propriétés de scattering sur des échelles spatiales fines, de l'ordre de 10-20 km. La diffusion est particulièrement marquée dans les bassins de la côte ouest, le sud et la chaîne côtière associée à la collision avec l'arc volcanique de Luçon à l'est. L'absorption augmente graduellement vers l'est et atteint son maximum sous l'arc volcanique. L'examen de l'accord entre données et modèles a posteriori montre sans ambiguïté les limites de l'hypothèse de diffusion isotrope dans un demi-espace sur un ensemble de stations situées le long des côtes. Ceci nous conduit à explorer les effets de la diffusion anisotrope dans un guide d'onde modélisant la croûte. La prise en compte de l'anisotropie améliore significativement l'accord du modèle aux données. En particulier, à basse fréquence (1-2 Hz), notre étude démontre la prédominance de rétro-diffusion. Ce résultat est compatible avec la présence de forts contrastes d'impédance dans la croûte et suggère la présence massive de fluides dans les zones de failles et de volcans à Taïwan. La mesure de l'anisotropie de la diffusion ouvre des perspectives nouvelles de caractérisation des hétérogénéités géophysiques de petite échelle. / During their propagation, seismic waves are attenuated by two phenomena: on one hand, absorption caused by the anelastic properties of the materials, and on the other hand, scattering caused by the presence of small-scale heterogeneities in the medium. The aim of this thesis is to map the properties of scattering and absorption of seismic waves in two extreme geophysical contexts with very different spatial scales. The first part of this memoir is devoted to the stratification of heterogeneities in the Moon. We use a new diffusion model in spherical geometry to invert measurements of the time of arrival of the maximum of energy and the seismic coda decay on data from the Apollo missions. Our inversions provide evidence of a very sharp contrast of scattering properties between the highly attenuating megaregolith and the transparent deep lunar mantle. Attenuation is largerly dominated by scattering and suggests the presence of fractures down to about 100 km depth, into the mantle. A new method for estimating the depth of shallow moonquakes based on diffusive signals was developped and confirms the existence of active faults around 50 km deep. The second part of this thesis is devoted to the attenuation structure of Taiwan, a region with a wide variety of geological structures in the context of two subduction zones. We use the MLTWA (Multiple Lapse Time Window Analysis) -a method based on the ratio between the coherent and incoherent energy of the seismic signal- to image the lateral variations of attenuation. We worked first with the classical hypothesis of isotropic scattering in a half-space. Our results provide evidence for a globally high level of attenuation with sharp contrasts of scattering properties across small spatial scales, of the order of 10-20 km. Scattering is particularly strong in the basins of the west coast, southern Taiwan and the eastern Coastal Range associated with the collision with the Luzon volcanic arc. Absorption increases gradually eastwards and reaches a maximum below the volcanic arc. A posteriori examination of the fit between data and model shows unambiguously the limits of the hypothesis of isotropic scattering in a half-space for a number of stations located along the coasts. This leads us to explore the effects of anisotropic scattering in a guide for seismic waves representing the crust. Taking anisotropy into account significantly improves the fitness of the model to the data. In particular, at low frequency (1-2 Hz), our study shows the prevalence of backscattering. This result is compatible with the presence of sharp contrasts of impedance in the crust and suggests the strong presence of fluids in fault zones and volcanoes in Taiwan. The measurement of scattering anisotropy opens new perspectives for characterizing small-scale geophysical heterogeneities.

Lateral perturbations in anelasticity (Q) and wave speed together provide important constraints on thermal and chemical structures in the mantle. In present-day tomography studies of global wave speed and anelasticity, the significance of 3-D wave speed and 3-D Q structures on surface wave travel times and amplitudes has not been well understood. In this dissertation, the effects of lateral perturbations in anelasticity (Q) and wave speed on surface wave observables are quantified based upon wave propagation simulations in 3-D earth models using a Spectral Element Method. Comparison between phase delays caused by 3-D wave speed structures and those caused by 3-D Q variations show that anelastic dispersion due to lateral perturbation in Q is important in long-period surface wave and can account for 15-20% observed phase delays. For amplitude perturbations, elastic focusing/defocusing effects associated with 3-D wave speed structures are dominant while energy dissipation is important in short-period (~ 50 s) surface waves but decreases quickly with increasing wave period. Anelastic focusing/defocusing associated with 3-D anelastic dispersion becomes more important than wave attenuation in longer period surface waves. In tomography studies, ray theory breaks down and finite frequency effects become important when the length scale of heterogenities are smaller than seismic wavelength. Finite frequency effects in 3-D earth models are investigated by comparing theoretical predictions of travel times and amplitudes with "ground truth" measurements made on synthetic seismograms generated in SEM simulations. The comparisons show that finite frequency effects are stronger in amplitudes than in phases, especially at long periods. / Ph. D.

Acquisition and interpretation of 3-D seismic data led DeAngelo and Hardage (2001) to describe the channel system in the south central Cut Bank area in Glacier County, Montana. The presence of a low velocity layer called Bentonite was also discovered in the area with the help of well-logs. Bentonite is a volcanic ash, which lies on both sides of the channel system and is absent within the channel. DeAngelo and Hardage (2001) shot a vertical seismic profiling (VSP) survey at well # 54-8 to analyze the formation structure in depth, since seismic signals around the reservoir area were unclear in the 3-D survey. This research attempts to estimate the attenuation properties of the Bentonite layer in the Cut Bank oil field. VSP data is processed for velocity information and estimation of seismic Q using the spectral ratios method (SRM). The SRM theoretically assumes that the propagating signal is a plane seismic wave traveling vertically from one point to another in a homogeneous model. The amplitudes at the start and end points are known and relate to each other with the attenuation coefficient in a frequency range. The relation between the seismic amplitudes at z distance from each other can be expressed as a linear function of frequency after a few modifications. SRM uses the linearity of the logarithmic ratio of the seismic amplitudes over a frequency range. In theory, ratios plotted against a frequency range must produce a flat line. However, in practice, the logarithmic ratios are expected to draw an approximate line (curve), where some of the data points deviate from the origin of the line. Thus fitting a line to the ratios curve and calculating the slope of this curve are necessary. Slope of the curve relates to the seismic attenuation coefficient and further to the seismic Q. The SRM results suggest that Bentonite may have a Q value as low as 5. This highly attenuative and thin (20 to 40 feet throughout the south central Cut Bank Unit) layer alters seismic signals propagating through it. A thorough analysis of the amplitude spectra suggests that seismic signals dramatically lose their energy when they pass through Bentonite. Low energy content of the signals below the Bentonite layer highlights that the recovery of the seismic energy is less likely despite the presence of multiples, which are known to affect the seismic signals constructively. Therefore, separation of reflected events is a greater challenge for the thin reservoir sand units lying underneath the Bentonite layer. Thus the Bentonite layer in the Cut Bank oil field has to be taken seriously and data processing should be done accordingly for better accuracy.

Cette thèse s’intéresse à divers aspects de l’atténuation sismique dans le manteau terrestre et aux implications de celle-ci quant à la structure de ce dernier. L’enjeu est de mieux comprendre les mécanismes d’atténuation ainsi que les mesures que l’on peut effectuer afin d’améliorer les modèles radiaux d’atténuation dont on dispose et in fine l’interprétation des modèles de tomographie. Je me suis concentrée sur deux exemples de mécanismes d’atténuation, appartenant à deux grands types d’atténuation : l’atténuation intrinsèque, liée à l’absorption par le milieu d’une partie de l’énergie sismique dissipée irréversiblement sous forme de chaleur, et l’atténuation extrinsèque, liée à la dispersion de cette énergie par le milieu. Dans le premier cas, j'ai regardé l’effet des transitions de phase sur l’atténuation des ondes sismiques. En appliquant un modèle thermomécanique développé par Ricard et al., 2009, pour prédire l’atténuation des ondes sismiques liée à la transition de phase uniquement et en comparant les valeurs obtenues aux mesures dont on dispose, j'ai pu contraindre la cinétique d’une transition de phase mantellique. Dans le second cas, j'ai testé l’effet de l’anisotropie comme mécanisme d'atténuation apparente, le but étant de trouver une distribution statistique d’orientation d’anisotropie pouvant reproduire la quasi-constance du facteur de qualité Q avec la fréquence, observée en sismologie et lors d’expériences de laboratoire (Knopoff, 1964), et aujourd’hui expliquée par un modèle ad-hoc seulement (Liu, 1976).Enfin, je me intéressée à mesurer cette atténuation sismique sur des enregistrements réels. Après avoir testé la méthode dite de la fréquence instantanée (Ford et al., 2012), je me suis concentrée sur deux régions, l’Amérique centrale et l’Alaska pour l'appliquer. Ces mesures sont ensuite interprétées en termes de modèle radial d’atténuation révélant un manteau inférieur hétérogène atténuant. Je montre aussi qu’une origine compositionnelle est la plus probable pour expliquer ces anomalies d’atténuation. / This thesis is devoted to various aspects of seismic attenuation in the Earth's mantle and the consequences on the mantle structure. The challenge is to better understand the attenuation mechanisms, as well as the measurements that can be done, in order to improve the published radial profiles of attenuation and in fine the interpretation of tomographic models.I focus on two examples of attenuation mechanisms, belonging to two kinds of attenuation: the intrinsic attenuation related to the absorption by the medium of a part of the seismic energy then irreversibly dissipated as heat, and the extrinsic attenuation related to the dispersion of the seismic energy by the medium. In the first case, I investigate the effect of phase transitions upon seismic attenuation. Applying the thermo-mechanical model developped by Ricard et al., 2012, to calculate the attenuation of seismic waves due to the phase transition only and comparing the obtained values to published measurements, I succeed in constraining the kinetics of a mantle phase transition. In the second case, I test the seismic anisotropy as a mechanism of extrinsic attenuation, the aim being to find a statistical distribution of anisotropy orientation and layer thicknesses that can reproduce the observed quasi-frequency independence of Q in seismology and laboratory experiments (Knopoff, 1964), and which is, today, only explained by an ad-hoc model (Liu, 1976).Finally, I was interested in measuring the seismic attenuation on real seismograms. After having tested the method of the instantaneous frequency (Ford et al., 2012), I applied it to seismic records sampling the mantle below Central America and Alaska. These measurements are then inverted for a radial profile of shear attenuation which reveals the existence of an attenuating zone in the lower mantle. I also show that these attenuation anomalies are likely to be of chemical origin.

L’évolution des assemblages pétrologiques avec l’augmentation de pression et de température est couramment perçue à l’équilibre et figée dans le temps. Le développement des sources synchrotron de rayons X permet à présent de mesurer in situ et en temps résolu les vitesses de transformations minéralogiques à haute pression (HP), haute température (HT). Cette thèse présente l’utilisation de cette technologie dans 2 contextes géologiques. (i) Le potentiel sismogène de la déstabilisation des minéraux hydratés dans les plaques en subduction est vérifié. Les cinétiques de déshydratation du talc, de la phase à 10Å et de l’antigorite ont été mesurées à HP-HT en système clos. Nous avons identifié que l’antigorite se déshydrate en passant par un stade intermédiaire. Toutes les vitesses de libération de fluides associées sont plus rapides que la déformation visqueuse des roches et sont donc compatibles avec le déclenchement de rupture. (ii) Les cinétiques de transition olivine-ringwoodite ont été déterminées dans la loupe de costabilité pour des compositions riches en Fe. Elles mettent en évidence une amorphisation partielle de l'olivine en début de transformation. Cela pourrait perturber de manière significative la vitesse des ondes sismiques lors de leur passage au niveau de la zone de transition mantellique. Par ailleurs, les temps caractéristiques de réaction et la réduction conséquente de la taille des grains, indiquent qu’une telle transition de phase induit une atténuation sismique importante. Ces résultats expérimentaux in situ HP-HT révèlent des mécanismes originaux de transition de phase et contribuent ainsi à une meilleure compréhension des modèles géodynamiques / The evolution of petrological assemblies with increasing pressure and temperature is commonly perceived at equilibrium and fixed within time. The development of X-ray synchrotron sources now enable to measure in situ, time-resolved rates of mineralogical transformations at high pressure (HP), high temperature (HT). This thesis presents the application of this technology in two geological settings. (i) The seismogenic ability of breakdown of hydrated minerals within the subducting slab is checked. The dehydration kinetics of talc, 10Å phase and antigorite were measured at HP-HT in a closed system. We have found that antigorite dehydrates through an intermediate stage. All associated rates of released fluids are faster than the viscous deformation of rocks and are therefore compatible with the trigger of rupture. (ii) The kinetics of olivineringwoodite transition were determined within the co-stability loop for Fe-rich compositions. They show a partial amorphization of olivine at the beginning of the transformation. This could significantly affect the velocities of seismic waves when crossing the mantle transition zone. Moreover, the characteristic times of this reaction and the substantial reduction in grain size, indicate that such a phase transition may induce a significant seismic attenuation. These in situ HP-HT experimental results reveal novel mechanisms of phase transition and thus contribute to a better understanding of geodynamic models

The primary focus of this thesis is to examine the detailed seismic structure of the northern Cascadia margin, including the Cascadia basin, the deformation front and the continental shelf. The results of this study are contributing towards understanding sediment deformation and tectonics on this margin. They also have important implications for exploration of hydrocarbons (oil and gas) and natural hazards (submarine landslides, earthquakes, tsunamis, and climate change). The first part of this thesis focuses on the role of gas hydrate in slope failure observed from multibeam bathymetry data on a frontal ridge near the deformation front off Vancouver Island margin using active-source ocean bottom seismometer (OBS) data collected in 2010. Volume estimates (∼ 0.33 km^3) of the slides observed on this margin indicate that these are capable of generating large (∼ 1 − 2 m) tsunamis. Velocity models from travel time inversion of wide angle reflections and refractions recorded on OBSs and vertical incidence single channel seismic (SCS) data were used to estimate gas hydrate concentrations using effective medium modeling. Results indicate a shallow high velocity hydrate layer with a velocity of 2.0 − 2.1 km/s that corresponds to a hydrate concentration of 40% at a depth of 100 m, and a bottom simulating reflector (BSR) at a depth of 265 − 275 m beneath the seafloor (mbsf). These are comparable to drilling results on an adjacent frontal ridge. Margin perpendicular normal faults that extend down to BSR depth were also observed on SCS and bathymetric data, two of which coincide with the sidewalls of the slump indicating that the lateral extent of the slump is controlled by these faults. Analysis of bathymetric data indicates, for the first time, that the glide plane occurs at the same depth as the shallow high velocity layer (100±10 mbsf). In contrast, the glide plane coincides with the depth of the BSR on an adjacent frontal ridge. In either case, our results suggest that the contrast in sediments strengthened by hydrates and overlying or underlying sediments where there is no hydrate is what causing the slope failure on this margin. The second part of this dissertation focuses on obtaining the detailed structure of the Cascadia basin and frontal ridge region using mirror imaging of few widely spaced OBS data. Using only a small airgun source (120 cu. in.), our results indicate structures that were previously not observed on the northern Cascadia margin. Specifically, OBS migration results show dual-vergence structure, which could be related to horizontal compression associated with subduction and low basal shear stress resulting from over-pressure. Understanding the physical and mechanical properties of the basal layer has important implications for understanding earthquakes on this margin. The OBS migrated image also clearly shows the continuity of reflectors which enabled the identification of thrust faults, and also shows the top of the igneous oceanic crust at 5−6 km beneath the seafloor, which were not possible to identify in single-channel and low-fold multi-channel seismic (MCS) data. The last part of this thesis focuses on obtaining detailed seismic structure of the Vancouver Island continental shelf from MCS data using frequency domain viscoacoustic full waveform inversion, which is first of its kind on this margin. Anelastic velocity and attenuation models, derived in this study to subseafloor depths of ∼ 2 km, are useful in understanding the deformation within the Tofino basin sediments, the nature of basement structures and their relationship with underlying accreted terranes such as the Crescent and the Pacific Rim terranes. Specifically, our results indicate a low-velocity zone (LVZ) with a contrast of 200 m/s within the Tofino basin sediment section at a depth 600 − 1000 mbsf over a lateral distance of 10 km. This LVZ is associated with high attenuation values (0.015 − 0.02) and could be a result of over pressured sediments or lithology changes associated with a high porosity layer in this potential hydrocarbon environment. Shallow high velocities of 4 − 5 km/s are observed in the mid-shelf region at depths > 1.5 km, which is interpreted as the shallowest occurrence of the Eocene volcanic Crescent terrane. The sediment velocities sharply increase about 10 km west of Vancouver Island, which probably corresponds to the underlying transition to the Mesozoic marine sedimentary Pacific Rim terrane. High attenuation values of 0.03 − 0.06 are observed at depths > 1 km, which probably corresponds to increased clay content and the presence of mineralized fluids. / Graduate / 0373 / 0372 / 0605 / subbarao@uvic.ca