SEARCH · Engineering Papers
Results for “crustal”
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.
Storage and Diffusion of Carbon and Nitrogen in Crustal Materials
Nitrogen is a primary constituent of Earth’s atmosphere and an essential component of many biological processes. Recent developments in understanding N incorporation into silicate minerals suggest that the solid earth and lithosphere could provide a substantial reservoir for the planetary N budget. Furthermore, determining the concentration and distribution of low abundances of N in silicates presents a significant challenge in microanalysis. We tested capabilities for high spatial resolution quantification of low levels of N using EPMA with the goal of developing standards to be used in coupled N and noble gas analyses of silicate minerals. We investigated the concentration and distribution of N in a suite of ammonium-bearing silicate minerals that reflect the dominant N-bearing phases of the lithosphere: tobelite, buddingtonite, ammonioleucite, tsaregorodtsevite, and hyalophane. These minerals form during NH 4 -rich alteration of K-bearing silicate minerals, and the bonding environment of N in each mineral is likely to differ based on crystal structure. Using an optimized analytical routine for N quantification enables sensitivity approaching 0.02 weight% N and high-resolution quantitative mapping of sub-micron variations in N contents.
Dynamic Rupture Modeling of the 1999 Chi-Chi, Taiwan Earthquake Using 3DFinite Element Method. Sensitivity Analysis of Slip Rate Function to Model Parameters of Crustal Weak-Zone
The main objective of this study was the dynamic rupture modeling of the M7.3 1999 Chi-Chi, Taiwan earthquake. The purpose of the numerical modeling was two folded. First, using models of rupture dynamics, derive kinematic rupture characteristics, including spatial and temporal variations of the slip rate time history and rupture speed, in the long-period ground motion generation area (LMGA), strong motion generation areas (SMGAs) and background fault area for a thrust shallow rupture. Second, investigate the sensitivity of source time function characteristics in the LMGA and SMGA areas to model parameterization of low-velocity weak zone, such as thickness and shear wave velocity.
Assessment of Crustal Leg Uncertainties for RSTT Path-dependent Uncertainty
The Regional Seismic Travel Time (RSTT) model and software were recently updated to allow for path-dependent uncertainty (PDU) estimates for travel-time (Begnaud et al. 2021; Begnaud et al. 2020). Begnaud et al. (2021) describes the process of improving the current distance-dependent uncertainty parameterization for RSTT using a random effects model (REM) to estimate slowness (inverse velocity) uncertainty as a mean squared error for each model parameter. The REM separates the error between observed slowness and model predicted slowness into bias and random components. The PDU is calculated by integrating these mean squared errors along a seismic-phase ray path.
Crustal Scale Travel Time Prediction with the SALSA3D Framework and Machine Learning
The SALSA3D project aims to improve our models used in travel time prediction. The current version uses tomographic modeling for propagation through the Earth’s mantle because of the large number (order of millions) of observations of seismic phases which primarily traverse the Earth’s mantle and the ability to pose the travel time problem as a set of linear equations. However, all seismic rays traverse the crust to reach receivers at the surface and therefore models of propagation through the crust are required. Therefore, the primary motivation for this study is to explore how to increase the scope of the SALSA3D project to phases which travel primarily through the crust. In this report, we evaluate new, machine learning based and physics-based methods to model these travel times for integration into the SALSA3D framework. Our results suggest that using our existing physics-based travel time tomography method is a viable approach for the regional to global scale, but better predictive capabilities can be achieved through a neural network trained on the region of interest for near-regional offsets. We suggest future iterations of SALSA3D should incorporate machine learning tools such as Physics-Informed Neural Networks or Bayesian Neural Networks.
Viscosity of the Mercurian Magma Ocean: Implications of Sulfur-Free and Sulfur-Bearing Magma Oceans for Differentiation and Crustal Petrogenesis
Mercury is believed to have had a magma ocean in its early evolution whose solidification dynamics contributed to the planet’s present day surface composition. The surface is characterized by two main geological provinces; the Northern Smooth Plains (NSP),and the Heavily Cratered Terrain (HCT) and Intercratered Terrian (IcT). These provinces are defined by the differences in morphology, composition, and age. The NSP are located in the northern hemisphere and have low Mg/Si, Ca/Si, and S/Si ratios and high Al/Si ratios, and are compositionally analogous to a magnesian alkali basalt. The HCT/IcT cover the majority of the planet and are older than the NSP, and are characterized by high Mg/Si, Ca/Si, and S/Si ratios and low Al/Si ratios. The dichotomy in surface composition motivates the study into the structure of the mantle as a product of magma ocean solidification. Two important and interrelated parameters that determine the efficiency of crystal fractionation during magma ocean solidification are convection vigor and viscosity of the magma ocean. Crystal suspension in a cooling magma ocean could be efficient or inefficient which suggests two end member possibilities for magma ocean solidification. The first is the formation of a compositionally stratified mantle reflecting a cumulate sequence produced along a fractional liquid line of descent. The second is the formation of a well-mixed, batch crystallized mantle by high viscosity and/or vigorously convecting magma ocean. Determining which of these end members, or where on the spectrum between them is most relevant to Mercury’s magma ocean can be experimentally constrained through high P-T falling-sphere viscometry experiments.
Dynamic Rupture Simulations of the Mw7.2 1992 Landers,California, Earthquake
An important part of the simulation technique implemented in the Recipe is the earthquake rupture characterization based on established empirical scaling relationships between source parameters and seismic moment. The scaling relationships used in the rupture characterization are derived from recorded data and are tested against dynamic and kinematic rupture models of crustal earthquakes. Currently, a revised three-stage source scaling relationship between the source rupture area and seismic moment of inland crustal earthquakes, initially proposed by Irikura and Miyake, is being used. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, including the PI, Arben Pitarka, indicate that in order to reproduce the near-fault long-period ground motions and the permanent displacements, observed during surface-rupturing earthquakes, the recipe needs modifications to the kinematic rupture parametrization in the shallow crustal weak zone. In particular, improvements are required in characterizing the slip rate function at shallow depths, above the sesimogenic zone. The main objective of this research project was the development of a dynamic rupture model for the Mw7.2 1992 Landers, California earthquake. The dynamic rupture model was used to derive kinematic rupture parameters, such as slip, slip-rate, rupture time distributions, and slip rate functions in the weak (0-4km) and seismogenic crustal zones. In addition to the knowledge gained from the dynamic rupture modeling of Chi-chi, Ridgecrest, Kumamoto and other crustal earthquakes, studied by the PI, the rupture characterization of the Landers earthquake can be helpful in updating the conventional recipe for broadband ground motion simulations of crustal earthquakes with shallow rupture. The project emphasizes the need for performing physics-based earthquake rupture modeling which ultimately will inform further improvements of rupture characterization. In particular, the modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, used in adjustments of empirical source scaling relationships for large crustal earthquakes.
Magmatic, hydrothermal and ore element transfer processes of the southeastern Archean Superior Province implied from electrical resistivity structure
Although magmatic- and metamorphic- derived fluids are widely recognized mineralizing agents, the role of crustal architecture in defining source and sink zones within the middle to lower crust and upper mantle of ancient orogens remains enigmatic. The globally largest and best-preserved Archean greenstone belts lie in the Superior Province, Canada. They provide an ideal location to investigate the influence of igneous construction and subsequent syn-deformational plutonism and metamorphism on the localization of metal-rich melts and fluids throughout the crustal column. Integration of three-dimensional magnetotelluric modelling and seismic reflection sections across the Abitibi subprovince reveals details of a 'whole-of-crust' magmatic and hydrothermal system. East-west low resistivity structures broadly underlie the surface traces of the major deformation zones that are host to significant gold endowment (>200 Moz), while mid-crustal cross trends suggest mineralized fluids flowed along 'pipes' within fault planes. Most low resistivity structures are inferred to represent domains containing interconnected zones of graphite and/or sulfide. Furthermore, these delineate relict mantle source/transit domains and crustal pathways enriched by the flow of magmas or metamorphic fluids genetically related to a late-stage pulse of ore-bearing magmatism, possibly as a result of slab break-off or delamination. Using the combined MT and seismic data, we develop a 3-D crustal-scale model which highlights how evolving orogenic architecture-controlled mass transfer and metallogenic processes developed.
First Synoptic Images of FUV Discrete Aurora and Discovery of Sinuous Aurora at Mars by EMM EMUS
Abstract We present the first measurements of Mars discrete aurora in the extreme ultraviolet (<110 nm) and the first synoptic aurora images in the far ultraviolet (110–180 nm). Auroral emission is detected in >75% of nightside images, with patterns shifting visibly over 15–20 min. Aurora is observed most frequently in regions of open magnetic topology (where crustal magnetic fields are very weak and/or vertical), with the brightest aurora where crustal fields are strongest. We present the first disk‐averaged spectrum of discrete aurora, with several O, C, and CO features as expected for electron impact primarily on CO 2 . We categorize discrete auroral morphology into three types: crustal field aurora, non‐crustal field patchy aurora, and a new type we call “sinuous” aurora, an elongated serpentine structure that stretches thousands of kilometers into the nightside from near midnight in the northern hemisphere. These observations point to a highly dynamic environment in Mars' magnetotail.
Whole-lithosphere shear during oblique rifting
Processes controlling the formation of continental whole-lithosphere shear zones are debated, but their existence requires that the lithosphere is mechanically coupled from base to top. We document the formation of a dextral, whole-lithosphere shear zone in the Death Valley region (DVR), southwest United States. Dextral deflections of depth gradients in the lithosphere-asthenosphere boundary and Moho are stacked vertically, defining a 20–50-km-wide, lower lithospheric shear zone with ~60 km of shear. These deflections underlie an upper-crustal fault zone that accrued ~60 km of dextral slip since ca. 8–7 Ma, when we infer that whole-lithosphere shear began. This dextral offset is less than net dextral offset on the upper-crustal fault zone (~90 km, ca. 13–0 Ma) and total upper-crustal extension (~250 km, ca. 16–0 Ma). We show that, before ca. 8–7 Ma, weak middle crust decoupled upper-crustal deformation from deformation in the lower crust and mantle lithosphere. Between 16 and 7 Ma, detachment slip thinned, uplifted, cooled, and thus strengthened the middle crust, which is exposed in metamorphic core complexes collocated with the whole-lithosphere shear zone. Midcrustal strengthening coupled the layered lithosphere vertically and therefore enabled whole-lithosphere dextral shear. Where thick crust exists (as in pre–16 Ma DVR), midcrustal strengthening is probably a necessary condition for whole-lithosphere shear.
Superthermal Electron Deposition on the Mars Nightside During ICMEs
Superthermal electron precipitation is one of the main sources supporting the Mars nightside ionosphere. It is expected that solar wind electron fluxes are to increase significantly during interplanetary coronal mass ejections (ICMEs) and therefore an enhanced nightside ionospheric density. This study is to quantify the variation of the precipitating and deposited electron fluxes during five extreme ICMEs encountered by Mars Global Surveyor (MGS). We find that energy fluxes correlate better with the upstream dynamic pressure proxy than number fluxes and electron fluxes increase more at high energies, which means that electrons tend to have a lower peak production altitude during storm times. The precipitating and net/deposited fluxes are increased up to an order of magnitude from low to high dynamic pressures. The estimated total electron content (TEC) is a few times of 10 14 m -2 for quiet times and on the order of 1015 m -2 for storm times, with an enhancement up to an order of magnitude locally near strong crustal fields. Crustal magnetic fields have an effect on the deposited fluxes with more prominent magnetic reflection over strong magnetic fields during quiet periods, which is significantly reduced during storm times. Lastly, we estimate a global energy input from downward fluxes of 1.1 × 10 8 and 5.5 × 10 8 W and the globally deposited energy from net fluxes of 2.0 × 10 7 and 1.6 × 10 8 W for quiet and storm time periods, a factor of 5 and 8 enhancement globally, respectively, but up to an order of magnitude locally near strong crustal fields.
The Role of Premagmatic Rifting in Shaping a Volcanic Continental Margin: An Example From the Eastern North American Margin
Both magmatic and tectonic processes contribute to the formation of volcanic continental margins. Such margins are thought to undergo extension across a narrow zone of lithospheric thinning (~100 km). New observations based on existing and reprocessed data from the Eastern North American Margin contradict this hypothesis. With ~64,000 km of 2-D seismic data tied to 40 wells combined with published refraction, deep reflection, receiver function, and onshore drilling efforts, we quantified along-strike variations in the distribution of rift structures, magmatism, crustal thickness, and early post-rift sedimentation under the shelf of Baltimore Canyon Trough (BCT), Long Island Platform, and Georges Bank Basin (GBB). Results indicate that BCT is narrow (80–120 km) with a sharp basement hinge and few rift basins. The seaward dipping reflectors (SDR) there extend ~50 km seaward of the hinge line. In contrast, the GBB is wide (~200 km), has many syn-rift structures, and the SDR there extend ~200 km seaward of the hinge line. Early post-rift depocenters at the GBB coincide with thinner crust suggesting “uniform” thinning of the entire lithosphere. Models for the formation of volcanic margins do not explain the wide structure of the GBB. We argue that crustal thinning of the BCT was closely associated with late syn-rift magmatism, whereas the broad thinning of the GBB segment predated magmatism. Correlation of these variations to crustal terranes of different compositions suggests that the inherited rheology determined the premagmatic response of the lithosphere to extension.
Does Magnetic Reconnection Occur in the Near Lunar Surface Environment?
The near lunar surface contains small-scale magnetic field structures that provide a natural test bed for observing plasmas with a non-zero Hall electric field, as well as potentially facilitating electron-only reconnection. This study presents observational evidence of magnetized electrons as well as demagnetized ions when THEMIS-ARTEMIS probe B reached an altitude of ~15 km above the lunar surface. Additionally, observations suggest the presence of a field line topology change and traversal of a closed magnetic field structure containing solar wind electrons, suggestive of magnetic reconnection having occurred at some point between the solar wind interplanetary magnetic field and a lunar crustal magnetic field. Thus, the observations presented here are consistent with previous studies that predict prominent Hall electric fields near lunar crustal magnetic fields and further suggest that the solar wind interplanetary magnetic field may reconnect with lunar crustal magnetic fields, most likely via electron-only reconnection.
Understanding relative metal prices and availability: Combining physical and economic perspectives
Abstract Scientists typically use physical indicators, such as average crustal abundance, and energy requirements in production, as measures of the availability of mineral resources. Economists, on the other hand, typically rely on measures such as prices, and extraction costs. This study investigates the role of crustal abundance, energy requirements, and other physical indicators in determining relative metal prices. The relationship between prices, and physical indicators is estimated statistically for a sample of 22 chemical elements over the period 1970–2013 using an economic market model. The results suggest that energy requirements in production explain 43% of observed variation in metal prices, crustal abundance 21%, and other physical indicators (toxicity, native metal status, and melting point) combine to explain 12% of the observed differences in metal prices.
Oxygen isotope ratios in zircon and garnet: A record of assimilation and fractional crystallization in the Dinkey Dome peraluminous granite, Sierra Nevada, California
Abstract The 119 Ma Dinkey Dome pluton in the central Sierra Nevada Batholith is a peraluminous granite and contains magmatic garnet and zircon that are complexly zoned with respect to oxygen isotope ratios. Intracrystalline SIMS analysis tests the relative importance of magmatic differentiation processes vs. partial melting of metasedimentary rocks. Whereas δ18O values of bulk zircon concentrates are uniform across the entire pluton (7.7‰ VSMOW), zircon crystals are zoned in δ18O by up to 1.8‰, and when compared to late garnet, show evidence of changing magma chemistry during multiple interactions of the magma with wall rock during crustal transit. The evolution from an early high-δ18O magma [δ18O(WR) = 9.8‰] toward lower values is shown by high-δ18O zircon cores (7.8‰) and lower δ 18O rims (6.8‰). Garnets from the northwest side of the pluton show a final increase in δ18O with rims reaching 8.1‰. In situ REE measurements show zircon is magmatic and grew before garnets. Additionally, δ18O in garnets from the western side of the pluton are consistently higher (avg = 7.3‰) relative to the west (avg = 5.9‰). These δ18O variations in zircon and garnet record different stages of assimilation and fractional crystallization whereby an initially high-δ18O magma partially melted low-δ18O wallrock and was subsequently contaminated near the current level of emplacement by higher δ18O melts. Collectively, the comparison of δ18O zoning in garnet and zircon shows how a peraluminous pluton can be constructed from multiple batches of variably contaminated melts, especially in early stages of arc magmatism where magmas encounter significant heterogeneity of wall-rock assemblages. Collectively, peraluminous magmas in the Sierran arc are limited to small <100 km2 plutons that are intimately associated with metasedimentary wall rocks and often surrounded by later and larger metaluminous tonalite and granodiorite plutons. The general associations suggest that early-stage arc magmas sample crustal heterogeneities in small melt batches, but that with progressive invigoration of the arc, such compositions are more effectively blended with mantle melts in source regions. Thus, peraluminous magmas provide important details of the nascent Sierran arc and pre-batholithic crustal structure.
Modeling Solids in Nuclear Astrophysics with Smoothed Particle Hydrodynamics
Abstract Smoothed particle hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust, which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches that can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations.