Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic sequences”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Geology and geochemistry of the Middle Proterozoic Eastern Ghat mobile belt and its comparison with the lower crust of the Southern Peninsular shield

Two prominent rock suites constitute the lithology of the Eastern Ghat mobile belt: (1) the khondalite suite - the metapelites, and (2) the charnockite suite. Later intrusives include ultramafic sequences, anorthosites and granitic gneisses. The chief structural element in the rocks of the Eastern Ghats is a planar fabric (gneissosity), defined by the alignment of platy minerals like flattened quartz, garnet, sillimanite, graphite, etc. The parallelism between the foliation and the lithological layering is related to isoclinal folding. The major structural trend (axial plane foliation trend) observed in the belt is NE-SW. Five major tectonic events have been delineated in the belt. A boundary fault along the western margin of the Eastern Ghats, bordering the low grade terrain has been substantiated by recent gravity and the deep seismic sounding studies. Field evidence shows that the pyroxene granulites (basic granulites) post-date the khondalite suite, but are older than the charnockites as well as the granitic gneisses. Polyphase metamorphism, probably correlatable with different periods of deformation is recorded. The field relations in the Eastern Ghats point to the intense deformation of the terrain, apparently both before, during and after metamorphism.

Rao, M. V. Subba↗
DOE OSTI2024

Making Phase-Picking Neural Networks More Consistent and Interpretable

Improving the interpretability of phase-picking neural networks remains an important task to facilitate their deployment to routine, real-time seismic monitoring. The popular phase-picking neural networks published in the literature lack interpretability because their output prediction scores do not necessarily correspond with the reliability of phase picks and can even be highly inconsistent depending on how we window the waveform data. Here, we show that systematically shifting the waveforms during training and using an antialiasing filter within the neural network architecture can substantially improve the consistency of the output prediction scores and can even make them scale with the signal-to-noise ratios of the waveforms. We demonstrate the improvements by applying these approaches to a commonly used phase-picking neural network architecture and using waveform data from the 2019 Ridgecrest earthquake sequence.

58 GEOSCIENCES↗
DOE OSTI2024

Seismic Monitoring near Ithaca, New York, Reveals Nonuniform Distribution of Microseismicity in an Intraplate Region

Abstract Cornell University intends to use a deep direct-use geothermal system to heat its Ithaca, New York, campus. In preparation for this project, the Cornell Seismic Network has been monitoring the background seismicity in this intraplate region since 2019. From January 2020 to June 2023, 95 events were detected within 20 km of the proposed geothermal well site, with local magnitudes ranging from −1.02 to 0.56. None of these events appear in regional or national catalogs. Events locate in a narrow geographic band, with one-fourth exhibiting multimodal hypocentral probability peaks both near the surface and at 1–4 km depth. We relocate events with a joint hypocenter and 1D velocity model inversion, in addition to a fully nonlinear method, and then compare observations with synthetic waveforms. Together, these approaches provide strong evidence for >95% of events locating at the surface or within the 3-km-thick sedimentary sequence. We explore how anthropogenic activity and regional topographic stress may contribute to frequent surficial events. This information is critical for characterizing the background microseismicity for comparison during future geothermal operations. Ithaca’s geology of Paleozoic sediments overlying Precambrian crystalline basement is typical of many continental interiors, so these results also provide insight into intraplate microseismicity patterns.

Geochemistry & Geophysics↗
NASA NTRS1973

Identification of Geostructures of continental crust, particularly as they relate to mineral-resource evaluation

