Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial constraints”

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 217 records · Page 12

Aerosol, Cloud and Trace Gas Observations Derived from Airborne Hyperspectral Radiance and Direct Beam Measurements in Recent Field Campaigns

The AERONET (AErosol RObotic NETwork) ground-based suite of sunphotometers provides measurements of spectral aerosol optical depth (AOD), precipitable water and spectral sky radiance, which can be inverted to retrieve aerosol microphysical properties that are critical to assessments of aerosol-climate interactions. Because of data quality criteria and sampling constraints, there are significant limitations to the temporal and spatial coverage of AERONET data and their representativeness for global aerosol conditions. The 4STAR (Spectrometer for Sky-Scanning, Sun-Tracking Atmospheric Research) instrument, jointly developed by NASA Ames and PNNL with NASA Goddard collaboration, combines airborne sun tracking and AERONET-like sky scanning with spectroscopic detection. Being an airborne instrument, 4STAR has the potential to fill gaps in the AERONET data set. Dunagan et al. [2013] present results establishing the performance of the instrument, along with calibration, engineering flight test, and preliminary scientific field data. The 4STAR instrument operated successfully in the SEAC4RS [Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys] experiment in Aug./Sep. 2013 aboard the NASA DC-8 and in the DoE [Department of Energy]-sponsored TCAP [Two Column Aerosol Project, July 2012 & Feb. 2013] experiment aboard the DoE G-1 aircraft (Shinozuka et al., 2013), and acquired a wealth of data in support of mission objectives on all SEAC4RS and TCAP research flights. 4STAR provided direct beam measurements of hyperspectral AOD, columnar trace gas retrievals (H2O, O3, NO2; Segal-Rosenheimer et al., 2014), and the first ever airborne hyperspectral sky radiance scans, which can be inverted to yield the same products as AERONET ground-based observations. In addition, 4STAR measured zenith radiances underneath cloud decks for retrievals of cloud optical depth and effective diameter. In this presentation, we provide an overview of the new 4STAR capabilities for airborne field campaigns, with an emphasis on comparisons between 4STAR and AERONET sky radiances, and retrievals of aerosol microphysical properties based on sky radiance measurements, column trace gas amounts from spectral direct beam measurements and cloud property retrievals from zenith mode observations for a few select case studies in the SEAC4RS and TCAP experiments. We summarize the aerosol, trace gas, cloud and airmass characterization studies made possible by the combined 4STAR direct beam, and sky/zenith radiance observations.

Aerosol-Climate↗

Union: A Unified HW-SW Co-Design Ecosystem in MLIR for Evaluating Tensor Operationson Spatial Accelerators

To meet the extreme compute demands for deep learning across commercial and scientific applications, dataflow accelerators are becoming increasingly popular. While these“domain-specific” accelerators are not fully programmable like CPUs and GPUs, they retain varying levels of flexibility with respect to data orchestration, i.e., dataflow and tiling optimizations to enhance efficiency. There are several challenges when designing new algorithms and mapping approaches to execute the algorithms for a target problem on new hardware. Previous works have addressed these challenges individually. To address this challenge as a whole, in this work, we present an HW-SW co-design ecosystem for spatial accelerators called Union within the popular MLIR compiler infrastructure. Our framework allows exploring different algorithms and their mappings on several accelerator cost models. Union also includes a plug-and-play library of accelerator cost models and mappers which can easily be extended. The algorithms and accelerator cost models are connected via a novel mapping abstraction that captures the map space of spatial accelerators which can be systematically pruned based on constraints from the hardware, workload, and mapper. We demonstrate the value of Union for the community with several case studies which examine offloading different tensor operations (CONV/GEMM/Tensor Contraction) on diverse accelerator architectures using different mapping schemes.

Jeong, Geonhwa↗

Soil Moisture and Sea Surface Salinity Derived from Satellite-Borne Sensors

