Esta dissertação apresenta o método de análise espectral de ondas superficiais (SASW) para a obtenção das variações do módulo cisalhante (Gmáx) com a profundidade, no domínio das deformações muito pequenas (abaixo de 0,001%). O SASW é um método sísmico in situ, não destrutivo, baseado na geração e detecção de ondas Rayleigh e na natureza dispersiva desta onda. Pela aplicação de um impacto na superfície do solo e detecção da onda em vários pontos, através de dois receptores, é construída uma curva de dispersão (velocidade de fase versus comprimento de onda). Esta curva de dispersão é, então, invertida. A inversão é um processo analítico para a reconstrução do perfil de velocidade de onda de cisalhamento (VS), partindo-se da curva de dispersão experimental de campo. O módulo de cisalhamento máximo de cada camada é facilmente obtido a partir do perfil de VS. No conteúdo teórico da dissertação discutem-se propriedades dinâmicas dos solos e descrevem-se as equações que dominam a propagação das ondas elásticas, tanto em meios homogêneos como em meios estratificados. A metodologia desenvolvida para a obtenção das curvas de dispersão, através da realização de ensaios SASW, apresenta os resultados obtidos em ensaios realizados na Cidade Universitária em São Paulo, sendo esses resultados comparados com estimativas feitas a partir de correlações baseadas em ensaios SPT existentes. Essas comparações permitem concluir que a metodologia SASW é uma boa alternativa para a determinação do perfil de rigidez (Gmáx) do solo, concordando com o nível de deformação envolvido nos ensaios. São desenvolvidos estudos de sensibilidade do método para verificar a influência na mudança dos parâmetros assumidos (peso específico, coeficiente de Poisson e espessuras das camadas) no processo de redução de dados (inversão) sobre o perfil final de VS, concluindo-se que o parâmetro que apresenta maior influência é o coeficiente de Poisson. / This dissertation presents the spectral-analysis-of-surface-waves (SASW) method as a tool for obtaining the variations in the modulus shear (Gmax) with depth in the field of very small strains (below 0,001%). The SASW method is a nondestructive in situ seismic method, based on the generation and measurement of Rayleigh wave and on its dispersive characteristic nature. Throughout the implementation of an impact on the soil surface and the detection of the wave at various points by two receptors a dispersion curve is constructed (phase velocity versus wave-length). This dispersion curve is then inverted. Inversion is an analytical process for reconstructing the shear wave velocity profile from the experimental field. The shear modulus of each layer is readily obtained from the shear wave velocity profile. The theoretical content of the dissertation presents dynamic properties of the soils and is described in the equations that dominate the propagation of elastic waves, both in homogeneous media and in stratified media. The methodology developed to obtain the dispersion curves through the implementation of SASW test is defined, and results from tests carried out at the University Campus in São Paulo are presented and compared with values obtained from correlations based on SPT tests. These comparisons indicate that the SASW method is a good alternative to determine the profile of stiffness (Gmax) of the soil, agreeing with the level of deformation involved in the tests. Studies on the methods sensitivity are developed to verify the influence on the changing of the parameters given (natural unit weight, Poisson coefficient and thickness of layers) in reduction of data (inversion) on the final profile of VS. The conclusion is that the Poisson coefficient is the parameter with greater influence.

Thesis advisor: John E. Ebel / Thesis advisor: Alan L. Kafka / As the best seismic indicator of shear modulus, shear-wave velocity is an important property in engineering problems in near-surface site characterization. Several surface-wave methods have been developed to obtain the subsurface shear-wave velocity structure. This thesis compared three surface-wave methods, Spectral Analysis of Surface Waves (SASW) (Nazarian et al., 1983), Multichannel Analysis of Surface Waves (MASW) (Park et al., 1999), and Refraction Microtremor (ReMi) (Louie, 2001), to determine which method gives the best estimation of the 1-D shear-wave velocity profile of near-surface soils. We collected seismic data at three sites in the greater Boston area where there are direct measurements of shear-wave velocities for comparison. The three methods were compared in terms of accuracy and precision. Overall, the MASW and the ReMi methods have comparable quality of accuracy, whereas the SASW method is the least accurate method with the highest percentage differences with direct measurements. The MASW method is the most precise method among the three methods with the smallest standard deviations. In general, the MASW method is concluded to be the best surface-wave method in determining the shear-wave velocities of the subsurface structure in the greater Boston area.

