Near-surface perennial water ice on Mars has been previously inferred down to latitudes of about 45{\deg} and could result from either water vapor diffusion through the regolith under current conditions or previous ice ages precipitations. In this paper we show that at latitudes as low as 25{\deg} in the southern hemisphere buried water ice in the shallow (< 1 m) subsurface is required to explain the observed surface distribution of seasonal CO2 frost on pole facing slopes. This result shows that possible remnants of the last ice age, as well as water that will be needed for the future exploration of Mars, are accessible significantly closer to the equator than previously thought, where mild conditions for both robotic and human exploration lie.

In this paper, we analyze water ice occurrences at the surface of Mars using near-infrared observations, and we study their distribution with a climate model. Latitudes between 45{\deg}S and 50{\deg}N are considered. Data from the Observatoire pour la Min\'eralogie, l'Eau, les Glaces et l'Actitit\'e and the Compact Reconnaissance Imaging Spectrometer for Mars are used to assess the presence of surface water ice as a function of location and season. A modeling approach combining the 1-D and 3-D versions of the General Circulation Model of the Laboratoire de M\'et\'eorologie Dynamique de Jussieu is developed and successfully compared to observations. Ice deposits 2-200 \mu m thick are observed during the day on pole facing slopes in local fall, winter and early spring. Ice extends down to 13{\deg} latitude in the Southern Hemisphere but is restricted to latitudes higher than 32{\deg} in the north. On a given slope, the pattern of ice observations at the surface is mainly controlled by the global variability of atmospheric water (precipitation and vapor), with local ground properties playing a lower role. Only seasonal surface ice is observed: no exposed patches of perennial ground ice have been detected. Surface seasonal ice is however sensitive to subsurface properties: the results presented in this study are consistent with the recent discovery of low latitude subsurface ice obtained through the analysis of CO2 frost.

We discuss the DISORT-based radiative transfer pipeline ('CRISM_LambertAlb') for atmospheric and thermal correction of MRO/CRISM data acquired in multispectral mapping mode (~200 m/pixel, 72 spectral channels). Currently, in this phase-one version of the system, we use aerosol optical depths, surface temperatures, and lower-atmospheric temperatures, all from climatology derived from Mars Global Surveyor Thermal Emission Spectrometer (MGS-TES) data, and surface altimetry derived from MGS Mars Orbiter Laser Altimeter (MOLA). The DISORT-based model takes as input the dust and ice aerosol optical depths (scaled to the CRISM wavelength range), the surface pressures (computed from MOLA altimetry, MGS-TES lower-atmospheric thermometry, and Viking-based pressure climatology), the surface temperatures, the reconstructed instrumental photometric angles, and the measured I/F spectrum, and then outputs a Lambertian albedo spectrum. The Lambertian albedo spectrum is valuable geologically since it allows the mineralogical composition to be estimated. Here, I/F is defined as the ratio of the radiance measured by CRISM to the solar irradiance at Mars divided by $\pi$. After discussing the capabilities and limitations of the pipeline software system, we demonstrate its application on several multispectral data cubes: the outer northern ice cap of Mars, Tyrrhena Terra, and near the landing site for the Phoenix mission. For the icy spectra near the northern polar cap, aerosols need to be included in order to properly correct for the CO_2 absorption in the H_{2}O ice bands at wavelengths near 2.0 $\mu$m. In future phases of software development, we intend to use CRISM data directly in order to retrieve the spatiotemporal maps of aerosol optical depths, surface pressure and surface temperature.

The observations of Mars by the CRISM and OMEGA hyperspectral imaging spectrometers require correction for photometric, atmospheric and thermal effects prior to the interpretation of possible mineralogical features in the spectra. Here, we report on a simple, yet non-trivial, adaptation to the commonly-used volcano-scan correction technique for atmospheric CO_2, which allows for the improved detection of minerals with intrinsic absorption bands at wavelengths between 1.9-2.1 $\mu$m. This volcano-scan technique removes the absorption bands of CO_2 by ensuring that the Lambert albedo is the same at two wavelengths: 1.890 $\mu$m and 2.011 $\mu$m, with the first wavelength outside the CO_2 gas bands and the second wavelength deep inside the CO_2 gas bands. Our adaptation to the volcano-scan technique moves the first wavelength from 1.890 $\mu$m to be instead within the gas bands at 1.980 $\mu$m, and for CRISM data, our adaptation shifts the second wavelength slightly, to 2.007 $\mu$m. We also report on our efforts to account for a slight ~0.001 $\mu$m shift in wavelengths due to thermal effects in the CRISM instrument.