Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic modeling”

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 73 records · Page 4

The Seismic Effect of Impacts on Asteroid Surface Morphology: Early Modeling Results

In recent years, spacecraft observations of asteroids 951 Gaspra, 243 Ida, 253 Mathilde, and 433 Eros have shown the overriding dominance of impact processes with regard to the structure and appearance of these small, irregular bodies. One currently unanswered question is: to what degree does seismic shaking from impacts affect the surface morphology of these bodies? Through detailed modeling, we begin to address this question.

Richardson, J. E.↗

The distribution of earthquakes with depth and stress in subducting slabs

The global variation of Benioff zone seismicity with depth and the orientation of stress axes of deep and intermediate earthquakes is explained using numerical models of subducting slabs. Models that match the seismicity and stress require a barrier to flow at the 670 km seismic discontinuity. The barrier may be a viscosity increase of at least an order of magnitude or a chemical discontinuity. Instantaneous flow is subparallel to the slabs for models with a viscosity increase but contorted for models with a chemical barrier. Log N (number of earthquakes) decreases linearly to 250-300 km depth and increases thereafter. Stress magnitude in the models shows the same pattern, in accord with experiments showing N proportional to e(k-sigma), with k a constant and sigma the stress magnitude. The models predict downdip compression in the slabs at depths below 300-400 km, as observed for earthquake stress axes.

Vassiliou, M. S.↗

Geophysical constraints on lunar bulk composition and structure - A reassessment

Theoretical lunar mantle seismic velocity and density profiles were calculated for a series of specific bulk composition models, thermal models, and evolutionary models. It was found that most of the calculated seismic profiles are in qualitative agreement with the Nakamura (1983) observational mantle seismic velocity model. With respect to the density modeling results, the most interesting finding is that none of the considered models yield mantle density increases that are sufficient to match the adapted upper bound on the lunar moment-of-inertia factor (0.3928). The results suggest that the moon possesses a dense metallic core with a mass in the range of 1 to 4 percent of the lunar mass and that the alumina-rich models of Morgan et al. (1978) and Taylor (1982) are best able to match the middle mantle velocity increase that characterizes the Nakamura model.

Hood, L. L.↗

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.↗

Differential shear wave attenuation and its lateral variation in the North Atlantic region

A digital data base of over 150 seismograms and a spectral radio technique are used to measure SS-S differential attenuation in the North Atlantic region. Differential attenuation is positively correlated with SS-S travel time residual, and both differential attentuation and travel time residual decrease with increasing seafloor age. Models are developed for seismic Q in which lateral variations include contributions from the asthenospheric low-Q zone as well as from lithospheric cooling. The Q models obtained under this assumption are in good agreement with those obtained from surface wave studies and are therefore preferred over those models with lateral variations confined to the upper 125 km. Systematic long-wavelength (1000-7000 km) variations in differential attenuation, corrected for seafloor age, are evident along the axis of the Mid-Atlantic Ridge. These variations can be qualitatively correlated with long-wavelength variations in SS-S differential travel time residuals and are attributed to along-axis differences in upper mantle temperature.

Sheehan, Anne F.↗

Testing Spatial Correlation of Subduction Interplate Coupling and Forearc Morpho-Tectonics

Subduction zones that are capable of generating great (Mw greater than 8) earthquakes appear to have a common assemblage of forearc morphologic elements. Although details vary, each have (from the trench landward), an accretionary prism, outer arc high, outer forearc basin, an inner forean: basin, and volcanic arc. This pattern is common in spite of great variation in forearc architecture. Because interseismic strain is known to be associated with a locked seismogenic plate interface, we infer that this common forearc morphology is related, in an unknown way, to the process of interseismic Strain accumulation and release in great earthquakes. To date, however, no clear relationship between the subduction process and the common elements of upper plate form has emerged. Whereas certain elements of the system, i.e. the outer arc high, are reasonably well- understood in a structural context, there is little understanding of the structural or topographic evolution of the other key elements like the inner arc and inner forearc basin, particularly with respect to the coupled zone of earthquake generation. This project developed a model of the seismologic, topographic, and uplift/denudation linkages between forearc topography and the subduction system by: 1) comparing geophysical, geodetic, and topographic data from subduction margins that generate large earthquakes; 2) using existing GPS, seismicity, and other data to model the relationship between seismic cycles involving a locked interface and upper-plate topographic development; and 3) using new GPS data and a range-scale topographic, uplift, and denudation analysis of the presently aseismic Cascadia margin to constrain topographic/plate coupling relationships at this poorly understood margin.

