Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earthquakes”

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 973 records · Page 54

Geodetic Observations of Interseismic Strain Segmentation at the Sumatra Subduction Zone

Deformation above the Sumatra subduction zone, revealed by Global Positioning System (GPS) geodetic surveys, shows nearly complete coupling of the forearc to the subducting plate south of 0.5 deg S and half as much to north. The abrupt change in plate coupling coincides with the boundary between the rupture zones of the 1833 and 1861 (Mw greater than 8) thrust earthquakes. The rupture boundary appears as an abrupt change in strain accumulation well into the interseismic cycle, suggesting that seismic segmentation is controlled by properties of the plate interface that persist occupied through more than one earthquake cycle. Structural evidence indicates that differences in basal shear stress may be related to elevated pore pressure in the north.

Prawirodirdjo, L.↗

Statistical Properties of SGR 1900+14 Bursts

We study the statistics of soft gamma repeater (SGR) bursts using a database of 187 events detected with BATSE and 837 events detected with the Rossi X-Ray Timing Explorer Proportional Counter Array: all events are from SGR 1900+14 during its 1998-1999 active phase. We find that the fluence or energy distribution of bursts is consistent with a power law of index 1.66, over 4 orders of magnitude. This scale-free distribution resembles the Gutenberg-Richter law for earthquakes and gives evidence for self-organized criticality in SGRS. The distribution of time intervals between successive bursts from SGR 1900+14 is consistent with a lognormal distribution. There is no correlation between burst intensity and the waiting times till the next burst, but there is some evidence for a correlation between burst intensity and the time elapsed since the previous burst. We also find a correlation between the duration and the energy of the bursts, but with significant scatter. In all these statistical properties, SGR bursts resemble earthquakes and solar flares more closely than they resemble any known accretion-powered or nuclear-powered phenomena. Thus, our analysis lends support to the hypothesis that the energy source for SGR bursts is internal to the neutron star and plausibly magnetic.

Gogus, Ersin↗

Ionospheric Signature of Surface Mine Blasts from Global Positioning System Measurements

Sources such as atmospheric or buried explosions and shallow earthquakes are known to produce infrasonic pressure waves in the atmosphere. Because of the coupling between neutral particles and electrons at ionospheric altitudes, these acoustic and gravity waves induce variations of the ionospheric electron density. The Global Positioning System (GPS) provides a way of directly measuring the total electron content in the ionosphere and, therefore, of detecting such perturbations in the upper atmosphere. In July and August 1996, three large surface mine blasts (1.5 Kt each) were detonated at the Black Thunder coal mine in eastern Wyoming. As part of a seismic and acoustic monitoring- experiment, we deployed five dual-frequency GPS receivers at distances ranging from 50 to 200 km from the mine and were able to detect the ionospheric perturbation caused by the blasts. The perturbation starts 10 to 15 min after the blast, lasts for about 30 min, and propagates with an apparent horizontal velocity of 1200 meters per second. Its amplitude reaches 3 x 10 (exp 14) el per square meters in the 7-3 min period band, a value close to the ionospheric perturbation caused by the M = 6.7 Northridge earthquake. The small signal-to-noise ratio of the perturbation can be improved by slant-stacking the electron content time-series recorded by the different GPS receivers taking into account the horizontal propagation of the perturbation. The energy of the perturbation is concentrated in the 200 to 300 second period band, a result consistent with previous observations and numerical model predictions. The 300 second band probably corresponds to gravity modes and shorter periods to acoustic modes, respectively. Using a 1-D stratified velocity model of the atmosphere we show that linear acoustic ray tracing fits arrival times at all GPS receivers. We interpret the perturbation as a direct acoustic wave caused by the explosion itself. This study shows that even relatively small subsurface events can produce ionospheric perturbations that are above the detection threshold of the GPS technique. By sensing derivative signals, which can be detected over a relatively broad region, it appears that GPS might be particularly useful for source characterization within the first acoustic quiet zone where infrasound would probably be ineffective. This suggests that dual-frequency GPS monitoring could contribute to Comprehensive Nuclear Test Ban Treaty verification.

