Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Response surface model”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Radiative Forcing by Well-Mixed Greenhouse Gases: Estimates from Climate Models in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4)

The radiative effects from increased concentrations of well-mixed greenhouse gases (WMGHGs) represent the most significant and best understood anthropogenic forcing of the climate system. The most comprehensive tools for simulating past and future climates influenced by WMGHGs are fully coupled atmosphere-ocean general circulation models (AOGCMs). Because of the importance of WMGHGs as forcing agents it is essential that AOGCMs compute the radiative forcing by these gases as accurately as possible. We present the results of a radiative transfer model intercomparison between the forcings computed by the radiative parameterizations of AOGCMs and by benchmark line-by-line (LBL) codes. The comparison is focused on forcing by CO2, CH4, N2O, CFC-11, CFC-12, and the increased H2O expected in warmer climates. The models included in the intercomparison include several LBL codes and most of the global models submitted to the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4). In general, the LBL models are in excellent agreement with each other. However, in many cases, there are substantial discrepancies among the AOGCMs and between the AOGCMs and LBL codes. In some cases this is because the AOGCMs neglect particular absorbers, in particular the near-infrared effects of CH4 and N2O, while in others it is due to the methods for modeling the radiative processes. The biases in the AOGCM forcings are generally largest at the surface level. We quantify these differences and discuss the implications for interpreting variations in forcing and response across the multimodel ensemble of AOGCM simulations assembled for the IPCC AR4.

Collins, W. D.↗
NASA NTRS2011

Stardust-NExT, Deep Impact, and the Accelerating Spin of 9P/Tempel 1

The evolution of the spin rate of Comet 9P/Tempel 1 through two perihelion passages (in 2000 and 2005) is determined from 1922 Earth-based observations taken over a period of 13 year as part of a World-Wide observing campaign and from 2888 observations taken over a period of 50 days from the Deep Impact spacecraft. We determine the following sidereal spin rates (periods): 209.023 +/- 0.025deg/dy (41.335 0.005 h) prior to the 2000 perihelion passage, 210.448 +/- 0.016deg/dy (41.055 +/- 0.003 h) for the interval between the 2000 and 2005 perihelion passages, 211.856 +/- 0.030deg/dy (40.783 +/- 0.006 h) from Deep Impact photometry just prior to the 2005 perihelion passage, and 211.625 +/- 0.012deg/dy (40.827 +/- 0.002 h) in the interval 2006-2010 following the 2005 perihelion passage. The period decreased by 16.8 +/- 0.3 min during the 2000 passage and by 13.7 +/- 0.2 min during the 2005 passage suggesting a secular decrease in the net torque. The change in spin rate is asymmetric with respect to perihelion with the maximum net torque being applied on approach to perihelion. The Deep Impact data alone show that the spin rate was increasing at a rate of 0.024 +/- 0.003deg/dy/dy at JD2453530.60510 (i.e., 25.134 dy before impact), which provides independent confirmation of the change seen in the Earth-based observations. The rotational phase of the nucleus at times before and after each perihelion and at the Deep Impact encounter is estimated based on the Thomas et al. (Thomas et al. [2007]. Icarus 187, 4-15) pole and longitude system. The possibility of a 180deg error in the rotational phase is assessed and found to be significant. Analytical and physical modeling of the behavior of the spin rate through of each perihelion is presented and used as a basis to predict the rotational state of the nucleus at the time of the nominal (i.e., prior to February 2010) Stardust-NExT encounter on 2011 February 14 at 20:42. We find that a net torque in the range of 0.3-2.5 x 10(exp 7) kg/sq m/sq s acts on the nucleus during perihelion passage. The spin rate initially slows down on approach to perihelion and then passes through a minimum. It then accelerates rapidly as it passes through perihelion eventually reaching a maximum post-perihelion. It then decreases to a stable value as the nucleus moves away from the Sun. We find that the pole direction is unlikely to precess by more than approx. 1deg per perihelion passage. The trend of the period with time and the fact that the modeled peak torque occurs before perihelion are in agreement with published accounts of trends in water production rate and suggests that widespread H2O out-gassing from the surface is largely responsible for the observed spin-up.