Goldfinger, Chris↗

Postseismic surface deformations due to lithospheric and asthenospheric viscoelasticity

A model is proposed for post seismic surface deformations attributing them to lithospheric and asthenospheric viscoelasticity. The model predicts that the deformations due to lithospheric viscoelasticity depend on the ratio of the effective shear modulus acting long after the lithospheric viscoelastic relaxation to that acting immediately following the earthquake. While such deformations are generally smaller than those associated with asthenospheric viscoelasticity, they occur on a shorter time scale and may be in opposite direction to both the motion occurring at the time of the earthquake and that occurring as the asthenospheric relaxation occurs.

Cohen, S. C.↗

Numerical determination of the transmissibility characteristics of a squeeze film damped forced vibration system

Numerical solutions of the governing equations of motion of a liquid squeeze film damped forced vibration system were carried out to examine the feasibility of using a liquid squeeze film to cushion and protect large structures, such as buildings, located in areas of high seismic activity. The mathematical model used was that for a single degree of freedom squeeze film damped spring mass system. The input disturbance was simulated by curve fitting actual seismic data with an eleventh order Lagranging polynomial technique. Only the normal component of the seismic input was considered. The nonlinear, nonhomogeneous governing differential equation of motion was solved numerically to determine the transmissibility over a wide range of physical parameters using a fourth-order Runge-Kutta technique. It is determined that a liquid squeeze film used as a damping agent in a spring-mass system can significantly reduce the response amplitude for a seismic input disturbance.

Sutton, M. A.↗

On the joint inversion of geophysical data for models of the coupled core-mantle system

Joint inversion of magnetic, earth rotation, geoid, and seismic data for a unified model of the coupled core-mantle system is proposed and shown to be possible. A sample objective function is offered and simplified by targeting results from independent inversions and summary travel time residuals instead of original observations. These data are parameterized in terms of a very simple, closed model of the topographically coupled core-mantle system. Minimization of the simplified objective function leads to a nonlinear inverse problem; an iterative method for solution is presented. Parameterization and method are emphasized; numerical results are not presented.

Voorhies, Coerte V.↗

A Frequency-Dependent Q Model Based Upon Transient Rheological Properties of Heterogeneous Materials

Seismic attenuation at low frequency has previously been modeled using two approaches: (1) Examination of all possible microstructural dislocation mechanisms that involve irreversible thermodynamics [Minster and Anderson 1981, P. Trans. R.S. Lond, A, 299]; (2) Parameterized superposition of one-dimensional analogue springs and dashpots, usually consisting of a three-parameter standard linear solid arrangement [Liu et al. 1976, G.J.R.A.S., 47]. The first approach is appealing in that a strong connection is made between micromechanical dissipation and short-term creep that can be studied in the laboratory. The second approach has the advantage of providing a robust 'fit' to actual seismic observations, but has little intrinsic physical content...

Ivins, Erik R.↗

Analysis and geological interpretation of gravity data from GEOS-3 altimeter

A number of detailed gravimetric geoids of portions of the world's oceans from marine gravity measurements were constructed. The geoids were constructed by computing 1 x 1 deg or 10 x 10 deg averages of free-air anomaly data and subtracting these values from currently used satellite derived Earth models. The resulting difference gravity anomalies are then integrated over a sphere using a simplified form of Stoke's equation to obtain a difference geoid. This difference geoid is added to the satellite derived model to obtain a 1 x 1 deg or 10 x 10 deg total gravimetric geoid. The geoid undulations are studied by comparison of the altimeter measurements with the morphology of the ocean floor. Utilizing a combination of altimetry data, gravity and seismic reflection data, geophysical models of the earth can be constructed.

Talwani, M.↗

Constraints on lunar structure