The author has identified the following significant results. Analysis of lineated lakes in the Umiat, Alaska area and comparison with known geology, gravity, and magnetic data in the the area suggest concealed structures exist at depth, possibly at or near basement, which may represent targets for petroleum exploration. Compilation of reconnaissance geologic data on 1:250,000 scale enlargements of ERTS-1 images near Corwin reveal structural and stratigraphic anomalies that suggest the Cretaceous sequence is less thick than supposed and is repeated in a series of plates superimposed by flat thrust faults. The structural style differs from that in coeval strata to the northeast, across the northwest-trending linear zone separating differing tectonic styles in older strata noted earlier. The regional extension of a fault known locally in the McCarthy area has been recognized; this fault appears to form the boundary of a significant terrane of mid-Paleozoic metamorphic rocks. ERTS-1 images are being used operationally, at 1:1,000,000 scale in the compilation of regional geologic maps, and at 1:250,000 scale in field mapping in the Brooks Range, in the study of faults in seismically active southern Alaska, in field-checking interpretations previously made from ERTS-1 imagery, and orthophoto base maps for geologic maps.

Gryc, G.↗
DOE OSTI2022

Safety Analysis for the Versatile Test Reactor Conceptual Design

Safety analysis simulations have been carried out using the SAS4A/SASSYS-1 safety analysis code system to update the baseline set of results for the Versatile Test Reactor conceptual design. Results from a previous revision of this work provided input for the Conceptual Safety Design Report summarizing the overall safety basis for the VTR. This work will also support future analyses that will be required for Chapter 13, “Accident Analysis,” of a preliminary safety analysis report. A number of transients have been evaluated based on categories of postulated accident sequences suggested for consideration in NUREG-1537. Ultimately, the full spectrum of event initiators and subsequent accident sequences will be defined through a probabilistic risk assessment of the plant design following the implementation of the Licensing Modernization Project guidelines as described in the VTR Safety Design Strategy. Since that process is ongoing, the present analysis addresses several transients that are expected to be bounding accident scenarios that represent the three key ways to perturb a reactor: through changes to the core inlet temperature, mass flow rate, or reactivity. These correspond to a loss of heat sink (LOHS), loss of flow (LOF) or station blackout (SBO), and transient overpower (TOP), respectively. Both Protected and Unprotected versions of these transients have been evaluated, along with additional transients that include transient overpower with one stuck rod, pump coast-down failures, seismic events, and overcooling events. At the current stage of design, transient simulation results for the Versatile Test Reactor indicate that large safety margins exist for many event initiators. However, several “enabling” assumptions have been made, in terms of both design features and design limits, in order to perform the analyses; and these assumptions will need to be revised as the design matures. In addition, a larger spectrum of events needs to be evaluated as part of the on-going probabilistic risk assessment process.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Managing Subsurface Pressure Buildup and Interference in Commercial-Scale CO 2 Storage Project with Proximal Injection Wells