Breast cancer is the most common cancer in women worldwide. In 2009, a novel imaging modality called Shear Wave Elastography (SWE), an ultrasound technique visualising the elasticity of tissue, was introduced to the field of clinical breast imaging. Because malignant tissues are generally stiffer than benign tissues, SWE supports the differentiation of benign / malignant solid breast lesions. However, no standard has yet been defined for the application and the evaluation of results. Furthermore, image evaluation has to be carried out directly from the ultrasound system, complicating long-term and multi-centre studies. This PhD thesis investigated the influences from the imaging process and image evaluation on SWE measurements. Various parameters were appraised with regard to their diagnostic performance, in order to define the best clinical standard. To define more complex image analysis, taking the parameters investigated into account, algorithms were devised to enable automatic assessment of B-mode and SWE images. In this work, influences from the imaging process and image evaluation on the SWE measurements were demonstrated. The influences investigated included: the impact from the region of interest and the imaging plane used; the individual variation in breast composition; the number of images considered and the pressure applied during imaging. The algorithms described within this work achieved a diagnostic accuracy similar to that of manual assessment by a radiology expert. This thesis demonstrated influences from the imaging process and image evaluation on the SWE measurements obtained. Taking these influences into consideration would complicate the clinical application of SWE imaging. However, automatic image evaluation as presented here would overcome this issue. Using the guidelines defined in this PhD thesis also allows for comparison of results taken from different imaging sites.

The upper mantle structure of Swedish part of Baltic Shield with its isotropic and anisotropic seismic velocity characteristics is investigated using telesesismic body waves (i.e. P waves and shear waves) recorded by the Swedish National Seismological Network (SNSN). Nonlinear high-resolution P and SV and SH wave isotropic tomographic inversions reveal velocity perturbations of ± 3 % down to at least 470 km below the network. Separate SV and SV models indicate several consistent major features, many of which are also consistent with P-wave results. A direct cell by cell comparison of SH and SV models reveals velocity differences of up to 4%. Numerical tests show that differences in the two S-wave models can only be partially caused by noise and limited resolution, and some features are attributed to the effect of large scale anisotropy. Shear-wave splitting and P-travel time residual analyses also detect anisotropic mantle structure. Distinct back-azimuth dependence of SKS splitting excludes single-layer anisotropy models with horizontal symmetry axes for the whole region. Joint inversion using both the P and S data reveals 3D self-consistent anisotropic models with well-defined mantle lithospheric domains. These domains of differently oriented anisotropy most probably retain fossil fabric since the domains' origin, supporting the idea of the existence of an early form of plate tectonics during formation of continental cratons already in the Archean. The possible disturbing effects of anisotropy on seismic tomography studies are investigated, and found to be potentially significant. P-wave arrival times were adjusted based on the estimates of mantle anisotropy, and re-inverted. The general pattern of the velocity-perturbation images was similar but changed significantly in some places, including the disappearance of a slab-like structure identified in the inversion with the original data. Thus the analysis demonstrates that anisotropy of quite plausible magnitude can have a significant effect on the tomographic images, and should not be ignored. If, as we believe, our estimates of anisotropy are reasonably correct, then the model based on the adjusted data should give a more robust and correct image of the mantle structure.

