Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic wave”

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 505 records · Page 28

Detection Limits and Near‐Field Ground Motions of Fast and Slow Earthquakes

Abstract We investigate theoretical limits to detection of fast and slow seismic events, and spatial variations of ground motion expected from M 6 earthquakes at short epicentral distances. The analyses are based on synthetic velocity seismograms calculated with the discrete wavenumber method assuming seismic velocities and attenuation properties of the crust in Southern California. The examined source properties include different magnitudes ( M −1.0 to M 6.0), static stress drops (0.1–10 MPa), and slow and fast ruptures (0.1–0.9 of shear wave velocity). For the M 6 events we also consider variations in rise times producing crack‐ and pulse‐type events and different rupture directivities. Slow events produce ground motion with considerably lower amplitude than corresponding regular fast earthquakes with the same magnitude, and hence are significantly more difficult to detect. The static stress drop and slip rise time also affect the maximum radiated seismic motion, and hence event detectability. Apart from geometrical factors, the saturation and depletion of seismic ground motion at short epicentral distances stem from radiation pattern, earthquake size (magnitude, stress drop), and rupture directivity. The rupture velocity, rise time, and directivity affect significantly the spatial pattern of the ground motions. The results can help optimizing detection of slow and fast small earthquakes and understand the spatial distribution of ground motion generated by large events.

58 GEOSCIENCES↗

CODAS Data from Oliktok Point, Beaufort Sea, Alaska

Cryosphere/Ocean Distributed Acoustic Sensing (CODAS) data collected from the Beaufort Sea, Alaska, using ~37.4 km of dark telecommunications fiber located at Oliktok Point, Alaska. Data were collected with a Silixa iDAS, using 10 m gauge length, 2 m spatial resolution, and 1000 Hz sample rate. Provided here are the DAS-recorded time series for the rapid refreeze event described in Baker & Abbott (2022) (see link below). This covers a date range of 2021-11-10 15:00 UTC to 2021-11-11 17:00 UTC. Data have been decimated to 100 Hz and 20 m (i.e., every 10th channel for 1831 channels, total), as used in Baker & Abbott (2022). Data have been extracted from raw format into 1-hour long .sac* files and organized into directories by channel number, spanning channels 100 to 18400. Time series units are nano-strainrate (nm/m/s). For distribution, data have been compressed into .zip files containing all time series files for 100 channels. *For information on the Seismic Analysis Code (SAC) file format: https://seiscode.iris.washington.edu/projects/sac

16 TIDAL AND WAVE POWER↗

Elastic wave velocities of Apollo 14, 15, and 16 rocks

Elastic wave velocities of two Apollo 14 rocks, 14053 and 14321, three Apollo 15 rocks, 15058, 15415, and 15545, and one Apollo 16 rock 60315 have been determined at pressures up to 10 kb. For sample 14321, the variation of the compressional wave velocities with temperature has been measured over the temperature range from 27 to 200 C. Overall elastic properties of these samples except sample 15415 are very similar to those of Apollo 11, 12, and 14 rocks and are concordant with Toksoz et al.'s (1972) interpretation that lunar upper crust is of basaltic composition. Temperature derivative of the P wave velocity for sample 14321 is a half to one order of magnitude larger than that for single crystalline minerals. This suggests that the seismic velocity in the lunar crust may be affected significantly by the temperature distribution.

Mizutani, H.↗

Implications of elastic wave velocities for Apollo 17 rock powders

Ultrasonic P- and S-wave velocities of lunar rock powders 172701, 172161, 170051, and 175081 were measured at room temperature and to 2.5 kb confining pressure. The results compare well with those of terrestrial volcanic ash and powdered basalt. P-wave velocity values up to pressures corresponding to a lunar depth of 1.4 km preclude cold compaction alone as an explanation for the observed seismic velocity structure at the Apollo 17 site. Application of small amounts of heat with simultaneous application of pressure causes rock powders to achieve equivalence of seismic velocities for competent rocks.

Talwani, P.↗

Science on the Moon: The Wailing Wall of Space Exploration

