Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geoscience”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

NRAP-open-IAM: A flexible open-source integrated-assessment-model for geologic carbon storage risk assessment and management

Large-scale implementation of geologic carbon storage (GCS) to help reduce atmospheric greenhouse gas emissions requires stakeholder confidence that injected CO2 will remain contained and that potential subsurface environmental risks are acceptably small and manageable. The U.S. Department of Energy’s National Risk Assessment Partnership (NRAP) has developed an open-source integrated assessment model (NRAP-Open-IAM) to help address questions about a potential GCS site’s ability to effectively contain injected CO 2 and protect groundwater and other overlying environmentally sensitive receptors. NRAP-Open-IAM allows a user to: (1) incorporate relevant site geologic and injection scenario data; (2) characterize important site features and events;(3) couple fast prediction models of various system components of the engineered geologic system; and (4) execute stochastic, dynamic simulation of whole GCS system performance, leakage risk assessment, and uncertainty quantification. NRAP-Open-IAM is available on GitLab (https://gitlab.com/NRAP/OpenIAM), and is accompanied by multiple application examples and detailed user and developer guides.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Toward more mechanistic representations of biogeochemical processes in river networks: Implementation and demonstration of a multiscale model

A recently introduced multiscale model for representing the effects of small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The model uses advection-dispersion-reaction equations for the channel network and one-dimensional advection-reaction subgrid models for the hyporheic zone. We summarize implementation in the integrated surface/subsurface hydrology modeling system Amanzi-ATS. The extension uses topologically defined meshes to represent stream/river networks and associates a hyporheic subgrid model with each channel grid cell. Biogeochemical reaction modeling capability residing in community software is accessed through an application programming interface. The implementation is verified against independent numerical solutions on a single reach. Mesh convergence studies show that commonly used semi-distributed representations can introduce significant spatial discretization error. Denitrification of farm runoff in a subbasin of the Portage River Basin in Ohio, USA is used to demonstrate the general-purpose reactive transport capability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Thermal and magnetic evolution of a crystallizing basal magma ocean in Earth's mantle

We present the thermochemical evolution of a downward crystallizing basal magma ocean (BMO) overlying the liquid outer core and probe its capability to dissipate enough power to generate and sustain an early dynamo. Furthermore, a total of 61 out of 112 scenarios for a BMO with imposed, present-day heat flux, Q BMO , values of 15, 18, and 21 TW and radiogenic heat, Q r , values of 4, 8, and 12 TW fully crystallized during the age of the Earth. Most of these models are energetically capable of inducing magnetic activity for the first 1.5 Gyrs, at least, with durations extending to 2.5 Gyrs; with final core-mantle boundary (CMB) temperatures of 4400 ±500 K-well within current best estimates for inferred temperatures. None of the models with Q BMO = 12 TW achieved a fully crystallized state, which may reflect a lower bound on the present-day heat flux across the CMB. BMO-powered dynamos exhibit strong dependence on the partition coefficient of iron into the liquid layer and its associated melting-point depression for a lower mantle composition at near-CMB conditions-parameters which are poorly constrained to date. Nonetheless, we show that a crystallizing BMO is a plausible mechanism to sustain an early magnetic field.

58 GEOSCIENCES↗
DOE OSTI2020

Large isotopic variability at the micron-scale in ‘Shuram’ excursion carbonates from South Australia

Ediacaran-aged (635–541 million years ago) marine sediments contain a large negative carbon isotope (δ 13 C) excursion, in which carbonate δ 13 C values reach -12‰ (VPDB). Known as the ‘Shuram’ excursion, many workers have interpreted this δ 13 C record as an unprecedented perturbation to the global carbon cycle, leading to speculation about a causal connection to the broadly contemporaneous rise of animal life. Others have interpreted the δ 13 C signal as a product of diagenesis, thereby minimizing its relevance for understanding the evolution of metazoans. In this work, we present SEM imaging and in-situ δ 13 C and δ 18 C values measured by secondary ion mass spectrometry (SIMS) to assess these competing hypotheses in the Wonoka Formation of South Australia. Our results from the minimum of the excursion show that rounded sedimentary grains of calcite have δ 13 C values between -12.8 to -10.6‰ and δ 18 C values between 17.8 to 15.5‰ (VPDB). Euhedral dolomite that appears to have grown unimpeded in open sedimentary pore spaces also is present. These early-stage dolomites are interpreted as early authigenic in origin and have δ 13 C values that reach +5‰, requiring a formation fluid with a substantially different δ 13 C composition from basin waters or bulk sediment. Together, these results provide little evidence for the hypothesis that a late diagenetic overprint has generated the ‘Shuram’ excursion in the Wonoka. Instead, they suggest the presence of a large carbon isotopic gradient in the surface environment, with shallow waters capable of precipitating carbonates with very low δ 13 C (down to 12‰) and deeper shelf and/or marine pore waters generating carbonates with positive carbon isotope values (up to +5‰). Because negative isotope excursions of similar magnitude are found in widely dispersed Ediacaran basins, it is likely that this gradient was characteristic of shelf environments of this period and that a still-unknown global process led to the ‘Shuram’ excursion in shallow water carbonates.

58 GEOSCIENCES↗
DOE OSTI2020

The timing of lunar solidification and mantle overturn recorded in ferroan anorthosite 62237

Ferroan anorthosite suite (FAS) rocks are widely interpreted to represent primordial lunar crust. Despite their importance in pinpointing the timing of lunar crust formation, robust chronological investigations for this rock type are scarce. Here, we report the Ar-Ar, Rb-Sr, and Sm-Nd isotopic systematics for the FAS troctolitic anorthosite 62237. The Ar-Ar isotopic system has been reset by a thermal event at 3710 ± 48 Ma, and the Rb-Sr isotopic systematics has been disturbed such that a Rb-Sr isochron age cannot be determined. However, an internal isochron for the Sm-Nd isotopic system has yielded an age of 4350 ± 73 Ma (MSWD = 2.0) with an initial ε 143 Nd CHUR of -0.53 ± 0.26. The mineral and whole-rock fractions of 62237 plot on the same internal isochron as FAS sample 60025. The combined datasets define an age of 4372 ± 35 Ma (MSWD = 4.0) with an initial ε 143 Nd CHUR of -0.17 ± 0.22. Literature Sm-Nd data for FAS and Mg-suite whole-rocks also plot on the 60025-62237 isochron. The coherence of data from both FAS and Mg-suite rocks examined thus far suggests that both rock suites formed contemporaneously from identical, or nearly identical, sources. In addition, the concordance of FAS and Mg-suite ages suggests that primordial crust solidification either involved both magmatic suites, or that Mg-suite magmatism was contemporaneous with FAS magmatism within resolution of the Sm-Nd chronometer. The ages for FAS and Mg-suite also coincide with the formation ages of the mare basalt source regions and urKREEP. Ferroan anorthosite suite rocks and urKREEP are thought to represent primordial LMO solidification products, whereas Mg-suite and the mare basalt source regions are argued to represent mixtures of various LMO crystallization products that were formed during density-driven overturn of the LMO. The concordance of ages implies that the 4372 ± 35 Ma Sm-Nd isochron records the age of mantle overturn, and that overturn occurred during, or shortly after, solidification of the LMO.

58 GEOSCIENCES↗
DOE OSTI2020

A numerical examination of the effect of sulfide dissolution on silicate weathering

Silicate weathering regulates climate on geological time scales as a net, climate-sensitive sink of atmospheric CO 2 . On the other hand, sulfuric acid produced through sulfide dissolution affects silicate weathering, diminishing its effectiveness as a climate regulator based on evidence from river chemistry. This study takes a theoretical approach to quantitatively examine the effect of sulfide dissolution, coupling a one-dimensional model of pyrite weathering to that of albite transformation to kaolinite. The coupled model reveals that when the reaction front of sulfide is deeper than that of silicates reacting with CO 2 , the silicate-weathering feedback on climate is not directly affected by sulfide dissolution, but only indirectly through oceanic processes such as reverse weathering and carbonate deposition. In turn, when sulfide dissolution occurs within zones of CO 2 -silicate reactions close to the surface, the feedback between climate and silicate weathering can be significantly weakened because sulfuric acid competes effectively with carbonic acid in mineral dissolution. As the depth of the sulfide reaction front can be changed by variations in water table depth, the latter feedback weakening could have occurred in concert with climate-change-driven water table fluctuations, amplifying atmospheric CO 2 variations during the Phanerozoic.

58 GEOSCIENCES↗
DOE OSTI2020

U-Pb zircon age constraints on the earliest eruptions of the Deccan Large Igneous Province, Malwa Plateau, India

Climate instability driven by emission of volatiles during emplacement of large igneous provinces (LIPs) is frequently invoked as a potential cause of mass extinctions. However, documenting this process in the geologic record requires a holistic understanding of eruption rates, the location of eruptive centers, and potential sources of climate-changing volatiles. We present new chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) U-Pb zircon geochronology from Malwa Plateau basalts on the northern margin of the Deccan LIP, India. These basalts have been previously interpreted as either an extension of the province's main volcanic stratigraphy or as an independent eruptive center active up to millions of years prior to the main eruptive phase. Our data instead demonstrate that the lower Malwa Plateau basalts are temporally correlative with the first pulse of Deccan volcanism and provide new constraints on its initiation and duration. Paleomagnetic data further indicate that upper Malwa Plateau basalts may be age-equivalent to the second, third, and fourth pulses of Deccan volcanism. The relative thicknesses of age-equivalent packages of basalt are consistent with eruption of the Deccan LIP from a southward-migrating eruptive center. The first eruptive pulse is coeval with a ~200 kyr Late Maastrichtian warming event preserved globally in contemporaneous stratigraphic sections. We propose that the first pulse of Deccan magmatism was more voluminous in the north, where it erupted through organic-rich sedimentary rocks of the Narmada-Tapti rift basin. Thermal metamorphism of these sediments could have been a source of sufficient CO 2 to drive the Late Maastrichtian warming event, which if true would reconcile the apparent dampened warming signals associated with later Deccan eruptive pulses.

