Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model species”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Validation of the Aeroplanets Model in a Present-Day Earth Simulation

We present a validation of the Aeroplanets/Aerochemistry model (Gronoff et al. 2012; Airapetian et al. 2016), a one-dimensional photochemistry model, with another one-dimensional photochemistry model implemented in Hu et al. 2016. Results from a simulation of the atmosphere of present-day Earth from both models were used in the validation. Realistic initial conditions and concentrations were used in the validation, including surface emissions of CO, CH4, NH3, N2O, NO, SO2, OCS, H2S, and H2SO4. As part of the validation, we compare the vertical profile of concentration of key chemical species, including O3, N2O, NO, NO2, HNO3, HO2, OH, and CH4. Our model was able to create vertical concentration profiles of these species that closely matched the results from Hu et al. 2016. With the successful validation, we plan to use the Aeroplanets model to explore the role of ions and aerosols, formed the interaction of solar energetic particles and the atmosphere of early Earth, in generating key greenhouse gas species and gas species that support the formation of molecules crucial for life.

Bradley Michael Hegyi↗
DOE OSTI2026

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗
DOE OSTI2023

Insights into Constraining Rate Coefficients in Fuel Oxidation Mechanisms Using Genetic Algorithm Optimization

Accurate fuel oxidation mechanisms can enable predictive capabilities that aid in advancing combustion technologies. High-level computational kinetics can yield reasonable rate coefficients with uncertainties, in some cases, below a factor of 2. Computed rate coefficients can be constrained further by optimizing against experimental data. Here, we explore the application of genetic algorithm (GA) optimization to constrain computed rate coefficients in complex fuel oxidation mechanisms in conjunction with temperature-dependent species mole fractions from jet-stirred reactor (JSR) measurements. Cyclohexane is a model candidate for understanding the reactivity of cyclic fuels. In this work, we optimize the rate coefficients of the most recent literature cyclohexane mechanism, which incorporates theoretically computed rate coefficients for the reaction networks stemming from the first and second O 2 addition pathways, against the experimental results of two separate literature JSR studies. Optimization consistency is evaluated by carrying out three GA optimizations: fitting to the temperature-dependent species mole fractions in each JSR experiment separately and simultaneously fitting the species mole fractions in both experiments. Local sensitivity analyses are used to identify five influential low-temperature oxidation reactions for optimization. Although the three optimizations do not yield identical rate coefficients, the direction of change in all five rate coefficients is consistent among the three optimizations. Performance of the models from the three optimizations is assessed against literature ignition delay times with differences in the level of agreement observed among the different optimizations. In conclusion, comparisons are made with our recent optimization work of a cyclopentane oxidation master-equation model against time-resolved species concentrations, and insights and improvements of the strategy for constraining rate coefficients using GA optimization are discussed.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

A mixture of grass–legume cover crop species may ameliorate water stress in a changing climate

Climate change models predict increasing precipitation variability in the mid-latitude regions of Earth, generating a need to reduce the negative impacts of these changes on crop production. Despite considerable research on how cover crops support agriculture in a changing climate, understanding is limited of how climate change influences the growth of cover crops. We investigated the early development of two common cover crop species—crimson clover (Trifolium incarnatum) and rye (Secale cereale)—and hypothesized that growing them in the mixture would ameliorate stress from drought or waterlogging. This hypothesis was tested in a 25-day greenhouse experiment, where the two factors (species number and water stress) were fully crossed in randomized blocks, and plant responses were quantified through survival, growth rate, biomass production and root morphology. Water stress negatively influenced the early growth of these two species in contrasting ways: crimson clover was susceptible to drought while rye performed poorly under waterlogging. Per-plant biomass in rye was always greater in mixture than in monoculture, while per-plant biomass of crimson clover was greater in mixture under drought. Both species grew longer roots in mixture than in monoculture under drought, and total biomass of mixtures did not differ significantly from the more-productive monoculture (rye) in any water condition. In the face of increasingly variable precipitation, growing crimson clover and rye together has potential to ameliorate water stress, a possibility that should be further investigated in field experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

