Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic refraction”

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 19 records

Landslide monitoring using seismic refraction tomography – The importance of incorporating topographic variations

Seismic refraction tomography provides images of the elastic properties of subsurface materials in landslide settings. Seismic velocities are sensitive to changes in moisture content, which is a triggering factor in the initiation of many landslides. However, the application of the method to long-term monitoring of landslides is rarely used, given the challenges in undertaking repeat surveys and in handling and minimizing the errors arising from processing time-lapse surveys. This work presents a simple method and workflow for producing a reliable time-series of inverted seismic velocity models. This method is tested using data acquired during a recent, novel, long-term seismic refraction monitoring campaign at an active landslide in the UK. Potential sources of error include those arising from inaccurate and inconsistent determination of first-arrival times, inaccurate receiver positioning, and selection of inappropriate inversion starting models. At our site, a comparative analysis of variations in seismic velocity to real-world variations in topography over time shows that topographic error alone can account for changes in seismic velocity of greater than ±10% in a significant proportion (23%) of the data acquired. The seismic velocity variations arising from real material property changes at the near-surface of the landslide, linked to other sources of environmental data, are demonstrated to be of a similar magnitude. Over the monitoring period we observe subtle variations in the bulk seismic velocity of the sliding layer that are demonstrably related to variations in moisture content. This highlights the need to incorporate accurate topographic information for each time-step in the monitoring time-series. The goal of the proposed workflow is to minimize the sources of potential errors, and to preserve the changes observed by real variations in the subsurface. Following the workflow produces spatially comparable, time-lapse velocity cross-sections formulated from disparate, discretely-acquired datasets. These practical steps aim to aid the use of the seismic refraction tomography method for the long-term monitoring of landslides prone to hydrological destabilization.

Whiteley, J. S.↗

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗

Crustal structure of the Southern Rio Grande rift determined from seismic refraction profiling

As part of a major cooperative seismic experiment, a series of seismic refraction profiles have been recorded in south-central New Mexico with the goal of determining the crustal structure in the southern Rio Grande rift. The data gathered greatly expand the seismic data base in the area, and consist of three interlocking regional profiles: a reversed E-W line across the rift, an unreversed N-S axial line, and an unreversed SW-SE line. The reversed E-W line shows no significant dip along the Moho (32 km thick crust) and a 7.7 km/s Pn velocity. Results from the N-S axial line and the NW-SE line indicate an apparent Pn velocity of 7.95 km/s and significant dip along the Moho with crustal thinning toward the south and southeast. When interpreted together, these data indicate a crustal thinning in the southern rift of 4-6 km with respect to the northern rift and the adjacent Basin and Range province, and establish the regional Pn velocity to be approximately 7.7 km/s. These results suggest that the Rio Grande rift can be identified as a crustal feature separate and distinct from the Basin and Range province.

Sinno, Y. A.↗

A seismic refraction technique used for subsurface investigations at Meteor Crater, Arizona

A seismic refraction technique for interpreting the subsurface shape and velocity distribution of an anomalous surface feature such as an impact crater is described. The method requires the existence of a relatively deep refracting horizon and combines data obtained from both standard shallow refraction spreads and distant offset shots by using the deep refractor as a source of initial arrivals. Results obtained from applying the technique to Meteor crater generally agree with the known structure of the crater deduced by other investigators and provide new data on an extensive fractured zone surrounding the crater. The breccia lens is computed to extend roughly 190 m below the crater floor, about 30 m less than the value deduced from early drilling data. Rocks around the crater are fractured as distant as 900 m from the rim crest and to a depth of at least 800 m beneath the crater floor.

Ackermann, H. D.↗

Early arrival waveform inversion using data uncertainties and matching filters with application to near-surface seismic refraction data

