Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismograms”

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 37 records · Page 2

Distinct meteoroid families identified on the lunar seismograms

The meteoroid impact-seismic activity data recorded by the Apollo lunar seismic network is examined. The study investigates the difference in temporal distribution between large and small impacts, clustering of impacts in a two-dimensional space of the time of the year and the time of the month, and the relationship of these observations with terrestrial observations. Several distinct families of meteoroids impacting the moon are identified. Most meteoroids producing small impact-seismic events, including ones associated with cometary showers, appear to approach from retrograde heliocentric orbits. In contrast, most meteoroids associated with large impact-seismic events appear to approach from prograde orbits; the observation is consistent with a hypothesis that many of them represent stony asteroidal material. It is suggested that the previously reported discrepancy between lunar and terrestrial meteoroid-flux estimates may be due to the differences in lunar and terrestrial detection efficiency among various families of meteoroids.

Oberst, Jurgen↗

Towards Simulating a Realistic Planetary Seismic Wavefield: The Contribution of the Megaregolith and Low-Velocity Waveguides

Lunar seismograms are distinctly different from their terrestrial counterparts. The Apollo lunar seismometers recorded moonquakes without distinct P- or S-wave arrivals; instead waves arrive as a diffuse coda that decays over several hours making the identification of body waves difficult. The unusual character of the lunar seismic wavefield is generally tied to properties of the megaregolith: it consists of highly fractured and broken crustal rock, the result of extensive bombardment of the Moon. The megaregolith extends several kilometers into the lunar crust, possibly into the mantle in some regions, and is covered by a thin coating of fine-scale dust. These materials possess very low seismic velocities that strongly scatter the seismic wavefield at high frequencies. Directly modeling the effects of the megaregolith to simulate an accurate lunar seismic wavefield is a challenging computational problem, owing to the inherent 3-D nature of the problem and the high frequencies (greater than 1 Hz) required. Here we focus on modeling the long duration code, studying the effects of the low velocities found in the megaregolith. We produce synthetic seismograms using 1-D slowness integration methodologies, GEMINI and reflectivity, and a 3-D Cartesian finite difference code, Wave Propagation Program, to study the effect of thin layers of low velocity on the surface of a planet. These codes allow us generate seismograms with dominant frequencies of approximately 1 Hz. For background lunar seismic structure we explore several models, including the recent model of Weber et al., Science, 2011. We also investigate variations in megaregolithic thickness, velocity, attenuation, and seismogram frequency content. Our results are compared to the Apollo seismic dataset, using both a cross correlation technique and integrated envelope approach to investigate coda decay. We find our new high frequency results strongly support the hypothesis that the long duration of the lunar seismic codes is generated by the presence of the low velocity megaregolith, and that the diffuse arrivals are a combination of scattered energy and multiple reverberations within this layer. The 3-D modeling indicates the extreme surface topography of the Moon adds only a small contribution to scattering effects, though local geology may play a larger role. We also study the effects of the megaregolith on core reflected and converted phases and other body waves. Our analysis indicates detection of core interacting arrivals with a polarization filter technique is robust and lends the possibility of detecting other body waves from the Moon.

Schmerr, Nicholas C.↗

The upper mantle beneath the Cascade Range: A comparison with the Gulf of California

Seismograms from 22 earthquakes along the northeast Pacific rim recorded in southern California form the data set for investigation of the upper mantle beneath the Cascade Range-Juan de Fuca region, a transitional area encompassing both very young ocean floor and a continental margin. These data consist of 853 seismograms (6 deg delta 42 deg) which produce 1068 travel times and 40 ray parameter estimates. These data are compared directly to another large suite of records representative of structure beneath the Gulf of California, an active spreading center. The spreading center model, GCA, was used as a starting point in WKBJ synthetic seismogram modeling and perturb GCA until the northeast Pacific data are matched. Application of wave field continuation to these two groups of data provides checks on model's consistency with the data as well as an estimate of the resolvability of differences between the two areas. Differences between the models derived from these two data sets are interpretable in terms of lateral structural variation beneath the two regimes.

Walck, M. C.↗

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↗

Core Angular Momentum and the IERS Sub-Centers Activity for Monitoring Global Geophysical Fluids: Core Angular Momentum and Earth Rotation - Part 1