LungMAP Portal Ecosystem: Systems-level Exploration of the Lung

An improved understanding of the human lung necessitates advanced systems models informed by an ever-increasing repertoire of molecular omics, cellular, imaging, and pathological datasets. To centralize and standardize information across broad lung research efforts we expanded the LungMAP.net website into a new gateway portal. This portal connects a broad spectrum of research networks, bulk and single-cell multi-omics data and a diverse collection of image data that span mammalian lung development, and disease. The data are standardized across species and technologies using harmonized data and metadata models that leverage recent advances including those from the Human Cell Atlas, diverse ontologies, and the LungMAP CellCards initiative. To cultivate future discoveries, we have aggregated a diverse collection of single-cell atlases for multiple species (human, rhesus, mouse), to enable consistent queries across technologies, cohorts, age, disease, and drug treatment. These atlases are provided as independent and integrated queryable datasets, with an emphasis on dynamic visualization, figure generation, re-analysis, cell-type curation, and automated reference-based classification of user-provided single-cell genomics datasets (Azimuth). Finally, as this resource grows, we intend to increase the breadth of available interactive interfaces, supported data types, data portals and datasets from LungMAP and external research efforts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗
DOE OSTI2021

Historical Contingency in Microbial Resilience to Hydrologic Perturbations

Development of reliable biogeochemical models requires a mechanistic consideration of microbial interactions with hydrology. Microbial response to and its recovery after hydrologic perturbations (i.e., resilience) is a critical component to understand in this regard, but generally difficult to predict because the impacts of future events can be dependent on the history of perturbations (i.e., historical contingency). Fundamental issues underlying this phenomenon include how microbial resilience to hydrologic perturbations is influenced by historical contingency and how their relationships vary depending on the characteristics of microbial functions. To answer these questions, we considered a simple microbial community composed of two species that redundantly consume a common substrate but specialize in producing distinct products and developed a continuous flow reactor model where the two species grow with trade-offs along the flow rate. Simulations of this model revealed that (1) the history of hydrologic perturbations can lead to the shifts in microbial populations, which consequently affect the community's functional dynamics, and (2) while historical contingency in resilience was consistently predicted for all microbial functions, it was more pronounced for specialized functions, compared to the redundant function. As a signature of historical contingency, our model also predicted the emergence of hysteresis in the transitions across conditions, a critical aspect that can affect transient formation of intermediate compounds in biogeochemistry. This work presents microbial growth traits and their functional redundancy or specialization as fundamental factors that control historical contingencies in resilience.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Scoping review of Culex mosquito life history trait heterogeneity in response to temperature

Abstract Background Mosquitoes in the genus Culex are primary vectors in the US for West Nile virus (WNV) and other arboviruses. Climatic drivers such as temperature have differential effects on species-specific changes in mosquito range, distribution, and abundance, posing challenges for population modeling, disease forecasting, and subsequent public health decisions. Understanding these differences in underlying biological dynamics is crucial in the face of climate change. Methods We collected empirical data on thermal response for immature development rate, egg viability, oviposition, survival to adulthood, and adult lifespan for Culex pipiens, Cx. quinquefasciatus, Cx. tarsalis , and Cx. restuans from existing literature according to the PRISMA scoping review guidelines. Results We observed linear relationships with temperature for development rate and lifespan, and nonlinear relationships for survival and egg viability, with underlying variation between species. Optimal ranges and critical minima and maxima also appeared varied. To illustrate how model output can change with experimental input data from individual Culex species, we applied a modified equation for temperature-dependent mosquito type reproduction number for endemic spread of WNV among mosquitoes and observed different effects. Conclusions Current models often input theoretical parameters estimated from a single vector species; we show the need to implement the real-world heterogeneity in thermal response between species and present a useful data resource for researchers working toward that goal. Graphical Abstract

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Deep learning of experimental electrochemistry for battery cathodes across diverse compositions