Calais, Eric↗

Environmental Perturbations Caused by the Impacts of Comets and Asteroids on Earth

The extinction mechanisms proposed at the Cretaceous-Tertiary geological boundary are reviewed and related to the impact of asteroids or comets in general. For impact energies below 10(exp 4) Megatons (less than 6 x 10(exp 4) years; asteroid diameter less than 650 m), blast, earthquake, and fire may destroy local areas up to 10(exp 5) square m. Tidal waves could flood a kilometer inland over entire ocean basins. The energy range from 105 to 106 Megatons (less than 2 x 10(exp 6) years; asteroid diameter less than 3 km) is transitional. Dust lifted, sulfur released from within impacting asteroids, and soot from fires started by comets can produce climatologically significant optical depths of 10. At energies beyond 10(exp 7) Megatons, blast and earthquake damage is regional (10(exp 6) square cm). Tsunami cresting to 100 m and flooding 20 km inland will sweep the coastal zones of the world's oceans. Fires will be set globally. Light levels may drop so low from the smoke, dust and sulfate that vision is not possible. At energies approaching 10(exp 9) Megatons the ocean surface waters may be acidified by sulfur. The combination of these effects would be devastating.

Toon, Owen B.↗

Feeling Gravity's Pull: Gravity Modeling. The Gravity Field of Mars

Most people take the constant presence of gravitys pull for granted. However, the Earth's gravitational strength actually varies from location to location. This variation occurs because mass, which influences an object's gravitational pull, is not evenly distributed within the planet. Changes in topography, such as glacial movement, an earthquake, or a rise in the ocean level, can subtly affect the gravity field. An accurate measurement of the Earth's gravity field helps us understand the distribution of mass beneath the surface. This insight can assist us in locating petroleum, mineral deposits, ground water, and other valuable substances. Gravity mapping can also help notice or verify changes in sea surface height and other ocean characteristics. Such changes may indicate climate change from polar ice melting and other phenomena. In addition, gravity mapping can indicate how land moves under the surface after earthquakes and other plate tectonic processes. Finally, changes in the Earth's gravity field might indicate a shift in water distribution that could affect agriculture, water supplies for population centers, and long-term weather prediction. Scientists can map out the Earth's gravity field by watching satellite orbits. When a satellite shifts in vertical position, it might be passing over an area where gravity changes in strength. Gravity is only one factor that may shape a satellite's orbital path. To derive a gravity measurement from satellite movement, scientists must remove other factors that might affect a satellite's position: 1. Drag from atmospheric friction. 2. Pressure from solar radiation as it heads toward Earth and. as it is reflected off the surface of the Earth 3. Gravitational pull from the Sun, the Moon, and other planets in the Solar System. 4. The effect of tides. 5. Relativistic effects. Scientists must also correct for the satellite tracking process. For example, the tracking signal must be corrected for refraction through the atmosphere of the Earth. Supercomputers can calculate the effect of gravity for specific locations in space following a mathematical process known as spherical harmonics, which quantifies the gravity field of a planetary body. The process is based on Laplace's fundamental differential equation of gravity. The accuracy of a spherical harmonic solution is rated by its degree and order. Minute variations in gravity are measured against the geoid, a surface of constant gravity acceleration at mean sea level. The geoid reference gravity model strength includes the central body gravitational attraction (9.8 m/sq s) and a geopotential variation in latitude partially caused by the rotation of the Earth. The rotational effect modifies the shape of the geoid to be more like an ellipsoid, rather than a perfect, circle. Variations of gravity strength from the ellipsoidal reference model are measured in units called milli-Galileos (mGals). One mGal equals 10(exp -5) m/sq s. Research projects have also measured the gravity fields of other planetary bodies, as noted in the user profile that follows. From this information, we may make inferences about our own planet's internal structure and evolution. Moreover, mapping the gravity fields of other planets can help scientists plot the most fuel-efficient course for spacecraft expeditions to those planets.