The monitoring of soil moisture and sea surface salinity over the Earth has been profoundly enhanced during the last thirteen years owing to new generation of satellite sensors. L-band radiometry is currently the only technology providing direct measurements of soil moisture, insensitive to surface roughness and distribution of elements in the soil, and the only technology that allows to measure sea surface salinity from space. The SMOS (Soil Moisture and Ocean Salinity) and SMAP (Soil Moisture Active Passive) satellite missions resolve global and local variability with a spatial resolution of approximately 43 km, a swath width close to 1000 km, and a sampling time, for each mission, of at least twice every 3 days. These resolutions and samplings can be increased by either merging data from the two sensors, and with complementary informations gathered with other passive or active sensors, or with in situ information at higher spatial resolution. Numerous scientific studies based on the use of this new type of measurement have led to a better understanding and constraint of the processes governing the variability of the water cycle, ocean circulation and the Earth's climate. The continuity of measurements, and the increased spatial and radiometric resolution is critical for fulfilling scientific needs. Future L-band radiometry missions currently being planned in Europe (the Copernicus Imaging Microwave Radiometer), and in China (the Ocean Salinity Mission) should provide better constraint on auxiliary parameters by combining multiple frequencies, but they will not have improved spatial resolution beyond SMOS and SMAP. The temporal continuity with SMOS and SMAP will most likely not be ensured. In parallel, new concepts are being developed to increase spatial resolution of both land and ocean parameters.

J. Boutin↗

Innovative laser-based methods for monitoring fuel retention in ITER

This paper addresses the challenge of tritium inventory management in ITER and future fusion reactors, highlighting the importance of accurate tritium measurement and its spatial distribution within the vacuum vessel. Given ITER’s operational constraints, especially the limit on tritium retention, precise measurement is essential for both safety and regulatory compliance. To tackle these questions, the paper presents the T-monitor diagnostic system developed by Forschungszentrum Jülich, which uses Laser-Induced Desorption (LID) in combination with Diagnostic Residual Gas Analysis (DRGA) to measure hydrogen isotope concentrations on the surface of divertor tiles. The system integrates a high-power laser, advanced optical components, and a Fast Scanning Mirror Unit (FSMU) for accurate laser spot positioning with rapid response. Designed to measure in situ tritium retention, the diagnostic provides high-resolution spatial mapping, vital for evaluating detritiation strategies. The laser heating process increases the divertor surface temperature to 1600 K within the laser spot, promoting hydrogen isotope desorption. Accurate measurements require the precise control of laser parameters, including pulse duration and spot size, with a target relative accuracy of 20%. The optical design includes both in-vessel and ex-vessel components, such as durable high-reflectivity mirrors made of gold and copper, selected not only for their infrared performance but also for their transmission of visible wavelengths for observation purposes. To protect optical components from contamination, a pneumatic shutter is used.

Hydrogen isotopes↗

Map-level baryonification: unified treatment of weak lensing two-point and higher-order statistics

Precision cosmology benefits from extracting maximal information from cosmic structures, motivating the use of higher-order statistics (HOS) at small spatial scales. However, predicting how baryonic processes modify matter statistics at these scales has been challenging. The baryonic correction model (BCM) addresses this by modifying dark-matter-only simulations to mimic baryonic effects, providing a flexible, simulation-based framework for predicting both two-point and HOS. We show that a 3-parameter version of the BCM can jointly fit weak lensing maps' two-point statistics, wavelet phase harmonics coefficients, scattering coefficients, and the third and fourth moments to within 2% accuracy across all scales ℓ < 2000 and tomographic bins for a DES-Y3-like redshift distribution ( z ≲ 2), using the FLAMINGO simulations. These results demonstrate the viability of BCM-assisted, simulation-based weak lensing inference of two-point and HOS, paving the way for robust cosmological constraints that fully exploit non-Gaussian information on small spatial scales.

79 ASTRONOMY AND ASTROPHYSICS↗

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗

Precipitable Water Variability Using SSM/I and GOES VAS Pathfinder Data Sets

