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 109 records · Page 6

Internal structure and properties of Mars

Theoretical physical models of the Martian interior are presented in the light of recent and revised data and constraints. These models include thermal evolution, densities, and seismic-wave velocities. The interior of Mars appears to be earthlike in many respects. Although thermal models indicate that Mars has passed its peak of evolution, it may still have an asthenosphere and may be moderately active tectonically. Mars has an Fe-FeS core with a radius of 1500-2000 km. The mantle is enriched in FeO with an olivine composition of about Fo75. Theoretically determined seismic-wave velocities are relatively well constrained in the mantle, with upper-mantle P-wave velocities ranging from 7.64 to 7.80 km/sec. However, there are wide variations in P-wave velocity in the core, dependent on composition. The shadow zone due to the core is larger than earth's.

Johnston, D. H.↗

The deep structure of lunar basins - Implications for basin formation and modification

Models for the crustal structure in the vicinity of nine impact basins, from an inversion of gravity and topographic data from the lunar nearside are presented. The models display a low-density nonmare crustal layer and a mare basalt layer, both of variable thicknesses. Assuming that topography in mare areas is isostatically compensated before the emplacement of mare basalts and that compensation of mare basalt units may be neglected, a decomposition of the gravity anomaly into contributions from Moho relief and mare fill is permitted. Minimum values for mare basalt thicknesses are obtained but because mare basalts and mantle material are similar in density, the thicknesses of the nonmare crust are estimated. An important constraint is the crustal thickness inferred from the Apollo 12 and 14 landing sites from seismic observations. The crustal thickness model indicates that the crust is thinner beneath each of the major nearside basins than in surrounding areas. New bounds on the volume of material ejected from each basin are derived. The geological implications of structural differences among basins for the processes of basin formation and modification are evaluated as functions of time on the moon.

Bratt, S. R.↗

The 26 December 2004 tsunami source estimated from satellite radar altimetry and seismic waves

The 26 December 2004 Indian Ocean tsunami was the first earthquake tsunami of its magnitude to occur since the advent of both digital seismometry and satellite radar altimetry. Both have independently recorded the event from different physical aspects. The seismic data has then been used to estimate the earthquake fault parameters, and a three-dimensional ocean-general-circulation-model (OGCM) coupled with the fault information has been used to simulate the satellite-observed tsunami waves. Here we show that these two datasets consistently provide the tsunami source using independent methodologies of seismic waveform inversion and ocean modeling. Cross-examining the two independent results confirms that the slip function is the most important condition controlling the tsunami strength, while the geometry and the rupture velocity of the tectonic plane determine the spatial patterns of the tsunami.

ocean general circulation model (OGCM)↗

Dynamically supported geoid highs over hotspots: Observation and theory

Hotspots are associated with long wavelength geoid highs, an association that is even stronger when the geoid highs associated with subduction zones are removed. These associations are quantified by expanding the hotspot distribution in spherical harmonics and calculating correlation coefficients as a function of harmonic degree. The hotspot distribution spectrum is essentially white, with peaks at degrees 2 and 6. It is correlated positively with the slab residual geoid for degrees 2 to 6, with low seismic velocity in the lower mantle at degree 2, and with low seismic velocity in the upper mantle at degree 6. A variety of fluid mechanical models were tested for hotspots, including lithospheric delamination and hot plumes, by calculating their predicted dynamic geoid responses and comparing them to the observations. These models include the effects of temperature dependent rheology. The preferred hotspot model, based on observations of the geoid and seismic tomography, has plumes preferentially occurring in regions of large scale background temperature highs in a mantle with substantial viscosity increase with depth, although other models are possible.

Richards, M. A.↗

Dynamically supported geoid highs over hotspots - Observation and theory

Hotspots are associated with long wavelength geoid highs, an association that is even stronger when the geoid highs associated with subduction zones are removed. These associations are quantified by expanding the hotspot distribution in spherical harmonics and calculating correlation coefficients as a function of harmonic degree. The hotspot distribution spectrum is essentially white, with peaks at degrees 2 and 6. It is correlated positively with the slab residual geoid for degrees 2 to 6, with low seismic velocity in the lower mantle at degree 2, and with low seismic velocity in the upper mantle at degree 6. A variety of fluid mechanical models were tested for hotspots, including lithospheric delamination and hot plumes, by calculating their predicted dynamic geoid responses and comparing them to the observations. These models include the effects of temperature dependent rheology. The preferred hotspot model, based on observations of the geoid and seismic tomography, has plumes preferentially occurring in regions of large scale background temperature highs in a mantle with substantial viscosity increase with depth, although other models are possible.