Artificial intelligence (AI) has emerged as a tool for discovering and optimizing novel battery materials. However, the adoption of AI in battery cathode representation and discovery is still limited due to the complexity of optimizing multiple performance properties and the scarcity of high-fidelity data. Here, we present a machine learning model (DRXNet) for battery informatics and demonstrate the application in the discovery and optimization of disordered rocksalt (DRX) cathode materials. We have compiled the electrochemistry data of DRX cathodes over the past 5 years, resulting in a dataset of more than 19,000 discharge voltage profiles on diverse chemistries spanning 14 different metal species. Learning from this extensive dataset, our DRXNet model can capture critical features in the cycling curves of DRX cathodes under various conditions. Our approach offers a data-driven solution to facilitate the rapid identification of novel cathode materials, accelerating the development of next-generation batteries for carbon neutralization.

25 ENERGY STORAGE↗
DOE OSTI2023

Comparing the Effect of Anthropogenically Amplified Halogen Natural Emissions on Tropospheric Ozone Chemistry Between Pre-Industrial and Present-Day

Reactive halogens (X + XO, X = I, Br or Cl) catalytically destroy a fraction of tropospheric ozone under present-day (PD) conditions, however, their distribution and potential impact on tropospheric ozone under pre-industrial (PI) conditions remain largely unexplored. This study uses the Community Atmosphere Model with Chemistry (CAM-Chem) to investigate the effect of anthropogenically amplified natural emissions of halogenated species and their subsequent chemistry on tropospheric ozone under PI and PD atmospheric conditions. Model results show that the global tropospheric ozone depletion due to natural halogens is slightly more sensitive in PI than PD, with percentage changes in tropospheric ozone burden (TOB) of -14.1 ± 0.6% for PI and -12.9 ± 0.6% for PD. Individually, the role of iodine and chlorine in ozone depletion is equivalent in both periods (ΔTOB I : ~-7% and ΔTOB Cl : ~-2.5%), while bromine plays a larger role in PI (ΔTOB Br : -5.5 ± 0.6%) versus PD (ΔTOB Br : -4.3 ± 0.7%). The increase in anthropogenic ozone precursor emissions from PI to PD has amplified the natural emission of inorganic halogens and led to a shift in the partitioning of inorganic halogens from reactive to reservoir species. Consequently, halogen-driven ozone depletion from the surface to the free troposphere is larger in PI than PD. In contrast, in the upper troposphere, the ozone depletion is larger in PD influenced mainly by stratospheric intrusion of reactive halogens from long-lived species. This study highlights the importance of including a complete chemical coupling of natural halogens and atmospheric pollutants in chemistry-climate models to adequately assess their effects on tropospheric ozone in a changing climate.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

Carbon chemistry in dense molecular clouds: Theory and observational constraints

For the most part, gas phase models of the chemistry of dense molecular clouds predict the abundances of simple species rather well. However, for larger molecules and even for small systems rich in carbon these models often fail spectacularly. Researchers present a brief review of the basic assumptions and results of large scale modeling of the carbon chemistry in dense molecular clouds. Particular attention is to the influence of the gas phase C/O ratio in molecular clouds, and the likely role grains play in maintaining this ratio as clouds evolve from initially diffuse objects to denser cores with associated stellar and planetary formation. Recent spectral line surveys at centimeter and millimeter wavelengths along with selected observations in the submillimeter have now produced an accurate inventory of the gas phase carbon budget in several different types of molecular clouds, though gaps in our knowledge clearly remain. The constraints these observations place on theoretical models of interstellar chemistry can be used to gain insights into why the models fail, and show also which neglected processes must be included in more complete analyses. Looking toward the future, larger molecules are especially difficult to study both experimentally and theoretically in such dense, cold regions, and some new methods are therefore outlined which may ultimately push the detectability of small carbon chains and rings to much heavier species.

