Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic sources”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Elastic-wave sensitivity-guided adaptive seismic survey design for cost-effective monitoring of geological carbon storage

Effective seismic monitoring is essential for verifying CO₂ containment, detecting potential leakage, and optimizing operational decisions in geologic carbon storage. Here, this study presents a time-adaptive, elastic-wave sensitivity-guided framework for designing cost-effective seismic monitoring layouts for tracking CO₂ plume migration. The method is based on elastic-wave sensitivity analysis, which quantifies how variations in subsurface properties impact seismic wavefields. Two complementary design strategies are developed: one based on selecting a fixed number of seismic sources (Method A), and the other based on selecting source–receiver pairs contributing to a fixed fraction of cumulative elastic-wave sensitivity energy (Method B). The optimization workflow to identify source–receiver configurations with the highest detection potential is demonstrated using a hypothetical GCS scenario at the Kimberlina site in California using simulations of elastic-wave sensitivity data at multiple post-injection timesteps. Results show that both strategies adapt to evolving plume geometries and wavefield sensitivities, with Method B offering broader spatial coverage and Method A ensuring simpler deployment. This framework enables site-specific, cost-effective, and risk-informed seismic survey designs, enhancing the ability to monitor CO₂ migration over time in evolving geological environments

58 GEOSCIENCES↗
DOE OSTI2023

High-resolution lithospheric structure of continental China from joint inversion of surface wave and gravity data

The lithospheric structure of continental China has been previously determined by seismic travel time tomography, surface wave tomography, and joint inversion of body wave and surface wave data. However, due to the inherent limitations of seismic data, the lithospheric structure of continental China is still not well resolved in the shallow part and in some regions where the station coverage is relatively sparse. In this study, we aim at improving the lithospheric structure by joint inversion of seismic surface wave data and satellite gravity data to take advantage of the uniform distribution and complementary strength of the gravity data. Here, the empirical relationship between velocity and density is used as a bridge for joint inversion of surface wave and gravity data. The joint inversion shear-wave velocity Vs (density) model, named as USTClitho1.0g, can fit both surface wave and gravity data well. This high-resolution Vs model can better fit the active airgun source seismic arrival times and better delineate some features in continental China, such as the velocity contrast across the north–south gravity lineament (NSGL), the lithosphere thinning in eastern China, the crustal footprint of the Hainan mantle plume, the likely magma chamber beneath volcanos in northeast China, the middle-lower crust low velocity layer beneath the Tibetan plateau, and the tearing of subducted Indian mantle lithosphere. Our joint inversion Vs model can provide a reference model for geosciences in continental China and surrounding areas.

58 GEOSCIENCES↗
DOE OSTI2023

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗
NASA NTRS1977

Faulting parameters derived from computer simulation of earthquakes

Seismic source parameters, average displacement, rupture length, and strain energy release are investigated by computer simulation using a coupled massive block model of the sliding along an active fault. Average displacements and energy release vary considerably with the degree of heterogeneity in the friction and elastic parameters used in the model. Strain energy release is determined primarily by the product of dynamic friction, rupture length, and average displacement. Interrelationships among the faulting parameters are consistent with theoretical arguments and experimental data. The variation in the frequency of occurrence of simulation events with strain energy release is different from the variation in the frequency of naturally occurring events with seismic energy.

Cohen, S. C.↗
NASA NTRS1993

Geologic controls on the West Antarctic Ice Sheet

The stability of the West Antarctic Ice Sheet is intertwined with its geologic history. The sub-ice geology and the possibility of active rifting and associated elevated heat-flux and volcanism might be determining factors in ice-sheet behavior. Seismic monitoring of natural events at the base of Ice Stream C reveals the presence of a young sedimentary basin beneath the ice stream. The sediments are presumed to be of glacio-marine origin, similar to those beneath Ice Stream B and in the Ross Sea. The young sediments are approximately 1/2 km thick at UpC camp, but thin abruptly southward to 100 m or less. We hypothesize the presence of a fault with a throw of 400 m to account for this (though we have not directly detected the fault), rather than invoking unrealistic basement dips. To extend these studies to critical inland regions, we suggest an expanded explosive-source seismic survey of the Byrd Subglacial Basin to determine the extent and character of the hypothesized rift basin. High-resolution seismic monitoring will detect layering in the sedimentary column, as well as possibly imaging faults directly.

Anandakrishnan, Sridhar↗
NASA NTRS1995

GPS detection of ionospheric perturbations following the January 17, 1994, northridge earthquake