We develop an early arrival waveform inversion (EAWI) technique for high-resolution near-surface velocity estimation by iteratively updating the P-wave velocity model to minimize the difference between the observed and calculated seismic refraction data. Traditional EAWI uses a least-squares penalty function and an acoustic forward-modeling engine. Conventional least-squares error is sensitive to data with low signal-to-noise ratio (S/N) and iterations of EAWI stop at a local-minimum data misfit or at the preassigned maximum number of iterations. These stopping criteria can result in overfitting the data. In addition, fitting the elastic field data with an acoustic modeling engine can introduce artifacts in velocity estimation, especially in land data with significant elastic effects. To overcome these challenges, we develop a robust EAWI (REAWI) method by (1) incorporating the data uncertainties into the penalty function and (2) mitigating the elastic effects using a matching filter workflow. The data uncertainties are estimated from waveform reciprocal errors. When full-waveform reciprocity is not available, trace interpolation is applied. The proposed method prevents closely fitting data with low S/N, avoids overall overfitting by stopping the iterations when a normalized chi-square ([Formula: see text]) waveform misfit of one is achieved, and is less affected by elastic effects. Numerical examples and application to near-surface refraction data at a groundwater contamination site suggest that the final REAWI models are more accurate than the corresponding EAWI models, at the same level of misfit. This is the first known application of a matching filter workflow to real land data. The final REAWI models satisfy an appropriate misfit between the real data and predicted elastic P-wave data, making this approach in this respect equivalent to elastic waveform inversion. We also develop a method to analyze model constraint by examining the energy of the wavefield Fréchet derivative thereby avoiding the influence of the data residuals in traditional Fréchet kernels.

Geochemistry & Geophysics↗

Seismic refraction tracks porosity generation and possible CO 2 production at depth under a headwater catchment

Significance To understand groundwater in weathered bedrock aquifers, we explore (bio)geochemical reactions using geophysical measurements. The result leads to a conceptual model for the three-dimensional distribution of geochemical alteration in the subsurface—well beyond the typical approach of inferring (bio)geochemical reactions in the subsurface from pits, boreholes, and outcrops. With this geophysical technique we discovered the prevalence of gas beneath the water table at this site. This geophysical viewing of the subsurface highlights patterns of water movement, reactions of clays, and microbial and abiotic gas-producing reactions recorded in the subsurface and thus illuminates storage, flow, and reactions of groundwater in bedrock systems. Protecting the groundwater resource beneath our feet will rely on such enhanced capacities to see into the subsurface.

58 GEOSCIENCES↗

Using Multiple Geophysical Methods to Refine a Stratigraphic Conceptual Site Model at a Nuclear Waste Site

At contaminated sites, there is a critical reliance on conceptual site models of which a key component is the geologic framework model (GFM), which is used to build predictive fate and transport modeling in support of remediation goals. Here, we detail an ongoing study at a nuclear waste site (Hanford Site, Washington, USA) where surface geophysical methods were used to complement existing borehole information and to help site new boreholes to update and refine the GFM. Ground-based electrical resistivity tomography (ERT), time-domain electromagnetics (TEM), and seismic data was acquired on the Hanford Site’s Central Plateau, where the presence of a hydraulically transmissive subsurface feature, or paleochannel, was suspected but uncertain based on contaminant concentrations at sparsely located boreholes. To assess and compare ERT, TEM, and seismic geophysical datasets, co-located data was acquired in an area where the existing GFM was more certain. Surrounding the suspected paleochannel, ERT and seismic refraction tomography showed similar subsurface structures consistent with a channelized feature. While TEM interpretation was limited in the suspected paleochannel area due to a thick resistive top layer and high electromagnetic noise, TEM showed greater utility in a different area where these factors were less prevalent. Based on the results of quasi-3D ERT inversions, we propose updates to the GFM, including identifying a paleochannel. Further, we evaluate the utility of geophysical methods and examine lessons learned that will be used as a basis for continued characterization efforts to support site management decisions and implementation efforts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The velocity structure of the lunar crust.

Seismic refraction data, obtained at the Apollo 14 and 16 sites, when combined with other lunar seismic data, allow a compressional wave velocity profile of the lunar near-surface and crust to be derived. The regolith, although variable in thickness over the lunar surface, possesses surprisingly similar seismic properties. Underlying the regolith at both the Apollo 14 Fra Mauro site and the Apollo 16 Descartes site is low-velocity brecciated material or impact derived debris. Key features of the lunar seismic velocity profile are: (1) velocity increases from 100 to 300 m/sec in the upper 100 m to about 4 km/sec at 5 km depth, (2) a more gradual increase from about 4 km/sec to about 6 km/sec at 25 km depth,(3) a discontinuity at a depth of 25 km, and (4) a constant value of about 7 km/sec at depths from 25 km to about 60 km.