Blake, Geoffrey A.↗
NASA NTRSDate not supplied

The Search for New Molecular Species in the Atmosphere of Titan

The atmosphere of Saturn’s moon, Titan, contains a vast array of minor molecular constituents originating from the dissociation of its primary compounds, N2and CH4. The extent of Titan’s photochemical network was revealed by ionospheric measurements from the Cassini orbiter, leaving substantial gaps in the corresponding neutral inventory typically probed through infrared spectroscopy. Utilizing the unprecedented sensitivity and high spectral resolution of the Atacama Large Millimeter/submillimeter Array (ALMA) at longer wavelengths, we searched for additional trace gases in Titan’s atmosphere during 2016, 2017,and 2019to further characterize Titan’s atmospheric and surface composition. These observations have resulted in the spectroscopic detections of the small cyclic molecule, cyclopropenylidene (c-C3H2), and the largest nitrile to date, methyl cyanoacetylene (CH3C3N), for the first time in Titan’s atmosphere –or anywhere else in the Solar System. We employed radiative transfer models to determine the volume mixing ratio of both species in Titan’s middle and upper atmosphere (altitudes > 350 km)and constrained the total column densities for the incorporation into future photochemical models of Titan’s atmosphere. Discrepancies between these results and current models necessitates the need for dedicated laboratory and modeling studies regarding the production pathways for trace complex organics in Titan’s atmosphere. Additionally, the measurement of previously undetected molecular species provides further incentive for the study of planetary atmospheres in the (sub)millimeter regime with facilities such as ALMA.

A E Thelen↗
DOE OSTI2026

Machine learning enables reconstruction of past fire regimes from charcoal-derived fire intensity and fuel composition

Background Fire is a foundational ecological process that shapes ecosystem structure, diversity, and resilience. Quantifying paleofire regime attributes such as frequency, severity, and intensity is essential for understanding the historical range of variability in fire behavior and its ecological effects. While frequency and severity are often reconstructed in paleofire studies, quantitative reconstructions of fire intensity remain limited. Recent work has shown that maximum pyrolysis temperature—a proxy for fire intensity—and plant species type can be inferred from charcoal using transmission Fourier-transform infrared (FTIR) spectroscopy. However, the sample preparation for transmission FTIR is destructive and time-consuming, limiting application and reuse of materials for other analyses. We evaluated reflectance FTIR spectroscopy as a non-destructive alternative for reconstructing combustion temperature and plant species from laboratory-generated charcoal. We also examined the influence of contrasting airflow environments (ambient air versus nitrogen-rich) on pyrolysis temperature and plant species reconstruction prediction accuracies and compared predictive performance between a novel, neural network–based deep learning model with the traditional modern analogue technique (MAT) using k-nearest neighbor functions. As proof of concept, we apply our enhanced methodology to ancient charcoal to demonstrate applicability at improving long-term fire regime reconstructions and the ability to link paleofire records with contemporary fire ecology. Results Our analysis shows that transmission and reflectance FTIR spectra yield comparable spectral profiles. However, sample preparation for reflectance FTIR is minimal and non-destructive, unlike transmission FTIR which is destructive. We demonstrate that oxygen environments improved reconstruction accuracy relative to nitrogen-rich conditions. Finally, our deep learning neural network (DL) achieved testing accuracies of 98.7% for temperature and 96.2% for species identification, outperforming MAT’s k-NN approach (89.8% and 65.9%, respectively). A Shapley importance analysis identified 5 key spectral regions that greatly influenced the model’s temperature or species categorization. When applied to ancient charcoal, our results show historic fires from the most recent past primarily burned at low intensities (400–500 °C), reflective of natural fire regimes in ponderosa pine forests. Our results corroborate charcoal morphology data that suggests all ancient charcoal originated from burned woody plant types. Conclusions By combining reflectance FTIR spectroscopy with a deep learning approach, we provide the first accuracies high enough to confidently identify both species and temperature from laboratory-produced charcoal, improving quantitative reconstructions of fire intensity and fuel composition from paleofire records. This opens a wide range of research into the link between fire and larger drivers (i.e., climate or human) and greater ecological understanding of fire regimes beyond that of burn scars or recent observations. These methodological improvements have direct relevance for fire management by improving interpretation of historical fire behavior, informing fuel–fire relationships, and providing a scalable analytical framework applicable to both long-term ecological studies and contemporary fire science.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

