Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic repository”

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 199 records · Page 11

Temperature effects on adsorption and capillarity water retention mechanisms in constrained unsaturated soils

Abstract This paper focuses on the impact of elevated temperatures on the adsorptive and capillarity water retention mechanisms of unsaturated soils under constrained (constant volume) conditions. This topic is critical for simulating the thermo-hydraulic behavior of soils in hydrogeological or geotechnical applications, including climate change effects on near surface soils, energy piles or soil borehole thermal energy storage systems in unsaturated soil layers, and buffers for geological nuclear waste repositories. A nonisothermal soil water retention curve (SWRC) that separately considers the temperature-dependency of the key parameters governing adsorptive and capillarity water retention mechanisms and soil physical parameters (e.g., surface tension, contact angle, adsorption capacity, cation exchange capacity, mean cavitation suction, air entry value and equilibrium film thickness) was developed to provide insights into the impact of temperature on water retention over the full suction range. The nonisothermal SWRC was validated using experimental data on high plasticity clays, with a good prediction of temperature effects on adsorption and capillarity water retention mechanisms in constrained unsaturated soils.

Lu, Yu↗

High-Temperature Thermodynamics of Cerium Silicates, A-Ce 2 Si 2 O 7 , and Ce 4.67 (SiO 4 ) 3 O

Lanthanide disilicates and oxyapatites have potential roles in high temperature applications as thermal (TBC) and environmental barrier coatings (EBC), or possible alteration phases in geological nuclear waste repositories. However, those Ce 3+ -bearing silicates have only been limitedly studied. In this work, we performed detailed structural and thermodynamic investigations on A-Ce 2 Si 2 O 7 (tetragonal, P4 1 ) and Ce 4.67 (SiO 4 ) 3 O (hexagonal, P6 3 /m). The high temperature structural behaviors and coefficients of thermal expansion were determined by in situ high temperature synchrotron X-ray diffraction (HT-XRD) implemented with Rietveld analysis and thermogravimetric analysis coupled with differential scanning calorimetry (TGA-DSC). A-Ce 2 Si 2 O 7 was found to be stable in N 2 and air up to ~1483 K with anisotropic thermal expansion along the a and c axes (α a = 12.3 × 10 -6 K -1 and α c = 12.4 × 10 -6 K -1 ). Ce 4.67 (SiO 4 ) 3 O had a slow partial oxidation between 533 K and 873 K to a new nonstoichiometric phase Ce 3+ 1.67- xCe4 + x Ce 3+ 3 (SiO 4 ) 3 O 1+0.5x , followed by a thermal decomposition to CeO 2 and SiO 2 at ~1000 K in air. By using high temperature oxide melt solution calorimetry at 973 K with lead borate as the solvent, the standard enthalpy of formation was determined for A-Ce 2 Si 2 O 7 (-3825.1 ± 6.0 kJ/mol) and Ce 4.67 (SiO 4 ) 3 O (-7391.3 ± 9.5 kJ/mol). Finally, these thermodynamic parameters were compared with those of CeO 2 , CeSiO 4 , and other silicate oxyapatites for examining their chemical stability in high temperature environments relevant for aeronautical applications, mineral formation, and nuclear fuel cycle.

36 MATERIALS SCIENCE↗

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Water transport pathway in clay interlayer upon dehydration [Slides]

Smectite (e.g., Montmorillonite): phyllosilicate minerals found in bentonites. Bentonites have been considered as key backfill barrier materials in deep geological nuclear waste repository concepts. Swelling/shrinking of montmorillonite (MMT) occurs with increasing/decreasing relative humidity. Our research question is, "Microscopically, how does the hydration/dehydration process occur?"

36 MATERIALS SCIENCE↗

Water transport pathway in clay interlayer upon dehydration [Slides]

Smectite (e.g., Montmorillonite): phyllosilicate minerals found in bentonites. Bentonites have been considered as key backfill barrier materials in deep geological nuclear waste repository concepts. Swelling/shrinking of montmorillonite (MMT) occurs with increasing/decreasing relative humidity. Microscopically, how does the hydration/dehydration process occur?

36 MATERIALS SCIENCE↗

DPC Direct Disposal Postclosure Thermal Modeling

Performance of geologic radioactive waste repositories depends on near-field and far-field processes, including km-scale flow and transport in engineered and natural barriers, that may require simulations of up to 1 M years of regulatory period. For a relatively short time span (less than 1000 years), the thermohydro-mechanical-chemical (THMC) coupled processes caused by heat from the waste package will influence near-field multiphase flow, chemical/reactive transport, and mechanical behaviors in the repository system. This study integrates the heat-driven perturbations in thermo-hydro-mechanical characteristics into thermo-hydro-chemical simulations using PFLOTRAN to reduce dimensionality and improve computational efficiency by implementing functions of stress-dependent permeability and saturation-temperature-dependent thermal conductivity. These process couplings are developed for spent nuclear fuel in dual-purpose canisters in two different hypothetical repositories: a shale repository and a salt repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Issues in Defining Software Architectures in a GIS Environment

