A study of the properties of a local dust storm with Mars Express OMEGA and PFS data

We present observations of a local dust storm performed by the OMEGA and PFS instruments aboard Mars Express. OMEGA observations are used to retrieve the dust single-scattering albedo in the spectral range 0.4-4.0 μm. The single-scattering albedo shows fairly constant values between 0.6 and 2.6 μm, and a sharp decrease at wavelengths shorter than 0.6 μm, in agreement with previous studies. It presents a small absorption feature due to ferric oxide at 0.9 μm, and a strong absorption feature due to hydrated minerals between 2.7 and 3.6 μm. We use a statistical method, the Independent Component Analysis, to determine that the dust spectral signature is decoupled from the surface albedo, proving that the retrieval of the single-scattering albedo is reliable, and we map the dust optical thickness with a conventional radiative transfer model. The effect of the dust storm on the atmospheric thermal structure is measured using PFS observations. We also simulate the thermal impact of the dust storm using a one-dimensional atmospheric model. A comparison of the retrieved and modeled temperature structures suggests that the dust in the storm should be confined to the 1-2 lowest scale heights of the atmosphere. However, the observed OMEGA reflectance in the CO₂ absorption bands does not support this suggestion.     

Constraints on the composition of the martian south polar cap from gravity and topography

The polar caps of Mars have long been acknowledged to be composed of unknown proportions of water ice, solid CO₂ (dry ice), and dust. Gravity and topography data are here analyzed over the southern cap to place constraints on its density, and hence composition. Using a localized spectral analysis combined with a lithospheric flexure model of ice cap loading, the best fit density of the volatile-rich south polar layered deposits is found to be 1271 kg m⁻³ with 1-σ limits of 1166 and 1391 kg m⁻³. The best fit elastic thickness of this geologically young deposit is 140 km, though any value greater than 102 km can fit the observations. The best fit density implies that about 55% dry ice by volume could be sequestered in these deposits if they were completely dust free. Alternatively, if these deposits were completely free of solid CO₂, the dust content would be constrained to lie between about 14 and 28% by volume. The bulk thermal conductivity of the polar cap is not significantly affected by these maximum allowable concentrations of dust. However, even if a moderate quantity of solid CO₂ were present as horizontal layers, the bulk thermal conductivity of the polar cap would be significantly reduced. Reasonable estimates of the present day heat flow of Mars predict that dry ice beneath the thicker portions of the south polar cap would have melted. Depending on the quantity of solid CO₂ in these deposits today, it is even possible that water ice could melt where the cap is thickest. If independent estimates for either the dust or CO₂ content of the south polar cap could be obtained, and if radar sounding data could determine whether this polar cap is presently experiencing basal melting or not, it would be possible to use these observations to place tight constraints on the present day heat flow of Mars.     

Analysis of condensates formed at the viking 2 lander site: The first winter

A thin light-colored ground covering appeared on the surface of Mars near the Viking 2 lander from L_s = 230° to L_s = 16°, a total of 249 Mars days, during the lander's first winter on the surface. This paper presents a reduction of applicable lander imagery during the period. Imaging sequences, relative surface albedo, spectral reflectance estimates, and limited photometric data are presented and compared with previous laboratory measurements. Photometric data are best fit by an average Minnaert k = 1.1 (blue), k = 1.0 (green), and k = 0.95 (red). Appearance and disappearance rates, spectral reflectance, and photometric data all tend to confirm an earlier proposal that the covering was a combination of H₂O and CO₂, which fell already condensed onto dust particles brought northward by the season's first major dust storm. Under this assumption, the covering thickness is estimated to be between 0.5 and a few millimeters.     

Dust aerosols above the south polar cap of Mars as seen by OMEGA

The time evolution of atmospheric dust at high southern latitudes on Mars has been determined using observations of the south seasonal cap acquired in the near infrared (1-2.65 μm) by OMEGA/Mars Express in 2005. Observations at different solar zenith angles and one EPF sequence demonstrate that the reflectance in the 2.64 μm saturated absorption band of the surface CO₂ ice is mainly due to the light scattered by aerosols above most places of the seasonal cap. We have mapped the total optical depth of dust aerosols in the near-IR above the south seasonal cap of Mars from mid-spring to early summer with a time resolution ranging from one day to one week and a spatial resolution of a few kilometers. The optical depth above the south perennial cap is determined on a longer time range covering southern spring and summer. A constant set of optical properties of dust aerosols is consistent with OMEGA observations during the analyzed period. Strong variations of the optical depth are observed over small horizontal and temporal scales, corresponding in part to moving dust clouds. The late summer peak in dust opacity observed by Opportunity in 2005 propagated to the south pole contrarily to that observed in mid spring. This may be linked to evidence for dust scavenging by water ice-rich clouds circulating at high southern latitudes at this season.     