<p>Prostate cancer (PCa) is the most common non-cutaneous cancer in men in the United States. Early detection of PCa is essential for improving treatment outcomes and survival rates. However, diagnosis of PCa at an early stage is challenged by the lack of an imaging method that can accurately visualize PCas. Because pathological processes change the mechanical properties of the tissue, elasticity imaging methods have the potential to differentiate PCas from other prostatic tissues. Acoustic radiation force impulse (ARFI) imaging is a relatively new elasticity imaging method that visualizes the local stiffness variations inside soft tissue.</p><p>The work presented in this dissertation investigates the feasibility of prostate ARFI imaging. Volumetric ARFI data acquisition and display methods were developed to visualize anatomic structures and pathologies in <italic>ex vivo </italic>human prostates. The characteristic appearances of various prostatic tissues in ARFI images were identified by correlating ARFI images with McNeal's zonal anatomy and the correlated histological slides, in which prostatic pathologies were delineated by a pathologist blinded to the ARFI images. The results suggest ARFI imaging is able to differentiate anatomic structures and identify suspicious PCa regions in the prostate.</p><p>To investigate the correlation between ARFI displacement amplitudes and the underlying tissue stiffness in the prostate ARFI images, the mechanical properties of prostatic tissues were characterized using a quantitative method, based upon shear wave elasticity imaging (SWEI). Co-registered ARFI and SWEI datasets were acquired in excised prostate specimens to reconstruct the shear moduli of prostatic tissues. The results demonstrated that variations in ARFI displacement amplitudes were inversely related to the underlying tissue stiffness; and the reconstructed shear moduli of prostatic tissues had good agreements with those reported in literature. The study suggests the matched ARFI and SWEI datasets provide complementary</p><p> information about tissue's elasticity. </p><p>To increase the efficiency of the data acquisition, a novel imaging sequence was developed to acquired matched ARFI-SWEI datasets without increasing the number of excitations compared to a conventional ARFI imaging sequence. Imaging parameters were analyzed both theoretically and experimentally. An analytical model was derived to quantify the fundamental accuracy limit in the reconstructed shear modulus, and demonstrated good agreement with the experimental data. The novel sequence was demonstrated in tissue-mimicking phantoms.</p><p>Finally, ARFI imaging sequences were developed in a transrectal probe, and ARFI images were presented from <italic>in vivo</italic> data acquired in patients under radical prostatectomy. The <italic>in vivo</italic> ARFI images demonstrated decreased contrast and resolution as compared to the matched <italic>ex vivo</italic> ARFI data. However, prostate anatomy and some PCa were successfully visualized in the <italic>in vivo</italic> ARFI images. Thus, we conclude that ARFI imaging has the potential to provide image guidance for locating cancerous regions during PCa diagnosis and treatment.</p> / Dissertation

<p>Acoustic radiation force impulses (ARFI) have been used to generated transverse-traveling mechanical waves in various biological tissues. The velocity of these waves is related to a medium's stiffness and thus can offer useful diagnostic information. Consequently, shear wave velocimetry has the potential to investigate cardiac disease states that manifest themselves as changes in tissue stiffness (e.g., ischemia).</p><p> The work contained herein focuses on employing ARFI-based shear wave velocimetry techniques, similar to those previously utilized on other organs (e.g., breast, liver), for the investigation of cardiac tissue. To this end, ARFI excitations were used to generate slow-moving (under 3 m/s) mechanical waves in exposed myocardium (with access granted through a thoracotomy); these waves were then tracked with ultrasonic methods. Imaging techniques to increase frame-rate, decrease transducer/tissue heating, and reduce the effects of physiological motion were developed. These techniques, along with two shear wave velocimetry methods (i.e., the Lateral Time-to-Peak and Radon sum transformation algorithms), were utilized to successfully track shear wave propagation through the mid-myocardial layer <italic>in vitro</italic> and <italic>in vivo</italic>. <italic>In vitro</italic> experiments focused on the investigation of a shear wave anisotropy through the myocardium. This experimentation suggests a moderate shear wave velocity anisotropy through regions of the mid-myocardial layer. <italic>In vivo</italic> experiments focused on shear wave anisotropy (which tend to corroborate the aforementioned <italic>in vitro</italic> results), temporal/spatial stability of shear wave velocity estimates, and estimation of wave velocity through the cardiac cycle. Shear wave velocity was found to cyclically vary through the cardiac cycle, with the largest estimates occurring during systole and the smallest occurring during diastole. This result suggests a cyclic stiffness variation of the myocardium through the cardiac cycle. A novel, on-axis technique, the displacement ratio rate (DRR) method, was developed and compared to conventional shear wave velocitmetry and ARFI imaging results; all three techniques suggest a similar cyclic stiffness variation.</p><p> Shear wave velocimetry shows promise in future investigations of myocardial elasticity. The DRR method may offer a means for transthoracic characterization of myocardial stiffness. Additionally, the future use of transesophageal and catheter-based transducers presents a way of generating and tracking shear waves in a clinical setting (i.e., when epicardial imaging is not feasible). Lastly, it is hoped that continued investigations into the physical basis of these ARFI-generated mechanical waves may further clarify the relationship between their velocity in myocardium and material stiffness.</p> / Dissertation