Determining moisture variability for all weather scenes is critical to understanding the earth's hydrologic cycle and global climate changes. Remote sensing from geostationary satellites provides the necessary temporal and spatial resolutions necessary for global change studies. Due to antenna size constraints imposed with the use of microwave radiometers, geostationary satellites have carried instruments passively measuring radiation at infrared wavelengths or shorter. The shortfall of using infrared instruments in moisture studies lies in its inability to sense terrestrial radiation through clouds. Microwave emissions, on the other hand, are mostly unaffected by cloudy atmospheres. Land surface emissivity at microwave frequencies exhibit both high temporal and spatial variability thus confining moisture retrievals at microwave frequencies to over marine atmospheres (a near uniform cold background). This study intercompares the total column integrated water content Precipitable Water, (PW) as derived from both the Special Sensor Microwave Imager (SSM/I) and the Geostationary Operational Environmental Satellite (GOES) VISSR Atmospheric Sounder (VAS) pathfinder data sets. PW is a bulk parameter often used to quantify moisture variability and is important to understanding the earth's hydrologic cycle and climate system. This research has been spawned in an effort to combine two different algorithms which together can lead to a more comprehensive quantification of global water vapor. The approach taken here is to intercompare two independent PW retrieval algorithms and to validate the resultant retrievals against an existing data set, namely the European Center for Medium range Weather Forecasts (ECMWF) model analysis data.

Lerner, Jeffrey A.↗

The Barberplaid illusion: Even Perceived Plaid Motion is Biased by Asymmetric Apertures

Our results show that the perceived direction of motion of plaids windowed by asymmetric spatial Gaussians is biased toward the long axis of the window. The bias increases as the relative angle between the plaid motion and the window increases, peaks at a relative angle of about 40 degrees, and then decreases. The peak bias was 14 degrees for a spatial frequency of 0.6 cpd and a window aspect ratio of 4.0. The biases increase as the window is elongated and decrease as the component spatial frequency increases. We tested the predictions of several models of human motion processing (cross correlation, motion energy, intersection of constraints, and vector sum), and show that none of these can predict our data. These results suggest that spatial integration of motion signals plays a crucial role in the perception of plaid motion.

Beutter, Brent↗

Optimizing GRACE/GRACE-FO Data and A priori Hydrological Knowledge for Improved Global Terrestial Water Storage Component Estimates

The comprehensive information of global terrestrial water storage (TWS) components (soil moisture, groundwater, snow, surface water) is essential for effective assessment of water resource availability, climate variation, and disaster mitigation measures. Observational data provided by the Gravity Recovery And Climate Experiment (GRACE) and GRACE Follow-On satellite missions offer global TWS variation (ΔTWS) in terms of an integrated water column. However, GRACE spatial resolution is relatively coarse (i.e., 3°), and the vertically integrated value cannot be separated into ΔTWS components directly. This study demonstrates the feasibility to estimate ΔTWS components at any desired spatial-vertical resolution by effectively maintaining the native resolution of the employed hydrological knowledge. It utilizes a least-squares with constraints (LSC) approach to rigorously incorporate GRACE and GRACE-FO data and a priori hydrological knowledge, with the aim to improve global ΔTWS components’ accuracy and spatial resolution. The 3°×3° GRACE mascon derived ΔTWS data is disaggregated into the 0.5°×0.5° anomalous soil moisture storage (ΔSMS), groundwater storage (ΔGWS), snow water equivalent (ΔSWE), and surface water storage (ΔSWS) based on the covariance information obtained from the Community Atmosphere Biosphere Land Exchange (CABLE) and the PCRaster Global Water Balance (PCR-GLOBWB) models. Evaluation with different ground measurements and satellite products between 2002 and 2019 exhibits significantly improved accuracy in all individual ΔTWS components. This improvement is of particular note in ΔGWS and ΔSWS, where the LSC approach increases the globally averaged correlation values by approximately 0.13 and 0.05, respectively. Reliable prior knowledge leads to a more accurate ΔTWS component estimate, and the use of ensemble-mean knowledge yields the best result.

