The effects and characteristics of atmospheric dust during martian global dust storm 2001A

We present retrieved trends in dust optical depth, dust effective radius and surface temperature from our analysis of Mars Global Surveyor Thermal Emission Spectrometer daytime data from global dust storm 2001A, and describe their significance for the martian dust cycle. The dust optical depth becomes correlated with surface pressure during southern spring and summer in years both with and without a global dust storm, indicating that global dust mixing processes are important at those seasons. The correlation is low at other times of the year. We found that the observed decay of optical depths at the later stages of the dust storm match, to first-order, theoretical values of clearing from Stokes–Cunningham fallout of the dust. Zonally averaged effective radius is constant within standard deviation of results (between 1.2 and 2.0 μm, with a global mean for all seasons of 1.7 μm), at all latitudes and seasons except at southern latitudes of 35° and higher around equinoxes in both martian years, where it is larger than average (2–3 μm). The emergence and disappearance of these larger particles correlates with observations of polar cap edge storms at those latitudes. Northern latitude observations under similar conditions did not yield a similar trend of larger average effective radii during the equinoxes. We also report on a linear correlation between daytime surface temperature drop and rise in optical depth during the global dust storm. Global dust storm 2001A produced a significant optical depth and surface temperature change.     

Retrieval algorithm for atmospheric dust properties from Mars Global Surveyor Thermal Emission Spectrometer data during global dust storm 2001A

We present a new parameter retrieval algorithm for Mars Global Surveyor Thermal Emission Spectrometer data. The algorithm uses Newtonian first-order sensitivity functions of the infrared spectrum in response to variations in physical parameters to fit a model spectrum to the data at 499, 1099, and 1301 cm⁻¹. The algorithm iteratively fits the model spectrum to data to simultaneously retrieve dust extinction optical depth, effective radius, and surface temperature. There are several sources of uncertainty in the results. The assumed dust vertical distribution can introduce errors in retrieved optical depth of a few tens of percent. The assumed dust optical constants can introduce errors in both optical depth and effective radius, although the systematic nature of these errors will not affect retrieval of trends in these parameters. The algorithm does not include the spectral signature of water ice, and hence data needs to be filtered against this parameter before the algorithm is applied. The algorithm also needs sufficient dust spectral signature, and hence surface-to-atmosphere temperature contrast, to successfully retrieve the parameters. After the application of data filters the algorithm is both relatively accurate and very fast, successfully retrieving parameters, as well as meaningful parameter variability and trends from tens of thousands of individual spectra on a global scale (Elteto, A., Toon, O.B. [2010]. Icarus, this issue). Our results for optical depth compare well with TES archive values when corrected by the single scattering albedo. Our results are on average 1–4 K higher in surface temperatures from the TES archive values, with greater differences at higher optical depths. Our retrieval of dust effective radii compare well with the retrievals of Wolff and Clancy (Wolff, M.J., Clancy, R.T. [2003]. J. Geophys. Res. 108 (E9), 5097) for the corresponding data selections from the same orbits.     

A microphysically-based approach to modeling emissivity and albedo of the martian seasonal caps

A new model of albedo and emissivity of the martian seasonal caps represented as porous CO₂ slabs containing spherical voids and dust particles is described. In the model, a radiative transfer model is coupled with a microphysical model in order to link changes in albedo and emissivity to changes in porosity caused by ice metamorphism. The coupled model is capable of reproducing temporal changes in the spectra of the caps taken by the Thermal Emission Spectrometer onboard the Mars Global Surveyor and it can be used as the forward model in the retrievals of the caps' physical properties (porosity, dust abundance, void and dust grain size) from the spectra. Preliminary results from such inversion studies are presented.     

Methane adsorption on a martian soil analog: An abiogenic explanation for methane variability in the martian atmosphere

Recent observations suggest methane in the martian atmosphere is variable on short spatial and temporal scales. However, to explain the variability by loss reactions requires production rates much larger than expected. Here, we report results of laboratory studies of methane adsorption onto JSC-Mars-1, a martian soil simulant, and suggest that this process could explain the observations. Uptake coefficient (γ) values were measured as a function of temperature using a high-vacuum Knudsen cell able to simulate martian temperature and pressure conditions. Values of γ were measured from 115 to 135 K, and the data were extrapolated to higher temperatures with more relevance to Mars. Adsorptive uptake was found to increase at lower temperatures and larger methane partial pressures. Although only sub-monolayer methane surface coverage is likely to exist under martian conditions, a very large mineral surface area is available for adsorption as atmospheric methane can diffuse meters into the regolith. As a result, significant methane may be temporarily lost to the regolith on a seasonal time scale. As this weak adsorption is fully reversible, methane will be re-released into the atmosphere when surface and subsurface temperatures rise and so no net loss of methane occurs. Heterogeneous interaction of methane with martian soil grains is the only process proposed thus far which contains both rapid methane loss and rapid methane production mechanisms and is thus fully consistent with the reported variability of methane on Mars.     