Sources such as atmospheric or buried explosions and shallow earthquakes producing strong vertical ground displacements produce pressure waves that propagate at infrasonic speeds in the atmosphere. At ionospheric altitudes low frequency acoustic waves are coupled to ionispheric gravity waves and induce variations in the ionoispheric electron density. Global Positioning System (GPS) data recorded in Southern California were used to compute ionospheric electron content time series for several days preceding and following the January 17, 1994, M(sub w) = 6.7 Northridge earthquake. An anomalous signal beginning several minutes after the earthquake with time delays that increase with distance from the epicenter was observed. The signal frequency and phase velocity are consistent with results from numerical models of atmospheric-ionospheric acoustic-gravity waves excited by seismic sources as well as previous electromagnetic sounding results. It is believed that these perturbations are caused by the ionospheric response to the strong ground displacement associated with the Northridge earthquake.

Calais, Eric↗
DOE OSTI2021

Data Analysis of the 2020 Central Idaho Mainshock-Aftershock Sequence

In an effort to inform the Senior Seismic Hazard Analysis Committee for the Idaho National Laboratory, we provide an improved aftershock catalog related to the March 31, 2020, Mw6.5 Stanley, Idaho earthquake from picks related to a temporary network of two real-time and 15 non-telemetered seismometers within the epicentral area. From the permanent and temporary (XP) real-time network, the USGS cataloged 1,946 aftershocks between April 1, 2020 and October 31, 2020. To improve aftershock location and magnitudes, we manually picked arrival times of P and S waves from off-line stations in the XP temporary network, generated a new crustal velocity model, and independently relocated each event using the HypoDD double-difference earthquake algorithm. We created our new velocity model from existing broadband and active source seismic campaign data that were acquired near the epicentral region prior to the 2020 earthquake. We compare arrival time differences, epicentral locations and depths between aftershocks recorded with the two catalogs. We find the addition of local stations provides tighter aftershock clustering that suggests an improved aftershock locations. To detect lower magnitude events, we employed deep learning. Our method solves common problems associated with detecting many events that have a low signal-to-noise ratio. From the machine learning database, we detected more than 74,000 aftershocks. Based on the number of identified earthquakes and Gutenberg-Richter relationships derived from the USGS catalog, we estimate that we have reduced the completion magnitude for the Stanley earthquake sequence to below M1 using this machine learning approach. We located each aftershock with our new velocity model. Our new velocity model and picks suggests aftershocks occurred mostly at shallower depths than assessed in the USGS catalog. These aftershocks align along two linear trends that suggest the activation of two unnamed primary faults.

58 GEOSCIENCES↗
DOE OSTI2020

Dense Seismic Array Study of a Legacy Underground Nuclear Test at the Nevada National Security Site

The complex postdetonation geologic structures that form after an underground nuclear explosion are difficult to constrain because increased heterogeneity around the damage zone affects seismic waves that propagate through the explosion site. Generally, a vertical rubble-filled structure known as a chimney is formed after an underground nuclear explosion that is composed of debris that falls into the subsurface cavity generated by the explosion. Compared with chimneys that collapse fully, leaving a surface crater, partially collapsed chimneys can have remnant subsurface cavities left in place above collapsed rubble. The 1964 nuclear test HADDOCK, conducted at the Nevada test site (now the Nevada National Security Site), formed a partially collapsed chimney with no surface crater. Understanding the subsurface structure of these features has significant national security applications, such as aiding the study of suspected underground nuclear explosions under a treaty verification. In this study, we investigated the subsurface architecture of the HADDOCK legacy nuclear test using hybrid 2D–3D active source seismic reflection and refraction data. The seismic data were acquired using 275 survey shots from the Seismic Hammer (a 13,000 kg weight drop) and 65 survey shots from a smaller accelerated weight drop, both recorded by ~1000 three-component 5 Hz geophones. First-arrival, P-wave tomographic modeling shows a low-velocity anomaly at ~200 m depth, likely an air-filled cavity caused by partial collapse of the rock column into the temporary postdetonation cavity. A high-velocity anomaly between 20 and 60 m depth represents spall-related compaction of the shallow alluvium. Hints of low velocities are also present near the burial depth (~364 m). The reflection seismic data show a prominent subhorizontal reflector at ~300 m depth, a short-curved reflector at ~200 m, and a high-amplitude reflector at ~50 m depth. Comparisons of the reflection sections to synthetic data and borehole stratigraphy suggest that these features correspond to the alluvium–tuff contact, the partial collapse cavity, and the spalled layer, respectively.