Lemoine, Frank↗

Fault Creep along the Southern San Andreas from Interferometric Synthetic Aperture Radar, Permanent Scatterers, and Stacking

Interferometric synthetic aperture radar (InSAR) provides a practical means of mapping creep along major strike-slip faults. The small amplitude of the creep signal (less than 10 mm/yr), combined with its short wavelength, makes it difficult to extract from long time span interferograms, especially in agricultural or heavily vegetated areas. We utilize two approaches to extract the fault creep signal from 37 ERS SAR images along the southem San Andreas Fault. First, amplitude stacking is utilized to identify permanent scatterers, which are then used to weight the interferogram prior to spatial filtering. This weighting improves correlation and also provides a mask for poorly correlated areas. Second, the unwrapped phase is stacked to reduce tropospheric and other short-wavelength noise. This combined processing enables us to recover the near-field (approximately 200 m) slip signal across the fault due to shallow creep. Displacement maps fiom 60 interferograms reveal a diffuse secular strain buildup, punctuated by localized interseismic creep of 4-6 mm/yr line of sight (LOS, 12-18 mm/yr horizontal). With the exception of Durmid Hill, this entire segment of the southern San Andreas experienced right-lateral triggered slip of up to 10 cm during the 3.5-year period spanning the 1992 Landers earthquake. The deformation change following the 1999 Hector Mine earthquake was much smaller (4 cm) and broader than for the Landers event. Profiles across the fault during the interseismic phase show peak-to-trough amplitude ranging from 15 to 25 mm/yr (horizontal component) and the minimum misfit models show a range of creeping/locking depth values that fit the data.

Lyons, Suzanne↗

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

Accuracy of High-Rate GPS for Seismology

We built a device for translating a GPS antenna on a positioning table to simulate the ground motions caused by an earthquake. The earthquake simulator is accurate to better than 0.1 mm in position, and provides the "ground truth" displacements for assessing the technique of high-rate GPS. We found that the root-mean-square error of the 1-Hz GPS position estimates over the 15-min duration of the simulated seismic event was 2.5 mm, with approximately 96% of the observations in error by less than 5 mm, and is independent of GPS antenna motion. The error spectrum of the GPS estimates is approximately flicker noise, with a 50% decorrelation time for the position error of approx.1.6 s. We that, for the particular event simulated, the spectrum of dependent error in the GPS measurements. surface deformations exceeds the GPS error spectrum within a finite band. More studies are required to determine whether a generally optimal bandwidth exists for a target group of seismic events.

Elosegui, P.↗

Tool for Viewing Faults Under Terrain

Multi Surface Light Table (MSLT) is an interactive software tool that was developed in support of the QuakeSim project, which has created an earthquake- fault database and a set of earthquake- simulation software tools. MSLT visualizes the three-dimensional geometries of faults embedded below the terrain and animates time-varying simulations of stress and slip. The fault segments, represented as rectangular surfaces at dip angles, are organized into collections, that is, faults. An interface built into MSLT queries and retrieves fault definitions from the QuakeSim fault database. MSLT also reads time-varying output from one of the QuakeSim simulation tools, called "Virtual California." Stress intensity is represented by variations in color. Slips are represented by directional indicators on the fault segments. The magnitudes of the slips are represented by the duration of the directional indicators in time. The interactive controls in MSLT provide a virtual track-ball, pan and zoom, translucency adjustment, simulation playback, and simulation movie capture. In addition, geographical information on the fault segments and faults is displayed on text windows. Because of the extensive viewing controls, faults can be seen in relation to one another, and to the terrain. These relations can be realized in simulations. Correlated slips in parallel faults are visible in the playback of Virtual California simulations.

Siegel, Herbert, L.↗

Adaptive Sampling of Time Series During Remote Exploration