Limits on the trapping of atmospheric CH4 in martian polar ice analogs

Recent detection of methane (CH₄) on Mars has generated interest in possible biological or geological sources, but the factors responsible for the reported variability are not understood. Here we explore one potential sink that might affect the seasonal cycling of CH₄ on Mars - trapping in ices deposited on the surface. Our apparatus consisted of a high-vacuum chamber in which three different Mars ice analogs (water, carbon dioxide, and carbon dioxide clathrate hydrates) were deposited in the presence of CH₄ gas. The ices were monitored for spectroscopic evidence of CH₄ trapping using transmission Fourier-Transform Infrared (FT-IR) spectroscopy, and during subsequent sublimation of the ice films the vapor composition was measured using mass spectrometry (MS). Trapping of CH₄ in water ice was confirmed at deposition temperatures <100 K which is consistent with previous work, thus validating the experimental methods. However, no trapping of CH₄ was observed in the ice analogs studied at warmer temperatures (140 K for H₂O and CO₂ clathrate, 90 K for CO₂ snow) with approximately 10 mTorr CH₄ in the chamber. From experimental detection limits these results provide an upper limit of 0.02 for the atmosphere/ice trapping ratio of CH₄. If it is assumed that the trapping mechanism is linear with CH₄ partial pressure and can be extrapolated to Mars, this upper limit would indicate that less than 1% is expected to be trapped from the largest reported CH₄ plume, and therefore does not represent a significant sink for CH₄.     

Extension of atmospheric dust loading to high altitudes during the 2001 Mars dust storm: MGS TES limb observations

Mars Global Surveyor (MGS) visible (solarband bolometer) and thermal infrared (IR) spectral limb observations from the Thermal Emission Spectrometer (TES) support quantitative profile retrievals for dust opacity and particle sizes during the 2001 global dust event on Mars. The current analysis considers the behavior of dust lifted to altitudes above 30 km during the course of this storm; in terms of dust vertical mixing, particle sizes, and global distribution. TES global maps of visible (solarband) limb brightness at 60 km altitude indicate a global-scale, seasonally evolving (over 190-240° solar longitudes, L_S) longitudinal corridor of vertically extended dust loading (which may be associated with a retrograde propagating, wavenumber 1 Rossby wave). Spherical radiative transfer analysis of selected limb profiles for TES visible and thermal IR radiances provide quantitative vertical profiles of dust opacity, indicating regional conditions of altitude-increasing dust mixing ratios. Observed infrared spectral dependences and visible-to-infrared opacity ratios of dust scattering over 30-60 km altitudes indicate particle sizes characteristic of lower altitudes (cross-section weighted effective radius, ), during conditions of significant dust transport to these altitudes. Conditions of reduced dust loading at 30-60 km altitudes present smaller dust particle sizes . These observations suggest rapid meridional transport at 30-80 km altitudes, with substantial longitudinal variation, of dust lifted to these altitudes over southern hemisphere atmospheric regions characterized by extraordinary (m/s) vertical advection velocities. By L_S=230° dust loading above 50 km altitudes decreased markedly at southern latitudes, with a high altitude (60-80 km) haze of fine (likely) water ice particles appearing over 10°S-40°N latitudes.     

Assessing spectral evidence of aqueous activity in two putative martian paleolakes

We evaluate the evidence for the presence of mineral spectral signatures indicative of the past presence of water at two putative paleolakes on Mars using observations by the Mars Reconnaissance Orbiter (MRO) Compact Reconnaissance Image Spectrometer for Mars (CRISM). CRISM spectra of both sites are consistent with laboratory spectra of Mg-rich phyllosilicates. Our analysis represents the first detailed evaluation of these locations. The spatial occurrence and association with topographic features within the craters is distinctly different for the two sites. The occurrence of these minerals supports the conclusion that water was once active in the areas sampled by these craters. The distribution of the phyllosilicates in Luqa does not provide distinctive evidence for the presence of a previous standing body of water and is consistent with either impact emplacement or post-impact alteration. For Cankuzo, the phyllosilicate distribution provides evidence of a layer in the crater wall indicative of aqueous activity, but does not require a paleolake.     

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.     