The primary mission of the Pan-American Center for Earth and Environmental Studies (PACES) is to advance the research areas that are relevant to NASA's Mission to Planet Earth program. One of the activities at PACES is the establishment of a repository for geographical, geological and environmental information that covers various regions of Mexico and the southwest region of the U.S. and that is acquired from NASA and other sources through remote sensing, ground studies or paper-based maps. The center will be providing access of this information to other government entities in the U.S. and Mexico, and research groups from universities, national laboratories and industry. Geographical Information Systems(GIS) provide the means to manage, manipulate, analyze and display geographically referenced information that will be managed by PACES. Excellent off-the-shelf software exists for a complete GIS as well as software for storing and managing spatial databases, processing images, networking and viewing maps with layered information. This allows the user flexibility in combining systems to create a GIS or to mix these software packages with custom-built application programs. Software architectural languages provide the ability to specify the computational components and interactions among these components, an important topic in the domain of GIS because of the need to integrate numerous software packages. This paper discusses the characteristics that architectural languages address with respect to the issues relating to the data that must be communicated between software systems and components when systems interact. The paper presents a background on GIS in section 2. Section 3 gives an overview of software architecture and architectural languages. Section 4 suggests issues that may be of concern when defining the software architecture of a GIS. The last section discusses the future research effort and finishes with a summary.

Acosta, Jesus↗

Evaluation of Techniques for Vegetation Removal from UAV-Based Photogrammetric DSM near WIPP Land Withdrawal Act Boundary, NM

The Waste Isolation Pilot Plant is located within New Mexico's karst landscape and is the United States' only deep geologic radioactive waste repository which isolates transuranic waste from defense activities underground in a bedded salt formation. Characteristic karst topography in this region is made up of material (such as limestone) that can be dissolved by water over a long period of time, and includes features such as springs, surface streams, sinking streams, caves, and sinkholes. Karst features can potentially impact the integrity of the waste repository in the future and compromise WIPP's performance, particularly when coupled with incompatible land-use activities within the Land Withdrawal Act (LWA) boundary as well as increased water withdrawals outside the LWA. In order to better capture the ground surface topography to support hydrology modeling efforts, the U.S. DOE has enlisted FIU to develop a high-resolution digital elevation model (DEM). In a previous field study, a high-resolution digital surface model (DSM) of an adjacent representative site, Basin 6, was developed by photogrammetry using aerial images captured with unmanned aerial vehicles. The DSM generated was based on point clouds representing the entire landscapes, including points of terrain, vegetation and infrastructure. These above-ground features need to be removed to create a DSM/DEM with good representation of the bare ground. In recent years, many techniques including both pre- and post-processing of DSMs, were developed for vegetation removal with varied success. This research focuses on applying and comparing various vegetation removal methods using point clouds generated of the Basin 6 pilot study area. The 3D RGB-based point cloud classification using Python scripts was a modification of the procedure by Themistocleous (2019). The method using Liblas classifies the point cloud three-dimensionally, whereas the method by Themistocleous (2019) is two-dimensional. The preliminary result of the method using Liblas, however, still had a large amount of vegetation remaining compared to the Pix4D Machine Learning method. The vegetation layer remains highly visible in the DEMs processed using the vegetation indices as opposed to the Pix4D Machine Learning method, which produces a DEM where mostly bare ground is seen. Further modifications to the input parameters are required to improve the vegetation vs. bare ground classification using this approach. The next step will be to improve the method using Liblas and to also test other vegetation removal methods found in the literature review to compare them against the two described here.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Panel Session 24B: Records, Knowledge, and Memory for Radioactive Waste Repositories: Generational Equity Focus

This panel focused on the latest thoughts, ideas, and methodologies being explored throughout the world on how to communicate with future generations regarding nuclear waste disposal. Scientists determined many years ago that geological disposal in a repository was the preferred solution for nuclear waste disposal given the longevity concerns of the waste(s). Future generations must be informed through records, memory keeping and permanent markers to ensure they are aware of, and knowledgeable of, the dangers associated with nuclear waste isolated from the biosphere. Panelists with presentations: You Want to Drill Where? Human Intrusion Messaging Considerations (Thomas Peake, Jonathan Major); Ethical Reflections on the Basic Reasons for RK and M Measures (Carl-Reinhold Brakenhielm); NEA Activities on Information, Data and Knowledge Management (Rebecca Tadesse); Records, Knowledge and Memory (RK and M) Across Generations: Recent Activities and Progress in Sweden (Claudio Pescatore)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DECOVALEX-2023: Task B Final Report