Large-scale decarbonization using carbon capture and storage (CCS) is likely to involve many commercial-scale CO 2 storage projects located in close proximity to each other. This close proximity raises concerns over pressure interference among the storage projects. Pressure interference between injection and storage efforts can reduce the practicable CO 2 storage resource and force wells to inject CO 2 at a lower rate to avoid the fracture pressure thresholds per United States Environmental Protection Agency (EPA) Class VI well regulations to preserve injection and confining zone integrity and potentially mitigate against inducing seismic activity. These analyses employ numerical full-physics reservoir modeling to evaluate how pressure buildup fronts and CO 2 plumes evolve under commercial-scale injection volumes of CO 2 in which multiple storage sites located in close proximity occur in tandem. The simulation models mimic injection at pseudo basin-scale and assume homogeneous saline formation(s) as storage targets with a pair of upper and lower homogeneous seal layer/s. These analyses specifically investigate the efficacy of two basin-wide reservoir pressure management strategies in addressing the technical challenges associated with pressure buildup and CO 2 plume commingling. The strategies explored include: 1) enlarging the area of injection well spacing (WS) and 2) storing CO 2 in a stacked sequence (SSS) of saline formations compared to a single formation. The storage and confining zones properties assumed were held common across the scenarios, unless specified otherwise. Analyses results show that after injecting 4 million tons per year for 30 years using 4 separate wells (each injecting 1 million metric tons per year), the radius of CO 2 plume extends to a mere 3 km or less from injection wells. Meanwhile, the radius of pressure buildup ranges on the order of tens to a few hundreds of kilometers, depending on the magnitude of pressure buildup threshold that one would use to define the front. CO 2 plume commingling from different injection wells appears to occur 50 years post-injection, especially under scenarios with narrowly spaced (i.e., < 5 km apart) injection well locations. Findings from sensitivity cases on the well spacing suggest that storage formations modeled would require different well spacing to avoid fracture pressure thresholds. For instance, modeled storage formations with high fracture gradients (i.e., 0.8 psi/ft) would need less than 5–km well spacing, whereas those with lower fracture gradients (i.e., 0.7 psi/ft) would need approximately 20–km well spacing, based on assumed modeling parameters. Under stacked injection, the pressure challenges (described above) still exist but are more alleviated due to distributing the same injection volume across more available reservoir volume. These analyses demonstrate that stacked-sequence storage can effectively address the challenges, while still providing the same target CO 2 storage volumes and allowing a large number of storage projects to be deployed in the same basin by better utilizing the available storage resource across different reservoir depths. Among cases modeled, the resulting pressure buildup front is most suppressed when each storage project distributes injection volumes over several wells, each of which injects a portion of the total CO 2 across the stacked sequence. This strategy results in the smallest CO 2 aerial footprint amongst scenarios evaluated but also shows the largest reduction in the pressure buildup at the top of perforation at the injection wells (upwards of approximately 42 percent compared to the commercial-scale single-formation storage), the result of which is crucial to maintain caprock integrity. The findings presented by this research draw attention to the importance of greater coordination among storage operators and regulatory stakeholders to foster the upscaling and deployment of CCS. These analyses provide insights into required decision-making when considering multi-project deployment in a shared basin. Because these analyses evaluate a very specific geologic situation, they bear further investigations across other geologic situations.

42 ENGINEERING↗
DOE OSTI2024

Graph-learning approach to combine multiresolution seismic velocity models

SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.

Geochemistry & Geophysics↗
DOE OSTI2023

Economic and operational investigation of CO 2 sequestration through enhanced oil recovery in unconventional reservoirs in Colorado, USA

The ongoing CCUS commercial projects are highly relied on the support of government incentives due to massive capital investment. Here, this study analyzes the economics of carbon capture utilization and sequestration (CCUS) projects, shows a state-wide CCUS deployment exercise, followed by simulation results of enhanced oil recovery (EOR) based CO 2 storage in unconventional reservoirs. The comprehensive economic analysis of capture, transportation, sequestration costs, enhanced 45Q tax credits, and EOR revenue implies the practicality of CO 2 -EOR to offset the high CCUS costs. With the economics understanding, we study the top CO 2 sources, existing CO 2 pipelines, and sequestration sinks in the state of Colorado, USA. This paper next presents results from EOR simulation in one section of the unconventional Denver-Julesburg (DJ) Basin Niobrara and Codell reservoirs. The simulation model is based on a geological static model, incorporated with hydraulic fracture stimulation, history matched to production, and calibrated to the microseismic and time-lapse surface seismic data. The CO 2 -EOR simulation results show that oil production can be increased and more CO 2 stored with: a longer primary production period; the presence of a shut-in period; higher injection rates; and multi-well injectors. The modeling results show that about 7–10 Mscf of CO 2 will be stored when recovering 1 stb of EOR oil. By adding the enhanced oil revenue and the carbon credits together, it is estimated that the most economic case can generate $\$13$ MM when oil price is assumed to be $\$80$/stb, and the EOR oil revenue is 3.4 times greater than that generated from 45Q incentives. It corresponds to the scenario that a five-year primary production is followed by CO 2 injection into four wells with the sequence of injection (4 MMscf/day for 6 months), shut-in (6 months) and production (12 months). The best practices in this study will provide valuable insights for similar CCUS projects in other unconventional fields. Furthermore, this study defines a term named “Carbon Neutrality Index (CNI)” by comparing the amount of CO 2 stored with that burned by EOR oil. The CNI value of 0 indicates the enhanced oil is carbon neutral; a negative CNI value implies there is a net reduction in carbon emission. The 4-year huff-n-puff (HnP) simulation leads to a positive CNI value, indicating that the EOR oil generated in this process is not carbon neutral yet.