58 GEOSCIENCES↗
DOE OSTI2020

HIMU geochemical signature originating from the transition zone

Plume volcanism may sample mantle sources deeper than mid-ocean ridge and arc volcanism. Ocean island basalts (OIBs) are commonly related to plume volcanism, and their diverse isotopic and elemental compositions can be described using a limited number of mantle endmembers. However, the origins and depths of these mantle endmembers are highly debated. Here we show that the HIMU (high μ, μ = 238 U/ 204 Pb) endmember may reside in the transition zone. Specifically, we report the geochemical signature of a high-pressure multiphase diamond inclusion, entrapped at 420–440 km depth and 1450 ±50 K, which matches exactly the geochemical patterns of the HIMU-rich OIBs. Since the HIMU component is variably sampled by almost all OIBs, our finding implies that the transition zone causes a major overprint of the geochemical features of mantle plumes. Furthermore, some mantle plumes, like those feeding Bermuda, St Helena, Tubuai and Mangaia, appear to be dominated by this source. Furthermore, our finding highlights the importance of the transition zone in highly incompatible element budget of the mantle.

58 GEOSCIENCES↗
DOE OSTI2020

Predicting the proximity to macroscopic failure using local strain populations from dynamic in situ X-ray tomography triaxial compression experiments on rocks