In all repository concepts for the geological disposal of radioactive waste, an engineered barrier system (EBS) is used to encapsulate the waste canister, or, to act as borehole or gallery seals. These systems are often based on bentonite clays due to their low permeability and high swelling capacity enabling the closure of engineering voids. However, in all repository concepts gases will be generated through the corrosion of metallic materials (under anoxic conditions), the radioactive decay of waste and the radiolysis of water. Thus, understanding the processes and mechanisms controlling the advective movement of gas (as a discrete phase) in clay-based materials is a key aspect when assessing the impact of gas flow in a repository safety case.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

GDSA Repository Systems Analysis Investigations in FY 2023

This report describes specific activities in the Fiscal Year (FY) 2023 associated with the Geologic Disposal Safety Assessment (GDSA) Repository Systems Analysis (RSA) work package funded by the Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy Office of Nuclear Energy (DOE-NE), Office of Spent Fuel and Waste Disposition (SFWD).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Numerical simulation of coupled THM behaviour of full-scale EBS in backfilled experimental gallery in the Horonobe URL

Bentonite-based engineered barrier system (EBS) is a key component of many repository designs for the geological disposal of high-level radioactive waste. Given the complexity and interaction of the phenomena affecting the barrier system, coupled thermo-hydro-mechanical (THM) numerical analyses are a potentially useful tool for a better understanding of their behaviour. In this context, a Task (the Horonobe EBS experiment) was undertaken to study, using numerical analyses, the thermo-hydro-mechanical (and thermo-hydro) interactions in bentonite based engineered barriers within the international cooperative project DECOVALEX 2023. One full-scale in-situ experiment and four laboratory experiments, largely complementary, were selected for modelling. The Horonobe EBS experiment is a temperature-controlled non-isothermal experiment combined with artificial groundwater injection. The Horonobe EBS experiment consists of the heating and cooling phases. Six research teams performed the THM or TH (depended on research team approach) numerical analyses using a variety of computer codes, formulations and constitutive laws. Finally, for each experiment, the basic features of the analyses are described and the comparison between calculations and laboratory experiments and field observations are presented and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Deep Isolation-Development of the Safety Case for Disposal of Radioactive Wastes in Horizontal Boreholes - 20028

Deep Isolation has developed a safe, secure, and permanent geological disposal method for high-level waste, including spent nuclear fuel as well as sealed sources and other highly radioactive materials. The method leverages well established directional drilling technology to create horizontal repositories deep underground. The combination of great depth (1-3 km) and the ability to precisely position repositories in a horizontal orientation provides access to geologic strata that are inaccessible for typical mined and deep vertical borehole repositories. A number of potential direct and indirect safety benefits accrue including:1) greater depth below regional freshwater aquifers; 2) increased flexibility and more geologic options for siting; 3) access to deep formations with persistent and sustained reducing conditions. In addition there are a number of safety elements related specifically to the horizontal repository geometry including: 1) passive direction of thermally driven fluid and radionuclide movement away from the vertical access hole and toward the 'dead end' portion of the repository; 2) mitigation of seismic hazards by orienting repositories parallel to local and regional fault structures. In this paper we explore and discuss some of the key technologic, geologic and hydrologic elements that support the deep horizontal borehole safety case. The stalemate seen across the globe on the disposal of spent nuclear fuel and high-level waste can be broken. Deep Isolation offers a novel option for safe, secure, and permanent deep geological disposal of nuclear waste that can be developed as centralized repositories or adapted to smaller regional or site-specific repositories located near waste sources. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗

Thermal loading effects on chalk hydromechanical behavior for nuclear waste disposal

Safe disposal of heat-generating nuclear waste depends on host rock stability under thermal, hydrological, and mechanical stresses. This study investigates the effect of thermal loading on mechanical behavior of the shallowly buried Ghareb formation chalk through triaxial and hydrostatic constant strain rate and creep tests at temperatures up to 100 ˚C and effective pressures up to 20.7 MPa. Experimental results show that thermal loading reduces the elastic moduli of chalk by 50–75%, and a transition occurs above 60 ˚C where creep rates increase rapidly. Water saturation nearly doubles the thermally induced strain compared to dry conditions and strongly decreases material rigidity. Thermal loading also leads to significant pore pressure increases under undrained conditions and leads to reductions in the apparent permeability during drained conditions. Laboratory experimental data were used to parameterize and develop a preliminary constitutive model for predicting future deformation during repository operations in the Ghareb. The strongly coupled effects – mechanical weakening, fluid pressure fluctuations, and permeability modification – demonstrate that elevated repository temperatures will have a pronounced effect on the near field Ghareb behavior during waste disposal operations. The findings indicate that the coupled interactions must be considered in predictive models and repository design to ensure long-term nuclear waste isolation and safety.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Corrosion of U233-Doped Uranium Oxide using Microfluidics Methods

