The uppermost crustal velocity structure of the Goldstream Valley, central Alaska is investigated using a series of five explosions that were detonated in schist bedrock and recorded at >120 local stations to develop 1-D and 3-D models of the upper crust. Simple refraction analyses reveal that both P- and S-wave arrival times are azimuth dependent, with the fastest velocities in the southeast and northeast directions. The Swave velocity structure of the upper crust is also determined through multiple filter analysis and a damped, least squares inversion of 0.2-2 sec period Rg waves. The shear wave velocity model from the surface-wave analysis is combined with the refraction analysis results to develop 1-D P- and S-wave models to a depth of 2 km. In order to better constrain P- and S-wave velocity variations both laterally and with depth throughout the Goldstream Valley, 3-D velocity models are produced using a numerical simulation model. / Thesis (MS) — Boston College, 2008. / Submitted to: Boston College. Graduate School of Arts and Sciences. / Discipline: Geology and Geophysics.

A new high-resolution 3-dimensional P-wave velocity model for Mt. Pinatubo volcano was developed by tomographic inversion of P-wave arrivals from 3,007 earthquakes recorded during a four month period from May to August, 1991. The arrivals were recorded by a network of seismic stations, consisting of seven pre-eruption stations and seven post-eruption stations. Two stations survived the June eruptions. First-arrival travel times were calculated using a finite-difference solution to the eikonal equation. An iterative, linearized approximation of the nonlinear tomography problem was used to solve separately for both velocity structure and hypocenter locations. Several inversions performed with different initial parameters and convergence schemes, and synthetic checkerboard reconstructions indicate a horizontal spatial resolution of velocity perturbations near 4 km. However, the network sparseness allows for a substantial trade-off between focal depth, origin time, and the vertical velocity profile. Many hypocenter clusters collapse from diffuse clouds into tighter features after 3-D relocation. These bands of earthquakes appear to represent fault-related structures. Three low-velocity (relative to the horizontal average) anomalies exist within the well-resolved portion of the velocity model. These anomalies are spatially associated with pre- and post-eruption earthquakes oriented along mapped surface fault zones. Similar anomalies observed at different volcanoes have been previously interpreted as magma related. The low-velocity anomalies at Pinatubo are interpreted as highly fractured, hot volumes of mostly competent rock, which may contain partial melt. / Master of Science

Title: Ray-based Born approximation Author: Libor Šachl Department: Department of Geophysics Supervisor: RNDr. Luděk Klimeš, DrSc. Supervisor's e-mail address: klimes@seis.karlov.mff.cuni.cz Abstract: One of the aims of this thesis was coding of program grdborn.for for computing the 2D and 3D ray-based Born approximation of the first order in an inhomoge- nous isotropic medium without attenuation. The computation of 3D amplitudes using the 2D Born approximation is based on the correction term, which is de- rived. The program is further used in computing the Born approximation in various models. We test its performance in three simple models. We study the effect of the discretization, the spurious waves introduced by the finite size of the grid etc. In the next step, we focus on the computations in more compli- cated models. We compute the Born seismograms in 2D heterogenous models. We study the diffracted waves, the effects of caustics etc. Keywords: Born approximation, ray theory, velocity model, perturbation 1