Belton, Michael J. S.↗
DOE OSTI2024

Fabrication of single-crystalline YFeO3 films with large antiferromagnetic domains

The antiferromagnetic orthoferrite YFeO3 possesses fascinating magnetic properties for spintronics, such as terahertz spin dynamics, ultrafast domain wall motion, and long magnon decay length. YFeO3 belongs to a special family of antiferromagnets that show an unusually strong non-trivial Kerr response due to its weak ferromagnetism. The highly stable antiferromagnetic domains without any spontaneous spin rotation transitions below the 645 K Néel temperature may be useful for nanoscale device applications. We report the successful fabrication of high-quality twinning-free (110)-oriented YFeO3 films by pulsed laser deposition. Detailed structural and magnetic characterization revealed that the crystal structure and magnetic properties of the YFeO3 films are comparable to bulk single crystals. We show that the spin rotation under high magnetic fields follows the two-sublattice approximation model. The film surface is atomically flat with step-terrace surface morphology. A longitudinal magneto-optic Kerr (MOKE) rotation of 10 mdeg was observed at room temperature, which is consistent with earlier reports on bulk single crystals. The in-plane anisotropy of the Kerr response corresponds to the obtained magnetic anisotropy from the SQUID measurement. The large MOKE signal enables the imaging of antiferromagnetic domains and their reversal. The domain size was found to be larger than 100 μm. These high-quality YFeO3 thin films facilitate the fabrication of antiferromagnetic spintronic devices and provide a convenient platform for studying various spin-related phenomena in thin films and at interfaces.

Physics↗
DOE OSTI2026

Data and Scripts Associated with "Modeling Ecohydrological Responses of Vegetation to Urban Microclimates Using the E3SM Land Model"

This dataset supports the study of vegetation ecohydrological responses to urban microclimates using the land component of the Energy Exascale Earth System Model (ELM) at four urban sites in Knoxville, Tennessee, USA. It includes the model inputs, simulation outputs, and associated scripts for running ELM simulations and analyzing the resulting data. The Model_Inputs folder includes static surface data, satellite-derived phenology (i.e., leaf area index), and atmospheric forcing data used to drive ELM simulations. Detailed descriptions of these datasets are provided in Section 2.3.2 of the associated manuscript. The Model_Outputs folder contains simulation results for the baseline, treatment, and ensemble experiments. Outputs from the baseline and treatment simulations are provided as raw ELM NetCDF files. Because the raw outputs from the 4,000-member ensemble are prohibitively large, the ensemble results are provided as summarized CSV files, which also serve as the source data for Figure 5 of the associated manuscript. The Scripts folder contains three components: E3SM, the core codebase of the Energy Exascale Earth System Model (E3SM); elm-olmt, the Offline Land Model Testbed (OLMT) used to perform the simulations; and knoxville_elm, which contains the analysis scripts used to process model outputs and generate the figures and results presented in the associated manuscript. Additional information is provided in Scripts_readme.txt within the Scripts directory.

Lu, Xiaoman [ORNL] (ORCID:0000000306698780)↗
DOE OSTI2020

Effective radiative forcing and adjustments in CMIP6 models

Abstract. The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W m−2, comprised of 1.81 (±0.09) W m−2 from CO2, 1.08 (± 0.21) W m−2 from other well-mixed greenhouse gases, −1.01 (± 0.23) W m−2 from aerosols and −0.09 (±0.13) W m−2 from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W m−2. The majority of the remaining 0.21 W m−2 is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W m−2, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

The volcanic signal in Goddard Institute for Space Studies three-dimensional model simulations