Kovach, R. L.↗

Apollo 14 and 16 Active Seismic Experiments, and Apollo 17 Lunar Seismic Profiling

Seismic refraction experiments were conducted on the moon by Apollo astronauts during missions 14, 16, and 17. Seismic velocities of 104, 108, 92, 114 and 100 m/sec were inferred for the lunar regolith at the Apollo 12, 14, 15, 16, and 17 landing sites, respectively. These data indicate that fragmentation and comminution caused by meteoroid impacts has produced a layer of remarkably uniform seismic properties moonwide. Brecciation and high porosity are the probable causes of the very low velocities observed in the lunar regolith. Apollo 17 seismic data revealed that the seismic velocity increases very rapidly with depth to 4.7 km/sec at a depth of 1.4 km. Such a large velocity change is suggestive of compositional and textural changes and is compatible with a model of fractured basaltic flows overlying anorthositic breccias. 'Thermal' moonquakes were also detected at the Apollo 17 site, becoming increasingly frequent after sunrise and reaching a maximum at sunset. The source of these quakes could possibly be landsliding.

Source record↗

Seismic properties of the crust and uppermost mantle of North America

Seismic refraction profiles for the North American continent were compiled. The crustal models compiled data on the upper mantle seismic velocity (P sub n), the crustal thickness (H sub c) and the average seismic velocity of the crystalline crust (V sub p). Compressional wave parameters were compared with shear wave data derived from surface wave dispersion models and indicate an average value for Poisson's ratio of 0.252 for the crust and of 0.273 for the uppermost mantle. Contour maps illustrate lateral variations in crustal thickness, upper mantle velocity and average seismic velocity of the crystalline crust. The distribution of seismic parameters are compared with a smoothed free air anomaly map of North America and indicate that a complidated mechanism of isostatic compensation exists for the North American continent. Several features on the seismic contour maps also correlate with regional magnetic anomalies.

Braile, L. W.↗

Age constraints for the present fault configuration in the Imperial Valley, California: Evidence for northwestward propagation of the Gulf of California rift system

Releveling and other geophysical data for the Imperial Valley of southern California suggest the northern section of the Imperial-Brawley fault system, which includes the Mesquite Basin and Brawley Seismic Zone, is much younger than the 4 to 5 million year age of the valley itself. A minimum age of 3000 years is calculated for the northern segment of the Imperial fault from correlations between surface topography and geodetically observed seismic/interseismic vertical movements. Calculations of a maximum age of 80,000 years is based upon displacements in the crystalline basement along the Imperial fault, inferred from seismic refraction surveys. This young age supports recent interpretations of heat flow measurements, which also suggest that the current patterns of seismicity and faults in the Imperial Valley are not long lived. The current fault geometry and basement morphology suggest northwestward growth of the Imperial fault and migration of the Brawley Seismic Zone. It is suggested that this migration is a manifestation of the propagation of the Gulf of California rift system into the North American continent.

Larsen, Shawn↗

Age constraints for the present fault configuration in the Imperial Valley, California - Evidence for northwestward propagation of the Gulf of California rift system

Releveling and other geophysical data for the Imperial Valley of southern California suggest the northern section of the Imperial-Brawley fault system, which includes the Mesquite Basin and Brawley Seismic Zone, is much younger than the 4 to 5 million year age of the valley itself. A minimum age of 3000 years is calculated for the northern segment of the Imperial fault from correlations between surface topography and geodetically observed seismic/interseismic vertical movements. Calculations of a maximum age of 80,000 years is based upon displacements in the crystalline basement along the Imperial fault, inferred from seismic refraction surveys. This young age supports recent interpretations of heat flow measurements, which also suggest that the current patterns of seismicity and faults in the Imperial Valley are not long lived. The current fault geometry and basement morphology suggest northwestward growth of the Imperial fault and migration of the Brawley Seismic Zone. It is suggested that this migration is a manifestation of the propagation of the Gulf of California rift system into the North American continent.

Larsen, Shawn↗