Dust deposition at the Mars Pathfinder landing site: observations and modeling of visible/near-infrared spectra

Temporal variations in the visible/near-infrared reflectance spectra of the radiometric calibration targets on the Mars Pathfinder (MPF) lander obtained by the Imager for Mars Pathfinder (IMP) camera reveal the effects of aeolian dust deposition at the MPF site throughout the mission. Sky brightness models in combination with two-layer radiative transfer models were used with these data to track changes in dust opacity on the radiometric calibration targets (RCTs) to constrain the dust deposition rate and the spectral properties of the deposited dust. Two-layer models were run assuming both linear and nonlinear dust accumulation rates, and suggest that RCT dust optical depth at the end of the 83-sol mission was 0.08 to 0.16, or on the order of 5- to 10-μm thickness for plausible values for dust porosity and grain size. These values correspond to dust fall rates of about 20-45 μm per Earth year, consistent with previous studies of dust deposition on Mars. The single scattering albedos of the dust derived from the models fall between those previously determined for atmospheric dust and bright soils. Comparisons of relative reflectance spectra calibrated using observed RCT radiances from late in the mission versus using radiances from modeled (dust-free) RCTs also reveal distinct spectral differences consistent with dust on the RCTs. Temporal variations in RCT dust opacity are not clearly linked to known passages of vortices at the MPF site, but overall suggest that deposition of dust onto the targets by local dust devils may be favored over erosion. Analyses of temporal changes in visible/near-infrared spectra will provide valuable information for future missions by constraining how dust deposition affects landed spacecraft operability (e.g., solar power availability), instrument calibration, and interpretations of surface mineralogy and composition.     

Constraints on martian lobate debris apron evolution and rheology from numerical modeling of ice flow

Radar observations in the Deuteronilus Mensae region by Mars Reconnaissance Orbiter have constrained the thickness and dust concentration found within mid-latitude ice deposits, providing an opportunity to more accurately estimate the rheology of ice responsible for the formation of lobate debris aprons based on their apparent age of ∼100 Myr. We developed a numerical model simulating ice flow under martian conditions using results from ice deformation experiments, theory of ice grain growth based on terrestrial ice cores, and observational constraints from radar profiles and laser altimetry. By varying the ice grain size, the ice temperature, the subsurface slope, and the initial ice volume we determine the combination of parameters that best reproduce the observed LDA lengths and thicknesses over a period of time comparable to the apparent ages of LDA surfaces (90-300 Myr). We find that an ice temperature of 205 K, an ice grain size of 5 mm, and a flat subsurface slope give reasonable ages for many LDAs in the northern mid-latitudes of Mars. Assuming that the ice grain size is limited by the grain boundary pinning effect of incorporated dust, these results limit the dust volume concentration to less than 4%. However, assuming all LDAs were emplaced by a single event, we find that there is no single combination of grain size, temperature, and subsurface slope which can give realistic ages for all LDAs, suggesting that some or all of these variables are spatially heterogeneous. Based on our model we conclude that the majority of northern mid-latitude LDAs are composed of clean (?4 vol%), coarse (?1 mm) grained ice, but regional differences in either the amount of dust mixed in with the ice, or in the presence of a basal slope below the LDA ice must be invoked. Alternatively, the ice temperature and/or timing of ice deposition may vary significantly between different mid-latitude regions. Either eventuality can be tested with future observations.     

High spectral resolution UV to near-IR observations of Mars using HST/STIS