The part of the grant was to use recordings of seismic waves travelling through the earth's core (PKP waves) to study the inner core rotation and constraints on possible density anomalies in the fluid core. The shapes and relative arrival times of such waves associated with a common source were used to reduce the uncertainties in source location and excitation and the effect of unknown mantle structure. The major effort of the project is to assemble historical seismograms with long observing base lines. We have found original paper records of SSI earthquakes at COL between 1951 and 1966 in a warehouse of the U.S. Geological Survey office in Golden, Colorado, extending the previous measurements at COL by Song and Richards [1996] further back 15 years. Also in Alaska, the University of Alaska, Fairbanks Geophysical Institute (UAFGI) has been operating the Alaskan Seismic Network with over 100 stations since the late 1960s. Virtually complete archives of seismograms are still available at UAFGI. Unfortunately, most of the archives are in microchip form (develocorders), for which the use of waveforms is impossible. Paper seismograms (helicorders) are available for a limited number of stations, and digital recordings of analog signals started around 1989. Of the paper records obtained, stations at Gilmore Dome (GLM, very close to COL), Yukon (FYU), McKinley (MCK), and Sheep Creek Mountain (SCM) have the most complete continuous recordings.

Song, Xia-Dong↗

GRAIL Refinements to Lunar Seismic Structure

The present ]day internal structure of the Moon provides insight not only into its own formation and evolution, but also that of all rocky planetary bodies. The most direct way to probe a planet fs interior structure is through seismology. As part of the Apollo lunar missions, four seismometers were deployed on the nearside surface of the Moon between the years 1969 and 1972. These instruments operated continuously until 1977, forming the only substantial extraterrestrial seismic data set in existence. These data have been used to constrain various aspects of the seismic velocity and density structure of the Moon. Typical 1-D models recognize a 30-60 km thick crust overlying a nearly constant ]velocity mantle, and extend to a depth of approximately 1000 km, below which the lack of penetrating moonquake ray ]paths precludes the seismic determination of deeper structure. Previously, the lack of observed moonquakes from the far side of the Moon has been used to infer the presence of a highly attenuating (possibly molten) core. Indirect geophysical measurements such as moment of inertia, magnetic induction, lunar laser ranging, and elemental abundances of mare basalts also place varying constraints on core size and state. In combination with seismic studies, these indirect measurements have been used to arrive at a commonly accepted model of the Moon's deepest interior that includes a solid inner and fluid outer core, overlain by a partial melt boundary layer. We recently applied modern array seismology techniques to the Apollo data and revealed detailed core structure, including the first direct confirmation of the presence of a solid inner core. Our study focused on the identification of core ]reflected phases in deep moonquake seismograms. The resulting model of the Moon fs innermost structure was found to be consistent with the commonly accepted model. However, the modeled layer radii may vary by tens of kilometers, as is expected when accounting for uncertainties such as moonquake location, timing errors, and potential seismic heterogeneities. In addition, the modeled velocities may vary with a 1-to-1 trade ]off with the modeled reflector depth. The GRAIL (Gravity Recovery and Interior Laboratory) mission, launched in Sept. 2011, placed two nearly identical spacecraft in lunar orbit. The two satellites make extremely high-resolution measurements of the lunar gravity field, which can be used to constrain the interior structure of the Moon using a "crust to core" approach. GRAIL fs constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten) will complement seismic investigations in several ways. Here we present a progress report on our efforts to advance our knowledge of the Moon fs internal structure using joint gravity and seismic analyses. We will focus on methodology, including 1) refinements to the seismic core constraint accomplished through array processing of Apollo seismic data, made by applying a set of travel time corrections based on GRAIL structure estimates local to each Apollo seismic station; 2) modeling deep lunar structure through synthetic seismograms, to test whether the seismic core model can reproduce the core reflections observed in the Apollo seismograms; and 3) a joint seismic and gravity inversion in which we attempt to fit a family of seismic structure models with the gravity constraints from GRAIL, resulting in maps of seismic velocities and densities that vary from a nominal model both laterally and with depth.

Weber, Renee↗

Probing the evolution of fault properties during the seismic cycle with deep learning

We use seismic waves that pass through the hypocentral region of the 2016 M6.5 Norcia earthquake together with Deep Learning (DL) to distinguish between foreshocks, aftershocks and time-to-failure (TTF). Binary and N-class models defined by TTF correctly identify seismograms in test with > 90% accuracy. We use raw seismic records as input to a 7 layer CNN model to perform the classification. Here we show that DL models successfully distinguish seismic waves pre/post mainshock in accord with lab and theoretical expectations of progressive changes in crack density prior to abrupt change at failure and gradual postseismic recovery. Performance is lower for band-pass filtered seismograms (below 10 Hz) suggesting that DL models learn from the evolution of subtle changes in elastic wave attenuation. Tests to verify that our results indeed provide a proxy for fault properties included DL models trained with the wrong mainshock time and those using seismic waves far from the Norcia mainshock; both show degraded performance. Our results demonstrate that DL models have the potential to track the evolution of fault zone properties during the seismic cycle. If this result is generalizable it could improve earthquake early warning and seismic hazard analysis.