58 GEOSCIENCES↗
DOE OSTI2021

SPiRaL: a multiresolution global tomography model of seismic wave speeds and radial anisotropy variations in the crust and mantle

SUMMARY SPiRaL is a joint global-scale model of wave speeds (P and S) and anisotropy (vertical transverse isotropy, VTI) variations in the crust and mantle. The model is comprised of >2.1 million nodes with five parameters at each node that capture velocity variations for P- and S-waves travelling at arbitrary directions in transversely isotropic media with a vertical symmetry axis (VTI). The crust (including ice, water, sediments and crystalline layers) is directly incorporated into the model. The default node spacing is approximately 2° in the lower mantle and 1° in the crust and upper mantle. The grid is refined with ∼0.25° minimum node spacing in highly sampled regions of the crust and upper mantle throughout North America and Eurasia. The data considered in the construction of SPiRaL includes millions of body wave traveltimes (crustal, regional and teleseismic phases with multiples) and surface wave (Rayleigh and Love) dispersion. A multiresolution inversion approach is employed to capture long-wavelength heterogeneities commonly depicted in global-scale tomography images as well as more localized details that are typically resolved in more focused regional-scale studies. Our previous work has demonstrated that such global-scale models with regional-scale detail can accurately predict both teleseismic and regional body wave traveltimes, which is necessary for more accurate location of small seismic events that may have limited signal at teleseismic distances. SPiRaL was constructed to predict traveltimes for event location and long-period waveform dispersion for seismic source inversion applications in regions without sufficiently tuned models. SPiRaL may also serve as a starting model for full-waveform inversion (FWI) with the goal of fitting waves with periods 10–50 s over multiple broad regions (thousands of kilometres) and potentially the globe. To gain insight to this possibility, we simulated waveforms for a small set of events using SPiRaL and independent waveform-based models for comparison. For the events tested, the performance of the traveltime-based SPiRaL model is shown to be generally on par with regional 3-D waveform-based models in three regions (western United States, Middle East, Korean Peninsula) suggesting SPiRaL may serve as a starting model for FWI over broad regions.

58 GEOSCIENCES↗
DOE OSTI2023

High Spatial Melt Rate Variability Near the Totten Glacier Grounding Zone Explained by New Bathymetry Inversion

Totten Glacier is a fast-moving East Antarctic outlet with the potential for significant future sea-level contributions. We deployed four autonomous phase-sensitive radars on its ice shelf to monitor ice-ocean interactions near its grounding zone and made active source seismic observations to constrain gravity-derived bathymetry models. We observe an asymmetry in basal melting with mean melt rates along the grounding zone differing by up to 20 m/a. Our new bathymetry model reveals that this melt rate asymmetry coincides with an asymmetry in water column thickness and that the low-melting ice-shelf portion is shielded from the main cavity circulation. A 2-year record yields year-to-year melt rate variability of 7–9 m/a with no seasonal cycle. Our results highlight the key role of bathymetry near grounding lines for accurate modeling of ice-shelf melt, and the importance of sustained multi-year monitoring, especially at ice-shelf cavities where the dominant melt rate drivers vary primarily inter-annually.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

In Situ Velocity‐Strain Sensitivity Near the San Jacinto Fault Zone Analyzed Through Train Tremors

Abstract We utilize train tremors as P‐wave seismic sources to investigate velocity‐strain sensitivity near the San Jacinto Fault Zone. A dense nodal array deployed at the Piñon Flat Observatory is used to detect and identify repeating train energy emitted from a railway in the Coachella valley. We construct P‐wave correlation functions across the fault zone and estimate the spatially averaged dt / t versus strain sensitivity to be 6.25 × 10 4 . Through numerical simulations, we explore how the sensitivity decays exponentially with depth. The optimal solution reveals a subsurface sensitivity of 1.2 × 10 5 and a depth decay rate of 0.05 km −1 . This sensitivity aligns with previous findings but is toward the higher end, likely due to the fractured fault‐zone rocks. The depth decay rate, previously unreported, is notably smaller than assumed in empirical models. This raises the necessity of further investigations of this parameter, which is crucial to study stress and velocity variations at seismogenic depth.

Geology↗
DOE OSTI2025

Control Mechanisms for Self‐Sealing in Activated Clay‐Rich Faults Through Controlled Hydraulic Injection Experiment