<p>The focus of this work is to use geologic and geophysical methods to better understand the faults and fracture systems at Puna, in southeastern Hawaii, and southern Montserrat, in the Lesser Antilles. The particular interest is understanding and locating the deep fracture networks that are necessary for fluid circulation in hydrothermal systems. The dissertation first presents a study in which identification of large scale faulting places Montserrat into a tectonic context. Then follow studies of Puna and Montserrat that focus on faults and fractures of the deep hydrothermal systems.</p><p>The first chapter consists of the results of the SEA-CALIPSO experiment seismic reflection data, recorded on a 48 channel streamer with the active source as a 2600 in3 airgun. This chapter discusses volcaniclastic debris fans off the east coast of Montserrat and faults off the west coast. The work places Montserrat in a transtensional environment (influenced by oblique subduction) as well as in a complex local stress regime. One conclusion is that the stress regime is inconsistent with the larger arc due to the influence of local magmatism and stress.</p><p>The second chapter is a seismic study of the Puna hydrothermal system (PHS) along the Kilauea Lower East Rift Zone. The PHS occurs at a left step in the rift, where a fracture network has been formed between fault segments. It is a productive geothermal field, extracting steam and reinjecting cooled, condensed fluids. A network of eight borehole seismometers recorded >6000 earthquakes. Most of the earthquakes are very small (< M.2), and shallow (1-3 km depth), likely the result of hydrothermal fluid reinjection. Deeper earthquakes occur along the rift as well as along the south-dipping fault plane that originates from the rift zone.</p><p>Seismic methods applied to the PHS data set, after the initial recording, picking, and locating earthquakes, include a tomographic inversion of the P-wave first arrival data. This model indicates a high seismic velocity under the field that is thought to be an intrusion and the heat source of the hydrothermal system. A shear wave splitting study suggested the PHS fracture system is largely oriented rift-parallel with some orthogonal fractures. Shear wave splitting data also were used in a tomographic inversion for fracture density. The fracture density is high in the PHS, which indicates high permeability and potential for extensive fluid circulation. This has been confirmed by high fluid flow and energy generation. The high fracture density is consistent with the interpretation of a transfer zone between the rift segments where a fracture mesh would be expected. In Puna the transfer zone is a relay ramp.</p><p>The results from the PHS are used as an example to examine the proposed hydrothermal system at St. George's Hill, Montserrat. In southern Montserrat, hot springs and fumaroles suggest a deep hydrothermal system heated by local magmatism. A magnetotelluric study obtained resistivity data that suggest focused alteration under southeastern Montserrat that is likely to be along fault segments. Several faults intersect under SGH, making it the probable center of the hydrothermal system. At Puna, and also Krafla, Iceland, where faults interact is an area of increased permeability, acting as a model to be applied to southern Montserrat. The conclusion is that in both Puna and Montserrat large faults interact to produce local areas of stress transfer that lead to fracturing and permeable networks; these networks allow for high-temperature hydrothermal circulation.</p> / Dissertation