Predicting the proximity of large-scale dynamic failure is a critical concern in the engineering and geophysical sciences. Here we use evolving contractive, dilatational, and shear strain deformation preceding failure in dynamic X-ray tomography experiments to examine which strain components best predict the proximity to failure. We develop machine learning models to predict the proximity to failure using time series of three-dimensional local incremental strain tensor fields acquired in rock deformation experiments under stress conditions of the upper crust. Three-dimensional scans acquired in situ throughout triaxial compression experiments provide a distribution of density contrasts from which we estimate the three-dimensional incremental strain that accumulates between each scan acquisition. Training machine learning models on multiple experiments of six rock types provides suites of feature importance that indicate the predictive power of each feature. Comparing the average importance of groups of features that include information about each strain component quantifies the ability of the contractive, dilatational and shear strain to predict the proximity of macroscopic failure. A total of 24 models of four machine learning algorithms with six rock types indicate that 1) the dilatational strain provides the best predictive power of the strain components, and 2) the intermediate values (25th-75th percentile) of the strain population provide the best predictive power of the statistics of the strain populations. In addition, the success of the predictions of models trained on one rock type and tested on other rock types quantifies the similarities and differences of the precursory strain accumulation process in the six rock types. These similarities suggest the potential existence of a unified theory of brittle rock deformation for a range of rock types.