The impact of martian mesoscale winds on surface temperature and on the determination of thermal inertia

Radiative control of surface temperature is a key characteristic of the martian environment and its low-density atmosphere. Here we show through meteorological modeling that surface temperature can be far from radiative equilibrium over numerous sloping terrains on Mars, where nighttime mesoscale katabatic winds impact the surface energy budget. Katabatic circulations induce both adiabatic atmospheric heating and enhancement of downward sensible heat flux, which then becomes comparable to radiative flux and acts to warm the ground. Through this mechanism, surface temperature can increase up to 20 K. One consequence is that warm signatures of surface temperature over slopes, observed through infrared spectrometry, cannot be systematically associated with contrasts of intrinsic soil thermal inertia. Apparent thermal inertia maps retrieved thus far possibly contain wind-induced structures. Another consequence is that surface temperature observations close to sloping terrains could allow the validation of model predictions for martian katabatic winds, provided contrasts in intrinsic thermal inertia can be ruled out. The thermal impact of winds is mostly discussed in this paper in the particular cases of Olympus Mons/Lycus Sulci and Terra Meridiani but is generally significant over any sloped terrains in low thermal inertia areas. It is even general enough to apply under daytime conditions, thereby providing a possible explanation for observed afternoon surface cooling, and to ice-covered terrains, thereby providing new insights on how winds could have shaped the present surface of Mars.     

The role of aqueous alteration in the formation of martian soils

Martian equatorial dark regions are dominated by unweathered materials and it has often been assumed that they have not been significantly altered from their source lithology. The suite of minerals present is consistent with a basaltic composition and there has been no need to invoke additional processes to explain the origin of these compositions. We have begun to question this result based on detailed observations using a variety of datasets. Locally derived dark soils have a mineralogy distinct from that of adjacent rocky surfaces; most notably a lower olivine content. This pattern is common for many surfaces across the planet. Previous work using detailed measurements acquired within the Gusev Plains has shown that olivine dissolution via acidic weathering may explain chemical trends observed between rock rinds and interiors. Mineralogical trends obtained from rocks and soils within the Gusev Plains are more prominent than the elemental trends and support previous results that indicate that dissolution of olivine has occurred. However, clear differences are also present in elemental abundances that indicate a variety of inputs and processes are likely responsible for the formation of martian dark soils. Despite the potential complexity of source materials and processes, it appears that most martian dark regions have likely experienced aqueous alteration and chemical weathering appears to be closely linked to the mechanical breakdown of materials. Regardless of the responsible mechanism, there appears to be a general, though not perfect, correlation between elevated olivine abundance and high-thermal inertia surfaces on Mars.     

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.     

Derivation of martian surface slope characteristics from directional thermal infrared radiometry

Directional thermal infrared measurements of the martian surface is one of a variety of methods that may be used to characterize surface roughness and slopes at scales smaller than can be obtained by orbital imagery. Thermal Emission Spectrometer (TES) emission phase function (EPF) observations show distinct apparent temperature variations with azimuth and emission angle that are consistent with the presence of warm, sunlit and cool, shaded slopes at typically ∼0.1 m scales. A surface model of a Gaussian distribution of azimuth independent slopes (described by θ-bar) is combined with a thermal model to predict surface temperature from each viewing angle and azimuth of the TES EPF observation. The models can be used to predict surface slopes using the difference in measured apparent temperature from 2 separate 60-70° emission angle observations taken ∼180° in azimuth relative to each other. Most martian surfaces are consistent with low to moderate slope distributions. The slope distributions display distinct correlations with latitude, longitude, and albedo. Exceptionally smooth surfaces are located at lower latitudes in both the southern highlands as well as in high albedo dusty terrains. High slopes are associated with southern high-latitude patterned ground and north polar sand dunes. There is little apparent correlation between high resolution imagery and the derived θ-bar, with exceptions such as duneforms. This method can be used to characterize potential landing sites by assuming fractal scaling behavior to meter scales. More precisely targeted thermal infrared observations from other spacecraft instruments are capable of significantly reducing uncertainty as well as reducing measurement spot size from 10s of kilometers to sub-kilometer scales.     

Modeling sublimation of ice exposed by new impacts in the martian mid-latitudes