GRACE↗

Detection of the SO2 atmosphere on Io with the Hubble Space Telescope

Observations of the trailing hemisphere of Io made with the Faint Object Spectrograph of the Hubble Space Telescope (HST) in March 1992 have resulted in the first detection of atmospheric SO2 absorption bands in the ultraviolet. These observations represent only the third positive means of detection of what is widely believed to be Io's primary atmospheric constituent. Below approximately 2130 A the geometric albedo of the satellite is dominated by SO2 gas absorption band signatures, which have been analyzed using models that include the effects of optical thickness, temperature, and spatial distribution. The disk-intergrated HST data cannot resolve the spatial distribution, but it is possible to define basic properties and set constraints on the atmosphere at the time of the observations. Hemispheric atmospheres with average column density N = 6 - 10 x 10(exp 15)/sq sm and T(gas) = 110 - 500 K fit the data, with preference for temperatures of approximately 200 - 250 K. Better fits are found as the atmosphere is spatially confined, with a limit of approximately 8% hemispheric areal coverage and N approximately equal to 3 x 10(exp 17)/sq cm with colder 110 - 250 K temepratures. A dense (N greater than or equal to 10(exp 16)/sq cm), localized component of SO2 gas, such as that possibly associated with active volcanoes, can generate the observed spectral constrast only when the atmosphere is cold (110 K) and an extended component such as Pele is included. The combination of a dense, localized atmosphere with a tenuous component (N less than 10(exp 16)/sq cm, either patchy or extended) also fits the data. In all cases the best fit models imply a disk-averaged column density larger than exospheric but approximately 10 - 30 times less than the previous upper limit from near-UV observations.

Ballester, G. E.↗

Sensitivity of CONUS Summer Rainfall to the Selection of Cumulus Parameterization Schemes in NU-WRF Seasonal Simulations

This study investigates the sensitivity of daily rainfall rates in regional seasonal simulations over the contiguous United States (CONUS) to different cumulus parameterization schemes. Daily rainfall fields were simulated at 24-km resolution using the NASA-Unified Weather Research and Forecasting (NU-WRF) Model for June-August 2000. Four cumulus parameterization schemes and two options for shallow cumulus components in a specific scheme were tested. The spread in the domain-mean rainfall rates across the parameterization schemes was generally consistent between the entire CONUS and most subregions. The selection of the shallow cumulus component in a specific scheme had more impact than that of the four cumulus parameterization schemes. Regional variability in the performance of each scheme was assessed by calculating optimally weighted ensembles that minimize full root-mean-square errors against reference datasets. The spatial pattern of the seasonally averaged rainfall was insensitive to the selection of cumulus parameterization over mountainous regions because of the topographical pattern constraint, so that the simulation errors were mostly attributed to the overall bias there. In contrast, the spatial patterns over the Great Plains regions as well as the temporal variation over most parts of the CONUS were relatively sensitive to cumulus parameterization selection. Overall, adopting a single simulation result was preferable to generating a better ensemble for the seasonally averaged daily rainfall simulation, as long as their overall biases had the same positive or negative sign. However, an ensemble of multiple simulation results was more effective in reducing errors in the case of also considering temporal variation.

Iguchi, Takamichi↗

A temperature profile diagnostic for radiation waves on OMEGA-60