We acquired high spectral and spatial resolution hyperspectral imaging spectrometer observations of Mars from near-UV to near-IR wavelengths (∼300 to 1020 nm) using the STIS instrument on the Hubble Space Telescope during the 1999, 2001, and 2003 oppositions. The data sets have been calibrated to radiance factor (I/F) and map-projected for comparison to each other and to other Mars remote sensing measurements. We searched for and (where detected) mapped a variety of iron-bearing mineral signatures within the data. The strong and smooth increase in I/F from the near-UV to the visible that gives Mars its distinctive reddish color indicates that poorly crystalline ferric oxides dominate the spectral properties of the high albedo regions (as well as many intermediate and low albedo regions), a result consistent with previous remote sensing studies of Mars at these wavelengths. In the near-IR, low albedo regions with a negative spectral slope and/or a distinctive ∼900 nm absorption feature are consistent with, but not unique indicators of, the presence of high-Ca pyroxene or possibly olivine. Mixed ferric-ferrous minerals could also be responsible for the ∼900 nm feature, especially in higher albedo regions with a stronger visible spectral slope. We searched for the presence of several known diagnostic absorption features from the hydrated ferric sulfate mineral jarosite, but did not find any unique evidence for its occurrence at the spatial scale of our observations. We identified a UV contrast reversal in some dark region spectra: at wavelengths shorter than about 340 nm these regions are actually brighter than classical bright regions. This contrast reversal may be indicative of extremely “clean” low albedo surfaces having very little ferric dust contamination. Ratios between the same regions observed during the planet-encircling dust storm of 2001 and during much clearer atmospheric conditions in 2003 provide a good direct estimate of the UV to visible spectral characteristics of airborne dust aerosols. These HST observations can help support the calibration of current and future Mars orbital UV to near-IR spectrometers, and they also provide a dramatic demonstration that even at the highest spatial resolution possible to achieve from the Earth, spectral variations on Mars at these wavelengths are subtle at best.     

A spherical Monte-Carlo model of aerosols: Validation and first applications to Mars and Titan

The atmospheres of Mars and Titan are loaded with aerosols that impact remote sensing observations of their surface. Here we present the algorithm and the first applications of a radiative transfer model in spherical geometry designed for planetary data analysis. We first describe a fast Monte-Carlo code that takes advantage of symmetries and geometric redundancies. We then apply this model to observations of the surface of Mars and Titan at the terminator as acquired by OMEGA/Mars Express and VIMS/Cassini. These observations are used to probe the vertical distribution of aerosols down to the surface. On Mars, we find the scale height of dust particles to vary between 6 km and 12 km depending on season. Temporal variations in the vertical size distribution of aerosols are also highlighted. On Titan, an aerosols scale height of 80 ± 10 km is inferred, and the total optical depth is found to decrease with wavelength as a power-law with an exponent of −2.0 ± 0.4 from a value of 2.3 ± 0.5 at 1.08 μm. Once the aerosols properties have been constrained, the model is used to retrieve surface reflectance properties at high solar zenith angles and just after sunset.     

Mariner 9 photometric observations of Mars from November 1971 through March 1972

The large quantity of Mariner 9 television pictures taken at phase angles ranging up to twice those accessible from Earth has been used to describe intergrated photometric properties of Mars. Although frame-to-frame variations emphasize the shortcomings of vidicon cameras when used in a photometric mode, statistical trends have yielded data comparable with Earth-based observations.Analysis of atmospheric parameters over a period of changing opacity, November 1971 through March 1972, has provided time-varying optical depths and light-scattering information. Linear function characterization of the Mars reflectance is clearly discrepant at large incidence or emission angles and at phase angles greater than 40°.     

Photometric imaging of Mars during the 1973 opposition

Images of Mars were obtained at 20 wavelengths between 0.35 and 1.10 μm during the 1973 opposition using the MIT MKIII silicon vidicon imaging system at the 88-in. telescope of the Mauna Kea Observatory, Hawaii. The observations covered a 3-hr interval on Day 2 of the 1973 Solis Lacus dust storm. At blue wavelengths two huge crescent-shaped anomalous brightenings surrounded the central dust storm in the early morning but faded during the morning; they were probably deposits of H₂O frost or low level haze that was degassed from the regolith in Solis Lacus. The north polar hood, which was nonsymmetrically oriented about the north pole, faded and receded during the observing period. Regions of known dust storm activity could be clearly distinguished, and the images revealed dust cloud activity in Chryse as well. Relative reflectance spectra generated from the images showed compositional variations within the Erythraeum M. dark area complex: Margaritifer S. and Eastern Erythraeum M. have similar ferrosilicate mineralogy, but have stronger ferrosilicate absorption bands than Pyrrhae R. and Eos. Western Erythraeum M. and northern Eos were contaminated by high albedo dust. Reflectance spectra obtained from orbit around Mars could provide higher spatial resolution and more complete coverage of the dark areas, yielding an opportunity to correlate geologic features with crustal petrology. An orbiting spectrometer could also provide spectra in the compositionally diagnostic wavelength range 1–4 μm, where the Earth's atmosphere contains intense absorption bands.     

Radiance factor calibration of near-infrared spectral images of Mars