Richards, Mark A.↗

Magnetization models for the source of the 'Kentucky anomaly' observed by Magsat

Both the aeromagnetic data and magnetic anomaly data obtained by Magsat indicate the presence of a very magnetic source region within the crust beneath Kentucky and Tennessee. A source model was previously developed to fit surface gravity and long-wavelength aeromagnetic data, using limited seismic constraint. For the present study the model was further developed, and it is demonstrated that the source region for the satellite anomaly is considerably more extensive than the Kentucky body sensu stricto. The extended source region is modeled using both prismatic model sources and dipole array sources. Magnetization directions for the source region found by inversion of various combinations of scalar and vector data are found to be close to the main field direction, implying the lack of a strong remanent component. It is shown by simulation that in a case (such as this) where the geometry of the source is known, if a strong remanent component is present its direction is determinable, but by scalar data as readily as vector data. Magnetization magnitude for the extended source region is about 3 A/m if the vertical extent of the source includes the whole of the crust.

Mayhew, M. A.↗

A kinematic model for the late Cenozoic development of southern California crust and upper mantle

A model is developed for the young and ongoing kinematic deformation of the southern California crust and upper mantle. The kinematic model qualitatively explains both the overall seismic structure of the upper mantle and much of the known geological history of the late Cenozoic as consequences of ongoing convection beneath southern California. In this model, the high-velocity upper-mantle anomaly of the Transverse ranges is created through the convergence and sinking of the entire thickness of subcrustal lihtosphere, and the low-velocity upper-mantle anomaly beneath the Salton Trough region is attributed to high temperatures and 1-4 percent partial melt related to adiabatic decompression during mantle upwelling.

Humphreys, Eugene D.↗

Origin of multi-ring basin ridge systems - An upper limit to elastic deformation based on a finite-element model

A finite-element model of a 500-km diameter lunar basin flooded with 10 km of basalt was used to represent the central portion of a multi-ring basin undergoing elastic deformation as a result of basin subsidence. Gravitational forces due to the weight of the basalt were set as the initial forces in the model; elastic constants computed from seismic velocities of lunar samples provided strength specifications. In addition, the crust was constrained by permitting no displacements at a depth of 60 km. Because of the strength of the lithosphere, surface subsidence amounted to only 20 m. This magnitude of elastic deformation does not equal that found on the lunar surface, but the orientation of principal stress and the position of circumferential compressional stress are compatible with the evolution of inner-ring and radial ridge systems by gravitational-induced downdropping.

Maxwell, T. A.↗

Information theory lateral density distribution for Earth inferred from global gravity field

Information Theory Inference, better known as the Maximum Entropy Method, was used to infer the lateral density distribution inside the Earth. The approach assumed that the Earth consists of indistinguishable Maxwell-Boltzmann particles populating infinitesimal volume elements, and followed the standard methods of statistical mechanics (maximizing the entropy function). The GEM 10B spherical harmonic gravity field coefficients, complete to degree and order 36, were used as constraints on the lateral density distribution. The spherically symmetric part of the density distribution was assumed to be known. The lateral density variation was assumed to be small compared to the spherically symmetric part. The resulting information theory density distribution for the cases of no crust removed, 30 km of compensated crust removed, and 30 km of uncompensated crust removed all gave broad density anomalies extending deep into the mantle, but with the density contrasts being the greatest towards the surface (typically + or 0.004 g cm 3 in the first two cases and + or - 0.04 g cm 3 in the third). None of the density distributions resemble classical organized convection cells. The information theory approach may have use in choosing Standard Earth Models, but, the inclusion of seismic data into the approach appears difficult.

Rubincam, D. P.↗

Terrestrial environment (climatic) criteria guidelines for use in aerospace vehicle development

Guidelines on terrestrial environment data specifically applicable for NASA aerospace vehicles and associated equipment development are provided. The general distribution of natural environmental extremes in the conterminous United States that may be needed to specify design criteria in the transportation of space vehicle subsystems and components is considered. Atmospheric attenuation was included, since certain Earth orbital experiment missions are influenced by the Earth's atmosphere. Climatic extremes for worldwide operational needs is also included. Atmospheric chemistry, seismic criteria, and a mathematical model to predict atmospheric dispersion of aerospace engine exhaust cloud rise and growth are discussed. Atmospheric cloud phenomena are considered.

Turner, R. E.↗

Estimating the Libration of Mercury by Remote Sensing of Gravity and Altimetry