Dicelis, G. E., Estudo Integrado da Camada de Basalto em Bebedouro SP, Bacia do Paraná: Relocalização dos sismos induzidos por poços tubulares, 2011. Dissertação de Mestrado Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, SP. A partir de 2004, a população do distrito de Andes, Bebedouro (SP), passou a sentir pequenos tremores de terra. Uma rede sismográfica instalada em 2005 mostrou que os sismos locais são induzidos por poços tubulares profundos usados para irrigação de lavouras (Assumpção et al., 2007). Neste projeto de pesquisa foi melhorada a localização desses sismos, usando análise da forma de onda e os tempos de chegada das ondas P e S, com o propósito de efetuar uma melhor análise da sismicidade induzida. Para isto foram aproveitados resultados de estudos anteriores como sísmica de refração, dispersão de ondas de superfície, sondagem elétrica vertical, função do receptor e TDEM. Para obter o melhor modelo de velocidades 1D, foram testados diferentes algoritmos de localização e técnicas de relocalização de hipocentros, para aproximadamente 3000 microtremores gravados pela rede sismográfica no período 2005 2010. Os resultados mostraram que as profundidades focais dos sismos relocalizados ficaram em sua maioria dentro da camada de basalto, a menos de 500m de profundidade. Este resultado confirma que os sismos estão relacionados com a perturbação do equilíbrio de tensões (pressões) geológicas, em zonas de fraturas já existentes dentro da camada de basalto. Foi usada uma combinação de correlação cruzada e tempos de chegada para analisar um conjunto de sismos, baseada na similaridade dos seus sismogramas. Foi melhorada a distribuição hipocentral o que permitiu determinar claramente uma solução do plano focal de um cluster, o que indicou a existência de uma falha normal de orientação WNW-ESSE e mergulho para o norte, com extensão NNE-SSW (eixo-T). Este mecanismo pode complementar os dados de esforços na região sudeste. / Dicelis, G. E., Integrated Study of the basalt layer in Bebedouro SP, Paraná Basin: Relocation of earthquakes induced by wells, 2011. Master dissertation - Institute of Astronomy, Geophysics and Atmospheric Sciences of the University of Sao Paulo, Brazil. Since 2004 the population of the district of Andes, Bebedouro SP, began to feel small earthquakes. A seismograph network installed in 2005 showed that earthquakes are caused by local deep wells used for irrigation of crops (Assumpção et al., 2007). In this research project we improved the location of these earthquakes, using analysis of waveform and the arrival times of P and S waves, in order to make a better analysis of the induced seismicity. To better relocate the events we determined a velocit model using results of previous studies such as seismic refraction, surface wave dispersion, vertical electrical sounding, receiver function and TDEM. The best fit 1D model was achieved testing different localization algorithms and techniques for relocation of hypocenters for approximately 3000 microtremors recorded by seismograph network to the period 2005 2010. The results showed that the focal depths of the relocated earthquakes are mostly within the basalt layer, less than 500m depth. This confirms that the earthquakes are related with disturbance of geological stresses (pressures) in pre-existing fracture zones within the basalt layer. We used a combination of cross-correlation and arrival times for analyze a set of earthquakes, based on the similarity of their recorded seismograms. The hypocentral distribution was improved allowing a clear determination of the best fit fault plane solution for one cluster, which indicates the existence of a normal fault with direction WNW-ESE and dip to the north, with NNE-SSW extension (T-axis). This mechanism may complement the stress data in the Southeast Brazil.

Dicelis, G. E., Estudo Integrado da Camada de Basalto em Bebedouro SP, Bacia do Paraná: Relocalização dos sismos induzidos por poços tubulares, 2011. Dissertação de Mestrado Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, SP. A partir de 2004, a população do distrito de Andes, Bebedouro (SP), passou a sentir pequenos tremores de terra. Uma rede sismográfica instalada em 2005 mostrou que os sismos locais são induzidos por poços tubulares profundos usados para irrigação de lavouras (Assumpção et al., 2007). Neste projeto de pesquisa foi melhorada a localização desses sismos, usando análise da forma de onda e os tempos de chegada das ondas P e S, com o propósito de efetuar uma melhor análise da sismicidade induzida. Para isto foram aproveitados resultados de estudos anteriores como sísmica de refração, dispersão de ondas de superfície, sondagem elétrica vertical, função do receptor e TDEM. Para obter o melhor modelo de velocidades 1D, foram testados diferentes algoritmos de localização e técnicas de relocalização de hipocentros, para aproximadamente 3000 microtremores gravados pela rede sismográfica no período 2005 2010. Os resultados mostraram que as profundidades focais dos sismos relocalizados ficaram em sua maioria dentro da camada de basalto, a menos de 500m de profundidade. Este resultado confirma que os sismos estão relacionados com a perturbação do equilíbrio de tensões (pressões) geológicas, em zonas de fraturas já existentes dentro da camada de basalto. Foi usada uma combinação de correlação cruzada e tempos de chegada para analisar um conjunto de sismos, baseada na similaridade dos seus sismogramas. Foi melhorada a distribuição hipocentral o que permitiu determinar claramente uma solução do plano focal de um cluster, o que indicou a existência de uma falha normal de orientação WNW-ESSE e mergulho para o norte, com extensão NNE-SSW (eixo-T). Este mecanismo pode complementar os dados de esforços na região sudeste. / Dicelis, G. E., Integrated Study of the basalt layer in Bebedouro SP, Paraná Basin: Relocation of earthquakes induced by wells, 2011. Master dissertation - Institute of Astronomy, Geophysics and Atmospheric Sciences of the University of Sao Paulo, Brazil. Since 2004 the population of the district of Andes, Bebedouro SP, began to feel small earthquakes. A seismograph network installed in 2005 showed that earthquakes are caused by local deep wells used for irrigation of crops (Assumpção et al., 2007). In this research project we improved the location of these earthquakes, using analysis of waveform and the arrival times of P and S waves, in order to make a better analysis of the induced seismicity. To better relocate the events we determined a velocit model using results of previous studies such as seismic refraction, surface wave dispersion, vertical electrical sounding, receiver function and TDEM. The best fit 1D model was achieved testing different localization algorithms and techniques for relocation of hypocenters for approximately 3000 microtremors recorded by seismograph network to the period 2005 2010. The results showed that the focal depths of the relocated earthquakes are mostly within the basalt layer, less than 500m depth. This confirms that the earthquakes are related with disturbance of geological stresses (pressures) in pre-existing fracture zones within the basalt layer. We used a combination of cross-correlation and arrival times for analyze a set of earthquakes, based on the similarity of their recorded seismograms. The hypocentral distribution was improved allowing a clear determination of the best fit fault plane solution for one cluster, which indicates the existence of a normal fault with direction WNW-ESE and dip to the north, with NNE-SSW extension (T-axis). This mechanism may complement the stress data in the Southeast Brazil.