A brief review is given of the constraints placed on lunar structure and composition by seismic data and density models. Bounds on the crustal velocity structure in Mare Cognitum are derived using travel-time data from artificial impacts, and a velocity model is determined on the basis of synthetic seismograms. It is shown that the velocities of P- and S-waves in the mantle can be fixed by a least-squares analysis of arrival times from meteor impacts and moonquakes, and that lunar density can be determined from the seismic structure, mean density, and moment of inertia. Olivine-pyroxene mixtures and certain olivine-rich compositions are found to be consistent with the seismic-velocity and density limits. Maximum radii are calculated for pure Fe and pure FeS cores, and it is concluded that the possibility of an ancient lunar magnetic dynamo may have to be reevaluated in the light of these figures.

Dainty, A. M.↗

The Use of Deep Moonquakes for Constraining the Internal Structure of the Moon

The installation of seismometers on the Moon s surface during the Apollo era provided a wealth of information that transformed our understanding of lunar formation and evolution. Seismic events detected by the nearside network were used to constrain the structure of the Moon s crust and mantle down to a depth of about 1000 km. The presence of an attenuating region in the deepest interior has been inferred from the paucity of farside events, as well as other indirect geophysical measurements. Recent re-analyses of the Apollo data have tentatively identified this region as a lunar core, although its properties are not yet constrained. Here we present new modeling in support of seismic missions that plan to build upon the knowledge of the Moon s interior gathered by Apollo. We have devised a method in which individual events can be linked to a known cluster using the observed S-P arrival time differences and azimuth to only two stations. Events can be further identified using each cluster's unique occurrence time signature

Weber, Renee↗

GRAIL Refinements to Lunar Seismic Structure

Gravity field measurements are perhaps the most numerous of the indirect observations relevant to the Moon's internal structure. Multiple recent missions have mapped the global lunar gravity field, each one improving upon the resolution of the last. The details of the Moon's deepest structure, including the parameters that define the lunar core, however, were still largely unaddressed by pre-GRAIL gravity measurements, which were not high enough resolution to resolve the tidal coefficients at a sufficient degree of accuracy. Current constraints on core size and state arise from other indirect measurements, including lunar laser ranging, magnetic induction studies, and analyses of elemental abundances in depth-derived mare basalts. These inferences vary widely, but when considered together with structure models derived from the seismic data gathered during the Apollo missions, a schematic of the lunar interior containing a partially molten deepest mantle layer overlying molten outer and solid inner core layers was obtained. Seismology provides the most direct constraints on the variables that govern the dynamic properties of the body. However, the GRAIL mission's high-resolution measurements of the lunar gravity field are being used to constrain the interior structure of the Moon using a "crust to core" approach. GRAIL's constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten) therefore complement seismic investigations. This work focuses on expanding our knowledge of the Moon's internal structure using joint gravity and seismic analyses, which will improve constraints on the deep lunar mantle and core.

Weber, Renee C.↗

The Lunar Geophysical Network Mission

Overarching Principles: Must be better than Apollo (coverage, duration, instrument performance); Learn from the Apollo experience. Lunar Geophysical Network (LGN) New Frontiers (NF)-class mission, as part of the NF-5 call. “This mission consists of several identical landers distributed across the lunar surface, each carrying geophysical instrumentation. The primary science objectives are to characterize the Moon’s internal structure, seismic activity, global heat flow budget, bulk composition, & magnetic field.” Global distribution of multiple stations. Each station should contain a seismometer, heat flow probe, electromagnetic sounder, laser retroreflector (lunar nearside). Each station must be long-lived (e.g., approximately10 years)to allow other stations (from other countries?) to be integrated with the anchor nodes to form the International Lunar Network. Why LGN? Planetary Science: Moon represents an end-member in planetary evolution (large small body, small rocky planet); Primary planetary differentiation preserved; Key to understanding terrestrial planet initial differentiation. Lunar Science: Heat flow probes yield crustal heat budget estimates; Combined with EMS (ElectroMagnetic Sounding), the temperature profile of the deep interior can be modeled along with mineralogy; Seismic and LLR (Lunar Laser Ranging) data also yield structure and compositional information of the lunar interior; High fidelity data from LGN would enhance the usefulness of the GRAIL (Gravity Recovery and Interior Laboratory) and SELENE (Selenological and Engineering Explorer) gravity data. Human Exploration: LGN must be established prior to renewed human lunar activity - we do not know the exact locations or causes of the shallow moonquakes (SMQs) - the largest magnitude seismic events recorded by Apollo (1 event per year of magnitude greater than or equal to 5); Establishing surface infrastructure near SMQ epicenters must be avoided.

Neal, C. R.↗