Transient calculations of the Goddard Institute for Space Studies general circulation model for the climatic signal of volcanic eruptions are analyzed. By compositing the output for two different volcanoes for scenario A and five different volcanoes for scenario B, the natural variability is suppressed and the volcanic signals are extracted. Significant global mean surface air temperature cooling and precipitation reduction are found for several years following the eruptions, with larger changes in the Northern Hemisphere (NH) than in the Southern Hemisphere. The global-average temperature response lasts for more than four years, but the precipitation response disappears after three years. The largest cooling in the model occurs in the NH summer of the year after spring eruptions. Significant zonal-average temperature reductions begin in the tropics immediately after the eruptions and extend to 45 deg S - 45 deg N in the year after the eruptions. In the second year, cooling is still seen from 30 deg S to 30 deg N. Because of the low variability in this model as compared to the real world, these signals may appear more significant here than they would by attempting to isolate them using real data. The results suggest that volcanoes can enhance the drought in the Sahel. No evidence was found that stratospheric aerosols from the low-atitude volcanic eruptions can trigger El Nino-Southern Oscillation (ENSO) events in this model.

Robock, Alan↗
DOE OSTI2023

Energy coupling and stoichiometry of Zn 2+ /H + antiport by the prokaryotic cation diffusion facilitator YiiP

YiiP from Shewanella oneidensis is a prokaryotic Zn 2+ /H + antiporter that serves as a model for the Cation Diffusion Facilitator (CDF) superfamily, members of which are generally responsible for homeostasis of transition metal ions. Previous studies of YiiP as well as related CDF transporters have established a homodimeric architecture and the presence of three distinct Zn 2+ binding sites named A, B, and C. In this study, we use cryo-EM, microscale thermophoresis and molecular dynamics simulations to address the structural and functional roles of individual sites as well as the interplay between Zn 2+ binding and protonation. Structural studies indicate that site C in the cytoplasmic domain is primarily responsible for stabilizing the dimer and that site B at the cytoplasmic membrane surface controls the structural transition from an inward facing conformation to an occluded conformation. Binding data show that intramembrane site A, which is directly responsible for transport, has a dramatic pH dependence consistent with coupling to the proton motive force. A comprehensive thermodynamic model encompassing Zn 2+ binding and protonation states of individual residues indicates a transport stoichiometry of 1 Zn 2+ to 2–3 H + depending on the external pH. This stoichiometry would be favorable in a physiological context, allowing the cell to use the proton gradient as well as the membrane potential to drive the export of Zn 2+ .

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2008

Dominance of ENSO-Like Variability in Controlling Tropical Ocean Surface Energy Fluxes in the Satellite Era

Ocean surface turbulent and radiative fluxes are critical links in the climate system since they mediate energy exchange between the two fluid systems (ocean and atmosphere) whose combined heat transport determines the basic character of Earth's climate. Moreover, interannual to decadal climate variability depends crucially on the nature of these exchange processes. For example, addressing the question of the degree to which the global hydrologic cycle is changing depends on our ability to observe and model these fluxes accurately. In this work we investigate the interannual to decadal variation of fluxes over the global tropics, especially the tropical oceans. Recent versions of satellite-derived fresh water flux estimates as well as some reanalyses (e.g. products from Remote Sensing Systems, the Woods Hole Oceanographic Institute, and Global Precipitation Climatology Project) suggest that increases in evaporation and precipitation over the past 20 years exceed those expected on the basis of climate model projected responses to greenhouse gas forcing. At the same time, it is well known that E1 Nino / Southern Oscillation behavior in the Pacific exhibits significant variability at scales longer than interannual. We examine here the degree to which surface fluxes attending these interannual to decadal fluctuations are related to ENSO. We examine consistency between these data sets and explore relationships between SST variations, flux changes and modulation of tropical Walker and Hadley circulations.

Robertson, F. R.↗
DOE OSTI2023

Helfrich-Hurault elastic instabilities driven by geometrical frustration