03 NATURAL GAS↗
DOE OSTI2025

A Deep Learning-Aided Workflow for Decoding the Stress Regime of Southern Nevada

The Rock Valley fault zone in southern Nevada has a notable history of seismic activity and is the site of a future direct comparison experiment of explosion and earthquake sources. This study aims to gain insight into regional tectonic processes by leveraging recent advances in seismic monitoring capabilities to elucidate the local stress regime. A crucial step in this investigation is the accurate determination of P-wave first-motion polarities, which play a vital role in resolving earthquake focal mechanisms of small earthquakes. Here, we deploy a deep learning-based method for automatic determination of first-motion polarities to vastly expand the polarity dataset beyond what has been reviewed by human analysts. By the integrating P-wave polarities with new measurements of S/P amplitude ratios, we obtain robust focal mechanism estimates for 1306 earthquakes with a local magnitude of 1 and above occurring between 2010 and 2023 in southern Nevada. We then use the focal mechanism catalog to examine the regional stress orientation, confirming an overall trans-tensional stress regime with smaller scale complexities illuminated by individual earthquake sequences. These findings demonstrate how detailed analyses of small earthquakes can provide fundamental information for understanding earthquake processes in the region and inform future experiments at the Nevada National Security Site.

58 GEOSCIENCES↗
NASA NTRS1988

A compilation of information and data on the Manson impact structure

A problem for the impact hypothesis for the Cretaceous-Tertiary (K-T) mass extinction is the apparent absence of an identifiable impact site. The Manson Impact Structure is a candidate because it is the largest recognized in the U.S.; it is relatively close to the largest and most abundant shocked quartz grains found at the K-T boundary; and its age is indistinguishable from that of the K-T boundary based on paleontological evidence, fission track dates, and preliminary Ar-40/Ar-39 measurements. The region of northwest central Iowa containing the Manson Impact Structure is covered by Quaternary glacial deposits underlain by Phanerozoic sedimentary rocks (mostly flat-lying carbonates) and Proterozoic red clastic, metamorphic, volcanic, and plutonic rocks. In a circular area about 22 miles (35 km) in diameter around Manson, Iowa, this normal sequence is absent or disturbed and near the center of the disturbed area granitic basement rocks have been uplifted some 20,000 ft (6000m). Attention was drawn to Manson initially by the unusual quality of the groundwater there. Within the structure three roughly concentric zones of rock associations have been identified: (1) displaced strata; (2) completely disrupted strata, and igneous and metamorphic rocks. Manson was established as an impact structure based on its circular shape, its central uplift, and the presence of shocked quartz within the granitic central uplift. A gravity survey identified locations of low-density brecciated rocks and high-density uplifted crystalline rocks, but the outer boundary of the structure could not be established. Aeromagnetic and ground magnetic surveys showed locations and depths of shallowly buried crystalline rock and the locations of faults. A refraction seismic survey identified the crystalline central uplift, determined that the average elevation of bedrock is 70 ft (20 m) higher outside the structure than within, and was used to map the bedrock topography within the structure. A connection between the Manson impact and the K-T boundary may be established or refuted through study of the impact energy, the impact time, and composition of host rock, possible impactors, and impact melts.

Hartung, Jack B.↗
NASA NTRS1993

Radiocarbon ages of lacustrine deposits in volcanic sequences of the Lomas Coloradas area, Socorro Island, Mexico