58 GEOSCIENCES↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2020

A new hydrous iron oxide phase stable at mid-mantle pressures

The amount of hydrogen stored in the Earth's interior is important for a range of issues, including the volatile incorporation during the Earth formation and the co-evolution of the atmosphere, the hydrosphere, and the interior. Recent experiments found titanium bearing ε-FeOOH in a hydrous basaltic system at 12–19 GPa and 1300 K. Pyrite-type FeOOH was found to be stable at pressures higher than 80 GPa. These discoveries suggest possible hydrogen storage in the mantle transition zone and in the mantle below 1800 km depths, respectively. However, it remains uncertain how the potential deep hydrogen storage can be connected to the shallower storage. Here, we report a new hydrous iron oxide (η-Fe 12 O 18+x / 2 H x , x ≈ 2) stable at pressures between the stability fields of the ε- and the pyrite-type FeOOH. Our experiment also shows that the new η phase can exist together with the major lower mantle minerals including bridgmanite and periclase, making it an important hydrogen-bearing phase in the Earth's deep interior. Because of its limited H 2 O storage capacity, which is less than 1/6 of the storage capacity of the pyrite-type phase and the ε phase, the stability of the η phase would result in H 2 O loss during water transport in the mid mantle and therefore limit the amount of H 2 O potentially stored in the Fe–O–H system of the lower mantle. As a result, the large channel in the crystal structure of the η phase could provide potential storage sites for other volatile elements in the deep mantle.

58 GEOSCIENCES↗
DOE OSTI2020

Pervasive low-velocity layer atop the 410-km discontinuity beneath the northwest Pacific subduction zone: Implications for rheology and geodynamics

Regional triplication waveforms of five intermediate-depth events are modeled to simultaneously obtain the compressional (P) and shear (SH) wave velocity structure beneath northwestern Pacific subduction zone. Both the P- and SH-wave velocity models for three different sub-regions show a low-velocity layer (LVL) with a thickness of ~55-80 km lying above the 410-km discontinuity with a ~900 km lateral extent from the Japan Sea to the northeastern Asian continental margin. With the dihedral angle approaching to zero around 400 km, a minute amount of melt atop the 410-km discontinuity caused by the hydrous slab might completely wet olivine grain boundaries and result in a low seismic velocity layer in this specific subduction zone. This mechanism suggests that the 410-LVL is a low viscosity zone that would partially decouple the upper mantle from the transition zone. We infer that the widespread 410-LVL provides evidence for a water-bearing mantle transition zone beneath the western Pacific subduction zone.

58 GEOSCIENCES↗
DOE OSTI2020