L'inversion des formes d'ondes (FWI) est une procédure d'imagerie sismique pour imager le sous-sol de la Terre. FWI est résolue comme un problème d'optimisation. En fonction du contenu en fréquence des données, la fonction objective de FWI peut être fortement non linéaire. Pour des données associées des réflexions, ce problème empêche notamment les méthodes basées sur le gradient de retrouver les grandes longueurs d'onde du modèle de vitesse. Dans cette thèse, nous proposons une variante de FWI basée sur la séparation des champs d'ondes, typiquement en champs montants et descendants, pour atténuer la non-linéarité du problème. Il consiste à décomposer le gradient de FWI en une partie de courte longueur d'onde et une partie de grande longueur d'onde après décomposition des champs d'ondes. L'inversion est effectuée d'une manière alternée entre ces deux parties. Nous appliquons cette méthode à plusieurs études de cas et montrons que la nouvelle approche est plus robuste en particulier pour la construction du modèle de grande longueur d'onde. / Full Waveform Inversion (FWI) is a seismic imaging procedure to image the subsurface of the Earth. FWI is resolved as an optimization problem . Depending on the frequency content of the data, the objective function of FWI may be highly nonlinear. If a data set mainly contains reflections, this problem particularly prevents the gradient-based methods from recovering the long wavelengths of the velocity model.In this thesis, I propose a variant of FWI based on the wavefield separation, typically between up- and down- going waves, to mitigate the nonlinearity of the problem. The new method consists of decomposing the gradient of FWI into a short-wavelength part and a long-wavelength part after wavefield decomposition. The inversion is performed in an alternating fashion between these two parts. We apply this method to several case studies and show that the new method is more robust especially for constructing the long-wavelength model.

Sedimentation is a fundamental operation in wastewater treatment works. A rational design of sedimentation tanks is currently achieved by plotting iso-percentile (iso-percentage) concentration removal profiles from flocculent settling data. A major drawback of the graphical iso-percentage method is that the iso-percentile lines are often manually interpolated and are mere hand drawn estimations. This is because the settling behaviour of sludge particles is highly non-linear. The manual analytical process is therefore very tedious, inaccurate and subjective. Hence, an optimised design of sedimentation tanks is necessary in order to eliminate the errors incurred during data analysis. In this study, a mechanistic iso-percentile flocculent model (referred to as the velocity flocculation model) is developed to simulate the behaviour of flocculating colloidal particles in turbid water. This model is based on the physical meanings of flocculent settling particles and on fractal theory. It is formulated to produce automated iso-percentile curves which are fundamental in the design of sedimentation tanks. The iso-percentile model was vertically integrated into a velocity model to produce a model expressing the velocity of particles as a function of removal rate. The velocity model has an obvious advantage over the iso-percentile model in that it is easy to contextualize. It can be reverted back to the iso-percentile trajectory analysis eliminating the need for extensive data interpolation and may in future eliminate the need for settling column analysis altogether. In the current study, the integrated velocity form is used to predict instantaneous flocculent settling velocity of fine suspended particles under near quiescent conditions. This is vital since it is difficult to obtain velocity values in-situ or directly from sedimentation tanks. Model validity and competency was tested by a direct comparison with existing literature models, such as Ozer’s model and Ramatsoma and Chirwa’s model. Model comparison was based on the goodness of fit, the least sum of square errors and mathematical consistency with known flocculent settling behaviour. The newly developed iso-percentile model achieved a more accurate simulation of physical experimental data, did not violate any of the mathematical constraints and yielded lower sum of square errors than originally achieved by Ozer and Ramatsoma and Chirwa. Notably, the proposed velocity model offers a distinctive advantage over conventional interpolated-iso-percentile based models which are prone to numerical errors during interpolation. Its performance (velocity model) was compared against Je and Chang’s velocity model. Higher velocity values were observed for the new model than for Je and Chang’s model implying that empirically based models would tend to under-predict the velocity values. The model developed in this study brings us one step closer to achieving full automation of the settling tank and clarifier design. / Dissertation (MEng)--University of Pretoria, 2013. / gm2014 / Chemical Engineering / unrestricted