One of the most intriguing scientific questions about Mercury is whether the planet possesses a liquid core. One way that this can be ascertained is through the measurement of the longitudinal libration of Mercury about its rotation axis. The MESSENGER (Mercury Surface, Space, Environment, Geochemistry, and Ranging) mission to Mercury will attempt this measurement by deriving the rotation of the planet from remotely sensed observations of the Mercury gravity field and the planet's shape. For a fully fluid core the libration is expected to be an approximately 350-m-amplitude oscillation (at the equator) about the mean rotation of the planet. Thus, the shape of the planet undergoes a librational oscillation in the same way that features on the surface participate in the motion. In this sense altimetric topography can be considered analogous to imaging, with the important exception that for altimetry, the long wavelengths corresponding to low degrees and orders in a spherical harmonic expansion, can be precisely determined from orbital measurements. The rotation of the planet thus introduces a time dependence of the topographic image in an inertial (i.e. fixed) reference frame and the libration is a variation in the planetary rotation rate, manifest as an oscillation of the topographic image. Similarly, gravity measures the distribution of mass within a planet and as such can be viewed as providing a vertically integrated 'image' of the internal density structure. For a librating planet the gravitational field will also reveal a variation in the rotation rate, manifest as an oscillation of the field with respect to an inertial reference frame. While the long-wavelength gravity field contains contributions from the radial distribution of mass with depth, the long wavelength terms are most sensitive to the mass distribution of the deep interior. Thus, both the topography and the gravity independently can be used to determine the (irregular) rotation of a planet, although the topographical method measures the rotation of the crust while the gravity measures the rotation of the distributed mass, and these are not necessarily identical. For example, for a planet that contains a fluid outer core that is effectively decoupled from the mantle the gravity field will not librate in the manner of a solid body. The libration obtained from the variable rotation of the gravity field will represent a combination of contributions due to a librating mantle and a differentially rotating fluid core. Differential rotation of the Earth's inner core had been theoretically predicted from three dimensional geodynamical models and subsequently reported from seismic observations. A difference between the libration of the lithosphere as determined from altimetry and the deep interior as determined from gravity, if it could be detected, might provide further insight into the nature of core-mantle coupling. Additional information is contained in the original extended abstract.

Smith, D. E.↗

Modeling surface deformations at complex strike-slip plate boundaries

Two new versions of the physical base-traction model of Li and Rice (1987) have been developed to study the loading processes and surface deformation near plate boundaries with geometric complexities, including a shallow creeping fault segment and two subparallel faults. Predictions of surface velocity and surface slip rate for the region between the San Andreas fault and the Pacific coast at the latitude of the Parkfield-Cholame segment using model parameters based on geologic, geodetic, and seismic considerations were found to be in good agreement with contemporary geodetic field data (not used in constrainning model parameters). Also, predicted surface velocity for the Salton Sea-Coachella Valley area were fit to recently derived contemporary geodetic data and very large baseline interferometry data, althoug here the model parameters are not well constrained.

Li, Victor C.↗

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

Lunar and Planetary Science XXXV: Asteroids, Meteors, and Comets

The session "Asteroids, Meteors, and Comets" contained the following reports:Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions; Macroscopic Voids in Small Asteroids: Effects of Cohesion; The Seismic Effect of Impacts on Asteroid Surface Morphology: Early Modeling Results; Photometric Studies of Eros from NEAR Data; Quantitative Aspects of Space Weathering: Implications for Regolith Breccia Meteorites and Asteroids; Diversity of Types of Hydrated Minerals on C-Class Asteroids; Mineralogical Variations Among High Albedo E-Type Asteroids: Implications for Asteroid Igneous Processes; Multi-Wavelength Observations of 2100 Ra-Shalom: Radar and Lightcurves; What are the P-type Asteroids Made Of?; Sodium Overabundance in Meteoroids from Meteor Spectroscopy; Migration Processes and Volatiles Inventory to the Inner Planets; Characterization of the Surface Properties of MUSES-C/Hayabusa Spacecraft Target; Asteroid 25143 Itokawa; and Sample Return Science by Hayabusa Near-Earth Asteroid Mission.

Source record↗

Crustal Structure of the Iceland Region from Spectrally Correlated Free-air and Terrain Gravity Data