Here, the Helfrich-Hurault (HH) elastic instability is a well-known mechanism behind patterns that form as a result of strain upon liquid crystal systems with periodic ground states. In the HH model, layered structures undulate and buckle in response to local, geometric incompatibilities in order to maintain the preferred layer spacing. Classic HH systems include cholesteric liquid crystals under electromagnetic field distortions and smectic liquid crystals under mechanical strains, where both materials are confined between rigid substrates. However, richer phenomena are observed when undulation instabilities occur in the presence of deformable interfaces and variable boundary conditions. Understanding how the HH instability is affected by deformable surfaces is imperative for applying the instability to a broader range of materials. In this review, the HH mechanism is reexamined and special focus is given to how the boundary conditions influence the response of lamellar systems to geometrical frustration. Lamellar liquid crystals confined within a spherical shell geometry are used as the model system. Made possible by the relatively recent advances in microfluidics within the past 15 years, liquid crystal shells are composed entirely of fluid interfaces and have boundary conditions that can be dynamically controlled at will. Past and recent work that exemplifies how topological constraints, molecular anchoring conditions, and boundary curvature can trigger the HH mechanism in liquid crystals with periodic ground states is examined. The review ends by identifying similar phenomena across a wide variety of materials, both biological and synthetic. The fact that the HH mechanism is a generic and often overlooked response of periodic materials to geometrical frustration is highlighted.

42 ENGINEERING↗
NASA NTRS1990

Venus geology and tectonics - Hotspot and crustal spreading models and questions for the Magellan mission

Spacecraft and ground-based observations of Venus have revealed a geologically young and active surface - with volcanoes, rift zones, orogenic belts and evidence for hotspots and crustal spreading - yet the processes responsible for these features cannot be identified from the available data. The Magellan spacecraft will acquire an unprecedented global data set which will provide a comprehensive and well resolved view of the planet. This will permit global geological mapping, an assessment of the style and relative importance of geological processes, and will help in the understanding of links between the surface geology and mantle dynamics of this earth-like planet.

Head, James W.↗
NASA NTRS1991

Phenomenological and experimental study of the thermal response of low density silica ablators to high enthalpy plasma flows

A complete analysis of the melting-ablator governing equations is carried out for porous ablators with the purpose of establishing an approximate analytical model capable of determining surface recession, temperature, and heat flux as well as in-depth temperature profiles. A detailed dimensional analysis of the melt-flow governing equations is performed for parameters applicable to high-enthalpy flows. Results of arc-jet flow experiments performed on NASA Reusable Surface Insulation materials of different densities and over a range of stagnation pressures and heat fluxes are discussed. A dimensionless correlation for melt run-offs is derived in terms of material and flow parameters; however, it is pointed out that the accuracy of the current data is sufficient to conclude the compliance with the model.

Henline, William D.↗
NASA NTRS2000

Evaluation of Bogus Vortex Techniques with Four-Dimensional Variational Data Assimilation

The effectiveness of techniques for creating "bogus" vortices in numerical simulations of hurricanes is examined by using the Penn State/NCAR nonhydrostatic mesoscale model (MM5) and its adjoint system. A series of four-dimensional variational data assimilation (4-D VAR) experiments is conducted to generate an initial vortex for Hurricane Georges (1998) in the Atlantic Ocean by assimilating bogus sea-level pressure and surface wind information into the mesoscale numerical model. Several different strategies are tested for improving the vortex representation. The initial vortices produced by the 4-D VAR technique are able to reproduce many of the structural features of mature hurricanes. The vortices also result in significant improvements to the hurricane forecasts in terms of both intensity and track. In particular, with assimilation of only bogus sea-level pressure information, the response in the wind field is contained largely within the divergent component, with strong convergence leading to strong upward motion near the center. Although the intensity of the initial vortex seems to be well represented, a dramatic spin down of the storm occurs within the first 6 h of the forecast. With assimilation of bogus surface wind data only, an expected dominance of the rotational component of the wind field is generated, but the minimum pressure is adjusted inadequately compared to the actual hurricane minimum pressure. Only when both the bogus surface pressure and wind information are assimilated together does the model produce a vortex that represents the actual intensity of the hurricane and results in significant improvements to forecasts of both hurricane intensity and track.