New impacts in the martian mid-latitudes have exposed near-surface ice. This ice is observed to slowly fade over timescales of months. In the present martian climate, exposed surface ice is unstable during summer months in the mid-latitudes and will sublimate. We model the sublimation of ice at five new impact sites and examine the implications of its persistence. Even with generally conservative assumptions, for most reasonable choices of parameters it is likely that over a millimeter of sublimation occurred in the period during which the ice was observed to fade. The persistence of visible ice through such sublimation suggests that the ice is relatively pure rather than pore-filling. Such ice could be analogous to the nearly pure ice observed by the Phoenix Lander in the “Dodo-Goldilocks” trench and suggests that the high ice contents reported by the Mars Odyssey Gamma Ray Spectrometer at high latitudes extend to the mid-latitudes. Our observations are consistent with a model of the martian ice table in which a layer with high volumetric ice content overlies pore-filling ice, although other structures are possible.     

Assessing the power law hypothesis for the size-frequency distribution of terrestrial and martian dust devils

Competing hypotheses for the diameter dependence of terrestrial and martian dust devil frequency are assessed using new field observations from two sites in the southwestern United States. We show that at diameters less than 12 m, our observed dust devil size-frequency distributions are better fit by an exponential function than by a power law formulation, and discuss the implications for larger dust devils on Earth and Mars.     

Principal components analysis of Mars in the near-infrared

Principal components analysis and target transformation are applied to near-infrared image cubes of Mars in a study to disentangle the spectra into a small number of spectral endmembers and characterize the spectral information. The image cubes are ground-based telescopic data from the NASA Infrared Telescope Facility during the 1995 and 1999 near-aphelion oppositions when ice clouds were plentiful [ [Clancy et al., 1996] and [56]], and the 2003 near-perihelion opposition when ice clouds are generally limited to topographically high regions (volcano cap clouds) but airborne dust is more common [Martin, L.J., Zurek, R.W., 1993. J. Geophys. Res. 98 (E2), 3221-3246]. The heart of the technique is to transform the data into a vector space along the dimensions of greatest spectral variance and then choose endmembers based on these new “trait” dimensions. This is done through a target transformation technique, comparing linear combinations of the principal components to a mineral spectral library. In general Mars can be modeled, on the whole, with only three spectral endmembers which account for almost 99% of the data variance. This is similar to results in the thermal infrared with Mars Global Surveyor Thermal Emission Spectrometer data [Bandfield, J.L., Hamilton, V.E., Christensen, P.R., 2000. Science 287, 1626-1630]. The globally recovered surface endmembers can be used as inputs to radiative transfer modeling in order to measure ice abundance in martian clouds [Klassen, D.R., Bell III, J.F., 2002. Bull. Am. Astron. Soc. 34, 865] and a preliminary test of this technique is also presented.     

Mineralogical constraints on the high-silica martian surface component observed by TES

The Thermal Emission Spectrometer (TES) has observed a high-silica material in the dark regions of Mars that is spectrally similar to obsidian glass and may have a volcanic origin. An alternate interpretation is that the spectrally amorphous material consists of clay minerals or some other secondary material, formed by chemical alteration of surface rocks. The regions where this material is observed (e.g., Acidalia Planitia) have relatively high spectral contrast, suggesting that the high-silica material exists as coarse particulates, indurated soils or cements, within rocks, or as indurated coatings on rock surfaces. The geologic interpretation of this spectral result has major implications for understanding magmatic evolution and weathering processes on Mars. One of the complications in interpreting spectral observations of glasses and clay minerals is that both are structurally and compositionally complex. In this study, we perform a detailed spectroscopic analysis of indurated smectite clay minerals and relate their thermal emission spectral features to structural and crystal chemical properties. We examine the spectral similarities and differences between smectite clay minerals and obsidian glass from a structural-chemical perspective, and make further mineralogical interpretations from previous TES results. The results suggest that neither smectite clays nor any clay mineral with similar structural and chemical properties can adequately explain TES observations of high-silica materials in some martian dark regions. If the spectrally amorphous materials observed by TES do represent an alteration product, then these materials are likely to be poorly crystalline aluminosilicates. While all clay minerals have Si/O ratios ?0.4, the position of the emissivity minimum at Mars suggests a Si/O ratio of 0.4-0.5. The spectral observation could be explained by the existence of a silica-rich alteration product, such as Al- or Fe-bearing opal, an intimate physical mixture of relatively pure silica and other aluminosilicates (such as clay minerals or clay precursors), or certain zeolites. The chemical alteration of basaltic rocks on Mars to phyllosilicate-poor, silica-rich alteration products provides a geologically reasonable and consistent explanation for the global TES surface mineralogical results.     

Ultraviolet dust aerosol properties as observed by MARCI