58 GEOSCIENCES↗

WaveDenoiser

We developed a robust deep learning model to effectively reduce background noise in the time domain from seismic waveforms. The deep learning model processes a 57-second three-component seismogram to predict and generate a denoised seismogram. The training was conducted using the benchmark STEAD dataset, which comprises globally distributed earthquake signals recorded at local distances ranging from 0 to 350 kilometers.

Chai, Chengping [Oak Ridge National Laboratory (OR↗

WaveDP

We developed a robust deep learning model designed to effectively reduce background noise and measure signal arrival times from seismic waveforms. This model processes a 57-second, three-component seismogram to predict the probability of Primary (P) and Secondary (S) waves for each timestamp, while also generating a denoised seismogram. Training was conducted using the benchmark STEAD dataset, which includes globally distributed earthquake signals recorded at local distances ranging from 0 to 350 kilometers.

Chai, Chengping [Oak Ridge National Laboratory (OR↗

Integration of Waveform Simulation Methods

The generation of synthetic seismograms through simulation is a fundamental tool of seismology required to run quantitative hypothesis tests. A variety of approaches have been developed throughout the seismological community and each has their own specific user interface based on their implementation. This causes a challenge to researchers who will need to learn new interfaces with each new software they wish to use and create substantial challenges when attempting to compare results from different tools. Here we provide a unified interface that facilitates interoperability amongst several simulation tools through a modern containerized Python package. Further, this package includes post-processing analysis modules designed to facilitate end-to-end analysis of synthetic seismograms. In this report we present the conceptual guidance and an example implementation of the new Waveform Simulation Framework.

58 GEOSCIENCES↗

Model Card for WaveDP

WaveDP was trained using both earthquake signals and pure noise waveforms from the Stanford Earthquake Dataset (STEAD). The model achieved a median signal-to-noise ratio (SNR) improvement of 7 dB. Because the model uses seismograms as direct input, the model is four times faster than an existing model for signal denoising that uses spectral representations of seismograms as input.

47 OTHER INSTRUMENTATION↗

Participation in the Apollo passive seismic experiment

Computer programs which were written to read digital tapes containing lunar seismic data were studied. Interpreting very early parts of the lunar seismogram as seismic body-wave phases enabled the determination of the structure of the outer part of the moon in the Fra Mauro region. The crust in the Fra Mauro region is 60 to 65 km-thick, overlaying a high velocity mantle. The crust is further divided into an upper part, 25 km thick, apparently made of material similar to the surficial basalts, and a lower part of seemingly different composition, possibly an anorthositic gabbro. The generation of the exceedingly long reverberating wave-train observed in lunar seismogram was also studied. This is believed to be due to an intense scattering layer with very high quality coefficient overlying a more homogeneous elastic medium. Titles and abstracts of related published papers are included.

Press, F.↗

Velocity structure and evolution of the moon

Seismic data from the Apollo Passive Seismic Network stations are analyzed to determine the velocity structure and to infer the composition and physical properties of the lunar interior. Data from artificial impacts (SIBV booster and LM-ascent stage) cover a distance range of 9 to 1750 km. Travel times and amplitudes, as well as theoretical seismograms, are used to derive a velocity model for the outer 150 km of the moon. The P-wave velocity model confirms an earlier report of a lunar crust in the eastern part of Oceanus Procellarum. The crust is about 60 km thick and may consist of two layers in the mare regions. Possible values for the P-wave velocity in the uppermost mantle are between 7.6 and 9.0 km/sec. The 9 km/sec velocity represents either a localized heterogeneous unit, or a thin layer less than about 40 km in thickness. The elastic properties of the deep interior, as inferred from the seismograms of natural events (meteoroid impacts and moonquakes) occurring at great distances, indicate that there is an increase in attenuation and a possible decrease of velocity at depths below about 1000 km.

Toksoz, M. N.↗

Lunar velocity structure and compositional and thermal inferences

Seismic data from Apollo Passive Seismic Network stations are analyzed to determine the velocity structure in the lunar crust and mantle. Since the publication of earlier results additional data has become available from the S-IVB impacts of the Apollo 16 and 17 missions and the LM ascent stage impact of the Apollo 17 mission. Data from these artificial impacts now cover a distance range of 9 to 1700 km. Travel times, amplitudes, and theoretical seismograms are used to derive a velocity structure for the outer 150 km of the moon. Preliminary results are given from the analysis of the complete set of seismograms.

Toksoz, M. N.↗

Seismic energy transmission in an intensively scattering environment

In order to account for some special features of lunar seismograms, namely, the gradual build-up of the signal, the extremely prolonged tail, and the lack of apparent coherence among three orthogonal components of ground motion, a statistical approach is proposed for describing transmission of seismic energy through a medium in which strong scattering takes place. A seismic diffusion theory is presented for a medium with randomly distributed scatterers of a given size distribution. A solution of the resulting diffusion equation for an impulsive energy source gives a curve which fairly closely reproduces the envelope of typical lunar impact seismograms. Since the model is based on constant diffusivity, long-range transmission will require a combination of diffusion and wave propagation treatments for accurate description.

Nakamura, Y.↗

The p-wave upper mantle structure beneath an active spreading centre - The Gulf of California

Over 1400 seismograms of earthquakes in Mexico are analyzed and data sets for the travel time, apparent phase velocity, and relative amplitude information are utilized to produce a tightly constrained, detailed model for depths to 900 km beneath an active oceanic ridge region, the Gulf of California. The data are combined by first inverting the travel times, perturbing that model to fit the p-delta data, and then performing trial and error synthetic seismogram modelling to fit the short-period waveforms. The final model satisfies all three data sets. The ridge model is similar to existing upper mantle models for shield, tectonic-continental, and arc-trench regimes below 400 km, but differs significantly in the upper 350 km. Ridge model velocities are very low in this depth range; the model 'catches up' with the others with a very large velocity gradient from 225 to 390 km.

Walck, M. C.↗

High-resolution seismic reflection survey at the Manson crater, Iowa

Approximately 17.4 km of high-resolution reflection data were acquired along an east-west radius of the Manson Impact Structure (MIS) to delineate the shallow (upper 300 m) subsurface structural configuration. The geometry of the shallow structure is poorly known due to a 30-90 m thick Pleistocene till cover. The resolution of the new seismic data is roughly 5-10 times that of existing Vibroseis data. Data quality varies rapidly along the line from exceptional to poor, due primarily to velocity variations associated with the geological complexity of the area. Preliminary results indicate subsurface structural blocks previously envisioned to be several hundreds of meters in size are actually an order of magnitude smaller and more complex. A seismogram-by-seismogram analysis is necessary to confidently identify intricate stratigraphic and structural relationships seen on preliminary CDP sections, as numerous faults, diffractions, and complicated reflection patterns create potential pitfalls.

Keiswetter, D. A.↗

Deep Moonquake Focal Mechanisms: Recovery and Implications

A defining characteristic of deep moonquakes is their tendency to occur with tidal periodicity, prompting previous studies to infer that they are related to the buildup and release of tidal stress within the Moon. In studies of tidal forcing, a key constraint is the focal mechanism: the fault parameters describing the type of failure moonquakes represent. The quality of the lunar seismic data and the limited source/receiver geometries of the Apollo seismic network prohibit the determination of deep moonquake fault parameters using first-motion polarities, as is typically done in terrestrial seismology. Without being able to resolve tidal stress onto a known failure plane, we can examine only gross qualities of the tidal stress tensor with respect to moonquake occurrence, so we cannot fully address the role of tidal stress in moonquake generation. We will examine the extent to which shear (S) and compression (P) wave amplitude ratios can constrain moonquake fault geometry by determining whether, for a given cluster, there exists a focal mechanism that can produce a radiation pattern consistent with the amplitudes measured by the Apollo instruments. Amplitudes are read in the ray coordinate frame, directly from seismograms for which the P and S arrivals are clearly identifiable on all long-period channels of the four Apollo stations. We apply an empirical station correction to account for site effects and the differences between P- and S-wave attenuation. Instead of focusing on the best fitting solution only, we formulate the inverse problem using a falsification criterion: all source orientations that do not reproduce the observed SV/P ratios within an error margin derived from the uncertainty of amplitude readings are rejected. All others are accepted as possible solutions. The inversion is carried out using an exhaustive grid search on a regular grid with predefined step size, encompassing all possible combinations of strike, dip and slip. To assess the sensitivity of the inversion for the uncertainty of the lunar interior structure, we carry out repeated inversions with different velocity structures. Our data set consist of a total of 106 events from 25 deep moonquake clusters. The largest contribution of 37 events originates from the most active cluster, A001, while other clusters are represented by 1 to 9 events. Since the definition of a cluster implies that all events share the same source orientation, a comparison of the inversion results of all events from one cluster will reduce ambiguities of the inversion. Once we obtain a suite of fault parameters for a given source, we can attempt to further constrain the focal mechanism with refined analyses of tidal stresses and predictions based on synthetic seismograms.

Knapmeyer, Martin↗