We present radiometrically calibrated spectrophotometric images of Mars taken at the NASA Infrared Telescope Facility (IRTF) near and during the 1995 and 1999 oppositions. Absolute intensity and radiance factor (r_F = I/F) values have been calculated for approximately 95% of the surface over all longitudes between −70° to 90° latitude in the 1.5- to 4.1-μm spectral region at a spectral resolution (Δλ/λ) of 1.5%. Values of radiance factor range from r_F = 0.4 to 0.6 at 2.2 μm for the bright regions such as Moab and Arabia to r_F = 0.12 to 0.3 at 2.2 μm around the dark regions Syrtis Major and Acidalia Planitia. Variations are seen due to seasonal dust and/or condensate cloud cover and viewing geometry. Our results are generally consistent with the few reported previous radiance factor determinations for Mars. These data are unique among ground-based data in their relatively high spatial resolution (?200 km/pixel at the sub-Earth point) and coverage combined with their spectral resolution and coverage. These radiometrically calibrated observations can be used as input to studies focusing on spectral unmixing of surface and atmospheric components, radiative transfer modeling of disk and limb radiances, and photometric modeling of the martian phase function.     

Effect of Mars analogue dust deposition on the automated detection of calcite in visible/near-infrared spectra

We have developed an artificial neural net detector for use on board Mars rovers that correctly identifies calcite under Mars analogue dust (JSC Mars-1 regolith simulant) layers up to ∼100 μm thickness and 80% aerial coverage. Both the detector output and the band depth of the ∼2300 nm CO⁼₃ absorption are linearly related to the surface area of exposed calcite. This detector provides a means for rapid and robust automated recognition of calcite on Mars in areas of active aeolian erosion.     

Statistical distribution of atmospheric dust devils

Theoretical predictions with regard to dust devil (apparent) angular size–frequency distribution are made and critically compared with Mars Exploration Rover (MER) Spirit optical observations. For an idealized horizontal viewing geometry one should expect that the number of dust devils having the apparent angular diameter greater than a given angle α is inversely proportional to α squared. The actual dependency for Spirit dust devils is in between the inverse-squared and simple inverse laws, and close to the latter one for small and moderate angles α. It is emphasized that such a comparison can be considered as a benchmark for completeness and adequateness of dust-devil optical observations and correctness of competing analytical formulations for dust devil size–frequency distribution.     

Modelling of atmospheric dust extinction and surface reflectance of Mars applying a radiative transfer simulation in the 2.0 and 2.7 μm CO2 bands

Infrared radiation spectra of Mars which can be measured by an orbiting Planetary Fourier Spectrometer (PFS) have been simulated in the spectral region from 1 to 50 μm. The radiative transfer simulation technique considers absorption, emission and multiple scattering by molecular (CO₂, H₂O, CO) and particulate (palagonite) species. It is shown that the contribution from atmospheric dust extinction and surface reflectance can be separated in the region of the CO₂ bands at 2.0 and 2.7 μm. Quantitative results of simultaneous surface pressure, reflectance and aerosol optical depth retrievals are discussed.     

Climate, weather, and north polar observations from the Mars Reconnaissance Orbiter Mars Color Imager

The Mars Reconnaissance Orbiter observes Mars from a nearly circular, polar orbit. From this vantage point, the Mars Color Imager extends the ∼5 Mars years record of Mars Global Surveyor global, visible-wavelength multi-color observations of meteorological events and adds measurements at three additional visible and two ultraviolet wavelengths. Observations of the global distribution of ozone (which anti-correlates with water vapor) and water ice and dust clouds allow tracking of atmospheric circulation. Regional and local observations emphasize smaller scale atmospheric dynamics, especially those related to dust lifting and subsequent motion. Polar observations detail variations related to the polar heat budget, including changes in polar frosts and ices, and storms generated at high thermal contrast boundaries.     

Scattering by nonspherical particles and optical properties of Martian dust

Ultraviolet spectral reflectance data obtained from the Mariner 9 during the 1971 Mars dust storm were analyzed at two wavelengths to deduce the complex index of refraction of Martian dust. Martian dust particles were considered to be irregularly shaped and randomly oriented nonspherical particles. The analysis shows that the Mariner 9 data are consistent with the scattering by irregular particles with an index of refraction m = 1.59 ? 0.0015i at λ = 0.305 μm and m = 1.55 ? 0.004i at λ = 0.268 μm. Both the real and the imaginary parts of refractive index found are considerably smaller than previously reported results using a spherical shape model for Martian dust particles.     

THEMIS observations of Mars aerosol optical depth from 2002-2008