Measuring multiple cosmogenic nuclides in glacial cobbles sheds light on Greenland Ice Sheet processes

The behavior of the Greenland Ice Sheet during the Pleistocene remains uncertain due to the paucity of evidence predating the Last Glacial Maximum. In this study, we employ a novel approach, cosmogenic nuclide analysis of individual subglacially-derived cobbles, which allows us to make inferences about ice sheet processes and subglacial erosion. From three locations in western Greenland, we collected 86 cobbles from the current ice sheet margin and nine cobbles exposed on the modern proglacial land surface. We measured the concentration of in situ 10 Be in all cobbles (n = 95) and 26 Al and 14 C in a subset (n = 14). Cobbles deposited during Holocene retreat have 10 Be exposure ages generally consistent with the timing of ice retreat determined by other methods. Conversely, most of the 86 subglacial cobbles contain very low concentrations of 10 Be (median 1.0×10 3 atoms g –1 ), although several have ~10 4 and one has ~10 5 atoms g –1 . The low concentrations of 10 Be in most subglacial cobbles imply that their source areas under the Greenland Ice Sheet are deeply eroded, preserving minimal evidence of surface or near-surface exposure. The presence of measurable 14 C in ten of the cobbles requires that they experienced cosmogenic nuclide production within the past ~30 ka; however, 14 C/ 10 Be ratios of ~6 suggest that nuclide production occurred during shielding by overlying material. Only two of the 86 subglacial cobbles definitively have cosmogenic nuclide concentrations consistent with prior surface exposure. Overall, isotopic analysis of subglacial cobbles indicates that much of western Greenland's subglacial landscape is characterized by deep erosion and minimal subaerial exposure.

58 GEOSCIENCES↗
DOE OSTI2021

Isotopic signatures as tools to reconstruct the primordial architecture of the Solar System

The significant reorganization of the early Solar System due to giant planet migration has hampered our understanding of where planetary bodies formed. Previously employed proxies for reconstructing the primordial planetary architecture, such as water content or oxidation state, are complicated by post-accretionary processes. Furthermore we investigate basaltic achondrites for their nucleosynthetic isotope signatures in the elements neodymium (Nd) and zirconium (Zr) and show that they are—similar to previously investigated chondritic meteorites—characterized by a relative deficit in isotopes produced by the s-process of nucleosynthesis. Importantly, these data are well correlated with nucleosynthetic signatures observed in other elements, demonstrating that s-process matter was heterogeneously distributed throughout the early Solar System. By comparing these isotopic signatures with potential proxies for Solar System reconstruction and computer modeling, we here argue that this isotopic heterogeneity in bulk meteoritic materials is linked to the original heliocentric distance of formation. Such scaling of nucleosynthetic signatures with heliocentric distance could permit reconstruction of the primordial architecture of the Solar System by ‘cosmolocating’ the accretion orbits of meteoritic parent bodies as a function of incorporated s-process matter.

58 GEOSCIENCES↗
DOE OSTI2021

Thermal expansion of liquid Fe-S alloy at high pressure

The local structure and density of liquid Fe-S alloys at high pressure have been determined in situ by combined angle and energy dispersive X-ray diffraction experiments in a multi-anvil apparatus, covering a large temperature and compositional range. Precise density measurements collected for increasing temperature allowed us to directly derive the thermal expansion coefficients for liquid Fe-S alloys as a function of composition. In turn, thermal expansion has been used to refine thermodynamic models and to address the crystallization regime of telluric planetary cores by comparing the adiabatic temperature gradient and the slope of the liquidus in the Fe-FeS system. For Fe-S cores of asteroids and small planetesimals, top-down solidification is the dominant scenario as the compositional domain for which the slope of the liquidus is greater than the adiabatic gradient is limited to a narrow portion on the Fe-rich side. However, bottom-up growth of the inner core is expected for S-poor cases, with this compositional domain expanding to more S-rich compositions with increasing pressure (size of the planetary body). In particular, bottom-up crystallization cannot be excluded for the Moon and Ganymede.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Origin of water masses in Floridan Aquifer System revealed by 81 Kr