Microearthquake (< M3.0) swarms occur frequently in volcanic environments, but do not always culminate in an eruption. Such non-eruptive swarms may be caused by stresses induced by magma intrusion, hydrothermal fluid circulation, or regional tectonic processes, such as slow-slip earthquakes. Strandline Lake, located 30 km northeast of Mount Spurr volcano in south-central Alaska, experienced an intense earthquake swarm between August 1996 and August 1998. The Alaska Volcano Observatory (AVO) catalog indicates that a total of 2,999 earthquakes were detected during the swarm period, with a maximum magnitude of Mw 3.1 and a depth range of 0-30 km below sea level (with the majority of catalog hypocenters located between 5-10 km BSL). The cumulative seismic moment of the swarm was 2.03e15 N-m, equivalent to a cumulative magnitude of Mw 4.2. Because of the swarm's distance from the nearest Holocene volcanic vent, seismic monitoring was poor and gas and GPS data do not exist for the swarm period. However, combined waveforms from a dense seismic network on Mount Spurr and from several regional seismic stations allow reanalysis of the swarm earthquakes. I first developed a new 1-D velocity model for the Strandline Lake region by re-picking and inverting precise arrival times for 27 large Strandline Lake earthquakes. The new velocity model reduced the average RMS for these earthquakes from 0.16 to 0.11s, and the average horizontal and vertical location errors from 3.3 to 2.5 km and 4.7 to 3.0 km, respectively. Depths of the 27 earthquakes ranged from 10.5 to 22.1 km with an average depth of 16.6 km. A moderately high b-value of 1.33 was determined for the swarm period, possibly indicative of magmatic activity. However, a similarly high b-value of 1.25 was calculated for the background period. 28 well-constrained fault plane solutions for both swarm and background earthquakes indicate a diverse mixture of strike-slip, dip-slip, and reverse faulting beneath Strandline Lake. Finally, five Interferometric Synthetic Aperture Radar (InSAR) images spanning the swarm period unambiguously show no evidence of surface deformation. While a shallow volcanic intrusion appears to be an unlikely cause of the Strandline Lake swarm based on the new well-constrained earthquake depths and the absence of strong surface deformation, the depth range of 10.5 to 22.1 km BSL for relocated earthquakes and the high degree of FPS heterogeneity for this swarm are similar to an earthquake swarm beneath Lake Tahoe, California in 2003 caused by a deep intrusion near the base of the crust (Smith et al, 2004). This similarity suggests that a deep crustal magmatic intrusion could have occurred beneath the Strandline Lake area in 1996-1998 and may have been responsible for the resulting microearthquake activity.

Although the petroleum industry has imaged reflections from the sides of salt domes, steeply dipping structures have not been imaged as reflectors outside of sedimentary basins; to do so requires appropriate data acquisition, prestack depth migration, and an excellent seismic velocity model. Poststack time migrated seismic images, normal moveout velocity analysis, well logs, and other geologic information are used to build the velocity model. In regions of interest outside of sedimentary basins, such as major strike-slip faults, seismic reflectivity is often sparse and little is known of detailed subsurface geology. Alternate methods of velocity model construction must be used. First arrival (refraction and turning ray) traveltime tomography is proposed to construct the preliminary velocity model for steep-dip prestack depth migration in settings with little a priori subsurface information. A densely spaced synthetic seismic data set with long-offset recording, modeled after a real survey across the San Andreas Fault, was constructed using a finite-difference algorithm. First arrival traveltimes were picked from the data and a velocity model was constructed using tomography. The velocity model was used to perform a Kirchhoff prestack depth migration of the synthetic shot gathers. The subsurface structure was sufficiently reconstructed that the velocity model could be refined through migration velocity analysis. A series of tomography tests was used to determine the spatial resolution limits of the velocity model. Isolated erroneous anomalies with sizes near the resolution limits were added to the velocity model derived from tomography and used as input for migration. This pessimistic test provided an adequate image and identifiable arrivals in migrated common image gathers, allowing the velocity model to be improved through migration moveout analysis. Data acquisition requirements for tomography include long recording offsets and times, larger sources, and dense spacings, very similar to the requirements for steep-dip reflection imaging. / Master of Science

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.