This work deals with the challenge of online adaptive data collection in a time series. A remote sensor or explorer agent adapts its rate of data collection in order to track anomalous events while obeying constraints on time and power. This problem is challenging because the agent has limited visibility (all its datapoints lie in the past) and limited control (it can only decide when to collect its next datapoint). This problem is treated from an information-theoretic perspective, fitting a probabilistic model to collected data and optimizing the future sampling strategy to maximize information gain. The performance characteristics of stationary and nonstationary Gaussian process models are compared. Self-throttling sensors could benefit environmental sensor networks and monitoring as well as robotic exploration. Explorer agents can improve performance by adjusting their data collection rate, preserving scarce power or bandwidth resources during uninteresting times while fully covering anomalous events of interest. For example, a remote earthquake sensor could conserve power by limiting its measurements during normal conditions and increasing its cadence during rare earthquake events. A similar capability could improve sensor platforms traversing a fixed trajectory, such as an exploration rover transect or a deep space flyby. These agents can adapt observation times to improve sample coverage during moments of rapid change. An adaptive sampling approach couples sensor autonomy, instrument interpretation, and sampling. The challenge is addressed as an active learning problem, which already has extensive theoretical treatment in the statistics and machine learning literature. A statistical Gaussian process (GP) model is employed to guide sample decisions that maximize information gain. Nonsta tion - ary (e.g., time-varying) covariance relationships permit the system to represent and track local anomalies, in contrast with current GP approaches. Most common GP models are stationary, e.g., the covariance relationships are time-invariant. In such cases, information gain is independent of previously collected data, and the optimal solution can always be computed in advance. Information-optimal sampling of a stationary GP time series thus reduces to even spacing, and such models are not appropriate for tracking localized anomalies. Additionally, GP model inference can be computationally expensive.

Thompson, David R.↗

Observing Traveling Ionospheric Disturbances Caused by Tsunamis Using GPS TEC Measurements

Ground-based Global Positioning System (GPS) measurements of ionospheric Total Electron Content (TEC) show variations consistent with atmospheric internal gravity waves caused by ocean tsunamis following two recent seismic events: the American Samoa earthquake of September 29, 2009, and the Chile earthquake of February 27, 2010. Fluctuations in TEC correlated in time, space, and wave properties with these tsunamis were observed in TEC estimates processed using JPL's Global Ionospheric Mapping Software. These TEC estimates were band-pass filtered to remove ionospheric TEC variations with wavelengths and periods outside the typical range of internal gravity waves caused by tsunamis. Observable variations in TEC appear correlated with the tsunamis in certain locations, but not in others. Where variations are observed, the typical amplitude tends to be on the order of 1% of the background TEC value. Variations with amplitudes ~ 0.1 - 0.2 TECU are observable with periods and timing affiliated with the tsunami. These observations are compared to estimates of expected tsunami-driven TEC variations produced by Embry Riddle Aeronautical University's Spectral Full Wave Model, an atmosphere-ionosphere coupling model, and found to be in good agreement in some locations, though there are cases when the model predicts an observable tsunami-driven signature and none is observed. These TEC variations are not always seen when a tsunami is present, but in these two events the regions where a strong ocean tsunami was observed did coincide with clear TEC observations, while a lack of clear TEC observations coincided with smaller tsunami amplitudes. There exists the potential to apply these detection techniques to real-time GPS TEC data, providing estimates of tsunami speed and amplitude that may be useful for early warning systems.

Galvan, David A.↗

Deep Learning Models for Planetary Seismicity Detection