Pu, Zhao-Xia↗
DOE OSTI2024

Crystal plasticity modeling and analysis for the transition from intergranular to transgranular failure in nickel-based alloy Inconel 740H at elevated temperature

The precipitation-strengthened Nickel alloy Inconel® 740H® (IN740H) exhibits increased ductility at higher applied strain rates during quasi-static tensile tests at an elevated temperature of 760°C. The examination of fracture surfaces in this context reveals a noteworthy transition of underlying fracture mechanisms from transgranular to intergranular fracture as the applied strain rate decreases from 1×10 -3 / s to 0.83×10 -4 / s . To thoroughly understand the mechanical response of IN740H under these conditions, this study develops a crystal plasticity finite element (CPFE) model. Further, this model incorporates various deformation mechanisms including dislocation slips, climb, and grain boundary sliding, which are relevant to the test conditions. The model is calibrated using data from both tensile tests at different strain rates and creep tests across a broad stress range at 760°C, enabling the accurate determination of model parameters for each mechanism. Simulation results well captured the experimental observations of different failure modes. At higher strain rates, the model shows a dominance of dislocation slip leading to heterogeneous plastic deformation and formation of transgranular shear bands causing the failure, while at lower strain rates, an increased activity of grain boundary sliding causes grain boundaries crack leading to intergranular failure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Characterization of Advanced Manufactured Surface Acoustic Wave Sensor for Temperature Measurement in Nuclear Applications

As part of a collaborative effort between Idaho National Laboratory and the Boise State University, printed SAW transducers were developed. The transducers were first modelled and evaluated for optimum performance, then fabricated using several piezoelectric substrates. These printed SAW-sensors were tested for temperature sensing ability and found to be highly sensitive and with a very linear response over the tested range.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2007

Seasonal Correlations of SST, Water Vapor, and Convective Activity in Tropical Oceans: A New Hyperspectral Data Set for Climate Model Testing

The analysis of the response of the Earth Climate System to the seasonal changes of solar forcing in the tropical oceans using four years of the Atmospheric Infrared Sounder (AIRS) and Advanced Microwave Sounding Unit (AMSU) data between 2002 and 2006 gives new insight into amplitude and phase relationships between surface and tropospheric temperatures, humidity, and convective activity. The intensity of the convective activity is measured by counting deep convective clouds. The peaks of convective activity, temperature in the mid-troposphere, and water vapor in the 0 - 30 N and 0 - 30 S tropical ocean zonal means occur about two months after solstice, all leading the peak of the sea surface temperature by several weeks. Phase is key to the evaluation of feedback. The evaluation of climate models in terms of zonal and annual means and annual mean deviations from zonal means can now be supplemented by evaluating the phase of key atmospheric and surface parameters relative to solstice. The ability of climate models to reproduce the statistical flavor of the observed amplitudes and relative phases for broad zonal means should lead to increased confidence in the realism of their water vapor and cloud feedback algorithms. AIRS and AMSU were launched into a 705 km altitude polar sun-synchronous orbit on the EOS Aqua spacecraft on May 4, 2002, and have been in routine data gathering mode since September 2002.

Atmospheric Infrared Sounder↗
NASA NTRS2007

Seasonal Correlations of SST, Water Vapor, and Convective Activity in Tropical Oceans: A New Hyperspectral Data Set for Climate Model Testing