Observations by the Mars Color Imager (MARCI) on board the Mars Reconnaissance Orbiter (MRO) in two ultraviolet (UV, Bands 6 and 7; 258 nm, and 320 nm, respectively) and one visible (Band 1, 436 nm) channels of the 2007 planet encircling dust storm are combined with those made by the two Mars Exploration Rovers (MERs) to better characterize the single scattering albedo (ω₀) of martian dust aerosols. Exploiting the low contrast of the surface in the UV (and blue) as well as the reduced importance of surface reflectance under very dusty conditions, we utilize the sampling of photometric angles by the MARCI cross-track geometry to synthesize an analog of the classical Emergence Phase Function (EPF). This so-called “pseudo-EPF”, used in conjunction with the “ground-truth” measurements provided by the MERs, is able to effectively isolate the effects of the dust ω₀. The motivation for this approach is the elimination of a significant portion of the type of uncertainty involved in many previous radiative transfer analyses. Furthermore, we produce a self-consistent set of complex refractive indices (m=n+ik) through our use of an explicit microphysical representation of the aerosol scattering properties. Because of uncertainty in the exact size of the dust particles during the epoch of the observations, we consider two effective particle radii (r_eff) to cover the range anticipated from the literature: 1.6 and 1.8 μm. The resulting set of model-data comparisons, ω₀, and m are presented along with an assessment of potential sources of error and uncertainty. Analysis of the Band 1 results is limited to ω₀ as a “proof-of-concept” for our approach through a comparison to contemporaneous CRISM EPF results at 440 nm. The derived ω₀ are: assuming , and 0.765, for Bands 6, 7, and 1, respectively; for , for the same band order. For either r_eff case, the total estimated error is 0.022, 0.019, and 0.010, again for Bands 6, 7, and 1. We briefly discuss our retrievals, including the asymmetry parameter (g) associated with our model phase functions, within the context of previous efforts, with an emphasis on the improved precision of our results compared to those in the literature. We also suggest several applications of our results, including an extension of the dust climatological record using MARCI Band 7 pseudo-EPFs outside of 2007 global dust event. Initial work on this particular application using a sample of 135 pseudo-EPFs near the MERs suggests that optical depth retrievals with a precision in the range 0.2-0.4 may be possible under moderate loading conditions (i.e., τ < 1.5).     

Multitemporal observations of identical active dust devils on Mars with the High Resolution Stereo Camera (HRSC) and Mars Orbiter Camera (MOC)

Active dust devils were observed in Syria Planum in Mars Observer Camera - Wide Angle (MOC-WA) and High Resolution Stereo Camera (HRSC) imagery acquired on the same day with a time delay of ∼26 min. The unique operating technique of the HRSC allowed the measurement of the traverse velocities and directions of motion. Large dust devils observed in the HRSC image could be retraced to their counterparts in the earlier acquired MOC-WA image. Minimum lifetimes of three large (avg. ∼700 m in diameter) dust devils are ∼26 min, as inferred from retracing. For one of these large dust devil (∼820 m in diameter) it was possible to calculate a minimum lifetime of ∼74 min based on the measured horizontal speed and the length of its associated dust devil track. The comparison of our minimum lifetimes with previous published results of minimum and average lifetimes of small (∼19 m in diameter, avg. min. lifetime of ∼2.83 min) and medium (∼185 m in diameter, avg. min. lifetime of ∼13 min) dust devils imply that larger dust devils on Mars are active for much longer periods of time than smaller ones, as it is the case for terrestrial dust devils. Knowledge of martian dust devil lifetimes is an important parameter for the calculation of dust lifting rates. Estimates of the contribution of large dust devils (>300-1000 m in diameter) indicate that they may contribute, at least regionally, to ∼50% of dust entrainment by dust devils into the atmosphere compared to the dust devils <300 m in diameter given that the size-frequency distribution follows a power-law. Although large dust devils occur relatively rarely and the sediment fluxes are probably lower compared to smaller dust devils, their contribution to the background dust opacity by dust devils on Mars could be at least regionally large due to their longer lifetimes and ability of dust lifting into high atmospheric layers.     

Fresnel modeling of hematite crystal surfaces and application to martian hematite spherules

The gray crystalline hematite at Meridiani Planum first discovered by the Mars Global Surveyor Thermal Emission Spectrometer (MGS-TES) instrument occurs as spherules that have been interpreted as concretions. Analysis of the TES and mini-TES spectra shows that no 390 cm⁻¹ feature is present in the characteristic martian hematite spectrum. Here, we incorporate the mid-IR optical constants of hematite into a simple Fresnel reflectance model to understand the effect of emission angle and crystal morphology on the presence or absence of the 390 cm⁻¹ feature in an IR hematite spectrum. Based on the results we offer two models for the internal structure of the martian hematite spherules.