Science on and from the Moon has important implications for expanding human knowledge and understanding, a prospect for the 21st Century that has been under discussion for at least the past 25 years. That having been said, however, there remain many issues of international versus national priorities, strategy, economy, and politics that come into play. The result is a very complex form of human behavior where science and exploration take center stage, but many other important human options are sacrificed. To renew this dialogue about the Moon, it seems we are already rushing pell-mell into it as has been done in the past. The U.S., Japan, China, India, and Russia either have sent or plan to send satellites and robotic landers there at this time. What does a return to the Moon mean, why are we doing this now, who should pay for it, and how? The only semblance of such a human enterprise seems to be the LHC currently coming online at CERN. Can it be used as a model of international collaboration rather than a sports or military event focused on national competition? Who decides and what is the human sacrifice? There are compelling arguments for establishing science on the Moon as one of the primary goals for returning to the Moon and venturing beyond. A number of science endeavors will be summarized, beyond lunar and planetary science per se. These include fundamental physics experiments that are background-limited by the Earth's magnetic dipole moment and noise produced by its atmosphere and seismic interior. The Moon is an excellent platform for some forms of astronomy. Other candidate Moon-based experiments vary from neutrino and gravitational wave astronomy, particle astrophysics, and cosmic-ray calorimeters, to space physics and fundamental physics such as proton decay. The list goes on and includes placing humans in a hostile environment to study the long-term effects of space weather. The list is long, and even newer ideas will come from this COSPAR conference. However, whatever the list the issue of cooperation and binding collaboration remains. As observers of Moon and other space enterprises, we all know that a room full of 60 scientists will not agree on much of anything and there will probably be 60! please for more funding. People have special interests and little common sense (e.g., conflict between NSF- and NASA-funding roadmaps). Scientists are no exception. Nevertheless, CERN has done it on Earth! Can we do the same on the Moon? Some of the present generation of proposals for science from and on the Moon, plus new ones, will witness a place in space exploration's future. It is clear, however, that the world has not thought this through adequately, except for talk about an international space federation whatever that is. An outpost on the Moon with humans permanently living there much like Antarctica on Earth may be in our future. However, such planning is our collective international responsibility and not that of special-interest investigators from individual nations unless they intend to pay for it.

Wilson, Thomas↗

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗

Seismic scattering and shallow structure of the moon in Oceanus Procellarum

Scattering in a high-Q medium has been the best hypothesis to date for explaining the observation of a long, reverberating train of waves in lunar seismographs. To test this hypothesis qualitatively, two experiments were devised which simulated this scattering and reproduced actual lunar seismographs. Pulses were propagated across a plate with grooves cut half-way through, and then were propagated along the edge of a plate with holes drilled within a skin depth of the edge. The seismographs of near impacts and moonquakes recorded by the Apollo 12 station in two frequency bands were studied. The impacts were those of S4B Saturn boosters and LM ascent stages. Interpretation of these data suggests the existence of a scattering layer 25 km thick with a Q of 5000. The density of the scatterers decreases with depth, suggesting that they are associated with cratering, or that they consist of cracks which anneal with depth.

Dainty, A. M.↗

Seismic investigation of the lunar interior

The velocity and attenuation structure of the moon below the crust is examined using surface events. The moon is divided into an upper mantle and a lower mantle, the division at a depth of about 500 km being marked by a reflector identified on polarization filtered record sections. The upper mantle has a P-wave velocity of about 8 km/sec, a Poisson's ratio of about 0.25 and a Q for P waves of about 5000. This region contains no partial melt and is depleted in volatiles, notably water. The lower mantle has a lower S-wave velocity and probably a lower P-wave velocity than the upper mantle, with a Poisson's ratio of about 0.34. The lower mantle has a Q for P waves of approximately 1500, substantially lower than the upper mantle but probably still high enough to preclude partial melting. The velocity structure and the current value of the moment of inertia factor indicate an increase of density below about 500 km, perhaps due to an increase in iron content. We do not have any information directly pertaining to seismic velocities below 1000 km depth.

Dainty, A. M.↗

An adaptive model-free robotic force control strategy for hydrodynamic real-time hybrid simulation of floating offshore wind turbines

Real-time hybrid simulation (RTHS) - a cyber-physical testing approach - promises to enhance the simulation fidelity of the model-scale experiments used to prototype floating offshore wind turbines (FOWTs). In hydrodynamic RTHS (hydro-RTHS), actuators emulate aerodynamic forces on model-scale FOWT specimens subjected to physical waves in a hydrodynamic laboratory. Robotic arms are promising candidates for actuation in hydro-RTHS due to their compact multi-degree-of-freedom (DOF) capabilities. Unlike classical RTHS for seismic applications, which typically relies on displacement control, hydro-RTHS requires 6-DOF force control on newly designed floating prototypes in a model-scale setting, which presents significant challenges, including modeling uncertainties, directional asymmetry, configuration drift, bandwidth limitations, and time-varying delays. To mitigate these constraints without extensive pre-test calibration, this study proposes an adaptive model-free robotic force control strategy that combines task-space explicit force control with a secondary joint-space pose-keeping task. The Adaptive Feedforward Compensator (AFC) is integrated into the force control loop to compensate for time-varying delay. Experimental testing was conducted using a Franka Emika Panda robotic arm with a 1:50 scale FOWT specimen under operational wind and wave conditions. Results demonstrate stable and consistent 6-DOF force tracking. Effective delay compensation was observed, with low-frequency delay reductions ranging from 71.4% to 91.8% and improvements in low-frequency surge force tracking of 25.0% to 52.1%. This study enhances robotic actuation performance in hydro-RTHS and introduces a force control strategy that supports reliable robotic operation in uncertain floating environments. Future work will explore disturbance-observer mechanisms to further enhance wave rejection capabilities under extreme wind and wave conditions.