Abstract In a high‐pressure injection fault activation experiment conducted at the Mont Terri underground research laboratory in Switzerland, the transmissivity of the Opalinus Clay fault significantly increased due to opening and shearing. The fluid injection, spanning a few hours, generated a 10 m radius fault activation patch. Subsequent pressure pulse tests conducted bi‐weekly for a year revealed the gradual return of fault transmissivity to its initial state. The study utilized fluid pressure decay analysis, optical fiber monitoring, continuous active source seismic measurements and borehole displacement sensors for measuring fault displacements. The fault zone exhibited a dilation of approximately 1.4 mm, associated with both normal and tangential movements during activation, resulting in a sudden transmissivity increase from 1 × 10 −12 to 3.2 × 10 −7 m 2 /s. Early post‐activation, transient compaction and the subsequent slow compaction were observed, transitioning to an extension regime. The pressure pulse tests demonstrated a rapid transmissivity drop by more than two orders of magnitude within the first 10 days, followed by a gradual and less pronounced decrease. Plastic shear and compaction dominated the transmissivity evolution until 70 days after injection ended, followed by a period where additional factors, such as clay mineral swelling, influenced the behavior. Extrapolation suggested a sealing process taking at least 50 years after the initial activation. Plain Language Summary A field‐scale fault activation experiment offers valuable insights into the elasto‐plastic processes governing the sealing of shale faults. The experiment reveals a rapid increase in the fault's transmissivity by approximately five orders of magnitude during activation. Subsequent observations show a gradual transmissivity decrease by about three orders of magnitude post‐activation, with slow long‐term plastic shear and compaction of the fault competing against secondary processes, notably clay mineral swelling. All conceptual models employed to interpret these field data converge on the estimation that the fault's return to its initial low transmissivity state would require a minimum of 50 years. Key Points High‐pressure injection fault activation experiment at the Mont Terri underground research laboratory Continuous transmissivity measurements record self‐sealing inside a clay‐rich fault zone Transmissivity undergoes a phase of domination by slow plastic compaction and shearing during the initial post‐activation period, with mineral swelling exerting its influence over the long term

Guglielmi, Yves↗
DOE OSTI2022

Validation of earthquake ground-motion models in southern California, USA, using precariously balanced rocks

Accurate estimates of earthquake ground shaking rely on uncertain ground-motion models derived from limited instrumental recordings of historical earthquakes. A critical issue is that there is currently no method to empirically validate the resultant ground-motion estimates of these models at the timescale of rare, large earthquakes; this lack of validation causes great uncertainty in ground-motion estimates. Here, we address this issue and validate ground-motion estimates for southern California utilizing the unexceeded ground motions recorded by 20 precariously balanced rocks. We used cosmogenic 10Be exposure dating to model the age of the precariously balanced rocks, which ranged from ca. 1 ka to ca. 50 ka, and calculated their probability of toppling at different ground-motion levels. With this rock data, we then validated the earthquake ground motions estimated by the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) seismic-source characterization and the Next Generation Attenuation (NGA)-West2 ground-motion models. We found that no ground-motion model estimated levels of earthquake ground shaking consistent with the observed continued existence of all 20 precariously balanced rocks. The ground-motion model I14 estimated ground-motion levels that were inconsistent with the most rocks; therefore, I14 was invalidated and removed. At a 2475 year mean return period, the removal of this invalid ground-motion model resulted in a 2−7% reduction in the mean and a 10−36% reduction in the 5th−95th fractile uncertainty of the ground-motion estimates. Our findings demonstrate the value of empirical data from precariously balanced rocks as a validation tool for removing invalid ground-motion models and, in turn, reducing the uncertainty in earthquake ground-motion estimates.

Geology↗
DOE OSTI2026

Improved Regional Moment Tensor Inversion for Moderately Large Earthquakes in the Western United States Using a 3D Earth Model Based on Full Waveform Tomography

The nature of seismic sources for moderately large (moment magnitude, M w 5.0–6.5) events are commonly characterized by their moment tensor (MT) solutions and obtained by inversion of regional distance (200–1600 km) long‐period (20–50 s) waveforms. Regional MT estimates are often calculated from average plane‐layered, one‐dimensional (1D) velocity models. However, 1D model calculations can produce misfits in the arrival times and waveform shapes that introduce errors, particularly at longer distances or for shorter periods, which are necessary for analyzing lower magnitude events. Approximate Earth models (e.g., 1D) representing broad areas may be inadequate, particularly in the crust and uppermost mantle of tectonically complex regions. In this study, we show how a three‐dimensional (3D) Earth model obtained from full waveform inversion tomography can improve waveform fits and decrease phase errors. We developed a platform and workflow to perform routine 3D MT inversions and inverted MTs for 25 earthquakes in the western United States and seven nuclear explosions using an average 1D and a recent 3D Earth model, WUS256 (Rodgers et al., 2022). Using the 3D model improves waveform fits (variance reduction and phase time shifts) compared with the 1D model, and the 3D MT solutions are stable across large distances. This study shows that 3D models obtained from full waveform tomography can improve MTs and source characterization especially at far regional distances (>800 km).