Predicting and matching radiation wave propagation with computational models has proven difficult. Information provided by experiments studying radiation flow has been limited when only radiation breakout is measured. We have developed the COAX (co-axial) diagnostic platform to provide spatial temperature profiles of a radiation wave through low density foams as a more detailed constraint for simulations. COAX uses a standard, laser-driven OMEGA-60 halfraum to drive radiation down a titanium-laden silicon oxide foam. Point-projection X-ray absorption spectroscopy perpendicular to the radiation flow measures the spatial profile of titanium ionization. The spectroscopic measurement utilizes a broadband capsule backlighter. Imaging and streak spectroscopy are used to characterize the size and spectrum of this source. Radiography provides an additional constraint by capturing the developing shock as the radiation flow becomes subsonic. The DANTE diagnostic is used to measure the halfraum temperature. Here, we provide a spectroscopic analysis of COAX data to determine temperature, and we describe experimental sources of uncertainty. The temperature is obtained by comparison to multi-temperature synthetic spectra post-processed from radiation-hydrodynamics simulations. Quantitative comparison between data and synthetic spectra generated from temperature profiles at relevant simulation times enable determination of a peak temperature of 114 ± 8 eV at 265 ± 22.4 μm from the halfraum. This represents an improvement over the temperature uncertainties of previous radiation flow experiments. Further refinements to the spectroscopic analysis could achieve ±4 eV. The combination between space-resolved spectroscopy and radiography enables us to determine the distance from the halfraum of both the radiation front and the shock front at the time of measurement. For the example shown in this paper the radiation front position is 600–630 μm at 3.43 ± 0.16 ns and the shock front position is 633 μm at 3.3 ± 0.24 ns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electron heat flux and propagating fronts in plasma thermal quench via ambipolar transport

The thermal collapse of a nearly collisionless plasma interacting with a cooling spot, in which the electron parallel heat flux plays an essential role, is both theoretically and numerically investigated. Here, we show that such thermal collapse, which is known as thermal quench in tokamaks, comes about in the form of propagating fronts, originating from the cooling spot, along magnetic field lines. The slow fronts, propagating with local ion sound speed, limit the aggressive cooling of plasma, which is accompanied by a plasma cooling flow toward the cooling spot. The extraordinary physics underlying such a cooling flow is that the fundamental constraint of ambipolar transport along the field line limits the spatial gradient of electron thermal conduction flux to the much weaker convective scaling, as opposed to the free-streaming scaling, so that a large electron temperature and, hence, pressure gradient can be sustained. The last ion front for a radiative cooling spot is a shock front where cold but flowing ions meet hot ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Significant DBSCAN+: Statistically Robust Density-based Clustering

Cluster detection is important and widely used in a variety of applications, including public health, public safety, transportation, and so on. Given a collection of data points, we aim to detect density-connected spatial clusters with varying geometric shapes and densities, under the constraint that the clusters are statistically significant. The problem is challenging, because many societal applications and domain science studies have low tolerance for spurious results, and clusters may have arbitrary shapes and varying densities. As a classical topic in data mining and learning, a myriad of techniques have been developed to detect clusters with both varying shapes and densities (e.g., density-based, hierarchical, spectral, or deep clustering methods). However, the vast majority of these techniques do not consider statistical rigor and are susceptible to detecting spurious clusters formed as a result of natural randomness. On the other hand, scan statistic approaches explicitly control the rate of spurious results, but they typically assume a single “hotspot” of over-density and many rely on further assumptions such as a tessellated input space. To unite the strengths of both lines of work, we propose a statistically robust formulation of a multi-scale DBSCAN, namely Significant DBSCAN+, to identify significant clusters that are density connected. As we will show, incorporation of statistical rigor is a powerful mechanism that allows the new Significant DBSCAN+ to outperform state-of-the-art clustering techniques in various scenarios. We also propose computational enhancements to speed-up the proposed approach. Experiment results show that Significant DBSCAN+ can simultaneously improve the success rate of true cluster detection (e.g., 10–20% increases in absolute F1 scores) and substantially reduce the rate of spurious results (e.g., from thousands/hundreds of spurious detections to none or just a few across 100 datasets), and the acceleration methods can improve the efficiency for both clustered and non-clustered data.

Computer Science↗

Probing Physics beyond the Standard Model through Combined Analyses of Next-generation Type Ia Supernova, Cosmic Microwave Background, and Baryon Acoustic Oscillation Surveys