We use infrared images obtained by the Thermal Emission Imaging System (THEMIS) instrument on-board Mars Odyssey to retrieve the optical depth of dust and water ice aerosols over more than 3.5 martian years between February 2002 (MY 25, Ls=330°) and December 2008 (MY 29, Ls=183°). These data provide an important bridge between earlier TES observations and recent observations from Mars Express and Mars Reconnaissance Orbiter. An improvement to our earlier retrieval [Smith, M.D., Bandfield, J.L., Christensen, P.R., Richardson, M.I., 2003. J. Geophys. Res. 108, doi:10.1029/2003JE002114] to include atmospheric temperature information from THEMIS Band 10 observations leads to much improved retrievals during the largest dust storms. The new retrievals show moderate dust storm activity during Mars Years 26 and 27, although details of the strength and timing of dust storms is different from year to year. A planet-encircling dust storm event was observed during Mars Year 28 near Southern Hemisphere Summer solstice. A belt of low-latitude water ice clouds was observed during the aphelion season during each year, Mars Years 26 through 29. The optical depth of water ice clouds is somewhat higher in the THEMIS retrievals at ∼5:00 PM local time than in the TES retrievals at ∼2:00 PM, suggestive of possible local time variation of clouds.     

Investigations of the variability of dust particle sizes in the martian atmosphere using the NASA Ames General Circulation Model

We present a Mars General Circulation Model (GCM) numerical investigation of the physical processes (i.e., wind stress and dust devil dust lifting and atmospheric transport) responsible for temporal and spatial variability of suspended dust particle sizes. Measurements of spatial and temporal variations in airborne dust particles sizes in the martian atmosphere have been derived from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) spectral and emission phase function data [Wolff, M.J., Clancy, R.T., 2003. J. Geophys. Res. (Planets) 108 (E9), doi:10.1029/2003JE002057. 1-1; Clancy, R.T., Wolff, M.J., Christensen, P.R., 2003. J. Geophys. Res. (Planets) 108 (E9), doi:10.1029/2003JE002058. 2-1]. The range of dust particle sizes simulated by the NASA Ames GCM is qualitatively consistent with TES-derived observations of effective dust particle size variability. Model results suggest that the wind stress dust lifting scheme (which produces regionally confined dust lifting) is the process responsible for the majority of the dust particle size variability in the martian atmosphere. Additionally, model results suggest that atmospheric transport processes play an important role in the evolution of atmospheric dust particles sizes during substantial dust storms on Mars. Finally, we show that including the radiative effects of a spatially variable particle size distribution significantly influences thermal and dynamical fields during the dissipation phase of the simulated global dust storm.     

Time variability of Martian bolometric albedo

A time-sequential set of bolometric albedo maps for Mars has been constructed from Viking Infrared Thermal Mapper data. The maps provide global coverage in longitude for latitudes -60° to +60° at 1° by 1° spatial resolution. Individual maps are constructed under strict geometric constraints for a narrow range of L_s. The set of albedo maps spans a Martian year and includes maps before, during, and after the global dust storms of 1977. Transient brightenings associated with local dust storms or condensate clouds are apparent in some of the maps. During dust-free periods, bolometric albedo maps are generally similar to classical, visual albedo map of Mars. The distribution of bolometric surface albedos is bimodal with typical, clear-sky, Lambert albedos of 0.27 and 0.16 for bright and dark areas, respectively. Atmospheric effects strongly influence apparent surface albedos, especially for dark areas. Neither bright nor dark regions show measurable, long-term variations of bolometric albedos during clear periods.     

Yearly and seasonal variations of low albedo surfaces on Mars in the OMEGA/MEx dataset: Constraints on aerosols properties and dust deposits

The time variations of spectral properties of dark martian surface features are investigated using the OMEGA near-IR dataset. The analyzed period covers two Mars years, spanning from early 2004 to early 2008 (includes the 2007 global dust event). Radiative transfer modeling indicates that the apparent albedo variations of low to mid-latitude dark regions are consistent with those produced by the varying optical depth of atmospheric dust as measured simultaneously from the ground by the Mars Exploration Rovers. We observe only a few significant albedo changes that can be attributed to surface phenomena. They are small-scaled and located at the boundaries between bright and dark regions. We then investigate the variations of the mean particle size of aerosols using the evolution of the observed dark region spectra between 1 and 2.5 μm. Overall, we find that the observed changes in the spectral slope are consistent with a mean particle size of aerosols varying with time between 1 and 2 μm. Observations with different solar zenith angles make it possible to characterize the aerosol layer at different altitudes, revealing a decrease of the particle size of aerosols as altitude increases.