Geosciences↗
DOE OSTI2022

NETL Plastic Pipes Project - June 2022 Progress Report

Plastic or composite pipelines have been the bane of the utility locating industry because conduits are neither conductive nor magnetic, which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipeline where the resistive cover allows for adequate penetration of radar energy. However, GPR cannot be used in areas where the soil cover is conductive due to significant clay and/or salt content. This study takes a comprehensive look at near-surface geophysical methods that potentially are useful for locating buried plastic/composite pipelines, either singly or in combination with other geophysical methods. Specifically, this modeling study uses computational numerical methods to forward model the response of GPR, resistivity, seismic, and gravity gradiometry methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft to 4 ft; (3) various degrees in contrast in physical properties (i.e., permittivity, elasticity, electrical resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar frequencies, seismic source frequency, electrode spacing).

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A Catalog of Temporally Localized Systematic Deviations in Global Body Wave Travel-Time Measurements

Accurate measurements of the arrival times of seismic waves are crucial for seismological analyses such as robust locations of earthquakes, characterization of seismic sources, and high-fidelity imaging of the Earth’s interior. However, these travel-time measurements can sometimes be contaminated by timing errors at the stations which record this data. In this study, we apply a classical approach, based on identifying time-dependence in measured body wave arrival times, to identify these timing errors in a dataset on the order of 107 individual measurements. We find timing deviations at a subset of the stations in our dataset and document the temporal location, extent, and severity of these errors, finding errors at 83 stations, and impacting ~100,000 measurements. This catalog of deviations may enable future investigators to obtain a more accurate dataset through the implementation of quality control measures to eliminate the contaminated data we have identified.

58 GEOSCIENCES↗
NASA NTRS1991

Waveform effects of a metastable olivine tongue in subducting slabs

Velocity models of subducting slabs with a kinetically-depressed olivine to beta- and gamma-spinel transition are constructed, and the effect that such structures would have on teleseismic P waveforms are examined using a full-wave finite-difference method. These 2D calculations yielded waveforms at a range of distances in the downdip direction. The slab models included a wedge-shaped, low-velocity metastable olivine tongue (MOTO) to a depth of 670 km, as well as a plausible thermal anomaly; one model further included a 10-km-thick fast layer on the surface of the slab. The principal effect of MOTO is to produce grazing reflections at wide angles off the phase boundary, generating a secondary arrival 0 to 4 seconds after the initial arrival depending on the take-off angle. The amplitude and timing of this feature vary with the lateral location of the seismic source within the slab cross-section.

Vidale, John E.↗
NASA NTRS1999

Using a Genetic Algorithm to Model Broadband Regional Waveforms for Crustal Structure in the Western United States

In this study, we analyze regional seismograms to obtain the crustal structure in the eastern Great Basin and western Colorado plateau. Adopting a for- ward-modeling approach, we develop a genetic algorithm (GA) based parameter search technique to constrain the one-dimensional crustal structure in these regions. The data are broadband three-component seismograms recorded at the 1994-95 IRIS PASSCAL Colorado Plateau to Great Basin experiment (CPGB) stations and supplemented by data from U.S. National Seismic Network (USNSN) stations in Utah and Nevada. We use the southwestern Wyoming mine collapse event (M(sub b) = 5.2) that occurred on 3 February 1995 as the seismic source. We model the regional seismograms using a four-layer crustal model with constant layer parameters. Timing of teleseismic receiver functions at CPGB stations are added as an additional constraint in the modeling. GA allows us to efficiently search the model space. A carefully chosen fitness function and a windowing scheme are added to the algorithm to prevent search stagnation. The technique is tested with synthetic data, both with and without random Gaussian noise added to it. Several separate model searches are carried out to estimate the variability of the model parameters. The average Colorado plateau crustal structure is characterized by a 40-km-thick crust with velocity increases at depths of about 10 and 25 km and a fast lower crust while the Great Basin has approximately 35- km-thick crust and a 2.9-km-thick sedimentary layer.

Bhattacharyya, Joydeep↗