Observations of Type Ia supernovae (SNe Ia), which probe the late Universe, together with baryon acoustic oscillations (BAO) and the cosmic microwave background (CMB), which probe the intermediate and early epochs, provide complementary constraints on the expansion history of the Universe. In this work, we forecast constraints on dark energy and other extensions to the standard cosmological model by combining the SN Ia sample expected from the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST), data from current and forthcoming CMB surveys, and BAO measurements from the Dark Energy Spectroscopic Instrument (DESI). For the CMB, we use temperature, polarization, and lensing power spectra (TT/EE/TE/ϕϕ) from the South Pole Telescope, the planned Advanced Simons Observatory, and a CMB-S4–like experiment. We derive constraints on ΛCDM and its extensions involving the dark energy equation-of-state parameters (w 0 , w a ) and the sum of neutrino masses ∑m ν using a Markov Chain Monte Carlo (MCMC) sampling framework. We find that the LSST Year 3 SN Ia sample can improve upon the DES Year 5 dark energy constraints by a factor of 2−2.5×, with the gains driven primarily by the significantly higher SN Ia density in the LSST sample. Similarly, DESI-DR3 shows up to a 1.8× improvement on dark energy parameters over DR2, driven largely by the substantial increase in the low-redshift sample. Combining CMB with LSST-Y3-SN Ia and DESI-DR3-BAO yields σ(w 0 ) = 0.028 and σ(w a ) = 0.11 for w 0 w a CDM cosmology with the results being largely independent of the CMB dataset. The constraints weaken by 10%–30% when freeing ∑m ν and spatial curvature. Moreover, the joint analysis of the three datasets can enable a 2σ–3σ detection of ∑m ν .

Raghunathan, Srinivasan [University of California;↗

On the stability of the boundary layer on a transonic swept wing

Both incompressible and compressible linear stability theory are applied to the three-dimensional compressible boundary layer on a particular transonic sweptback wing of infinite span. A spatial stability theory is used which identifies the growth direction with the real part of the complex angle of the group velocity. It is found that in the forward, but not the rear, crossflow instability region, the maximum amplification rates of the steady disturbances may be calculated to within about 10% by the incompressible stability theory. There is little difference between the sixth and eighth-order compressible theories. The maximum amplification rate of the steady disturbances at any chordwise station is closely related to the maximum crossflow at that station independent of the Reynolds number. For other than crossflow instability, there can be large differences between the incompressible and compressible theories, both as to the amplification rate and the angle of the wavenumber vector for maximum instability. Amplitude ratios of individual wave components are obtained by integrating the spatial amplification rate along the growth direction subject to the constraint that the wavenumber vector is irrotational. This procedure yields steady disturbances aligned with the local potential flow direction whose wavelengths are nearly independent of downstream distance.

Mack, L. M.↗

The autocorrelation properties of fluctuations in the cosmic X-ray background

The HEAO 1 A-2 all-sky survey data have been analyzed to investigate the autocorrelation function of the surface brightness fluctuations of the extragalactic 2-10 keV X-ray background on angular scales ranging from 3 to 27 deg. The derived upper limits are compatible with optical estimates of the spatial correlation function of Abell's clusters and set interesting constraints on possible AGN clustering.

Persic, M.↗

Temporal changes in the geographic distribution, elevation, and potential origin of the Martian outflow channels

Observational evidence of outflow channel activity on Mars suggests that water was abundant in the planet's early crust. However, with the decline in the planet's internal heat flow, a freezing front developed within the regolith that propagated downward with time and acted as a thermodynamic sink for crustal H2O. One result of this thermal evolution is that, if the initial inventory of water on Mars was small, the cryosphere may have grown to the point where all the available water was taken up as ground ice. Alternatively, if the inventory of H2O exceeds the current pore volume of the cryosphere, then Mars has always possessed extensive bodies of subpermafrost groundwater. We have investigated the relative age, geographic distribution, elevation, and geologic setting of the outflow channels in an effort to accomplish the following: (1) identify possible modes of origin and evolutionary trends in their formation; (2) gain evidence regarding the duration and spatial distribution of groundwater in the crust; and (3) better constraint estimates of the planetary inventory of H2O.

Tribe, S.↗