Extensive eruptions of alkalic basalt from low-elevation fissures and vents on the southern flank of the dormant volcano, Cerro Evermann, accompanied the most recent phase of volcanic activity on Socorro Island, and created the Lomas Coloradas, a broad, gently sloping terrain comprising the southern part of the island. We obtained 14C ages of 4690 +/- 270 BP (5000-5700 cal BP) and 5040 +/- 460 BP (5300-6300 cal BP) from lacustrine deposits that occur within volcanic sequences of the lower Lomas Coloradas. Apparently, the sediments accumulated within a topographic depression between two scoria cones shortly after they formed. The lacrustine environment was destroyed when the cones were breached by headward erosion of adjacent stream drainages. This was followed by the eruption of a thin basaltic flow from fissures near the base of the northernmost cone. The flow moved downslope for a short distance and into the drainages that presently bound the study area on the east and west. The flow postdates development of the present drainage system and may be very recent. Our 14C data, along with historical accounts of volcanic activity over the last century, including submarine eruptions that occurred a few km west of Socorro in early 1993, underscore the high risk for explosive volcanism in this region and the need for a detailed volcanic hazards plan and seismic monitoring.

NASA Center ARC↗
DOE OSTI2023

Global Characteristics of Observable Foreshocks for Large Earthquakes

Abstract Foreshocks are the only currently widely identified precursory seismic behavior, yet their utility and even identifiability are problematic, in part because of extreme variation in behavior. Here, we establish some global trends that help identify the expected frequency of foreshocks as well the type of earthquake most prone to foreshocks. We establish these tendencies using the global earthquake catalog of the U.S. Geological Survey National Earthquake Information Center with a completeness level of magnitude 5 and mainshocks with Mw≥7.0. Foreshocks are identified using three clustering algorithms to address the challenge of distinguishing foreshocks from background activity. The methods give a range of 15%–43% of large mainshocks having at least one foreshock but a narrower range of 13%–26% having at least one foreshock with magnitude within two units of the mainshock magnitude. These observed global foreshock rates are similar to regional values for a completeness level of magnitude 3 using the same detection conditions. The foreshock sequences have distinctive characteristics with the global composite population b-values being lower for foreshocks than for aftershocks, an attribute that is also manifested in synthetic catalogs computed by epidemic-type aftershock sequences, which intrinsically involves only cascading processes. Focal mechanism similarity of foreshocks relative to mainshocks is more pronounced than for aftershocks. Despite these distinguishing characteristics of foreshock sequences, the conditions that promote high foreshock productivity are similar to those that promote high aftershock productivity. For instance, a modestly higher percentage of interplate mainshocks have foreshocks than intraplate mainshocks, and reverse faulting events slightly more commonly have foreshocks than normal or strike-slip-faulting mainshocks. The western circum-Pacific is prone to having slightly more foreshock activity than the eastern circum-Pacific.

Geochemistry & Geophysics↗
DOE OSTI2021

PLATO hare-and-hounds exercise: asteroseismic model fitting of main-sequence solar-like pulsators

ABSTRACT Asteroseismology is a powerful tool to infer fundamental stellar properties. The use of these asteroseismic-inferred properties in a growing number of astrophysical contexts makes it vital to understand their accuracy. Consequently, we performed a hare-and-hounds exercise where the hares simulated data for six artificial main-sequence stars and the hounds inferred their properties based on different inference procedures. To mimic a pipeline such as that planned for the PLATO mission, all hounds used the same model grid. Some stars were simulated using the physics adopted in the grid, others a different one. The maximum relative differences found (in absolute value) between the inferred and true values of the mass, radius, and age were 4.32, 1.33, and 11.25 per cent, respectively. The largest systematic differences in radius and age were found for a star simulated assuming gravitational settling, not accounted for in the model grid, with biases of −0.88 per cent (radius) and 8.66 per cent (age). For the mass, the most significant bias (−3.16 per cent) was found for a star with a helium enrichment ratio outside the grid range. Moreover, an ∼7 per cent dispersion in age was found when adopting different prescriptions for the surface corrections or shifting the classical observations by ±1σ. The choice of the relative weight given to the classical and seismic constraints also impacted significantly the accuracy and precision of the results. Interestingly, only a few frequencies were required to achieve accurate results on the mass and radius. For the age the same was true when at least one l = 2 mode was considered.

Cunha, M. S. (ORCID:0000000182377343)↗
DOE OSTI2021

Stratigraphically Controlled Stress Variations at the Hydraulic Fracture Test Site-1 in the Midland Basin, TX