Research in planetary seismology is fundamentally constrained by a lack of data. Seismo-logical science products of future missions can typically only be informed by theoretical signal/noise characteristics of the environment or likely Earth-analogues. Although objectives can be re-assessed after some initial data-collection upon lander arrival, transfer of high-resolution data back to Earth is costly on lander power usage. Over the last several years, development of GPU computing techniques and open-source high-level APIs have led to rapid advances in deep learning within the fields of computer vision, natural language processing, and collaborative filtering. These techniques are actively being adapted in seismology for a variety of tasks, including: earthquake detection, seismic phase discrimination, and ground-motion prediction. Until the recent detection of mars quakes during the Mars InSight mission, the only other measurements of seismicity recorded outside of Earth was on the Moon during the Apollo missions between 1969 to 1977. These unique data sets have been periodically revisited using new seismological methods, including ambient noise interferometry and Hidden Markov Models. Our objective is to develop a deep learning seismic detector and use it to catalog moonquakes from the Apollo 17 Lunar Seismic Profiling Experiment (LSPE) and compare the results with those obtained by other methods. Additionally, we will assess the accuracy tradeoff between using a training set of lunar data and one composed of Earth seismicity. In this document, we present preliminary results using a prototype classifier trained on a small set of earthquakes that was able to obtain detections for LSPE moonquakes with a greater accuracy than a recent study using Hidden Markov Models.

Civilini, F.↗

TPSAS-NF1676L-13211-DND

On March 11, 2011 the Great Tohoku Earthquake occurred approximately 70 kilometers off the coast of Japan. This magnitude 9.0 earthquake was closely followed by a massive tsunami that reached 7 meters in height. Using NASA Earth Observation Systems (EOS), we will assess the damage to the area impacted by this disaster. Landsat 5 Thematic Mapper (TM) will be used to create normalized difference vegetation index (NDVI) maps and MODIS on Aqua and Terra will be used to produce enhanced vegetation index (EVI) maps. CALIPSO and HYSPLIT will also be used to track smoke released from the Fukushima Daiichi Nuclear Power Plant to assess where the radioactive steam traveled to. An analysis of the economic impact of this disaster will be conducted to help policy makers in the United States and in Japan. This research will highlight the ability of NASA EOS in providing support for disaster response around the world.

Malcom Jones, Jr.↗

A preliminary review and evaluation of historical unrest and eruptive activity at Misti volcano, Arequipa, Peru

Since the Spanish arrived in Arequipa in the mid-16th century, 11 potential small-volume eruptions or periods of unrest are reported at Misti. Yet, documentation of these events is scant, often based upon secondary or tertiary sources, and has never been assessed using modern volcanological methods. As a result, Misti?s historical record contains significant uncertainty and iterated errors. We assess these events with modern volcanological knowledge and when possible return to primary published and archival sources. Additionally, we scrutinized historical newspapers, photos, and contemporary travelers? accounts to better understand the volcano?s activity since the early-19th century. Misti?s earliest reported activity in 1542CE illustrates the problems associated with the volcano?s historical record. The 1542CE event first appears in the literature in the late-19th century. Despite some uncertainty of its validity in these reports, it remains listed as a potential eruption. Numerous contemporaneous accounts document a major earthquake in Arequipa on 22 January 1582CE, causing major damage and widespread fear of an eruption. Upon investigation of the volcano, such fear was quickly dispelled. We suggest that Misti?s 1542CE event did not occur, rather its origin is the fear induced by the 1582CE tectonic earthquake and an iterated typographic error in the date. Historical documents also sometimes record when Misti was explicitly not active, such as in a letter from a resident of Arequipa published in 1872CE. Events are also missing from Misti?s reported history, such as a period of enhanced degassing in 1971CE that is recorded in newspaper articles and eye-witness accounts. In summary, we consider three of Misti?s reported events to be false, four to be possible but unlikely, three are likely, and the 1985CE degassing crisis certainly occurred. We further uncovered preliminary evidence of up to four additional degassing crises at the volcano during the 19th and 20th centuries.

Christopher Harpel↗

Coastal Land Change Due to Tectonic Processes and Implications for Relative Sea-Level Rise in the Samoan Islands