17 WIND ENERGY↗

Characterization of damage processes in Montney siltstone under triaxial compression using acoustic emission and diagnostic imaging

SUMMARY Crack nucleation and rock failure processes in a fine-grained siltstone (Montney Formation) under triaxial compression are investigated using combined diagnostic techniques, including ultrasonic-wave measurement, acoustic-emission (AE) monitoring, computed tomography (CT) scanning, and thin-section imaging. The sample displays a weak-to-moderate inherent seismic anisotropy and noticeable stress-induced anisotropy prior to failure. No AE event was detected until the applied axial stress reached 95 per cent of the peak value. The signal-to-noise ratio is relatively low, however, and detectable AE events are more diffuse than those observed in highly brittle rocks. The AE locations correlate with a shear fracture zone imaged by CT scanning. AE moment-tensor analysis reveals that events with larger relative magnitudes are characterized by high volumetric (tensile or compressive) components, and the initiation of the failure zone is dominated by combined shear-tensile (ST) failure. Stress inversion of the AE events with high tensile components is in good agreement with the known applied stress. Microscopic imaging of thin sections from the failed sample shows that the failure zone is an en echelon structure consisting of a major fracture with branching micro and minor cracks. This failure mechanism is consistent with a ST source mechanism and is interpreted to be associated with the fine granular structure and mineral composition of Montney siltstone.

Jia, Suzie Q. (ORCID:0000000326496190)↗

Strength and compressibility of returned lunar soil.

Two oedometer and three direct shear tests have been performed in vacuum on a 200 g sample of lunar soil from Apollo 12 (12001, 119). The compressibility data have been used to calculate bulk density and shear wave velocity versus depth on the lunar surface. The shear wave velocity was found to increase approximately with the one-fourth power of the depth, and the results suggest that the Apollo 14 Active Seismic Experiment may not have detected the Fra Mauro formation at a depth of 8.5 m, but only naturally consolidated lunar soil. The shear data indicate that the strength of the lunar soil sample is about 65% that of a ground basalt simulant at the same void ratio.

Carrier, W. D., III↗

Image-domain DAS 3D VSP elastic transmission tomography

SUMMARY Full-wavefield elastic imaging of active-source seismic data acquired by downhole receivers commonly offers higher-resolution subsurface images in the vicinity of a borehole compared to conventional surface seismic data sets, which can lack higher-frequency wavefield components due to longer travel paths and increased attenuation. An increasingly used approach for downhole acquisition is vertical seismic profiling (VSP), which has become more attractive when coupled with distributed acoustic sensing (DAS) using optical fibres installed in wells. The main difficulty for generating high-quality images with full-wavefield imaging tools for DAS VSP data, though, is the need for an accurate velocity model. To build plausible velocity models using active-source DAS VSP data, we adopt a 3-D image-domain elastic transmission tomography technique, originally developed for surface-recorded passive (microseismic) data, by exchanging the source and receiver positions (i.e. reciprocity) to mimic a passive-seismic surface monitoring scenario. The inversion approach exploits various images for each source constructed through time-reverse imaging (TRI) of downgoing P- and S-wave first-arrival waveforms. The TRI process uses the kinetic term of the (extended) PS energy imaging condition that exhibits sufficient sensitivity to velocity model errors. The method automatically updates the P- and S-wave velocity models to optimize image focusing via adjoint-state inversion. We illustrate the efficacy of the adopted elastic inversion technique using an active-source DAS 3-D VSP field data set acquired in the North Slope of Alaska. The numerical experiments demonstrate that the inverted elastic velocity models can be further used in full-wavefield acoustic/elastic imaging algorithms to obtain accurate subsurface images.

Oren, Can (ORCID:000000019534265X)↗

Seismic Elastic Double-Beam Characterization of Faults and Fractures for CO₂ Storage Site Selection

Site characterization for underground injection and storage of gigatonne-scale CO₂ requires reliable and cost-effective methods to detect and characterize faults and fractures and to assess their stress state and fault activation potential. This is critical, as wastewater injection and disposal have been shown to activate faults and induce earthquakes, and CO₂ leakage remains a key concern for long-term storage. In this project, we developed seismic methods to detect and characterize large-scale sedimentary and crystalline basement faults and associated small-scale fractures below conventional seismic imaging resolution using multicomponent (9C) surface seismic data. Machine learning was used to automatically interpret large-scale faults, providing key information for estimating the maximum magnitude of potential induced earthquakes. High-fidelity imaging was achieved by exploiting redundancy across multiple elastic wave modes, where independent images from different modes and frequencies cross-validate each other. We also used our nonlinear signal comparison (NLSC) method for ground roll removal, improving data quality in complex near-surface conditions. The methods were validated using field data acquired in central Montana. Results show that basement faults extend into the sedimentary section and that small-scale fractures are widespread above the basement. The inferred stress orientation is consistent with regional stress data, and the estimated maximum induced earthquake magnitude is small (Mw ~2.3). The developed workflow provides a practical approach for fault and fracture characterization and for assessing induced seismicity and leakage risk. It is directly applicable to CO₂ storage site selection and to other subsurface systems.