The aim of work this year has been to investigate the role of alpha (α)-radiation and hydrogen (H 2 ) on the corrosion of uranium oxide (UO 2 ) using a microfluidic device. The microfluidic device, termed the Particle-Attached Microfluidic Electrochemical Cell, (PAMEC), enables monitoring of the UO 2 electrochemical corrosion potential (E corr ) under different environments, including de-aerated conditions and in the presence of dissolved H 2 . The Si 3 N 4 window allows us to study morphological and chemical changes under an electron microscope. We have previously demonstrated that the PAMEC matches the results from bulk electrochemical tests with UO 2 [1]. We used the high specific activity uranium (U) isotope, 233 U, (t 1/2 = 160,000 years) incorporated into UO 2 , to generate a localized a-field. The objective of the experiments were to mimic the radiation environment that would be experience at the surface of aged spent nuclear fuel (SNF) during long-term geologic disposal under anoxic conditions. Wittman et al. [2] predicted that in the presence of a pure a-radiation field and under H 2 conditions, the concentration of the radiolytic oxidant H 2 O 2 would be suppressed or even eliminated. In a UO 2 corrosion experiment this would be exhibited through a lowering of the measured UO 2 corrosion potential compared to identical conditions in the absence of dissolved H 2 . We found that the predictions of Wittman and co-workers were supported and that the corrosion potential of the 233 U-doped UO 2 in solution lowered with presence of H 2 gas and increased in the absence of H 2 , when under anoxic conditions. The PAMEC experiments indicate that H 2 O 2 has been eliminated in a solution sparged with H 2 while exposed to an a-radiation field. The corrosion potential of the a-doped 233 U(10%)- 238 UO 2 in a solution sparged with Ar/H 2 matched the corrosion potential of 238 UO 2 in solution sparged with air. This clearly demonstrated that the H 2 O 2 had been eliminated and that the only oxidant present was O 2 in this system in complete agreement with the modeling results of Wittman et al. These results point to the need to improve the Fuel Matrix Degradation (FMD) Process Model training data set that is being used in the FMD surrogate model that is being developed for the repository program. The incorporation of realistic radiation chemistry will improve the scientific basis for the FMD Model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Genomic reconstruction of fossil and living microorganisms in ancient Siberian permafrost

Abstract Background Total DNA (intracellular, iDNA and extracellular, eDNA) from ancient permafrost records the mixed genetic repository of the past and present microbial populations through geological time. Given the exceptional preservation of eDNA under perennial frozen conditions, typical metagenomic sequencing of total DNA precludes the discrimination between fossil and living microorganisms in ancient cryogenic environments. DNA repair protocols were combined with high throughput sequencing (HTS) of separate iDNA and eDNA fraction to reconstruct metagenome-assembled genomes (MAGs) from ancient microbial DNA entrapped in Siberian coastal permafrost. Results Despite the severe DNA damage in ancient permafrost, the coupling of DNA repair and HTS resulted in a total of 52 MAGs from sediments across a chronosequence (26–120 kyr). These MAGs were compared with those derived from the same samples but without utilizing DNA repair protocols. The MAGs from the youngest stratum showed minimal DNA damage and thus likely originated from viable, active microbial species. Many MAGs from the older and deeper sediment appear related to past aerobic microbial populations that had died upon freezing. MAGs from anaerobic lineages, including Asgard archaea, however exhibited minimal DNA damage and likely represent extant living microorganisms that have become adapted to the cryogenic and anoxic environments. The integration of aspartic acid racemization modeling and metaproteomics further constrained the metabolic status of the living microbial populations. Collectively, combining DNA repair protocols with HTS unveiled the adaptive strategies of microbes to long-term survivability in ancient permafrost. Conclusions Our results indicated that coupling of DNA repair protocols with simultaneous sequencing of iDNA and eDNA fractions enabled the assembly of MAGs from past and living microorganisms in ancient permafrost. The genomic reconstruction from the past and extant microbial populations expanded our understanding about the microbial successions and biogeochemical alterations from the past paleoenvironment to the present-day frozen state. Furthermore, we provided genomic insights into long-term survival mechanisms of microorganisms under cryogenic conditions through geological time. The combined strategies in this study can be extrapolated to examine other ancient non-permafrost environments and constrain the search for past and extant extraterrestrial life in permafrost and ice deposits on Mars.

54 ENVIRONMENTAL SCIENCES↗