Seismic refraction studies have provided critical, but spatially restricted constraints on the structure of the Icelandic crust. To obtain a more comprehensive regional view of this tectonically complicated area, we spectrally correlated free-air gravity anomalies against computed gravity effects of the terrain for a crustal thickness model that also conforms to regional seismic and thermal constraints. Our regional crustal thickness estimates suggest thickened crust extends up to 500 km on either side of the Greenland-Scotland Ridge with the Iceland-Faeroe Ridge crust being less extended and on average 3-5 km thinner than the crust of the Greenland-Iceland Ridge. Crustal thickness estimates for Iceland range from 25-35 km in conformity with seismic predictions of a cooler, thicker crust. However, the deepening of our gravity-inferred Moho relative to seismic estimates at the thermal plume and rift zones of Iceland suggests partial melting. The amount of partial melting may range from about 8% beneath the rift zones to perhaps 20% above the plume core where mantle temperatures may be 200-400 C above normal. Beneath Iceland, areally limited regions of partial melting may also be compositionally and mechanically layered

Leftwich, T. E.↗

Crustal Structure of the Iceland Region from Spectrally Correlated Free-air and Terrain Gravity Data

Seismic refraction studies have provided critical, but spatially restricted constraints on the structure of the Icelandic crust. To obtain a more comprehensive regional view of this tectonically complicated area, we spectrally correlated free-air gravity anomalies against computed gravity effects of the terrain for a crustal thickness model that also conforms to regional seismic and thermal constraints. Our regional crustal thickness estimates suggest thickened crust extends up to 500 km on either side of the Greenland-Scotland Ridge with the Iceland-Faeroe Ridge crust being less extended and on average 3-5 km thinner than the crust of the Greenland-Iceland Ridge. Crustal thickness estimates for Iceland range from 25-35 km in conformity with seismic predictions of a cooler, thicker crust. However, the deepening of our gravity-inferred Moho relative to seismic estimates at the thermal plume and rift zones of Iceland suggests partial melting. The amount of partial melting may range from about 8% beneath the rift zones to perhaps 20% above the plume core where mantle temperatures may be 200-400 C above normal. Beneath Iceland, areally limited regions of partial melting may also be compositionally and mechanically layered and intruded. The mantle plume appears to be centered at (64.6 deg N, 17.4 deg W) near the Vatnajokull Glacier and the central Icelandic neovolcanic zones.

Leftwich, T. E.↗

Crustal Deformation and the Seismic Cycle across the Kodiak Islands, Alaska

The Kodiak Islands are located approximately 130 to 250 km from the Alaska-Aleutian Trench where the Pacific plate is underthrusting the North American plate at a rate of about 57 mm/yr. The southern extent of the 1964 Prince William Sound (${M-w}$ = 9.2) earthquake rupture occurred offshore and beneath the eastern portion of the Kodiak Islands. Here we report GPS results (1993-2001) from northern Kodiak Island that span the transition between the 1964 uplift region along the eastern coast and the region of coseismic subsidence further inland. The horizontal velocity vectors range from 22.9 $\pm$ 2.2 mm/yr at N26.3$\deg$W $\pm$ 2.5$\deg$, about 150 km from the trench, to 5.9 $\pm$ 1.3 mm/yr at N65.9$\deg$W $\pm$ 6.6$\deg$, about 190 km from the trench. Near the northeastern coast of Kodiak the velocity vector above the shallow, locked main thrust zone is between the orientation of PCFC-NOAM plate motion (N22$/deg$W) and the trench-normal (N3O$\deg$W). Further west, our geodetic results suggest the accumulation of shear strain that will be released eventually as left-lateral motion on upper plate faults such as the Kodiak Island fault. These results are consistent with the hypothesis that the difference between the Pacific-North American plate motion and the orientation of the down going slab would lead to 4-8 mm/yr of left-lateral slip. Short-term geodetic uplift rates range from 2 - 14 mm/yr, with the maximum uplift located near the axis of maximum subsidence during the 1964 earthquake. We evaluated alternate interseismic models for Kodiak to test the importance of various mechanisms responsible for crustal deformation rates. These models are based on the plate interface slip history inferred from earlier modeling of coseismic and post-seismic geodetic results. The horizontal (trench perpendicular) and vertical deformation rates across Kodiak are consistent with a model that includes the viscoelastic response to : (1) a downgoing Pacific plate interface that is locked at shallow depths,(2) coseismic slip in the 1964 and (3) interseismic creep below the seismogenic zone. The change in orientation of the horizontal velocity vector occurs down-dip from the locked main thrust zone. In southern Kodiak, the coseismic slip in the 1964 earthquake was smaller than in the northern Kodiak region; yet, the horizontal, interseismic velocities as a function of distance from the trench are comparable to those in northern Kodiak. Based on the earthquake history prior to, and following the 1964 earthquake, we hypothesize that the plate interface in southern Kodiak slips in more frequent large earthquakes than in northern Kodiak.

Sauber, Jeanne↗