SEARCH · Search NASA
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.
For treaty monitoring, the ability to accurately locate seismic events is essential. Since many event location techniques rely on the comparison of observed to predicted seismic wave travel times through the Earth, it is important to be able to accurately estimate travel times and their uncertainties.
The Earth’s Coriolis force has been well–known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear–wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear–wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well–isolated phases. Here, we recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.
The “light mantle” deposit at the base of South Massif in the Moon's Taurus‐Littrow Valley was a primary science target for the Apollo 17 exploration. The possibility that it was a landslide triggered by ejecta from Tycho Crater is critical for establishing the age of Tycho and constraining recent lunar impact chronology; however, the mechanism of emplacement of the deposit has recently been questioned. The newly opened 73001/73002 double drive tube from Station 3 sampled 70.6 cm deep into the regolith and represents the first stratigraphic section of an extraterrestrial landslide deposit returned to Earth. Here we provide an overview of the stratigraphy of the 73001/73002 core based on top to bottom variations revealed by coordinated laboratory analyses and explore constraints on the emplacement of the light mantle deposit. Briefly, the upper ∼10 cm of 73002 contains a disturbed zone from space weathering and emplacement of ejecta from a nearby crater that excavated and ejected basaltic material. Below 10 cm is a nearly uniform unit of immature regolith. These data support a single event for the emplacement of the deposit at this location, followed by weathering and mixing of materials from nearby crater ejecta in the upper 10 cm. Slight variations in chemistry and clast components may reflect the relative stratigraphy of the South Massif slope, with material toward the bottom of 73001 originating from lower slopes and material from higher up in the core representing regolith from higher up the South Massif slopes.
The San Juan Basin (SJB), located in southwestern Colorado northwestern New Mexico, containing Paleozoic through middle Cenozoic strata, formed as a partitioned basin approximately 80 Ma in response to the Laramide orogeny. The SJB is a commercially mature, petroleum-producing basin, and currently being explored for CCUS and geothermal resources. Modern heat flow within the SJB is spatially variable with higher geothermal gradients in its northern and eastern portions. The temporal history of the basin thermal history is essential for understanding capacity for carbon storage and geothermal exploration. We present thermochronometric analyses of two surface and two subsurface samples to constrain spatial-temporal thermal evolution of the northwestern San Juan Basin. Apatite (U-Th)/He thermochronometric analyses were conducted on Cliff House and Kirtland Formation outcrops and subsurface samples of the Ojo Alamo and Pictured Cliffs Formations from a DOE funded, CCUS project pilot-well. In addition, vitrinite reflectance data and 1D basin modelling constrain possible time-temperature pathways of the samples in an area where no previous studies have constrained the uplift/exhumation of the Hogback Monocline using thermochronometry. Modeling suggests that elevated subsurface temperatures developed simultaneously with regional late Oligiocene volcanism in the San Juan Volcanic field and timing of maximum burial in the SJB. Late Miocene to Pliocene cooling/exhumation through the apatite (U-Th)/He partial retention zone support stratigraphic evidence of the Colorado Plateau uplift, possibly regional epeirogenic uplift driven by mantle processes.
Global seismographic networks (GSNs) emerged during the late nineteenth and early twentieth centuries, facilitated by seminal international developments in theory, technology, instrumentation, and data exchange. The mid- to late-twentieth century saw the creation of the World-Wide Standardized Seismographic Network (1961) and International Deployment of Accelerometers (1976), which advanced global geographic coverage as seismometer bandwidth increased greatly allowing for the recording of the Earth's principal seismic spectrum. The modern era of global observations and rapid data access began during the 1980s, and notably included the inception of the GEOSCOPE initiative (1982) and GSN (1988). Through continual improvements, GEOSCOPE and the GSN have realized near-real time recording of ground motion with state-of-art data quality, dynamic range, and timing precision to encompass 180 seismic stations, many in very remote locations. Data from GSNs are increasingly integrated with other geophysical data (e.g., space geodesy, infrasound and Interferometric Synthetic Aperture Radar). Globally distributed seismic data are critical to resolving crust, mantle, and core structure; illuminating features of the plate tectonic and mantle convection system; rapid characterization of earthquakes; identification of potential tsunamis; global nuclear test verification; and provide sensitive proxies for environmental changes. As the global geosciences community continues to advance our understanding of Earth structure and processes controlling elastic wave propagation, GSN infrastructure offers a springboard to realize increasingly multi-instrument geophysical observatories. Here, we review the historical, scientific, and monitoring heritage of GSNs, summarize key discoveries, and discuss future associated opportunities for Earth Science.
Fluids generated from subducted slabs participate in the cycling of deep carbon in the crust and upper mantle. In these fluids, aqueous carbon species vary in oxidation state between +IV and -IV depending on whether the fluids are oxidizing or reducing, respectively. Most studies of subduction-zone fluids have focused on oxidized carbon species. However, recent studies of natural samples have demonstrated the occurrence of deep, reducing fluids, generated in both subducted oceanic upper mantle and crustal rocks. CH 4 -H 2 -rich fluid inclusions in subducted carbonated serpentinites have demonstrated the existence of abiotic, immiscible, hydrocarbon fluids at upper mantle conditions. To investigate the formation of such immiscible hydrocarbon fluids during the evolution of subducted carbonated serpentinites, we used equilibrium constants from the Deep Earth Water model to carry out predictive chemical mass transfer modeling to simulate the alteration reactions. A novel feature of the models was the inclusion of an immiscible hydrocarbon fluid containing six components (CH 4,f , C 2 H 6,f , C 3 H 8,f , isoC 4 H 10,f , CO 2,f , H 2,f ). This feature enabled prediction of the formation of a separate immiscible fluid in equilibrium with aqueous species and minerals. We developed a predictive reaction path model of invasive H 2,f reacting with carbonated serpentinites and interstitial aqueous fluids for comparison with the natural samples from the Lanzo Massif, western Italian Alps. Over a range of temperatures and pressures, immiscible hydrocarbon fluids formed in association with altered mineral assemblages. Reaction progress caused the transformation of carbonated serpentinites and the formation of clinopyroxene, brucite, graphite, and hydrocarbon fluids, along with changes of pH, logfO 2 , and aqueous species. CH 4,f was the most abundant hydrocarbon species in all the models. The overall results at 2.0 GPa and 400 to 450 °C were consistent with the natural samples from the Lanzo Massif. Interestingly, large amounts of H 2 O formed due to oxidation of H 2 . More hydrocarbons and H 2 O formed in models with lower fluid/rock mass ratios or with more reactant H 2 . Models at different pressure and temperature conditions showed similar results with some variation in the relative stabilities of aragonite, graphite and olivine solid solution, and associated differences in mineral sequences, hydrocarbon fluids, values of aqueous species, and the final log fO 2 and pH. As a result, our models strongly support the laboratory and field evidence that reduction of carbonated serpentinites by infiltrating H 2 fluids can cause the formation of immiscible, abiotic hydrocarbon fluids in subduction zones.
SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.
The chemical compositions of island arc basalts (IAB) reflect contributions from the mantle as well as fluids and melts from the subducting slab. Addition of radiogenic seawater Sr to oceanic crust through hydrothermal alteration and subsequent subduction is often invoked to explain elevated 87 Sr/ 86 Sr signatures in modern IAB. However, changes in the 87 Sr/ 86 Sr of island arc magmatic rocks through time has not been investigated, limiting our understanding of the factors influencing the Sr budgets of arcs throughout Earth’s history. To address this, we compiled 87 Sr/ 86 Sr values from island arc magmatic rocks ranging in age from modern to Paleoproterozoic, only including data from island arc localities that best preserve initial magmatic 87 Sr/ 86 Sr. Median initial 87 Sr/ 86 Sr values are consistently elevated compared to depleted mantle 87 Sr/ 86 Sr over this period, indicating persistent enrichment in radiogenic Sr in island arcs. Moreover, the elevation in island arc 87 Sr/ 86 Sr relative to the depleted mantle is variable. A notable rise in island arc 87 Sr/ 86 Sr during the late Neoproterozoic coincides with a steep increase in seawater 87 Sr/ 86 Sr and Sr concentration. To investigate this potential connectivity, we modeled the 87 Sr/ 86 Sr of island arc magmas between 0 and 830 Ma with inputs of depleted mantle 87 Sr/ 86 Sr, seawater 87 Sr/ 86 Sr, and seawater Sr concentration. The model reproduces the overall trajectory of the compiled data. We interpret the observed temporal variation in island arc 87 Sr/ 86 Sr values and its close association with fluctuations in seawater chemistry as evidence that changes in marine geochemistry have strongly influenced the Sr isotopic record of island arc magmas over time.
Iron(Fe)-water reactions in a magma ocean can influence water storage and density of planets. These reactions can form Fe-O-H phases, whose density, melting, and electronic properties at planetary interior conditions are important for informing planetary models. Here, we study natural goethite (α-FeOOH) that is shock-compressed along its principal Hugoniot. Analysis of our velocity interferometer system for any reflector (VISAR) results extends the equation of state to over 800 GPa. X-ray diffraction and VISAR reflectivity results indicate the onset of melting occurs at ~95 GPa with complete melting by 166 GPa, which may be relevant to low seismic velocity anomalies observed above the core-mantle boundary. Analysis of X-ray emission spectroscopy results up to 285 GPa shows the spin crossover of Fe, with dominantly low spin Fe above ~265 GPa in the melt, supporting formation of dense basal magma oceans in terrestrial planets. Using our measured FeOOH densities, we model planetary interiors up to 10 Earth masses. Assuming FeOOH forms via iron-water reactions, the radius decreases by up to 28%, while the density increases by up to 165% compared to the unreacted cases, providing an avenue to investigate water storage and evolution in super-Earths and sub-Neptunes.
We report the detection of an atmosphere on a rocky exoplanet, GJ 1132 b, which is similar to Earth in terms of size and density. The atmospheric transmission spectrum was detected using Hubble WFC3 measurements and shows spectral signatures of aerosol scattering, HCN, and CH{sub 4} in a low mean molecular weight atmosphere. We model the atmospheric loss process and conclude that GJ 1132 b likely lost the original H/He envelope, suggesting that the atmosphere that we detect has been reestablished. We explore the possibility of H{sub 2} mantle degassing, previously identified as a possibility for this planet by theoretical studies, and find that outgassing from ultra-reduced magma could produce the observed atmosphere. In this way we use the observed exoplanet transmission spectrum to gain insights into magma composition for a terrestrial planet. The detection of an atmosphere on this rocky planet raises the possibility that the numerous powerfully irradiated super-Earth planets, believed to be the evaporated cores of sub-Neptunes, may, under favorable circumstances, host detectable atmospheres.
Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.
Subduction of oceanic lithosphere transports surface H 2 O into the mantle. Recent studies show that dense SiO 2 in the form of stishovite, an abundant mineral in subducted oceanic crust at depths greater than ~270 km, has the potential to host and transport a considerable amount of H 2 O into the lower mantle, but the H 2 O storage capacity of SiO 2 phases at high pressure and temperature remains uncertain. We investigate the hydration of stishovite and its higher-pressure polymorphs, β-stishovite and seifertite, with in situ X-ray diffraction experiments at high pressures and temperatures. The H 2 O contents in SiO 2 phases are quantified based on observed increases in unit cell volume relative to the anhydrous SiO 2 system. Density functional theory (DFT) computations permit calibration of water content as a function of volume change based on interstitial substitution of H 2 O. Regression of our experimental data indicates an H 2 O storage capacity in stishovite of ~3.5 wt% in the transition zone and shallow lower mantle, decreasing to about 0.8 wt% at the base of the mantle. We find that SiO 2 -bearing subducted oceanic crust can accommodate all the H 2 O in slab lithosphere that survives sub-arc dehydration. Hydration of silica phases in subducted oceanic crust and their unparalleled capacity to host significant amounts of H 2 O even at high mantle temperatures provides a unique mechanism for transport and storage of water in the deepest mantle.
Abstract A function of global monitoring of nuclear explosions is the development of Earth models for predicting seismic travel times for more accurate calculation of event locations. Most monitoring agencies rely on fast, distance-dependent one-dimensional (1D) Earth models to calculate seismic event locations quickly and in near real-time. RSTT (Regional Seismic Travel Time) is a seismic velocity model and computer software package that captures the major effects of three-dimensional crust and upper mantle structure on regional seismic travel times, while still allowing for fast prediction speed (milliseconds). We describe updates to the RSTT model using a refined data set of regional phases (i.e., Pn, Pg, Sn, Lg) using the Bayesloc relative relocation algorithm. The tomographic inversion shown here acts to refine the previous RSTT public model ( rstt201404um ) and displays significant features related to areas of global tectonic complexity as well as further reduction in arrival residual values. Validation of the updated RSTT model demonstrates significant reduction in median epicenter mislocation (15.3 km) using all regional phases compared to the iasp91 1D model (22.1 km) as well as to the current station correction approach used at the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (18.9 km).
We report Northwest Africa (NWA) 7042 is an intermediate, permafic shergottite consisting of two generations of zoned olivine (early olivine Fo 41-76 , and late-stage olivine Fo 46-56 ), complexly zoned pyroxene (En 35-64 Fs 22-46 Wo 5-34 ), maskelynitized feldspar and shock melt (An 5-30 Ab 16-61 Or 1-47 ), and accessory merrillite, apatite, ilmenite, titanomagnetite, Fe-Cr-Ti spinels, pyrrhotite, and baddeleyite. The olivine grains have been pervasively altered resulting in orange-brown stained cores with dark-brown margins, surrounded by significantly less altered overgrowth rims. This textural relationship suggests that the alteration of the cores took place at magmatic temperatures prior to crystallization of the rims on Mars. While this type of olivine alteration is rare in shergottites, it is similar to deuterically altered olivine in basalts and gabbros on Earth, caused by residual water-rich magmatic fluids. Launch-generated shock veins in NWA 7042 crosscut and displace several of the altered olivine grains indicating that the alteration occurred before ejection of the meteorite from Mars. However, NWA 7042 has also been subjected to extensive terrestrial weathering, with terrestrial calcite/barite veins occurring throughout the sample. Transmission electron microscopy analysis of the olivine alteration also reveals an absence of high-temperature phases typically associated with hydrous deuteric alteration on Earth, complicating our interpretation. The potential presence of deuterically altered olivine in NWA 7042 has significant implications, as it is the third martian meteorite where deuteric alteration of olivine has been observed (the others being NWA 10416, ALH 77005). The different mantle sources for the parental melts of these three meteorites would suggest many, if not all martian mantle reservoirs have the potential to produce water-rich magmas.
Here, we report for the first time the formation of a HP-PdF 2 -type FeCl 2 phase (space group Pa$\bar{3}$), through high pressure-temperature (P-T) reactions in the hydrous systems (Mg 0.6 Fe 0.4 )SiO 3 –H 2 O–NaCl and FeO2H–NaCl in a laser-heated diamond-anvil cell up to 108 GPa and 2000 K. Applying single-crystal X-ray diffraction (XRD) analysis to individual submicrometer-sized grains, we have successfully determined the crystal structure of the as-synthesized FeCl 2 phase, in agreement with our theoretical structure search results. In situ high P-T XRD data revealed the substitution of Cl for OH(O) in such a cubic Pa$\bar{3}$ structure, demonstrating that this topology is a potential host for both H and Cl in the deep Earth. The chemical analysis of the recovered sample showed that the post-perovskite phase contains considerable amounts of Na 2 O and Fe 2 O 3 . The coexistence of the cubic FeCl 2 phase and post-perovskite suggests that the lowermost mantle could be a potential reservoir of Cl. The possible presence of volatiles such as H and Cl in the deep lower mantle would impact the composition and iron valence state of the post-perovskite phase.
The D″ region of the lower mantle, which lies just above the core–mantle boundary, is distinct from the bulk of the lower mantle in that it exhibits complex seismic heterogeneity and seismic anisotropy. Seismic anisotropy in this region is likely to be largely due to the deformation-induced texture (crystallographic preferred orientation) development of the constituent mineral phases. Thus, seismic anisotropy can provide a marker for deformation processes occurring in this dynamic region of the Earth. Post-perovskite-structured (Mg,Fe)SiO3 is believed to be the dominant mineral phase in many regions of the D”. As such, understanding deformation mechanisms and texture development in post-perovskite is important for the interpretation of observed seismic anisotropy. Here, we report on high-pressure diamond anvil cell deformation experiments on NaMgF3 neighborite (perovskite structure) and post-perovskite. During deformation, neighborite develops a 100 texture, as has been previously observed, both in NaMgF3 and MgSiO3 perovskite. Upon transformation to the post-perovskite phase, an initial texture of {130} at high angles to compression is observed, indicating that the {100} planes of perovskite become the ~{130} planes of post-perovskite. Further compression results in the development of a shoulder towards (001) in the inverse pole figure. Plasticity modeling using the elasto-viscoplastic self-consistent code shows this texture evolution to be most consistent with deformation on (001)[100] with some contribution of glide on (100)[010] and (001)<110> in NaMgF3 post-perovskite. The transformation and deformation mechanisms observed in this study in the NaMgF3 system are consistent with the behavior generally observed in other perovskite–post-perovskite systems, including the MgSiO3 system. This shows that NaMgF3 is a good analog for the mantle bridgmanite and MgSiO3 post-perovskite.
Rock strength has long been linked to lithospheric deformation and seismicity. However, independent constraints on the related elastic heterogeneity are missing, yet could provide key information for solid Earth dynamics. Using coseismic Global Navigation Satellite Systems (GNSS) data for the 2011 M9 Tohoku-oki earthquake in Japan, we apply an inverse method to infer elastic structure and fault slip simultaneously. We find compliant material beneath the volcanic arc and in the mantle wedge within the partial melt generation zone inferred to lie above ~100 km slab depth. We also identify low-rigidity material closer to the trench matching seismicity patterns, likely associated with accretionary wedge structure. Along with traditional seismic and electromagnetic methods, our approach opens up avenues for multiphysics inversions. Those have the potential to advance earthquake and volcano science, and in particular once expanded to InSAR type constraints, may lead to a better understanding of transient lithospheric deformation across scales.
The Mid-Pliocene represents the most recent interval in Earth history with climatic conditions similar to those expected in the coming decades. Mid-Pliocene sea level estimates therefore provide important constraints on projections of future ice sheet behavior and sea level change but differ by tens of meters due to local distortion of paleoshorelines caused by mantle dynamics. We combine an Australian sea level marker compilation with geodynamic simulations and probabilistic inversions to quantify and remove these post-Pliocene vertical motions at continental scale. Dynamic topography accounts for most of the observed sea level marker deflection, and correcting for this effect and glacial isostatic adjustment yields a Mid-Pliocene global mean sea level of +16.0 (+10.4 to +21.5) m (50th/16th to 84th percentiles). Recalibration of recent high-end sea level projections using this revised estimate implies a more stable Antarctic Ice Sheet under future warming scenarios, consistent with midrange forecasts of sea level rise that do not incorporate a marine ice cliff instability.