Of the major coastal land change mechanisms responsible for relative sea-level change, tectonic subsidence is generally quoted as ranging from < mm/yr to 1 cm/yr. However, we documented coseismic and ongoing post-earthquake surface displacements from continuous GPS and tide gauge/altimetry data that indicated rapid subsidence on two of the major Samoan Islands of 12 - 20 cm during and following the 8.1 2009 Tonga-Samoa earthquake. Earlier results and our modeling of GRACE-derived gravimetric data provided a preliminary forecast of future relative sea-level rise through rapid land subsidence [Han et al., 2019]. Of course these numerical forecasts of time-dependent deformation are only as good as our input observations and our assumed rheological models. As part of our current NASA Earth Surface and Interior study, we are obtaining a wider range of data to constrain and test alternate models of ongoing postseismic deformation across American Samoa and Upolu, Samoa: (1) times series of altimetry plus tide gauge data processed to complement the cGPS data available to provide high-temporal resolution, point measurements of uplift/subsidence, (2) InSAR derived observations of surface deformation across the highly vegetated Samoan Islands, (3) evaluating and using NASA satellite lidar data (ICESat-I & ICESat-II, GEDI) for fusion with multi-source topographic data sets and for estimating topographic change on the decadal time scale. We are evaluating and using these new observations to better understand and separate out local, island-wide, and multi-island subsidence patterns and to evaluate the high impact of rising sea-level in a tectonically active region.

Jeanne Sauber-Rosenberg↗

Finite Element Modeling of Spent Nuclear Fuel in Dry Storage under Seismic Excitation

Spent nuclear fuel (SNF) is currently stored at independent spent fuel storage installations (ISFSIs) at nuclear power plant sites across the US. The US Department of Energy (DOE) Spent Fuel and Waste Science and Technology (SFWST) program is currently working to understand the response of SNF during a seismic event. This investigation is part of a long-term effort by DOE to quantify the external mechanical loads on SNF. The SFWST program has recently completed campaigns evaluating SNF behavior during normal conditions of transport and in 30 cm package drop scenarios; both campaigns included a combination of testing and modeling. Likewise, the investigation into seismic response of SNF will include full-size shake table tests of a vertical storage overpack containing instrumented surrogate fuel assemblies, as well as pre- and post-test finite element modeling. This work focuses on the development of finite element models simulating the planned test configuration. The seismic response is calculated at two scales: at the package level, and at the fuel level. The package level model evaluates the behavior of the storage overpack and its contents in response to the ground motions of a representative earthquake. The motion of the package is input to the fuel level model, which consists of a detailed finite element model of a 17x17 or 16x16 pressurized water reactor (PWR) fuel assembly or of a 10x10 boiling water reactor (BWR) fuel assembly. The stress and strain on individual components such as fuel rods and guide tubes are extracted from this model. This work details the SNF response to various input ground motions, considering representative earthquakes in the Central and Eastern US under various site conditions. Also discussed is the development of a “reduced scope” finite element model of the SNF, which consists of a single fuel rod while retaining the dynamic response of the entire fuel assembly. This reduced scope model addresses the significant computational expense of modeling the SNF under long-duration input ground motions and facilitates simulating additional parameter variations of the model beyond the as-tested configuration. The modeling results provide important insight that will help to inform future testing and modeling efforts by the SFWST program to close the knowledge gap related to external loads on SNF.

Kadooka, Kevin↗

A Comprehensive Analysis of Real-World Accelerometer Data Quality in a Global Smartphone-based Seismic Network

The proliferation of low-cost sensors in smartphones has facilitated numerous applications; however, large-scale deployments often encounter performance issues. Sensing heterogeneity, which refers to varying data quality due to factors such as device differences and user behaviors, presents a significant challenge. In this research, we perform an extensive analysis of 3-axis accelerometer data from the MyShake system, a global seismic network utilizing smartphones. We systematically evaluate the quality of approximately 22 million 3-axis acceleration waveforms from over 81 thousand smartphone devices worldwide, using metrics that represent sampling rate and noise level. We explore a broad range of factors influencing accelerometer data quality, including smartphone and accelerometer manufacturers, phone specifications (release year, RAM, battery), geolocation, and time. Our findings indicate that multiple factors affect data quality, with accelerometer model and smartphone specifications being the most critical. In addition, we examine the influence of data quality on earthquake parameter estimation and show that removing low-quality accelerometer data enhances the accuracy of earthquake magnitude estimation.

58 GEOSCIENCES↗