The Impact of Anisotropic Error Correlation Modelling

Most data assimilation systems assume isotropic forecast error correlations, but results from two dimensional Kalman Filter experiments indicate that the correlations can be far from isotropic. In this paper we use a simple two-dimensional data assimilation system, which analyses trace chemical species such as ozone, to assess different approaches to modelling the error correlations. We compare assimilation results using isotropic correlations with results obtained using different approaches to modelling anisotropic correlations: first, using correlations based on the concentrations of the trace chemicals, and secondly using an advective correlation model. We show that these relatively cheap ways of modelling anisotropic correlations give objectively better results than using isotropic correlations. We discuss the possible extension of these approaches to a full 3-D meteorological data assimilation system.

Swinbank, R.↗
DOE OSTI2021

Modeling the performance and faradaic efficiency of solid oxide electrolysis cells using doped barium zirconate perovskite electrolytes

Y-doped BaZrO 3 (BaZr 1 x Y x O 3 δ , or “BZY”), a proton-conducting ceramic featuring high bulk conductivity and good chemical stability, is a promising electrolyte material for solid oxide electrolysis cells. Further doping with Ce and/or Yb (creating materials “BCZY” and “BCZYYb”) can improve conductivity and sintering properties, but at significant penalty to cells’ faradaic efficiency (FE). Studies have proposed that reduction of lattice Ce can occur in the hydrogen electrode, which consumes some hydrogen produced by the hydrogen evolution reaction, leading to decreased FE. Despite studies suggesting this phenomenon, the mechanism is largely unknown. In this work, we developed a multiphysics model to study the transport of multiple defect species and the performance of BZY, BCZY, and BCZYYb, capturing the tradeoff between enhanced performance at the cost of FE for BCZY and BCZYYb electrolytes compared to BZY. We also found that increasing the water content of the anode gas supply lowers the current output of the cell but results in better FE. The model, which uses several parameters previously unavailable in the literature, was validated to experiments varying temperature, steam water content, and electrolyte material, as well as two performance metrics (performance curves and FE). Results verify and explain observed trends, informing future work on Ce-doped BZY electrolytes.

08 HYDROGEN↗
NASA NTRS1995

Modeling of nonequilibrium space plasma flows

Godunov-type numerical solution of the 20 moment plasma transport equations. One of the centerpieces of our proposal was the development of a higher order Godunov-type numerical scheme to solve the gyration dominated 20 moment transport equations. In the first step we explored some fundamental analytic properties of the 20 moment transport equations for a low b plasma, including the eigenvectors and eigenvalues of propagating disturbances. The eigenvalues correspond to wave speeds, while the eigenvectors characterize the transported physical quantities. In this paper we also explored the physically meaningful parameter range of the normalized heat flow components. In the second step a new Godunov scheme type numerical method was developed to solve the coupled set of 20 moment transport equations for a quasineutral single-ion plasma. The numerical method and the first results were presented at several national and international meetings and a paper describing the method has been published in the Journal of Computational Physics. To our knowledge this is the first numerical method which is capable of producing stable time-dependent solutions to the full 20 (or 16) moment set of transport equations, including the full heat flow equation. Previous attempts resulted in unstable (oscillating) solutions of the heat flow equations. Our group invested over two man-years into the development and implementation of the new method. The present model solves the 20 moment transport equations for an ion species and thermal electrons in 8 domain extending from a collision dominated to a collisionless region (200 km to 12,000 km). This model has been applied to study O+ acceleration due to Joule heating in the lower ionosphere.