Groundwater in coastal aquifers serves as an essential resource in densely populated areas throughout the world. However, hydrological connection to the ocean leaves it at risk of salinization via changes in climate and hydrological cycles. Therefore, an accurate hydrological characterization of coastal aquifers is of the utmost importance. Although localized salinization of shallow aquifers has commonly been studied, geochemical constraints on the time scale of regional freshwater and seawater in coastal aquifers are limited especially in carbonate aquifers. Such information is not only crucial for water resource management, but also for understanding solute flux over the land-ocean boundary and for interpreting the geochemical signatures that each water mass carries. Toward this goal, we report tracer-based, subsurface residence times of groundwater in the Floridan Aquifer System in southern Florida, which is one of the most productive aquifers on Earth. The first application of Kr 81-an emerging tracer of old groundwater-to this aquifer has identified freshwater recharge during the last glacial period as anticipated from previous studies using other geochemical tools including radiocarbon, stable isotopes and noble gas concentrations. Moreover, freshwater of Holocene age was also detected at 100 km-distance from the recharge area, suggesting the possibility of an active flow system in the upgradient region. A contribution of fossil seawater that predates the last glacial maximum was locally identified, suggesting slow seawater circulation and limited but significant solute flux from the ocean into the aquifer for dolomitization. This also implies that the deep saline coastal groundwater can potentially serve as a paleo-seawater archive.

58 GEOSCIENCES↗
DOE OSTI2021

The gallium isotopic composition of the Moon

Here, in this study, we present new Ga isotope data from a suite of 28 mare basalts and lunar highland rocks. The δ 71 Ga values of these samples range from -0.10 to +0.66‰ (where δ 71 Ga is the relative difference between the 71 Ga/ 69 Ga ratio of a sample and the Ga-IPGP standard), which is an order of magnitude more heterogeneous than δ 71 Ga values in terrestrial magmatic rocks. The cause of this isotopic heterogeneity must be established to estimate the bulk δ 71 Ga value of the Moon. In general, low-Ti basalts and ferroan anorthosite suite (FAS) rocks have δ 71 Ga values that are lower than high-Ti basalts and KREEP-rich rocks. The observation that rocks derived from later forming LMO cumulates have higher δ 71 Ga values suggests that Ga isotopes are fractionated by processes that operate within the chemically evolving LMO, rather than localized degassing or volatile redistribution. Correlations between indices of plagioclase removal from the LMO (e.g. Eu/Eu*) with Ga isotope ratios suggest that a Δ 71 Ga plagioclase-melt of -0.3‰, (where Δ 71 Ga plagioclase-melt is the isotopic fractionation associated with crystallization of plagioclase from a melt), could drive the observed isotopic fractionation in high-Ti mare basalts and KREEP-rich rocks. This would be consistent with the observation that FAS rocks have δ 71 Ga values that are lower than mare basalts. However, the addition of KREEP-like material into the mare basalt source regions would not contribute enough Ga to perturb the isotopic composition outside of analytical uncertainty. Thus, basalts derived from early formed LMO cumulates such as those from Apollo 15, would preserve light Ga isotopic compositions despite containing modest amounts of urKREEP. We estimate that the δ 71 Ga value of the LMO was ~0.14‰ prior to the onset of plagioclase crystallization and extraction. Whether this δ 71 Ga value is representative of the initial BSM cannot be ascertained from the current dataset. It remains plausible that the Moon accreted with a heavier Ga isotopic composition than the Earth. Alternatively, the Moon and Earth could have accreted with similar isotopic compositions (BSE = 0.00 ± 0.06‰, Kato et al., 2017) and volatile loss drove the LMO to higher δ 71 Ga values prior to formation of the lunar crust.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