We investigated the relationship between stratigraphy, stress, and microseismicity at the Hydraulic Fracture Test Site-1. The site comprises two sets of horizontal wells in the Wolfcamp shale and a deviated well drilled after hydraulic fracturing. Regional stresses indicate normal/strike-slip faulting with E-W compression. Stress measurements in vertical and horizontal wells show that the minimum principal stress varies with depth. Strata with high clay and organic content show high values of the least compressive stress, consistent with the theory of viscous stress relaxation. By integrating data from core, logs, and the hydraulic fracturing stages, we constructed a stress profile for the Wolfcamp sequence, which predicts how much pressure is required for hydraulic fracture growth. We applied the results to fracture orientation data from image logs to determine the population of pre-existing faults that are expected to slip during stimulation. We also determined microseismic focal plane mechanisms and found slip on steeply dipping planes striking NW, consistent with the orientations of potentially active faults predicted by the stress model. This case study represents a general approach for integrating stress measurements and rock properties to predict hydraulic fracture growth and the characteristics of injection-induced microseismicity.

microseismicity↗
NASA NTRS1992

Impact-generated winds on Venus: Causes and effects

The pressure of the dense atmosphere of Venus significantly changes the appearance of ejecta deposits relative to craters on the Moon and Mercury. Conversely, specific styles and sequences of ejecta emplacement can be inferred to represent different intensities of atmospheric response winds acting over different timescales. Three characteristic timescales can be inferred from the geologic record: surface scouring and impactor-controlled (angle and direction) initiation of the long fluidized run-out flows; nonballistic emplacement of inner, radar-bright ejecta facies and radar-dark outer facies; and very late reworking of surface materials. These three timescales roughly correspond to processes observed in laboratory experiments that can be scaled to conditions on Venus (with appropriate assumptions): coupling between the atmosphere and earlytime vapor/melt (target and impactor) that produces an intense shock that subsequently evolves into blast/response winds; less energetic dynamic response of the atmosphere to the outward-moving ballistic ejecta curtain that generates nonthermal turbulent eddies; and late recovery of the atmosphere to impact-generated thermal and pressure gradients expressed as low-energy but long-lived winds. These different timescales and processes can be viewed as the atmosphere equivalent of shock melting, material motion, and far-field seismic response in the target. The three processes (early Processes, Atmospheric Processes, and Late Recovery Winds) are discussed at length.

Schultz, Pater H.↗
NASA NTRS2002

The Late Pliocene Eltanin Impact: Documentation From Sediment Core Analyses

The expeditions ANT-XII/4 (1995) and ANT-XVIII/5a (2001) of the RV POLARSTERN collected extensive bathymetric and seismic data sets as well as sediment cores from an area in the Bellingshausen Sea (eastern Pacific Southern Ocean) that allow the first comprehensive geoscientific documentation of an asteroid impact into a deep ocean (approx. 5 km) basin, named the Eltanin impact. Impact deposits have now been recovered from a total of more than 20 sediment cores collected in an area covering about 80,000 km2. Combined biomagnetostratigraphic dating places the impact event into the earliest Matuyama Chron, a period of enhanced climate variability. Sediment texture analyses and studies of sediment composition including grain size and microfossil distribution reveal the pattern of impact- related sediment disturbance and the sedimentary processes immediately following the impact event. The pattern is complicated by the San Martin Seamounts (approx. 57.5 S, 91 W), a large topographic elevation that rises up to 3000 m above the surrounding abyssal plain in the area affected by the Eltanin impact. The impact ripped up sediments as old as Eocene and probably Paleocene that have been redeposited in a chaotic assemblage. This is followed by a sequence sedimented from a turbulent flow at the sea floor, overprinted by fall-out of airborne meteoritic ejecta that settled trough the water column. Grain size distribution reveals the timing and interaction of the different sedimentary processes. The gathered estimate of ejecta mass deposited over the studied area, composed of shock-melted asteroidal material and unmelted meteorites including fragments up to 2.5 cm in diameter, point to an Eltanin asteroid larger than the 1 km in diameter size originally suggested as a minimum based on the ANT-XII/4 results. This places the energy released by the impact at the threshold of those considered to cause environmental disturbance at a global scale and it makes the impact a likely transport mechanism explaining the presence of extinct Cenozoic microfossils in the transantarctic Sirius Unit. Although a crater structure representing Eltanin ground zero has not been discovered, the distribution pattern of sediment disturbance and ejecta deposits now allows to better determine the central target area north of the San Martin Seamounts.

