Browsing by Author "Chapman, Martin C."
Now showing 1 - 20 of 51
Results Per Page
Sort Options
- Analysis of Seismic Signatures Generated from Controlled Methane and Coal Dust Explosions in an Underground MineMurphy, Michael M. (Virginia Tech, 2008-11-07)Examination of seismic records during the time interval of the Sago Mine disaster in 2006 revealed a small amplitude signal possibly associated with an event in the mine. Although the epicenter of the signature was located in the vicinity where the explosion occurred, it could not be unequivocally attributed to the explosion. More needs to be understood about the seismicity from mine explosions in order to properly interpret critical seismic information. A seismic monitoring system located at NIOSH's Lake Lynn Experimental Mine has monitored nineteen experimental methane and dust based explosions. The objective of the study was to analyze seismic signatures generated by the methane and dust explosions to begin understanding their characteristics at different distances away from the source. The seismic signatures from these different events were analyzed using standard waveform analysis procedures in order to estimate the moment magnitude and radiated seismic energy. The procedures used to analyze the data were conducted using self-produced programs not available with existing commercial software. The signatures of the explosions were found to be extremely complex due a combination of mine geometry and experimental design, both of which could not be controlled for the purposes of the study. Geophones located approximately 600 m (1970 ft) and over from the source collected limited data because of the attenuation of the seismic waves generated by the methane explosion. A combination of the methods used to characterize the seismic signatures allowed for differentiation between experimental designs and the size of the explosion. The factors having the largest impact on the signatures were the mine geometry, size of the methane explosion, construction of the mine seal and location of the mine seal. A relationship was derived to correlate the radiated seismic energy to the size of the explosion. Recommendations were made, based upon the limitations of this study, on methods for better collection of seismic data in the future.
- Analyzing the Effect of Moving Resonance on Seismic Response of Structures using Wavelet TransformsNaga, Pradeep (Virginia Tech, 2011-08-05)Nonlinear structures, when subjected to multiple ground motion records that are scaled to consistent ground motion intensity show significant variation in their response. This effect of ground motion randomness on the variation of structural response is defined as Record-to-Record (RTR) Variability. Ground motion characteristics that contribute to this variability in response includes the variation of signal composition (frequency content) with time (spectral nonstationarity).The phenomenon of moving resonance which occurs when the frequency content of the ground motion shifts in a similar manner as the natural frequencies of the structural response, is likely a contributor to variability. This brings the need to further understand the sources of variability due to moving resonance. The present study was carried out to develop a method to analyze the time-frequency content of a ground motion to assess the occurrence of moving resonance and to quantify its potential in effecting the structural systems. Bilinear elastic and elastoplastic hysteretic behavior was considered. Detailed analysis is done to quantify the effect of moving resonance on structural systems due to 22 far field ground motion records. The wavelet coefficient plots gave very good detail of the characteristics of the ground motions that were not clear from the acceleration time histories and response spectra plots. Instances of moving resonance were found out to be significant. Amplification due to moving resonance was found to be quite large. One instance studied in detail (accelerogram of Northridge earthquake at Beverly Hills) had peak displacement amplified by 6 times compared to the amount of peak displacement expected if the system did not exhibit moving resonance. Based on the analyses results, the characteristics of the ground motion records that don't cause significant moving resonance effect on structural systems were observed. Similarly, the characteristics of the ground motions that do cause moving resonance effect on structural systems were examined.
- Apparent Stress Estimates of Seismic Events Associated with Underground Mining Operations in Sedimentary Rock Units: Implications for Improving Underground Roof Collapse Warning SystemsGrant, Elige Buckhanan (Virginia Tech, 2006-07-07)We analyze recordings of seismic events induced by underground mining operations at the Moonee Colliery, located in Australia, and at the Springfield Pike Limestone Quarry, located in the United States. The data were recorded underground in the vicinity of active mining operations and were initially used by the monitoring systems at the two mines to evaluate the potential for roof failure based on the temporal and spatial distribution of the seismic activity. In an earlier study, the roof failures at the two underground mine locations were found to be the result of two distinct failure processes, both of which were correlated with escalated seismic activity before the roof collapsed. In this study, we reexamine the recordings of these seismic emissions for a further assessment of the state of instability in the roof. We estimate the static seismic moment and radiated seismic energy for each recorded seismic event induced by mining operations at the two underground mine locations. These seismic source parameters are estimated from source spectra have been corrected for the instrument response, propagation effects and bandwidth limitations. The apparent stress, which provides an estimate of the stress drop (or stress release) associated with a seismic event, is then determined from the product between the modulus of rigidity and the ratio between the radiated seismic energy and static seismic moment. The validity of constant stress drop scaling for the seismic events at the two underground mine locations is tested. Estimation of the seismic source parameters indicate that the stress drop of the mining induced events increases over three orders of magnitude of increasing seismic moment (106 N·m ≤ M0 ≤ 109 N·m) and indicate a divergence from constant stress drop scaling. When these results are compiled with the results from seven other independent studies, which analyzed the seismicity associated with a variety of seismogenic environments, this trend is found to span over ten orders of magnitude of seismic moment (106 N·m ≤ M0 ≤ 1016 N·m). The observation that the mining induced events do not conform to constant stress drop scaling may assist in gaining a better understanding of the evolution of the roof failure process. We have found that the stress drop at one of the studied mines appears to increase through time prior to a roof collapse. More data are necessary to test this hypothesis. If this hypothesis is validated, it would have important implications for monitoring roof stability. Incorporation of near-real-time estimates of the stress drop into the existing seismic monitoring protocol may provide improved warning of imminent roof collapse hazards.
- Application of Fatigue Theories to Seismic Compression Estimation and the Evaluation of Liquefaction PotentialLasley, Samuel James (Virginia Tech, 2015-08-21)Earthquake-induced liquefaction of saturated soils and seismic compression of unsaturated soils are major sources of hazard to infrastructure, as attested by the wholesale condemnation of neighborhoods surrounding Christchurch, New Zealand. The hazard continues to grow as cities expand into liquefaction- and seismic compression-susceptible areas hence accurate evaluation of both hazards is essential. The liquefaction evaluation procedure presented herein is based on dissipated energy and an SPT liquefaction/no-liquefaction case history database. It is as easy to implement as existing stress-based simplified procedures. Moreover, by using the dissipated energy of the entire loading time history to represent the demand, the proposed procedure melds the existing stress-based and strain-based liquefaction procedures in to a new, robust method that is capable of evaluating liquefaction susceptibility from both earthquake and non-earthquake sources of ground motion. New relationships for stress reduction coefficient (r_d) and number of equivalent cycles ($n_{eq}$) are also presented herein. The r_d relationship has less bias and uncertainty than other common stress reduction coefficient relationships, and both the $n_{eq}$ and $r_d$ relationships are proposed for use in active tectonic and stable continental regimes. The $n_{eq}$ relationship proposed herein is based on an alternative application of the Palmgren-Miner damage hypothesis, shifting from the existing high-cycle, low-damage fatigue framework to a low-cycle framework more applicable to liquefaction analyses. Seismic compression is the accrual of volumetric strains caused by cyclic loading, and presented herein is a "non-simplified" model to estimate seismic compression. The proposed model is based on a modified version of the Richart-Newmark non-linear cumulative damage hypothesis, and was calibrated from the results of drained cyclic simple shear tests. The proposed model can estimate seismic compression from any arbitrary strain time history. It is more accurate than other "non-simplified" seismic compression estimation models over a greater range of volumetric strains and can be used to compute number-of-equivalent shear strain cycles for use in "simplified" seismic compression models, in a manner consistent with seismic compression phenomenon.
- Applications of Double-Difference Tomography for a Deep Hard Rock MineKerr, Jeffrey Bryan (Virginia Tech, 2011-11-16)Seismicity at a deep hard rock mine can be a precursor to ground failure events. Seismicity data can be used in double-difference tomography, which produce tomograms showing velocity distributions in the rock mass that can be used to infer relative stress of the rock mass. The data set used for the double-difference tomography inversion was from Creighton Mine in Sudbury, Ontario, Canada, and consisted of two months of data averaging 150 microseismic events per day. Three separate studies were conducted to evaluate the applications of double-difference tomography on a deep hard rock mine. These studies produced mine scale tomograms, stope scale tomograms of two active stopes, and stope scale tomograms for a cluster of events. TomoDD was used for the tomographic inversion, with other commercial programs used to view the results. All three studies produced results consistent with prior mine knowledge and basic concepts of rock mass stress redistribution. Mine scale tomograms accurately displayed a low velocity where the mined ore body is known to be with adjacent high velocity, stope scale tomograms of the two stopes both correctly demonstrated a low velocity relaxed zone near the stope following a production blast, and stope scale tomograms of an event cluster displayed consistency in results for two clusters in periods before, during, and after each cluster. The three studies show that double-difference tomography is a promising tool for observing rock mass stress redistribution that provides a baseline evaluation for the potential uses of the technology in a deep hard rock mine.
- An Assessment of Hypocenter Errors Associated with the Seismic Monitoring of Induced Hydro-fracturing in Hydrocarbon ReservoirsGilliland, Ellen (Virginia Tech, 2009-10-14)Expanding the standard, single-well recording geometry used to monitor seismicity during hydro-fracture treatments could provide more accurate hypocenter locations and seismic velocities, improving general reservoir characterization. However, for the real, two-well data set obtained for this project, only S-wave picks were available, and testing resulted in anomalous hypocenter location behavior. This study uses a hypocenter location algorithm and both real and synthetic data sets to investigate how the accuracy of the velocity model, starting hypocenter location, recording geometry, and arrival-time picking error affect final hypocenter locations. Hypocenter locations improved using a velocity model that closely matched the observed sonic log rather than a smoothed version of this model. The starting hypocenter location did not affect the final location solution if both starting and final locations were between the wells. Two solutions were possible when the true solution was not directly between the wells. Adding realistic random picking errors to synthetic data closely modeled the dispersed hypocenter error pattern observed in the real data results. Adding data from a third well to synthetic tests dramatically reduced location error and removed horizontal geometric bias observed in the two-well case. Seismic event data recorded during hydro-fracture treatments could potentially be used for three-dimensional joint hypocenter-velocity tomography. This would require observation wells close enough to earthquakes to record P- and S-wave arrivals or wells at orientations sufficient to properly triangulate hypocenter locations. Simulating results with synthetic tests before drilling could optimize survey design to collect data more effectively and make analysis more useful.
- The Automation of Numerical Models of Coseismic TsunamisWiersma, Codi Allen (Virginia Tech, 2019-08-26)The use of tsunami models for applications of 'now-casting', which is the prediction of the present and near future behavior, has limited exploration, and could potentially be of significant usefulness. Tsunamis are most often caused by earthquakes in subduction zones, which generates coupled uplift and subsidence, and displaces the water column. The behavior of the fault failure is difficult to describe in the short term, often requiring seismic waveform inversion, which takes a length of time on the order of weeks to months to properly model, and is much too late for any use in a now-casting sense. To expedite this length of time, a series of source models are created with variable fault geometry behaviors, using fault parameters from Northern Oceanic and Atmospheric Administration's Short-term Inundation and Forecasting of Tsunamis (SIFT) database, in order to model a series of potential tsunami behaviors using the numerical modelling package, GeoClaw. The implementation of modeling could identify areas of interest for further study that are sensitive to fault failure geometry. Initial results show that by varying the geometry of sub-faults of a given earthquake, the resulting tsunami models behave fairly differently with different wave dispersion behavior, both in pattern and magnitude. While there are shortcomings of the potential geometries the code created in this study, and there are significant improvements that can be made, this study provides a good starting point into now-casting of tsunami models, with future iterations likely involving statistical probability in the fault failure geometries.
- Body and surface wave ambient noise seismic interferometry across the Salton Sea Geothermal Field, CaliforniaSabey, Lindsay Erin (Virginia Tech, 2015-01-13)Virtual source gathers were generated using the principles of seismic interferometry from 135 hours of ambient noise recorded during a controlled-source survey across the Salton Sea Geothermal Field in southern California. The non-uniform nature of the noise sources violated a primary assumption of the method and generated artifacts in the data. The artifacts generated by the high-energy impulsive sources (e.g. earthquakes, shots) were removable using traditional methods of amplitude normalization prior to cross-correlation. The continuous source artifacts generated by the geothermal wells and highways required an unconventional approach of utilizing only normalized impulsive sources to successfully reduce the artifacts. Virtual source gathers were produced successfully that contained strong surface waves at 0.4-2.5 Hz, an order of magnitude below the corner frequency of the geophones, and modest body waves at 22-30 Hz, which are generally more difficult to obtain due to the need for many large, well-distributed subsurface sources. The virtual source gathers compare well to nearby explosive shots and are more densely spaced, but have a much lower signal-to-noise ratio. Analysis of the surface waves was complicated by strong higher-order modes. Spectral analysis of virtual source gathers required utilization of the geothermal plant energy, which produced usable signal at offsets required for mode separation. The virtual source dispersion curve compared well to a dispersion curve from a nearby explosive shot. P-waves were observed on the virtual source gathers. Creation of a low-quality multichannel reflection stack revealed two weak reflectors in the upper 2 km.
- Comparison of 4.5 Hz Geophones and a Broadband Seismometer in a Real Field DeploymentRasmussen, Tyler Wyatt (Virginia Tech, 2019-06-18)An analysis of waveforms, power spectral density and array responses was performed using geophones and broadband seismometers, co-deployed as part of a geologically motivated study. Broadband seismometers record excellent waveforms but, due to cost and deployment effort, wavefields are usually spatially aliased above ~0.1 Hz. Industry rapidly deploys many thousands of inexpensive, passive geophones to record full, unaliased seismic wavefields; however, waveform quality is limited below the instrument's natural frequency of ≥2 Hz. In 2012, coincident passive and controlled-source seismic surveys were deployed to investigate tectonics in Idaho and Oregon. Broadband stations were deployed at quiet sites every 15 km, taking experienced professionals >1 person-days per station. Fifty 4.5 Hz geophones and "Texan" seismographs at 200-m spacing were deployed per person-day by inexperienced students. Geophone data were continuously recorded for 3 nights and 1 day, while broadband seismometers were deployed for ~2 years. The spectral and array responses of these real deployments were compared. For a M7.7 teleseismic event, the broadband seismometer and geophone recorded nearly identical waveforms down to <0.03 Hz (32 s) and matching power spectral density down to 0.02 Hz (50 s). For quiet ambient noise, the waveforms strongly correlate down to <0.25 Hz (4 s) and the power spectral density match to the low-frequency side of the microseismic peak at ~0.15 Hz (~7 s). By deploying a much larger number of geophones, waveforms can be stacked to reduce instrument self-noise and beamforming can be used to identify wavefield azimuth and apparent velocity. Geophones can be an effective tool in ambient noise seismology down to ~7 seconds and can be used to record large seismic events effectively down to tens of seconds, well below the natural frequency of the instruments. A well-designed deployment of broadbands and geophones can enable full wavefield studies from long period to short period. Scientific and societal applications that could benefit from the improved unaliased wavefield bandwidth include local to regional seismicity, strong ground motion, magma migration, nuclear source discrimination, and crustal studies.
- Continental Tectonics from Dense Array Seismic Imaging: Intraplate Seismicity in Virginia and a Steep Cratonic Margin in IdahoDavenport, Kathy (Virginia Tech, 2016-09-21)Dense array seismic techniques can be applied to multiple types of seismic data to understand regional tectonic processes via analysis of crustal velocity structure, imaging reflection surfaces, and calculating high-resolution hypocenter locations. The two regions presented here include an intraplate seismogenic fault zone in Virginia and a steep cratonic margin in eastern Oregon and Idaho. The intraplate seismicity study in Virginia consisted of using 201 short-period vertical-component seismographs, which recorded events as low as magnitude -2 during a period of 12 days. Dense array analysis revealed almost no variation in the seismic velocity within the hypocentral zone, indicating that the aftershock zone is confined to a single crystalline-rock terrane. The 1-2 km wide cloud of hypocenters is characterized by a 29° strike and 53° dip consistent with the focal mechanism of the main shock. A 5° bend along strike and a shallower dip angle below 6 km points toward a more complex concave shaped fault zone. The seismic study in Idaho and Oregon was centered on the inversion of controlled-source wide-angle reflection and refraction seismic P- and S-wave traveltimes to determine a seismic velocity model of the crust beneath this part of the U.S. Cordillera. We imaged a narrow, steep velocity boundary within the crust that juxtaposes the Blue Mountains accreted terranes and the North American craton at the western Idaho shear zone. We found a 7 km offset in Moho depth, separating crust with different seismic velocities and Poisson's ratios. The crust beneath the Blue Mountains terranes is consistent with an intermediate lithology dominated by diorite. In the lower crust there is evidence of magmatic underplating which is consistent with the location of the feeder system of the Columbia River Basalts. The cratonic crust east of the WISZ is thicker and characterized by a felsic composition dominated by granite through most of the crust, with an intermediate composition layer in the lower crust. This sharp lithologic and rheologic boundary strongly influenced subsequent deformation and magmatic events in the region.
- Crustal Structure in a Mesozoic Extensional Terrane: The South Georgia Rift and the Epicentral Area of the 1886 Charleston, South Carolina, EarthquakeBuckner, Jesse Conard (Virginia Tech, 2011-01-28)On August 31, 1886 a large scale earthquake occurred in Summerville, S.C. causing severe damage in the coastal city of Charleston. Although intensive geological and geophysical studies have been conducted in the area, uncertainty remains about the details of the event. Recently evidence from seismic reflection profiles have shed light on the tectonic environment of the area. The epicentral area of the 1886 event lies within the South Georgia Rift, a Mesozoic rift terrane. Previous studies have revealed clues to the geologic structure and evolution of this feature. SEISDATA4 is the largest seismic reflection profile recorded in the area. By re-processing the line, information about the tectonic structure of the area was revealed. The early Mesozoic extensional basin that hosted the 1886 earthquake and is host to the modern seismicity recorded in the area, extends several kilometers to the south and west of Charleston, along SEISDATA4. Cenozoic and Mesozoic faults were resolved within the basin and along its northwestern boundary that is distinguished by a strong gradient in the magnetic field. However, the question as to which fault was responsible for the rupture of 1886 still remains. The refraction analysis provides better resolution of the lithology in Lower Mesozoic section. The termination of the strong reflection at the base of the Atlantic Coastal Plain occurs in a section of the profile that shows major disruption of the underlying reflections, and suggests that the termination of a lower Mesozoic basalt flow responsible for the reflection may be related to tectonic deformation.
- Determination of Heterogeneity by High-Resolution Seismic Reservoir Characterization in the Heavy Oil Temblor Reservoir of Coalinga Field, CaliforniaMahapatra, Sailendra Nath (Virginia Tech, 2005-08-12)The research focuses on analysis and subsurface imaging of siliciclastics rocks on steam-affected 3D poststack seismic data, merged from different vintages, from the Temblor Formation in the Coalinga heavy oil reservoir in the San Joaquin basin, California. The objective was identification, delineation, and demarcation of reservoir heterogeneities by seismostratigraphic and seismogeomorphic analysis. The proximity of the San Andreas Transforms greatly controlled basin evolution and caused substantial reservoir heterogeneity by changing the depositional environment from shallow marine to near-shore fluvial. Moreover, two unconformities dissect the reservoir interval. The seismic dataset exhibits erratic, distorted reflection strengths and amplitudes caused by steam-injection-aided production. A petrophysical analysis based on Gassmann fluid substitution suggests a 27% P-wave velocity decrease in steam-saturated intervals. Seismic to well log ties were problematic and vexing due to the resulting statics, wavelet changes, and line mismatches. Mapping and flattening on a deeper horizon, however, allowed mapping of the internal unconformities and well ties which were crucial for seismostratigraphic sequence identification. Visualization of seismic attributes brought out stratification patterns and two distinct, laterally and vertically extensive, porous, and interconnected facies tracts interpreted as incised valley fills and tidal-to-subtidal deposits as evidenced by bright, steam related amplitudes. Seismic attribute analysis, Geobody Visualization and Interpretation, and structure and isochron maps brought out two prominent channel-systems, recut and restacked in the central part of the area. These deposits were identified on seismic data and correlated to high-gamma coarsening-upward sands on logs and cores. The deeper one, shifting towards SSE with depth, lies between the Base Temblor and Buttonbed unconformities both in the southwestern and northwestern parts of the study area and is scattered in the western-central portion. The shallower one originates in southwestern corner below the Top Temblor unconformity shifts towards ESE-SE with depth, and runs nearly parallel to the Top Temblor unconformity. It cuts across the Valv unconformity in central part creating a channel incision, and follows the Buttonbed unconformity towards the north. The investigation segmented the reservoir into channels, non-channel bearing, and unconformity-bounded subunits which will allow the operator to improve steam injection and optimize placement of oil producing infill wells.
- Disaggregated Seismic Hazard and the Elastic Input Energy Spectrum: An Approach to Design Earthquake SelectionChapman, Martin C. (Virginia Tech, 1998-06-25)The design earthquake selection problem is fundamentally probabilistic. Disaggregation of a probabilistic model of the seismic hazard offers a rational and objective approach that can identify the most likely earthquake scenario(s) contributing to hazard. An ensemble of time series can be selected on the basis of the modal earthquakes derived from the disaggregation. This gives a useful time-domain realization of the seismic hazard, to the extent that a single motion parameter captures the important time-domain characteristics. A possible limitation to this approach arises because most currently available motion prediction models for peak ground motion or oscillator response are essentially independent of duration, and modal events derived using the peak motions for the analysis may not represent the optimal characterization of the hazard. The elastic input energy spectrum is an alternative to the elastic response spectrum for these types of analyses. The input energy combines the elements of amplitude and duration into a single parameter description of the ground motion that can be readily incorporated into standard probabilistic seismic hazard analysis methodology. This use of the elastic input energy spectrum is examined. Regression analysis is performed using strong motion data from Western North America and consistent data processing procedures for both the absolute input energy equivalent velocity, (Vea), and the elastic pseudo-relative velocity response (PSV) in the frequency range 0.5 to 10 Hz. The results show that the two parameters can be successfully fit with identical functional forms. The dependence of Vea and PSV upon (NEHRP) site classification is virtually identical. The variance of Vea is uniformly less than that of PSV, indicating that Vea can be predicted with slightly less uncertainty as a function of magnitude, distance and site classification. The effects of site class are important at frequencies less than a few Hertz. The regression modeling does not resolve significant effects due to site class at frequencies greater than approximately 5 Hz. Disaggregation of general seismic hazard models using Vea indicates that the modal magnitudes for the higher frequency oscillators tend to be larger, and vary less with oscillator frequency, than those derived using PSV. Insofar as the elastic input energy may be a better parameter for quantifying the damage potential of ground motion, its use in probabilistic seismic hazard analysis could provide an improved means for selecting earthquake scenarios and establishing design earthquakes for many types of engineering analyses.
- Dynamics of the North American Plate: Numerical Development, Mantle Flow Modeling, and Receiver Function AnalysisLiu, Shangxin (Virginia Tech, 2021-06-15)With only approximately one quarter of plate margins composed of subduction zones, North American plate is an unique continental plate featured with a western active continental margin atop widespread slow seismic velocity anomalies in the asthenosphere, an eastern passive continental margin covering several localized regions of slow seismic velocity, and a strong central cratonic root (Laurentia). The coexistence of the prominent thermal and compositional structures beneath the North American plate complicates the construction of numerical models needed to investigate the dynamics of the whole plate. Recently, a new generation mantle convection code, ASPECT (Advanced Solver for Problems in Earth ConvecTion) equipped with fully adaptive mesh refinement (AMR) technology opens up the potential to build a multi-scale global mantle flow model with a local high-resolution focus beneath the North America plate. Given the immature state of this new code for mantle flow modeling in 3-D spherical shell geometry at the beginning of my doctoral study, I first developed a new geoid algorithm for the 3-D spherical AMR numerical modeling based on ASPECT. Then I systematically benchmarked the velocity, dynamic topography, and geoid solutions from ASPECT through analytical kernel approach in the uniform mesh. I further verified the accuracy of the AMR mantle flow computation in the 3-D spherical shell geometry. Based on the improved ASPECT code, I construct global mantle flow models to investigate the driving forces for the North American plate motion. I focus on the comparison between the effects of near-field slabs (Aleutian, central American, and Caribbean slabs) and far-field slabs (primarily those around western Pacific subduction margins) and find that the far-field slabs provide the dominant driving forces for the North American plate. I further identified that interpreting the extremely slow seismic anomalies associated with the partial melt in the uppermost mantle around southwestern U.S. as purely thermal in origin results in considerably excessive resistance to North American plate motion. My numerical experiments prove that a significantly reduced velocity-to-density scaling (0.05 or smaller in our models) from the original thermal scaling coefficients (0.25 in our models) for these negative seismic shear-velocity anomalies must be incorporated into the construction of the buoyancy field to predict North American plate motion. I also examine the role of the lower mantle buoyancy including the ancient descending Kula-Farallon plates and the active upwelling below the Pacific margin of North American plate. Lower mantle buoyancy primarily affects the amplitudes, as opposed to the patterns of both North American and global plate motions. Another part of this dissertation reports the receiver function analysis along a recent dense seismic array across the eastern U.S from the western border of Ohio to the Atlantic coast of Virginia. 3D stacking yields shallowing trends of 410-km and 660-km discontinuities and thinning transition zone thickness from the inland to the coast. These results are hard to reconcile with any of the three existing hypotheses regarding the vertical mantle flow patterns beneath the eastern U.S., including edge-driven convection excited by the craton edge, hydrous upwelling from the dehydration of the deep Farallon slab, and the sinking of the delaminated or dripped mantle lithospheric block below the central West Virginia/Virginia border. A hydro-thermal upwelling beneath the eastern U.S. coastal plain due to hydrated transition zone and the neighboring passive hot upwelling induced by the descending Farallon slab in the lower mantle is consistent with the results from 3D stacking. The hydro-thermal upwelling hypothesis is also able to reconcile the shallower tectonic processes and deeper mantle dynamics below the eastern U.S. through its dehydration melting atop 410-km discontinuity. Overall, this dissertation documents the technical details on the improvements of the ASPECT code in mantle flow modeling and provides new insights into the dynamics and evolution of the North American continent.
- Evidence of Hydro-Seismicity in the Tennessee Seismic ZoneKing, Gregory Lester (Virginia Tech, 2023-11-17)Reservoir level and discharge data were collected for three Tennessee Valley Authority (TVA) dams (Chickamauga, Watts Bar and Douglas) on the Tennessee river in the Eastern Tennessee Seismic Zone (ETSZ). Earthquake catalog data was also collected for the ETSZ from 1980-2018. Well levels from 2 U.S. Geological Survey groundwater monitoring wells were also collected. Rainfall data were collected for the ETSZ. Reservoir, well and rainfall data were examined for correlation with earthquake occurrence rates. High Pearson correlation and anti-correlation coefficients (.7-.9) were obtained for the reservoir level and dam discharge volume rates vs monthly earthquake counts. Pearson correlation coefficients for rainfall and well level vs. monthly earthquake counts were small (.0-.2). A large difference in earthquake rates between the months of February (low rate) and April (high rate) over the 39 year study period was observed. The difference in earthquakes rates for February and April are statistically significant at the 95% confidence interval using the two sample Poisson rate test. The high correlations for the reservoir level and flow discharge vs. earthquake counts provides strong evidence of hydro-seismicity occurring in the ETSZ from seasonal fluctuations in reservoir level and discharge flow.
- An Examination of Site Response in Columbia, South Carolina: Sensitivity of Site Response to "Rock" Input Motion and the Utility of Vs(30)Lester, Alanna Paige (Virginia Tech, 2005-06-01)This study examines the sensitivity of calculated site response in connection with alternative assumptions regarding input motions and procedures prescribed in the IBC 2000 building code, particularly the use of average shear wave velocity in the upper 30 meters as an index for engineering design response spectra. Site specific subsurface models are developed for four sites in and near Columbia, South Carolina using shear wave velocity measurements from cone penetrometer tests. The four sites are underlain by thin coastal plain sedimentary deposits, overlying high velocity Paleozoic crystalline rock. An equivalent-linear algorithm is used to estimate site response for vertically incident shear waves in a horizontally layered Earth model. Non-linear mechanical behavior of the soils is analyzed using previously published strain-dependent shear modulus and damping degradation models. Two models for material beneath the investigated near-surface deposits are used: B-C outcrop conditions and hard rock outcrop conditions. The rock outcrop model is considered a geologically realistic model where a velocity gradient, representing a transition zone of partially weathered rock and fractured rock, overlies a rock half-space. Synthetic earthquake input motions are generated using the deaggregations from the 2002 National Seismic Hazard Maps, representing the characteristic Charleston source. The U. S. Geological Survey (2002) uniform hazard spectra are used to develop 2% in 50 year probability of exceedance input ground motions for both B-C boundary and hard rock outcrop conditions. An initial analysis was made for all sites using an 8 meter thick velocity gradient for the rock input model. Sensitivity of the models to uncertainty of the weathered zone thickness was assessed by randomizing the thickness of the velocity gradient. The effect of the velocity gradient representing the weathered rock zone increases site response at high frequencies. Both models (B-C outcrop conditions and rock outcrop conditions) are compared with the International Building Code (IBC 2000) maximum credible earthquake spectra. The results for both models exceed the IBC 2000 spectra at some frequencies, between 3 and 10 Hz at all four sites. However, site 2, which classifies as a C site and is therefore assumed to be the most competent of the four sites according to IBC 2000 design procedures, has the highest calculated spectral acceleration of the four sites analyzed. Site 2 has the highest response because a low velocity zone exists at the bottom of the geotechnical profile in immediate contact with the higher velocity rock material, producing a very large impedance contrast. An important shortcoming of the IBC 2000 building code results from the fact that it does not account for cases in which there is a strong rock-soil velocity contrast at depth less than 30 meters. It is suggested that other site-specific parameters, specifically, depth to bedrock and near-surface impedance ratio, should be included in the IBC design procedures.
- Geophysical Imaging of Earth Processes: Electromagnetic Induction in Rough Geologic Media, and Back-Projection Imaging of Earthquake AftershocksBeskardes, Gungor Didem (Virginia Tech, 2017-06-04)This dissertation focuses on two different types of responses of Earth; that is, seismic and electromagnetic, and aims to better understand Earth processes at a wider range of scales than those conventional approaches offer. Electromagnetic responses resulting from the subsurface diffusion of applied electromagnetic fields through heterogeneous geoelectrical structures are utilized to characterize the underlying geology. Geology exhibits multiscale hierarchical structure which brought about by almost all geological processes operating across multiple length scales and the relationship between multiscale electrical properties of underlying geology and the observed electromagnetic response has not yet been fully understood. To quantify this relationship, the electromagnetic responses of textured and spatially correlated, stochastic geologic media are herein presented. The modelling results demonstrate that the resulting electromagnetic responses present a power law distribution, rather than a smooth response polluted with random, incoherent noise as commonly assumed; moreover, they are examples of fractional Brownian motion. Furthermore, the results indicate that the fractal behavior of electromagnetic responses is correlated with the degree of the spatial correlation, the contrasts in ground electrical conductivity, and the preferred orientation of small-scale heterogeneity. In addition, these inferences are also supported by the observed electromagnetic responses from a fault zone comprising different lithological units and varying wavelengths of geologic heterogeneity. Seismic signals generated by aftershocks are generally recorded by local aftershock networks consisted of insufficient number of stations which result in strongly spatially-aliased aftershock data. This limits aftershock detections and locations at smaller magnitudes. Following the 23 August 2011 Mineral, Virginia earthquake, to drastically reduce spatial aliasing, a temporary dense array (AIDA) consisting of ~200 stations at 200-400 m spacing was deployed near the epicenter to record the 12 days of the aftershocks. The backprojection imaging method is applied to the entire AIDA dataset to detect and locate aftershocks. The method takes advantage of staking of many seismograms and improves the signal-to-noise ratio for detection. The catalog obtained from the co-deployed, unusually large temporal traditional network of 36 stations enabled a quantitative comparison. The aftershock catalog derived from the dense AIDA array and the backprojection indicates event detection an order of magnitude smaller including events as small as M–1.8. The catalog is complete to magnitude –1.0 while the traditional network catalog was complete to M–0.27 for the same time period. The AIDA backprojection catalog indicate the same major patterns of seismicity in the epicentral region, but additional details are revealed indicating a more complex fault zone and a new shallow cluster. The b-value or the temporal decay constant were not changed by inclusion of the small events; however, they are different for two completeness periods and are different at shallow depth than greater depth.
- Global Structure of the Mantle Transition Zone Discontinuities and Site Response Effects in the Atlantic and Gulf Coastal PlainGuo, Zhen (Virginia Tech, 2019-09-03)This thesis focuses on two different topics in seismology: imaging the global structures of the mantle transition zone discontinuities and studying the site response effects in the Atlantic and Gulf Coastal Plain. Global structures of the mantle transition zone discontinuities provide important constraints on thermal structures and dynamic processes in the mid mantle. In this thesis, global topographic structures of the 410- and 660-km discontinuities are obtained from finite-frequency tomography of SS precursors. The finite-frequency sensitivities of SS waves and precursors are calculated based on a single-scattering (Born) approximation and can be used for data selection. The new global models show a number of smaller-scale features that were absent in back-projection models. Good correlation between the mantle transition zone thickness and wave speed variations suggests dominantly thermal origins for the lateral variations in the transition zone. The high-resolution global models of the 410- and 660-km discontinuities in this thesis show strong positive correlation beneath western North America and eastern Asia subduction zones with both discontinuities occurring at greater depths. Wavespeed and anisotropy models support vertical variations in thermal structure in the mid mantle, suggesting return flows from the lower mantle occur predominantly in the vicinity of stagnant slabs and the region overlying the stagnant slabs. In oceanic regions, the two discontinuities show a weak anti-correlation, indicating the existence of a secondary global far-field return flow. The Atlantic and Gulf Coastal Plain is covered by extensive Cretaceous and Cenozoic marine sediments. In this thesis, the site response effects of sediments in the Coastal Plain region relative to the reference condition outside that region are investigated using Lg and coda spectral ratios. The high-frequency attenuation factors (kappa) in the Coastal Plain are strongly correlated with the sediment thickness. At frequencies between 0.1-2.86 Hz, the Lg spectral ratio amplitudes are modeled as functions of frequency and thickness of the sediments in the Coastal Plain. Analysis of the residuals from the stochastic ground motion prediction method suggests that incorporating the site response effects as functions of sediment thickness may improve ground motion prediction models for the Coastal Plain region.
- High density oilfield wastewater disposal causes deeper, stronger, and more persistent earthquakesPollyea, Ryan M.; Chapman, Martin C.; Jayne, Richard S.; Wu, Hao (Nature, 2019)Oilfield wastewater disposal causes fluid pressure transients that induce earthquakes. Here we show that, in addition to pressure transients related to pumping, there are pressure transients caused by density differences between the wastewater and host rock fluids. In northern Oklahoma, this effect caused earthquakes to migrate downward at ~0.5 km per year during a period of high-rate injections. Following substantial injection rate reductions, the downward earthquake migration rate slowed to ~0.1 km per year. Our model of this scenario shows that the density-driven pressure front migrates downward at comparable rates. This effect may locally increase fluid pressure below injection wells for 10+ years after substantial injection rate reductions. We also show that in north-central Oklahoma the relative proportion of high-magnitude earthquakes increases at 8+ km depth. Thus, our study implies that, following injection rate reductions, the frequency of high-magnitude earthquakes may decay more slowly than the overall earthquake rate.
- Hybrid Steel FramesAtlayan, Ozgur (Virginia Tech, 2013-04-22)The buildings that are designed according to the building codes generally perform well at severe performance objectives (like life safety) under high earthquake hazard levels. However, the building performance at low performance objectives (like immediate occupancy) under low earthquake hazards is uncertain. The motivation of this research is to modify the design and detailing rules to make the traditional systems perform better at multi-level hazards. This research introduces two new structural steel systems: hybrid Buckling Restrained Braced Frames (BRBF) and hybrid steel Moment Frames (MF). The "hybrid" term for the BRBF system comes from the use of different steel material including carbon steel (A36), high-performance steel (HPS) and low yield point (LYP) steel. The hybridity of the moment frames is related to the sequence in the plastification of the system which is provided by using weaker and stronger girder sections. Alternative moment frame connections incorporating the use of LYP steel plates are also investigated. The hybrid BRBF approach was evaluated on seventeen regular (standard) frames with different story heights, seismic design categories and building plans. By varying the steel areas and materials in the BRB cores, three hybrid BRBFs were developed for each regular (standard) frame and their behavior was compared against each other through pushover and incremental dynamic analyses. The benefits of the hybridity were presented using different damage measures such as story accelerations, interstory drifts, and residual displacements. Collapse performance evaluation was also provided. The performance of hybrid moment frames was investigated on a design space including forty-two moment frame archetypes. Two different hybrid combinations were implemented in the designs with different column sections and different strong column-weak beam (SC/WB) ratios. The efficiency of the hybrid moment frame in which only the girder sizes were changed to control the plastification was compared with regular moment frame designs with higher SC/WB ratios. As side studies, the effect of shallow and deep column sections and SC/WB ratios on the moment frame behavior were also investigated. In order to provide adequate ductility in the reduced capacity bays with special detailing, alternative hybrid moment frame connections adapting the use of low strength steel were also studied.
- «
- 1 (current)
- 2
- 3
- »