02 PETROLEUM↗

Validation of (not-historical) large-event near-fault ground-motion simulations for use in civil engineering applications

Ground-motion simulations generated from physics-based wave propagation models are gaining increasing interest in the engineering community for their potential to inform the performance-based design and assessment of infrastructure residing in active seismic areas. A key prerequisite before the ground-motion simulations can be used with confidence for application in engineering domains is their comprehensive and rigorous investigation and validation. This article provides a four-step methodology and acceptance criteria to assess the reliability of simulated ground motions of not historical events, which includes (1) the selection of a population of real records consistent with the simulated scenarios, (2) the comparison of the distribution of Intensity Measures (IMs) from the simulated records, real records, and Ground-Motion Prediction Equations (GMPEs), (3) the comparison of the distribution of simple proxies for building response, and (4) the comparison of the distribution of Engineering Demand Parameters (EDPs) for a realistic model of a structure. Specific focus is laid on near-field ground motions (<10km) from large earthquakes (M w 7), for which the database of real records for potential use in engineering applications is severely limited. The methodology is demonstrated through comparison of (2490) near-field synthetic records with 5 Hz resolution generated from the Pitarka et al (2019) kinematic rupture model with a population of (38) pulse-like near-field real records from multiple events and, when applicable, with NGA-W2 GMPEs. Lastly, the proposed procedure provides an effective method for informing and advancing the science needed to generate realistic ground-motion simulations, and for building confidence in their use in engineering domains.

42 ENGINEERING↗

Apollo 14 active seismic experiment.

Explosion seismic refraction data indicate that the lunar near-surface rocks at the Apollo 14 site consist of a regolith 8.5 meters thick and characterized by a compressional wave velocity of 104 meters per second. The regolith is underlain by a layer with a compressional wave velocity of 299 meters per second. The thickness of this layer, which we interpret to be the Fra Mauro Formation, is between 16 and 76 meters. The layer immediately beneath this has a velocity greater than 370 meters per second. We found no evidence of permafrost.

Watkins, J. S.↗

Lunar seismicity, structure, and tectonics

Seismic data is used to develop a lunar model consisting of five zones. These include a 50-60 km thick crust characterized by seismic velocities appropriate to plagioclase rich materials, the 250 km thick upper mantle characterized by seismic velocities consistent with an olivine-pyroxene composition, the 500 km thick middle mantle characterized by a high Poisson ratio, the lower mantle characterized by high shear wave attenuation, and a core of radius between 170 and 360 km characterized by a greatly reduced compressional wave velocity.

Lammlein, D. R.↗

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↗

Host-region parameters for an adjustable model for crustal earthquakes to facilitate the implementation of the backbone approach to building ground-motion logic trees in probabilistic seismic hazard analysis

The backbone approach to constructing a ground-motion logic tree for probabilistic seismic hazard analysis (PSHA) can address shortcomings in the traditional approach of populating the branches with multiple existing, or potentially modified, ground-motion models (GMMs) by rendering more transparent the relationship between branch weights and the resulting distribution of predicted accelerations. To capture epistemic uncertainty in a tractable manner, there are benefits in building the logic tree through the application of successive adjustments for differences in source, path, and site characteristics between the host region of the selected backbone GMM and the target region for which the PSHA is being conducted. The implementation of this approach is facilitated by selecting a backbone GMM that is amenable to such host-to-target adjustments for individual source, path, and site characteristics. The NGA-West2 GMM of Chiou and Youngs (CY14) has been identified as a highly adaptable model for crustal seismicity that is well suited to such adjustments. Rather than using generic source, path, and site characteristics assumed appropriate for the host region, the final suite of adjusted GMMs for the target region will be better constrained if the host-region parameters are defined specifically on the basis of their compatibility with the CY14 backbone GMM. To this end, making use of a recently developed crustal shear-wave velocity profile consistent with CY14, we present an inversion of the model to estimate the key source and path parameters, namely the stress parameter and the anelastic attenuation. With these outputs, the effort in constructing a ground-motion logic tree for any PSHA dealing with crustal seismicity can be focused primarily on the estimation of the target-region characteristics and their associated uncertainties. The inversion procedure can also be adapted for any application in which different constraints might be relevant.

Stafford, Peter J.↗