Gombosi, Tamas↗
DOE OSTI2022

Copper availability governs nitrous oxide accumulation in wetland soils and stream sediments

Denitrification is microbially-mediated through enzymes containing metal cofactors. Laboratory studies of pure cultures have highlighted that the availability of Cu, required for the multicopper enzyme nitrous oxide reductase, can limit N2O reduction. However, in natural aquatic systems, such as wetlands and hyporheic zones in stream beds, the role of Cu in controlling denitrification remains incompletely understood. In this study, we collected soils and sediments from three natural environments -- riparian wetlands, marsh wetlands, and a stream -- to investigate their nitrogen species transformation activity at background Cu levels and different supplemented Cu loadings. All of the systems contained solid-phase associated Cu below or around geological levels (40 - 280 nmol g-1) and exhibited low dissolved Cu (3-50 nM), which made them appropriate sites for evaluating the effect of limited Cu availability on denitrification. In laboratory incubation experiments, high concentrations of N2O accumulated in all microcosms lacking Cu amendment except for one stream sediment sample. With Cu added to provide dissolved concentrations at trace levels (10-300 nM), reduction rate of N2O to N 2 in the wetland soils and stream sediments was enhanced. A kinetic model could account for the trends in nitrogen species by combining the reactions for microbial reduction of NO 3 - to NO 2 -/N 2 O/N 2 and abiotic reduction of NO 2 - to N 2 . The model revealed that the rate of N 2 O to N 2 conversion increased significantly in the presence of Cu. For riparian wetland soils and stream sediments, the kinetic model also suggested that overall denitrification is driven by abiotic reduction of NO2- in the presence of inorganic electron donors. This study demonstrated that natural aquatic systems containing Cu at concentrations less than or equal to crustal abundances may display incomplete reduction of N 2 O to N 2 that would cause N 2 O accumulation and release to the atmosphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Study of Martian Polar Caps with GISS ROCKE-3D GCM

Martian polar caps consist of both H 2 O and CO 2 ice. While H 2 O ice is mainly passive on modern Mars, it may have not been the case in recent Martian history, when its obliquity was higher, or when it was changing rapidly. The distribution of ice species in the snowpack affects its physical and thermodynamic properties. In the upper layers, it determines its albedo and thermal emissivity. Thus understanding the mutual effect between these ices and their interaction with the atmosphere is crucial for understanding the evolution of Martian polar regions. In this study, we employ a newly-developed Exotic Ices snow model coupled to the NASA Goddard Institute for Space Studies (GISS) ROCKE-3D planetary General Circulation Model (GCM) [1] to study the behavior of Martian polar caps. ROCKE-3D is a planetary GCM developed at NASA GISS as an extension of its Earth climate model, modelE [2]. It has been extensively used to simulate climate of various planets, including Mars (e.g. [3,4]). The Exotic Ices snow model was specially developed for planetary applications which involve more than one condensable in the atmosphere, in which case snow can contain multiple species of ice (CO2 and H 2 O in the Mars case). For each species of ice, the model uses their proper physical properties and phase diagram, but otherwise it treats all species of ice on an equal footing. The combined effects on albedo, thermal inertia and mutual insulation are treated accordingly. The snowpack interacts with the atmospheric dust cycle, and can accumulate a prognostic amount of dust, though the effect of dust on snow properties is not currently treated explicitly, and is prescribed. In this study, we first validate our model against the modern Martial climate, for which we use mission results from Mars Climate Sounder (atmospheric temperature and dust optical depth), SPICAM on Mars Express (atmospheric water), and Viking 2 (surface pressure). We investigate the effect of snow radiative properties on CO 2 and water cycles and the ability of our model to accurately reproduce those with minimal model tuning. We then perform simulations for several obliquities from a recent Martian past, and investigate the behavior of the Martian polar caps in such conditions.

Mars↗