Lunar seismic exploration.
Seismic experimentation methods to be performed on moon, giving P and S wave travel time curves for models of lunar interior
SEISMIC WAVE↗
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.
A critical assessment of viscous models of trench topography and corner flow
Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.
Zhang, J.↗
A critical assessment of viscous models of trench topography and corner flow
Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.
Zhang, J.↗
Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology
We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.
Reches, Ze'ev↗
The mechanics of the South Iceland Seismic Zone
The mechanics of the South Iceland Seismic Zone is examined by means of boundary element modeling. The differences between the observed geometry of the major tectonically active features and that of a geometrically orthogonal ridge-transform system (which assumes that the upper part of the earth's crust is an elastic plate containing vertical cuts) are examined. It is suggested that north-south faults can accommodate transform deformation only if the faults are longer or more numerous than those observed so far. This is considered reasonable because earthquake surface rupture lengths are commonly less than the inferred fault length at depth. The South Iceland Seismic Zone is subject to sequences of large earthquake every 45-112 years. In comparing the seismic moment release derived from earthquake magnitudes with that predicted by the models, it is confirmed that the system of north-south faults can act as a transform fault.
Hackman, M. Christine↗
Crustal structure, geophysical models and contemporary tectonism of the Colorado Plateau
A regional analysis of the crust and upper mantle of the Colorado Plateau is presented, using existing geophysical and geological data combined with new surface wave dispersion and groundwater geothermometry data; the tectonic implications of these models are also investigated. Surface wave and seismic refraction data indicate that the crust of the interior of the Colorado Plateau is 44 + or - 3 km thick, and its crustal structure is typical of stable continental areas. Pn velocities, however, appear to be lower (7.8 km/s) than would be expected in a stable region, while silica geothermometry indicates that the average heat flow for the plateau is 55 mW per sq m (1.3 HFU).
Keller, G. R.↗
Absorption band Q model for the Earth
Attenuation in solids and liquids, as measured by the quality factor Q, is typically frequency dependent. In seismology, however, Q is usually assumed to be independent of frequency. Body wave, surface wave, and normal mode data are used to place constraints on the frequency dependence of Q in the mantle. Specific features of the absorption band model are: low-Q in the seismic band at both the top and the base of the mantle, low-Q for long-period body waves in the outer core, an inner core Q sub s that increases with period, and low Q sub p/Q sub s at short periods in the middle mantle.
Anderson, D. L.↗
The formation of antipodal-impact terrains on Mars
The regions antipodal to Mars' three largest impact basins, Hellas, Isidis, and Argyre, were assessed for evidence of impact induced seismically disrupted terrains. Photogeology and preliminary computer modeling using the Simplified Arbitrary Lagrangian-Eulerian computer program suggest such terrain could have been formed by the Hellas and Isidis impacts. Previously, the only assessment of potential surface disruptions on Mars antipodal to large impacts was by Peterson (1978). Herein, the results give credence to Peterson's assertions. For Hellas, the computer model clearly suggests that antipodal pressures were probably strong enough to have fractured the crust and disrupted the surface, and may account for aspects of volcanism at Alba Patera. It is suggested that the Hellas impact produced sufficient focused energy at its antipode to produce deep fracture in the Martian crust centered below the current caldera for Alba Petera. With the generation and evolution of magma in the Tharsis area, the fracture system provided a ready conduit for early stage eruptions of Martians lavas derived from deep in the crust and/or upper mantle, perhaps of komatiitic composition.
Williams, D. A.↗
Charge Generation and Propagation in Igneous Rocks
Resistivity changes, ground potentials, electromagnetic (EM) and luminous signals prior to or during earthquakes have been reported, in addition to ground uplift and tilt, and to changes in the seismic wave propagation parameters. However, no physical model exists that ties these diverse phenomena together. Through time-resolved impacts experiments it has been observed that, when igneous rocks (gabbro, diorite, granite) are impacted at low velocities (approx. 100 m/sec), highly mobile electronic charge carriers are generated, spreading from a small volume near the impact point, causing electric potentials, EM and light emission. The rock becomes momentarily conductive. When impacted at higher velocities (approx. 1.5 km/sec), the propagation of the P and S waves is registered through the transient piezoelectric response of quartz. At the same time, the rock volume is filled with mobile charge carriers, and a positive surface potential is registered. During the next 1-2 msec the surface potential oscillates, due to electron injection from ground. These observations are consistent with positive holes, e.g. defect electrons in the O(2-) sublattice, that can travel via the O 2p-dominated valence band of the silicate minerals at the speed of a phonon-mediated charge transfer. Before activation, the positive hole charge carriers lay dormant in form of positive hole pairs, PHP, electrically inactive, chemically equivalent to peroxy links in the structures of constituent minerals. PHPs are introduced by way of hydroxyl (O3Si-OH) incorporated into nominally anhydrous minerals when they crystallize in water-laden environments. Given that sound waves of even relatively low intensity appear to cause PHPs dissociation, thus generating mobile positive holes, it is proposed that microfracturing during rock deformation cause PHP dissociation. Depending on where and how much the rock volume is stressed, the positive holes are expected to form fluctuating charge clouds in the earthquake source region that may account for earthquake-related electrical phenomena and the reported low frequency EM signals.
Freund, Friedemann↗
Determining the Shallow Surface Velocity at the Apollo 17 Landing Site
Many studies have been performed to determine the shallow surface velocity model at the Apollo 17 landing site. The Lunar Seismic Profiling Experiment (LSPE) had both an active component with eight explosive packages (EPs) and a passive experiment collecting data at various time intervals. Using the eight EPs, the initial shallow surface velocity model was determined to be 250 m/s in the first layer of depth 248 m, 1200 m/s with a depth of 927 m in the second layer, and 4000 m/s down to a depth of 2 km in the third layer. Have performed variations on this study to produce new velocity models shown. Recent studies have also been reanalyzing the passive LSPE data and have found three different thermal moonquake event types occurring at different times within the lunar day. The current goal of the project is to collocate the thermal moonquakes to physical surface features to determine the breakdown of lunar rocks. However, to locate shallow surface events, an accurate velocity model is needed. Presented a thermal moonquake location algorithm using first order approximation, including surface events only. To improve these approximations, a shallow surface velocity is needed.
Phillips, D.↗
Planetary Seismology
Of the many geophysical means that can be used to probe a planet's interior, seismology remains the most direct. In addition to Earth, seismometers have been installed on Venus, Mars, and the Moon. Given that the seismic data gathered on the Moon (now over 40 years ago) revolutionized our understanding of the Moon and are still being used today to produce new insight into the state of the lunar interior, it is no wonder that many future missions, both real and conceptual, plan to take seismometers to other planets. To best facilitate the return of high-quality data from these instruments, as well as to further our understanding of the dynamic processes that modify a planet's interior, various modeling approaches are used to quantify parameters such as the amount and distribution of seismicity, tidal deformation, and seismic structure of the terrestrial planets. In addition, recent advances in wavefield modeling have permitted a renewed look at seismic energy transmission and the effects of attenuation and scattering, as well as the presence and effect of a core, on recorded seismograms. In this talk I will discuss some of these methods and review the history of planetary seismology.
planetary seismology↗
Axial focusing of impact energy in the Earth's interior: Proof-of-principle tests of a new hypothesis
A causal link between major impact events and global processes would probably require a significant change in the thermal state of the Earth's interior, presumably brought about by coupling of impact energy. One possible mechanism for such energy coupling from the surface to the deep interior would be through focusing due to axial symmetry. Antipodal focusing of surface and body waves from earthquakes is a well-known phenomenon which has previously been exploited by seismologists in studies of the Earth's deep interior. Antipodal focusing from impacts on the Moon, Mercury, and icy satellites has also been invoked by planetary scientists to explain unusual surface features opposite some of the large impact structures on these bodies. For example, 'disrupted' terrains have been observed antipodal to the Caloris impact basis on Mercury and Imbrium Basin on the Moon. Very recently there have been speculations that antipodal focusing of impact energy within the mantle may lead to flood basalt and hotspot activity, but there has not yet been an attempt at a rigorous model. A new hypothesis was proposed and preliminary proof-of-principle tests for the coupling of energy from major impacts to the mantle by axial focusing of seismic waves was performed. Because of the axial symmetry of the explosive source, the phases and amplitudes are dependent only on ray parameter (or takeoff angle) and are independent of azimuthal angle. For a symmetric and homogeneous Earth, all the seismic energy radiated by the impact at a given takeoff angle will be refocused (minus attenuation) on the axis of symmetry, regardless of the number of reflections and refractions it has experienced. Mantle material near the axis of symmetry will experience more strain cycles with much greater amplitude than elsewhere and will therefore experience more irreversible heating. The situation is very different than for a giant earthquake, which in addition to having less energy, has an asymmetric focal mechanism and a larger area. Two independent proof-of-principle approaches were used. The first makes use of seismic simulations, which are being performed with a realistic Earth model to determine the degree of focusing along the axis and to estimate the volume of material, if any, that experiences significant irreversible heating. The second involves two-dimensional hydrodynamic code simulations to determine the stress history, internal energy, and temperature rise as a function of radius along the axis.
Boslough, M. B.↗
Seismic Assessment of a Long-Term Lunar Habitat
The establishment of secure earth-independent long-term lunar habitats has been envisioned by numerous government agencies and private companies. Recent advancements in assessing seismic hazards caused by shallow moonquakes have highlighted the importance of incorporating this phenomenon into the design of robust and resilient lunar structures. However, further research is required to explore lunar habitat design that considers seismic loads. This paper proposes assessing the structural response of a lunar habitat made of sulfur concrete covered with a regolith layer. The numerical model of the structure is subjected to gravitational, internal pressure and seismic loads. The seismic analysis of the structure is carried out using spectral and nonlinear time history methods. Conditional mean spectra for shallow moonquakes with return periods of 75, 475, 970 and 2475 years are used in the seismic analysis. The records used for the temporal analyses were ground motions that agree with a preliminary seismic hazard on the Moon. The results of temporal analyses reveal that shallow moonquakes with return periods greater than 475 years can lead to the loss of the global stability of the structure. Consequently, the findings imply that seismic loads have the potential to impose unacceptable demands on lunar structures constructed from in-situ materials like sulfur concrete. Hence, it is imperative to incorporate seismic considerations in the design process for developing resilient and long-term lunar habitats.
Lunar habitat↗
Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) Model - An Unified Concept for Earthquake Precursors Validation
The paper presents a conception of complex multidisciplinary approach to the problem of clarification the nature of short-term earthquake precursors observed in atmosphere, atmospheric electricity and in ionosphere and magnetosphere. Our approach is based on the most fundamental principles of tectonics giving understanding that earthquake is an ultimate result of relative movement of tectonic plates and blocks of different sizes. Different kind of gases: methane, helium, hydrogen, and carbon dioxide leaking from the crust can serve as carrier gases for radon including underwater seismically active faults. Radon action on atmospheric gases is similar to the cosmic rays effects in upper layers of atmosphere: it is the air ionization and formation by ions the nucleus of water condensation. Condensation of water vapor is accompanied by the latent heat exhalation is the main cause for observing atmospheric thermal anomalies. Formation of large ion clusters changes the conductivity of boundary layer of atmosphere and parameters of the global electric circuit over the active tectonic faults. Variations of atmospheric electricity are the main source of ionospheric anomalies over seismically active areas. Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model can explain most of these events as a synergy between different ground surface, atmosphere and ionosphere processes and anomalous variations which are usually named as short-term earthquake precursors. A newly developed approach of Interdisciplinary Space-Terrestrial Framework (ISTF) can provide also a verification of these precursory processes in seismically active regions. The main outcome of this paper is the unified concept for systematic validation of different types of earthquake precursors united by physical basis in one common theory.
Pulinets, S.↗
Broadscale Postseismic Gravity Change Following the 2011 Tohoku-Oki Earthquake and Implication for Deformation by Viscoelastic Relaxation and Afterslip
The analysis of GRACE gravity data revealed post-seismic gravity increase by 6 micro-Gal over a 500 km scale within a couple of years after the 2011 Tohoku-Oki earthquake, which is nearly 40-50% of the co-seismic gravity change. It originates mostly from changes in the isotropic component corresponding to the M(sub rr) moment tensor element. The exponential decay with rapid change in a year and gradual change afterward is a characteristic temporal pattern. Both viscoelastic relaxation and afterslip models produce reasonable agreement with the GRACE free-air gravity observation, while their Bouguer gravity patterns and seafloor vertical deformations are distinctly different. The post-seismic gravity variation is best modeled by the bi-viscous relaxation with a transient and steady state viscosity of 10(exp 18) and 10(exp 19) Pa s, respectively, for the asthenosphere. Our calculated higher-resolution viscoelastic relaxation model, underlying the partially ruptured elastic lithosphere, yields the localized post-seismic subsidence above the hypocenter reported from the GPS-acoustic seafloor surveying.
GRACE↗
Comparison of Magnetic Fields from a New Satellite Magnetic Model of the Lithosphere with Magnetic Fields Predicted from Crust 1.0
Both magnetic and seismic techniques can provide information about the Moho (Mohorovicic discontinuity). We develop a new technique that provides a better estimate of the magnetic thickness of the crust, as compared with previous approaches. It uses prior knowledge from seismology (Crust 1.0), a new high-degree model from CHAMP (CHAllenging Mini-satellite Payload) and Swarm (LCS-1 - a model of Earth's lithospheric field) and a newly developed technique. The technique is appropriate for regions where induced magnetization dominates over remanent magnetization. We compare the predictions from LCS-1 with those from Crust 1.0, with some simple assumptions, and find that the correlations increase until about spherical harmonic degree 30, and then decrease globally. Spatially, the correlations between the seismic and magnetic techniques are strongest over North America and Australia, and weakest over South America and northern Africa. Strong correlations also exist between the two approaches over the Antarctic, northern Europe, and Greenland. While we might expect the seismic and magnetic approaches to correlate over well-characterized regions (i.e. North America), and show weaker correlations over poorly-characterized regions (i.e. South America and north Africa), the strong correlation in the Antarctic and Greenland is puzzling, because both of these regions are poorly-characterized. We discuss some possible explanations, and implications, of this attempt to correlate seismic and magnetic approaches to characterizing the lithosphere.
Seismic↗
Glacio-Seismotectonics: Ice Sheets, Crustal Deformation and Seismicity
The last decade has witnessed a significant growth in our understanding of the past and continuing effects of ice sheets and glaciers on contemporary crustal deformation and seismicity. This growth has been driven largely by the emergence of postglacial rebound models (PGM) constrained by new field observations that incorporate increasingly realistic rheological, mechanical, and glacial parameters. In this paper, we highlight some of these recent field-based investigations and new PGMs, and examine their implications for understanding crustal deformation and seismicity during glaciation and following deglaciation. The emerging glacial rebound models outlined in the paper support the view that both tectonic stresses and glacial rebound stresses are needed to explain the distribution and style of contemporary earthquake activity in former glaciated shields of eastern Canada and Fennoscandia. However, many of these models neglect important parameters, such as topography, lateral variations in lithospheric strength and tectonic strain built up during glaciation. In glaciated mountainous terrains, glacial erosion may directly modulate tectonic deformation by resetting the orogenic topography and thereby providing an additional compensatory uplift mechanism. Such effects are likely to be important both in tectonically active orogens and in the mountainous regions of glaciated shields.
Sauber, Jeanne↗
Model history of the lunar surface
Lunar surface history model, discussing mascons, chemical composition, isochron ages and seismic and electrical properties of samples
Hawkins, J. W.↗