Gersonde, R.↗
DOE OSTI2022

Earthquake Phase Association Using a Bayesian Gaussian Mixture Model

Earthquake phase association algorithms aggregate picked seismic phases from a network of seismometers into individual seismic events and play an important role in earthquake monitoring and research. Dense seismic networks and improved phase picking methods produce massive seismic phase datasets, particularly for earthquake swarms and aftershocks occurring closely in time and space, making phase association a challenging problem. Here, we present a new association method, the Gaussian Mixture Model Association (GaMMA), that combines the Gaussian mixture model with earthquake location, origin time, and magnitude estimation. We treat earthquake phase association as an unsupervised clustering problem in a probabilistic framework, where each earthquake corresponds to a cluster of P and S phases with a hyperbolic moveout of arrival times and a decay of amplitude with distance. We use the multivariate Gaussian distribution to model the collection of phase picks of an event; and the mean of the multivariate Gaussian distribution is given by the predicted arrival time and amplitude from the causative event. We carry out the pick assignment to each earthquake and determine earthquake source parameters (i.e., earthquake location, origin time, and magnitude) under the maximum likelihood criterion using the Expectation-Maximization algorithm. The GaMMA method does not require typical association steps of other algorithms, such as grid-search or supervised training. The results for both synthetic tests and for the 2019 Ridgecrest earthquake sequence show that GaMMA effectively associates phases from a temporally and spatially dense earthquake sequence while producing useful estimates of earthquake location and magnitude.

58 GEOSCIENCES↗
NASA NTRS1996

Transient Creep of a Composite Lower Crust: A Polymineralic Basis for Rapidly Evolving Postseismic Deformation Modes - 2

Postseismic horizontal strain and displacement following the June 28, 1992, Landers, California, earthquake (M(sub W) 7.3) is broad scale and cannot be explained solely by delayed afterslip located at the rupturing fault trace. Both the observed strain at Pifion Flat Observatory (PFO) and observed Global Positioning System receiver velocities evolve rapidly after the Landers-Big Bear earthquake sequence. The observed exponential decay of these motions, with timescales of 4-34 days, may reflect a soft creep rheology in the lower crust and brittle-ductile transition zone or even within the seismogenic crust itself. Here a simple model of a two-dimensional screw dislocation in a layered Maxwell viscoelastic Earth is employed in conjunction with a composite rheology to demonstrate that the short timescale transient response modes (approx. = 4-34 days) are consistent with the behavior of a biviscous lower crust. The lowest viscosity of this system is derivable from laboratory experimental data on the long-term creep of natural quarztites, and the highest viscosity is consistent with isostasy-related lower crustal flow in a continental extensional tectonic environment. The model predicts significant stress relaxation at the base of the seismogenic crust. Near the base of the seismogenic zone, and about 4 km away from the mainshock, the rate of predicted relaxation is of the order of 0.01 MPa/ d during the first 20 days of postseismic flow. Oblate spheroidal inclusions at 5% concentration levels that are both aligned and fairly flat in shape and that have a viscosity of 3-4 x 10(exp 15) Pa s are consistent with both the amplitude and decay time of horizontal crustal strain observed at PFO after the Landers mainshock. It is speculated that the structures exposed in cross sections and in seismic reflection profiles of the lower crust that have mylonitic associations are, in part, the cause of such rapid postseismic evolution in southeastern California. Unmylonitized quartz-rich rock at sufficiently elevated temperatures could also contribute to the rapid decay modes.

Ivins, Erik R.↗