The analysis of the response of the Earth Climate System to the seasonal changes of solar forcing in the tropical oceans using four years of the Atmospheric Infrared Sounder (AIRS) and Advanced Microwave Sounding Unit (AMSU) data between 2002 and 2006 gives new insight into amplitude and phase relationships between surface and tropospheric temperatures, humidity, and convective activity. The intensity of the convective activity is measured by counting deep convective clouds. The peaks of convective activity, temperature in the mid-troposphere, and water vapor in the 0-30 N and 0-30 S tropical ocean zonal means occur about two months after solstice, all leading the peak of the sea surface temperature by several weeks. Phase is key to the evaluation of feedback. The evaluation of climate models in terms of zonal and annual means and annual mean deviations from zonal means can now be supplemented by evaluating the phase of key atmospheric and surface parameters relative to solstice. The ability of climate models to reproduce the statistical flavor of the observed amplitudes and relative phases for broad zonal means should lead to increased confidence in the realism of their water vapor and cloud feedback algorithms. AIRS and AMSU were launched into a 705 km altitude polar sun-synchronous orbit on the EOS Aqua spacecraft on May 4, 2002, and have been in routine data gathering mode since September 2002.

Atmospheric Infrared Sounder (AIRS)↗
NASA NTRS1995

Use of midlatitude soil moisture and meteorological observations to validate soil moisture simulations with biosphere and bucket models

Soil moisture observations in sites with natural vegetation were made for several decades in the former Soviet Union at hundreds of stations. In this paper, the authors use data from six of these stations from different climatic regimes, along with ancillary meteorological and actinometric data, to demonstrate a method to validate soil moisture simulations with biosphere and bucket models. Some early and current general circulation models (GCMs) use bucket models for soil hydrology calculations. More recently, the Simple Biosphere Model (SiB) was developed to incorporate the effects of vegetation on fluxes of moisture, momentum, and energy at the earth's surface into soil hydrology models. Until now, the bucket and SiB have been verified by comparison with actual soil moisture data only on a limited basis. In this study, a Simplified SiB (SSiB) soil hydrology model and a 15-cm bucket model are forced by observed meteorological and actinometric data every 3 h for 6-yr simulations at the six stations. The model calculations of soil moisture are compared to observations of soil moisture, literally 'ground truth,' snow cover, surface albedo, and net radiation, and with each other. For three of the stations, the SSiB and 15-cm bucket models produce good simulations of seasonal cycles and interannual variations of soil moisture. For the other three stations, there are large errors in the simulations by both models. Inconsistencies in specification of field capacity may be partly responsible. There is no evidence that the SSiB simulations are superior in simulating soil moisture variations. In fact, the models are quite similar since SSiB implicitly has a bucket embedded in it. One of the main differences between the models is in the treatment of runoff due to melting snow in the spring -- SSiB incorrectly puts all the snowmelt into runoff. While producing similar soil moisture simulations, the models produce very different surface latent and sensible heat fluxes, which would have large effects on GCM simulations.

Robock, Alan↗
NASA NTRS1993

Simultaneous laboratory measurements of CO2 and H2O adsorption on palagonite: Implications for the Martian climate and volatile reservoir

We are measuring the simultaneous adsorption of H2O and CO2 on palagonite materials in order to improve the formulation of climate models for Mars. We report on the initial co-adsorption data. Models of the Martian climate and volatile inventory indicate that the regolith serves as one of the primary reservoirs of outgassed volatiles and that it exchanges H2O and CO2 with the atmosphere in response to changes in insolation associated with astronomical cycles. Physical adsorbate must exist on the surfaces of the cold particulates that constitute the regolith, and the size of that reservoir can be assessed through laboratory measurements of adsorption on terrestrial analogs. Many studies of the independent adsorption of H2O and CO2 on Mars analog were made and appear in the literature. Empirical expressions that relate the adsorptive coverage of each gas to the temperature of the soil and partial pressure have been derived based on the laboratory data. Numerical models incorporate these adsorption isotherms into climatic models, which predict how the adsorptive coverage of the regolith and hence, the pressure of each gas in the atmosphere will vary as the planet moves through its orbit. These models suggest that the regolith holds several tens to hundreds of millibars of CO2 and that during periods of high obliquity warming of the high-latitude regolith will result in desorption of the CO2, and a consequent increase in atmospheric pressure. At lower obliquities, the caps cool and the equator warms forcing the desorption of several tens of millibars of CO2, which is trapped into quasipermanent CO2 caps.

Zent, A. P.↗