The development of a site-specific seismic hazard curve for a soil site requires the incorporation of site effects into the hazard calculation through the use of a site-specific amplification function. This study investigates the effect on the resulting soil hazard curves of different approaches to compute the site-specific amplification function. Amplification functions and their standard deviations can be developed using equivalent linear site response analyses. This study investigates the amplification function predictions of one-dimensional (1D) and two-dimensional (2D) site response analyses. For 1D analysis, one set of analyses are performed using time series (TS) input motions while a second set is performed using random vibration theory (RVT). One-dimensional site response analyses are performed for a shallow and a deep soil site and the results are compared for seismic hazard predictions. The influence of spatial variability introduced through randomization of site shear wave velocity (V[subscript S]) is also investigated. Shear wave velocity profile randomization does not significantly change the predicted amplification functions, except for the RVT analysis near the site period. At these periods, (V[subscript S]) randomization reduces the amplification function predicted by RVT making it more similar to the TS analysis prediction. The surface hazard at a site is dependent on the median amplification factor and its associated standard deviation. Spatial variability and uncertainties in soil properties across a site are often taken into account by modeling multiple 1D profiles in 1D site response analyses. However, this approach assumes that analyzing multiple 1D profiles captures accurately the effects of the true multi-dimensional spatial variability of the soil properties. This study investigates the results of two-dimensional (2D) site response analyses that incorporate spatial variability in the (V[subscript S]) profile through Monte Carlo simulation. Two-dimensional site response analyses are performed for 2D random fields generated with various statistical parameters (i.e. vertical and horizontal correlation distances) to investigate the effect of different levels of spatial variability on surface response across a region of interest (ROI). Two-dimensional site response analyses are performed for a shallow site. Results indicate that horizontal correlation distance has more influence on the analyses results than the vertical correlation distance. As the horizontal correlation distance increases, the median surface response spectrum across the ROI decreases. This reduction in median surface response is more pronounced around the site period. The influence of the vertical correlation distance is more pronounced when the horizontal correlation distance is large. As the vertical correlation distance increases, the median surface response spectrum across the ROI increases, which is more pronounced around the period of the motion. The predictions of 1D and 2D site response analyses modeling the (V[subscript S]) variability are compared. 1D analyses are performed on separately generated 1D (V[subscript S]) profiles (infinite horizontal correlation) and on the (V[subscript S]) profiles across the ROI of each 2D (V[subscript S]) field realization generated for 2D analysis (finite horizontal correlation). The results indicate that both sets of 1D analyses predict lower median response than 2D analyses. The 1D analyses with finite horizontal correlation display comparable levels of variability in the site response, however 1D analyses with infinite horizontal correlation display higher variability. / text

Evaluation of the effect of particle cementation on the stiffness of uniform sand was carried out by measuring compression wave velocities (VP) and shear wave velocities (VS) on both clean and artificially cemented specimens. Piezoelectric transducers (PT) were used to perform the majority of the measurements. Shear wave velocity (VS), shear moduli (G) and material damping ratio (D) of clean and artificially cemented specimens were also determined using resonant column (RC) testing. Linear (shearing strains ≤ 0.001%) and nonlinear (shearing strains > 0.001%) behavior of the specimens were evaluated in the resonant column tests. The sand selected for this investigation is commonly known as Hickory sand, from the Hickory formation, western Llano uplift, Texas. This material was selected for its grain geometry and gradation; it consists of uniformly graded sand with rounded particles. The sand specimens were artificially cemented with a solution of hydrated sodium silicate and water. Sodium silicate is an alkaline compound obtained from the reaction of sodium hydroxide and silica. All artificially cemented specimens and uncemented hickory sand specimens were formed by pluviation through air. The microstructure of the specimens was visually assessed with images obtained from both optical and scanning electron (SEM) microscopes. These images confirmed that the procedure used to form artificially cemented specimens provides cementation around the contacts while some grain-to-grain contact appears to be preserved. Seismic and drained strength measurements on Hickory sand specimens were obtained from different cement concentrations and compared with results from clean sand specimens. Among the findings of this investigation are: (1) the procedure to artificially cement sand specimens in the laboratory was successful, (2) the slopes (nP and nS) obtained from the relationships between compression and shear wave velocities with effective isotropic confining pressure in log-log scale decrease as the cement content increases, and (3) as increasing amounts of cement are added to the sand particles, the nonlinearity of the specimens increases up to certain amount of cement, after which the nonlinearity of the specimen decreases and tends towards rock-like behavior. / text

The behavior of unsaturated soils depends heavily on material properties and soil conditions. In Geotechnical Engineering, compacted soils are frequently used as fill material, and quality control is vital to the construction process. There are few methods available to estimate the parameters associated with unsaturated soils based on field measurements, and a relationship between these factors could reduce testing time and lower construction costs. Undrained triaxial tests were performed on four clays representing a range of material properties in an effort to reach the maximum dry density, which provides the highest bearing capacity. Each clay was compacted at optimum moisture content, as well as wet and dry of optimum. Measurements were taken using the GeoGauge and shear wave velocities. An empirical approach was used to estimate the effect of a density gradient on soil suction. A relationship between the normal stress and matric suction produced a strong trend when plotted against a function of stiffness and the void ratio, which represents a density gradient. Another relationship between the GeoGauge and shear wave stiffness measurements was found, but no relationship with the material properties of the samples was observed, indicating that more in-depth research is needed to find a stronger relationship.