Flow patterns of lobate debris aprons and lineated valley fill north of Ismeniae Fossae, Mars: Evidence for extensive mid-latitude glaciation in the Late Amazonian

A variety of Late Amazonian landforms on Mars have been attributed to the dynamics of ice-related processes. Evidence for large-scale, mid-latitude glacial episodes existing within the last 100 million to 1 billion years on Mars has been presented from analyses of lobate debris aprons (LDA) and lineated valley fill (LVF) in the northern and southern mid-latitudes. We test the glacial hypothesis for LDA and LVF along the dichotomy boundary in the northern mid-latitudes by examining the morphological characteristics of LDA and LVF surrounding two large plateaus, proximal massifs, and the dichotomy boundary escarpment north of Ismeniae Fossae (centered at 45.3°N and 39.2°E). Lineations and flow directions within LDA and LVF were mapped using images from the Context (CTX) camera, the Thermal Emission Imaging Spectrometer (THEMIS), and the High Resolution Stereo Camera (HRSC). Flow directions were then compared to topographic contours derived from the Mars Orbiter Laser Altimeter (MOLA) to determine the down-gradient components of LDA and LVF flow. Observations indicate that flow patterns emerge from numerous alcoves within the plateau walls, are integrated over distances of up to tens of kilometers, and have down-gradient flow directions. Smaller lobes confined within alcoves and superposed on the main LDA and LVF represent a later, less extensive glacial phase. Crater size-frequency distributions of LDA and LVF suggest a minimum (youngest) age of 100 Ma. The presence of ring-mold crater morphologies is suggestive that LDA and LVF are formed of near-surface ice-rich bodies. From these observations, we interpret LDA and LVF within our study region to result from formerly active debris-covered glacial flow, consistent with similar observations in the northern mid-latitudes of Mars. Glacial flow was likely initiated from the accumulation and compaction of snow and ice on plateaus and in alcoves within the plateau walls as volatiles were mobilized to the mid-latitudes during higher obliquity excursions. Together with similar analyses elsewhere along the dichotomy boundary, these observations suggest that multiple glacial episodes occurred in the Late Amazonian and that LDA and LVF represent significant reservoirs of non-polar ice sequestered below a surface lag for hundreds of millions of years.     

Concentric crater fill in the northern mid-latitudes of Mars: Formation processes and relationships to similar landforms of glacial origin

Hypotheses accounting for the formation of concentric crater fill (CCF) on Mars range from ice-free processes (e.g., aeolian fill), to ice-assisted talus creep, to debris-covered glaciers. Based on analysis of new CTX and HiRISE data, we find that concentric crater fill (CCF) is a significant component of Amazonian-aged glacial landsystems on Mars. We present mapping results documenting the nature and extent of CCF along the martian dichotomy boundary over −30 to 90°E latitude and 20-80°N longitude. On the basis of morphological analysis we classify CCF landforms into “classic” CCF and “low-definition” CCF. Classic CCF is most typical in the middle latitudes of the analysis area (∼30-50°N), while a range of degradation processes results in the presence of low-definition CCF landforms at higher and lower latitudes. We evaluate formation mechanisms for CCF on the basis of morphological and topographic analyses, and interpret the landforms to be relict debris-covered glaciers, rather than ice-mobilized talus or aeolian units. We examine filled crater depth-diameter ratios and conclude that in many locations, hundreds of meters of ice may still be present under desiccated surficial debris. This conclusion is consistent with the abundance of “ring-mold craters” on CCF surfaces that suggest the presence of near-surface ice. Analysis of breached craters and distal glacial deposits suggests that in some locations, CCF-related ice was once several hundred meters higher than its current level, and has sublimated significantly during the most recent Amazonian. Crater counts on ejecta blankets of filled and unfilled craters suggests that CCF formed most recently between ∼60 and 300 Ma, consistent with the formation ages of other martian debris-covered glacial landforms such as lineated valley fill (LVF) and lobate debris aprons (LDA). Morphological analysis of CCF in the vicinity of LVF and LDA suggests that CCF is a part of an integrated LVF/LDA/CCF glacial landsystem. Instances of morphological continuity between CCF, LVF, and LDA are abundant. The presence of formerly more abundant CCF ice, coupled with the integration of CCF into LVF and LDA, suggests the possibility that CCF represents one component of the significant Amazonian mid-latitude glaciation(s) on Mars.     

Amazonian northern mid-latitude glaciation on Mars: A proposed climate scenario

Recent geological observations in the northern mid-latitudes of Mars show evidence for past glacial activity during the late Amazonian, similar to the integrated glacial landsystems in the Dry Valleys of Antarctica. The large accumulation of ice (many hundreds of meters) required to create the observed glacial deposits points to significant atmospheric precipitation, snow and ice accumulation, and glacial flow. In order to understand the climate scenario required for these conditions, we used the LMD (Laboratoire de Météorologie Dynamique) Mars GCM (General Circulation Model), which is able to reproduce the present-day water cycle, and to predict past deposition of ice consistent with geological observations in many cases. Prior to this analysis, however, significant mid-latitude glaciation had not been simulated by the model, run under a range of parameters.In this analysis, we studied the response of the GCM to a wider range of orbital configurations and water ice reservoirs, and show that during periods of moderate obliquity (? = 25-35°) and high dust opacity (τ_dust = 1.5-2.5), broad-scale glaciation in the northern mid-latitudes occurs if water ice deposited on the flanks of the Tharsis volcanoes at higher obliquity is available for sublimation. We find that high dust contents of the atmosphere increase its water vapor holding capacity, thereby moving the saturation region to the northern mid-latitudes. Precipitation events are then controlled by topographic forcing of stationary planetary waves and transient weather systems, producing surface ice distribution and amounts that are consistent with the geological record. Ice accumulation rates of ∼10 mm yr⁻¹ lead to the formation of a 500-1000 m thick regional ice sheet that will produce glacial flow patterns consistent with the geological observations.     

MARCI and MOC observations of the atmosphere and surface cap in the north polar region of Mars

We used MGS-MOC and MRO-MARCI daily mapping images of the North Polar Region of Mars from 16 August 2005 (L_s = 270°) to 21 May 2009 (L_s = 270°), covering portions of three consecutive martian years (MY 27-MY 29), to observe the seasonal behavior of the polar ice cap and atmospheric phenomena. The rate of cap regression was similar in MY 28 and MY 29, but was advanced by 3.5° of L_s (∼7-8 sols) in MY 29. The spatial and temporal behaviors of dust and condensate clouds were similar in the two years and generally in accord with prior years. Dust storms (>100 km²) were observed in all seasons, with peak activity occurring at L_s = 10-20° from 50°N to 70°N and at L_s = 135-140° from 70°N to 90°N. The most active quadrant was 0-90°W in MY 28, shifting to 180-270°W in MY 29. The majority of regional storms in both years developed in longitudes from 10°W to 60°W. During late summer the larger storms obscure the North Polar Region in a cloud of dust that transitions to north polar hood condensate clouds around autumnal equinox.Changes in the distribution of perennial ice deposits, especially in Olympia Planum, were observed between the 2 years, with the MY 29 ice distribution being the most extensive observed to date. Modeling suggests that the small, bright ice patches on the residual cap are not the result of slope or elevation effects. Rather we suggest that they are the result of local meteorological effects on ice deposition. The annual darkening and brightening of peripheral areas of the residual cap around summer solstice can be explained by the sublimation of a brighter frost layer revealing an underlying darker, ice rich layer that itself either sublimes to reveal brighter material below or acts as a cold trap, attracting condensation of water vapor that brightens the surface. An alternative explanation invokes transport and deposition of dust on the surface from the cap interior, and later removal of that dust. The decrease in cap albedo and accompanying increase in near surface atmospheric stability may be related to the annual minimum of polar storm activity near northern summer solstice.     

A geomorphic analysis of Hale crater, Mars: The effects of impact into ice-rich crust

Hale crater, a 125 × 150 km impact crater located near the intersection of Uzboi Vallis and the northern rim of Argyre basin at 35.7°S, 323.6°E, is surrounded by channels that radiate from, incise, and transport material within Hale’s ejecta. The spatial and temporal relationship between the channels and Hale’s ejecta strongly suggests the impact event created or modified the channels and emplaced fluidized debris flow lobes over an extensive area (>200,000 km²). We estimate ∼10¹⁰ m³ of liquid water was required to form some of Hale’s smaller channels, a volume we propose was supplied by subsurface ice melted and mobilized by the Hale-forming impact. If 10% of the subsurface volume was ice, based on a conservative porosity estimate for the upper martian crust, 10¹² m³ of liquid water could have been present in the ejecta. We determine a crater-retention age of 1 Ga inside the primary cavity, providing a minimum age for Hale and a time at which we propose the subsurface was volatile-rich. Hale crater demonstrates the important role impacts may play in supplying liquid water to the martian surface: they are capable of producing fluvially-modified terrains that may be analogous to some landforms of Noachian Mars.     

Convective instability in the martian middle atmosphere

Dry convective instabilities in Mars’s middle atmosphere are detected and mapped using temperature retrievals from Mars Climate Sounder observations spanning 1.5 martian years. The instabilities are moderately frequent in the winter extratropics. The frequency and strength of middle atmospheric convective instability in the northern extratropics is significantly higher in MY 28 than in MY 29. This may have coupled with changes to the northern hemisphere mid-latitude and tropical middle atmospheric temperatures and contributed to the development of the 2007 global dust storm. We interpret these instabilities to be the result of gravity waves saturating within regions of low stability created by the thermal tides. Gravity wave saturation in the winter extratropics has been proposed to provide the momentum lacking in general circulation models to produce the strong dynamically-maintained temperature maximum at 1-2 Pa over the winter pole, so these observations could be a partial control on modeling experiments.     

Adsorptive fractionation of HDO on JSC MARS-1 during sublimation with implications for the regolith of Mars

A chamber was constructed to simulate the boundary between the ice table, regolith and atmosphere of Mars and to examine fractionation between H₂O and HDO during sublimation under realistic martian conditions of temperature and pressure. Thirteen experimental runs were conducted with regolith overlying the ice. The thickness and characteristic grain size of the regolith layer as well as the temperature of the underlying ice was varied. From these runs, values for the effective diffusivity, taking into account the effects of adsorption, of the regolith were derived. These effective diffusivities ranged from 1.8 × 10⁻⁴ m² s⁻¹ to 2.2 × 10⁻³ m² s⁻¹ for bare ice and from 2.4 × 10⁻¹¹ m² s⁻¹ to 2.0 × 10⁻⁹ m² s⁻¹ with an adsorptive layer present. From these, latent heats of adsorption of 8.6 ± 2.6 kJ mol⁻¹ and 9.3 ± 2.8 kJ mol⁻¹ were derived at ice-surface temperatures above 223 ± 8 K and 96 ± 28 kJ mol⁻¹ and 104 ± 31 kJ mol⁻¹ respectively for H₂O and HDO were derived at colder temperatures. For temperatures below 223 K, the effective diffusivity of HDO was found to be lower than the diffusivity of H₂O by 40% on average, suggesting that the regolith was adsorptively fractionating the sublimating gas with a fractionation factor of 1.96 ± 0.74. Applying these values to Mars predicts that adsorbed water on the regolith is enriched in HDO compared to the atmosphere, particularly where the regolith is colder. Based on current observations, the D/H ratio of the regolith may be as high as 21 ± 8 times VSMOW at 12°S and L_S = 357° if the regolith is hydrated primarily by the atmosphere, neglecting any hydration from subsurface ice.     

The evolution of exposed ice in a fresh mid-latitude crater on Mars

Recent observations of the surface of Mars have shown several fresh mid-latitude craters. Some of these craters show exposed ice (Byrne, S. et al. [2009]. Science 325, 1674-1676.). In some craters, albedo of ice slowly decreases, while in others, it remains nearly constant. We attempt to determine influence of the regolith structure on the rate of sublimation of ice. For this purpose we performed numerical simulations describing evolution of the exposed ice in model craters located at middle latitudes.We consider a new model for the structure and evolution of the material at- and beneath the crater floors. In contrast to the previous study by Dundas and Byrne (Dundas, C.M., Byrne, S. [2010]. Icarus 206, 716-728.) we do not investigate sublimation of dirty ice, and the related formation of a sublimation lag. Instead, we consider sublimation of a pure ice layer on top of layered regolith. In our model the observed reflectivity decreases due to the sublimation-driven changes of the optical properties of thinning clean ice. This offers an alternative to the deposition of the dust embedded in ice (sublimation lag).We have shown that in our model among many parameters affecting ice sublimation rate, volumetric fraction of water ice in the subsurface beneath the crater has the strongest influence. Hence observed darkening of the ice patch on the crater floor might be sufficient to determine the content of water ice in the subsurface. Our calculations show that an albedo decrease of fresh ice patches in mid-latitude craters can be explained by either strong dust sedimentation or, if this is excluded, by sublimation of a thin layer of water ice from the regolith with large thermal inertia. This is consistent with a large volumetric fraction of water ice beneath the crater floor and contributes to evidence for an extended subsurface water reservoir on Mars.The overall conclusion of our work is that a thin post-impact surface ice coating over ice-rich ground beneath the crater floors is consistent with the observations.     

Residual south polar cap of Mars: Stratigraphy, history, and implications of recent changes

The residual south polar cap (RSPC) of Mars includes a group of different depositional units of CO₂ ice undergoing a variety of erosional processes. Complete summer coverage of the RSPC by ∼6-m/pixel data of the Context Imager (CTX) on Mars Reconnaissance Orbiter (MRO) has allowed mapping and inventory of the units in the RSPC. Unit maps and estimated thicknesses indicate the total volume of the RSPC is currently <380 km³, and represents less than 3% of the total mass of the current Mars atmosphere. Scarp retreat rates in the CO₂ ice derived from comparison of High Resolution Imaging Science Experiment (HiRISE) data with earlier images are comparable to those obtained for periods up to 3 Mars years earlier. These rates, combined with sizes of depressions suggest that the oldest materials were deposited more than 125 Mars years ago. Most current erosion is by backwasting of scarps 1-12 m in height. This backwasting is initiated by a series of scarp-parallel fractures. In the older, thicker unit these fractures form about every Mars year; in thinner, younger materials they form less frequently. Some areas of the older, thicker unit are lost by downwasting rather than by the scarp retreat. A surprising finding from the HiRISE data is the scarcity of visible layering of RSPC materials, a result quite distinct from previous interpretations of layers in lower resolution images. Layers ∼0.1 m thick are exposed on the upper surfaces of some areas, but their timescale of deposition is not known. Late summer albedo changes mapped by the CTX images indicate local recycling of ice, although the amounts may be morphologically insignificant. Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) data show that the primary material of all the different forms of the RSPC is CO₂ ice with only small admixtures of water ice and dust.     

Groundwater discharge and gully formation on martian slopes

Young gullies and gully deposits on walls of martian craters have been cited as evidence that liquid water flowed on the surface of Mars relatively recently. Effects of variable environmental conditions at the surface of Mars are modeled and applied to the case of groundwater emergence from shallow aquifers to investigate whether groundwater is a viable source to enable the erosion of these gullies. The model includes detailed treatment of ice growth in the aquifer. Model results indicate that groundwater discharge can be maintained under the current environmental conditions if the aquifer permeability is like that of terrestrial gravel or higher, if the aquifer is 350 K or warmer, or if the aquifer is a brine with a freezing point depressed to 250 K or below. Groundwater discharge cannot be maintained for the conservative case of a cold, pure water, semi-pervious aquifer. Cold (275 K) pure water pervious (gravel) aquifers, warm (350 K) pure water semi-pervious aquifers, and cold (275 K) CaCl₂ brine semi-pervious aquifers all exhibit a dependence of discharge on season, latitude and slope orientation in our modeling. Seasonal, latitudinal and azimuthal discharge variations are strongest for cold CaCl₂ brine semi-pervious aquifers, with discharges from this aquifer type favoring equator-facing slopes at mid and high southern latitudes. At all latitudes and slope azimuths under our nominal conditions, the cold pure water pervious aquifer, the cold pure water semi-pervious aquifer and the cold CaCl₂ brine semi-pervious aquifer all freeze completely shortly after the simulations are started. Discharge restarts in the summer for the cold pure water pervious aquifer and the cold brine aquifer, but discharge does not restart for the cold pure water semi-pervious aquifer. The warm pure water semi-pervious aquifer maintains daily seeps throughout the year at all but high latitudes. In the case of the cold pure water pervious aquifer, approximately 500,000 m³ of water could be discharged from a mid-latitude, 150-m thick aquifer with a 20-m wide seepage face orientated towards the equator or the pole after a single undermining-induced event before ice growth seals the seepage face. For a brine semi-pervious aquifer with the same dimensions, 200-300 m³ of water could be released from a mid-latitude 20-m wide equator-facing seepage face before the fresh exposure is sealed for the fall and winter seasons. Our results do not rule out groundwater emergence as a means of creating some recent gullies, but they indicate that rather special and perhaps unusual conditions would be required.     

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.     

Dust devil sediment flux on Earth and Mars: Laboratory simulations

Laboratory simulations using the Arizona State University Vortex Generator (ASUVG) were run to simulate sediment flux in dust devils in terrestrial ambient and Mars-analog conditions. The objective of this study was to measure vortex sediment flux in the laboratory to yield estimations of natural dust devils on Earth and Mars, where all parameters may not be measured. These tests used particles ranging from 2 to 2000 μm in diameter and 1300 to 4800 kg m⁻³ in density, and the results were compared with data from natural dust devils on Earth and Mars. Typically, the cores of dust devils (regardless of planetary environment) have a pressure decrease of ∼0.1-1.5% of ambient atmospheric pressure, which enhances the lifting of particles from the surface. Core pressure decreases in our experiments ranged from ∼0.01% to 5.00% of ambient pressure (10 mbar Mars cases and 1000 mbar for Earth cases) corresponding to a few tenths of a millibar for Mars cases and a few millibars for Earth cases. Sediment flux experiments were run at vortex tangential wind velocities of 1-45 m s⁻¹, which typically correspond to ∼30-70% above vortex threshold values for the test particle sizes and densities. Sediment flux was determined by time-averaged measurements of mass loss for a given vortex size. Sediment fluxes of ∼10⁻⁶-10⁰ kg m⁻² s⁻¹ were obtained, similar to estimates and measurements for fluxes in dust devils on Earth and Mars. Sediment flux is closely related to the vortex intensity, which depends on the strength of the pressure decrease in the core (ΔP). This study found vortex size is less important for lifting materials because many different diameters can have the same ΔP. This finding is critical in scaling the laboratory results to natural dust devils that can be several orders of magnitude larger than the laboratory counterparts.     

The effects of the martian regolith on GCM water cycle simulations

This paper describes General Circulation Model (GCM) simulations of the martian water cycle focusing on the effects of an adsorbing regolith. We describe the 10-layer regolith model used in this study which has been adapted from the 1-D model developed by Zent, A.P., Haberle, R.M., Houben, H.C., Jakosky, B.M. [1993. A coupled subsurface-boundary layer model of water on Mars. J. Geophys. Res. 98 (E2), 3319-3337, February]. Even with a 30-min timestep and taking into account the effect of surface water ice, our fully implicit scheme compares well with the results obtained by Zent, A.P., Haberle, R.M., Houben, H.C., Jakosky, B.M. [1993. A coupled subsurface-boundary layer model of water on Mars. J. Geophys. Res. 98 (E2), 3319-3337, February]. This means, however, that the regolith is not able to reproduce the diurnal variations in column water vapour abundance of up to a factor of 2-3 as seen in some observations, with only about 10% of the atmospheric water vapour column exchanging with the subsurface on a daily basis. In 3-D simulations we find that the regolith adsorbs water preferentially in high latitudes. This is especially true in the northern hemisphere, where perennial subsurface water ice builds up poleward of 60° N at depths which are comparable to the Odyssey observations. Much less ice forms in the southern high latitudes, which suggests that the water ice currently present in the martian subsurface is not stable under present conditions and is slowly subliming and being deposited in the northern hemisphere. When initialising the model with an Odyssey-like subsurface water ice distribution the model is capable of forcing the simulated water cycle from an arbitrary state close to the Mars Global Surveyor Thermal Emission Spectrometer observations. Without the actions of the adsorbing regolith the equilibrated water cycle is found to be a factor of 2-4 too wet. The process by which this occurs is by adsorption of water during northern hemisphere summer in northern mid and high latitudes where it remains locked in until northern spring when the seasonal CO₂ ice cap retreats. At this time the water diffuses out of the regolith in response to increased temperature and is returned to the residual water ice cap by eddie transport.     

Lidar atmospheric measurements on Mars and Earth

The LIDAR instrument operating from the surface of Mars on the Phoenix Mission measured vertical profiles of atmospheric dust and water ice clouds at temperatures around −65 °C. An equivalent lidar system was utilized for measurements in the atmosphere of Earth where dust and cloud conditions are similar to Mars. Coordinated aircraft in situ sampling provided a verification of lidar measurement and analysis methods and also insight for interpretation of lidar derived optical parameters in terms of the dust and cloud microphysical properties. It was found that the vertical distribution of airborne dust above the Australian desert is quite similar to what is observed in the planetary boundary layer above Mars. Comparison with the in situ sampling is used to demonstrate how the lidar derived optical extinction coefficient is related to the dust particle size distribution. The lidar measurement placed a constraint on the model size distribution that has been used for Mars. Airborne lidar measurements were also conducted to study cirrus clouds that form in the Earth’s atmosphere at a similar temperature and humidity as the clouds observed with the lidar on Mars. Comparison with the in situ sampling provides a method to derive the cloud ice water content (IWC) from the Mars lidar measurements.     

Comparison of PFS and TES observations of temperature and water vapor in the martian atmosphere

The interval from L_s = 330° in Mars Year (MY) 26 until L_s = 84° in MY 27 has been used to compare and validate measurements from the Mars Global Surveyor Thermal Emission Spectrometer (TES) and the Mars Express Planetary Fourier Spectrometer (PFS). We studied differences between atmospheric temperatures observed by the two instruments. The best agreement between atmospheric temperatures was found at 50 Pa between 40°S and 40°N latitude, where differences were within ±5 K. For other atmospheric levels, differences as large as ∼25 K were observed between the two instruments at some locations. The largest temperature differences occurred mainly over the Hellas Planitia, Argyre Planitia, Tharsis and Valles Marineris regions.On this basis we report on the variability of the martian atmosphere during the 5.5 martian years of Mars climatology obtained by combining the two data sets from TES and PFS. Atmospheric temperatures at 50 Pa responded to the global-scale dust storms of MY 25 and in MY 28 raising temperatures from ∼220 K to ∼250 K during the daytime. An atmospheric temperature of ∼140 K at 50 Pa was observed poleward of 70°N during northern winter and poleward of 60°S during southern winter each year in both the PFS and TES results. Water vapor observed by the two spectrometers showed consistent seasonal and latitudinal variations.     

Water vapor variability in the north polar region of Mars from Viking MAWD and MGS TES datasets

Atmospheric water vapor abundances in Mars’ north polar region (NPR, from 60° to 90°N) are mapped as function of latitude and longitude for spring and summer seasons, and their spatial, seasonal, and interannual variability is discussed. Water vapor data are from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) and the Viking Orbiter (VO) Mars Atmospheric Water Detector (MAWD). The data cover three complete northern spring-summer seasons in 1977-1978, 2000-2001 and 2002-2003, and shorter periods of spring-summer seasons during 1975, 1999 and 2004. Long term interannual variability in the averaged NPR abundances may exist, with Viking MAWD observations showing twice as much water vapor during summer as the MGS TES observations more than 10 martian years (MY) later. While the averaged abundances are very similar in TES observations for the same season in different years, the spatial distributions in the early summer season do vary significantly year over year. Spatial and temporal variabilities increase between L_s ∼ 80-140°, which may be related to vapor sublimation from the North Polar Residual Cap (NPRC), or to changes in circulation. Spatial variability is observed on scales of ∼100 km and temporal variability is observed on scales of <10 sols during summer. During late spring the TES water vapor spatial distribution is seen to correlate with the low topography/low albedo region of northern Acidalia Planitia (270-360°E), and with the dust spatial distribution across the NPR during late spring-early summer. Non-uniform vertical distribution of water vapor, a regolith source or atmospheric circulation ‘pooling’ of water vapor from the NPRC into the topographic depression may be behind the correlation with low topography/low albedo. Sublimation winds carrying water vapor off the NPRC and lifting surface dust in the areas surrounding the NPRC may explain the correlation between the water vapor and dust spatial distributions. Correlation between water vapor and dust in MAWD data are only observed over low topography/low albedo area. Maximum water vapor abundances are observed at L_s = 105-115° and outside of the NPRC at 75-80°N; the TES data, however, do not extend over the NPRC and thus, this conclusion may be biased. Some water vapor appears to be released in plumes or ‘outbursts’ in the MAWD and TES datasets during late spring and early summer. We propose that the sublimation rate of ice varies across the NPRC with varying surface winds, giving rise to the observed ‘outbursts’ at some seasons.     

Aram Chaos and its constraints on the surface heat flux of Mars

The surface heat flux of a planet is an important parameter to characterize its internal activity and to determine its thermal evolution. Here we report on a new method to constrain the surface heat flux of Mars during the Hesperian. For this, we explore the consequences for the martian surface heat flux from a recently presented new hypothesis for the formation of Aram Chaos (Zegers, T.E., Oosthoek, J.H.P., Rossi, A.P., Blom, J.K., Schumacher, S. [2010]. Earth Planet. Sci. Lett. 297, 496-504. doi:10.1016/j.epsl.2010.06.049.). In this hypothesis the chaotic terrain is thought to have formed by melting of a buried ice sheet. The slow sedimentation and burial of the ice sheet led to an increased thermal insulation of the ice and subsequently to a temperature increase high enough to trigger melting and the formation of the subsurface lake. As these processes highly depend on the thermal properties of the subsurface and especially on the surface heat flux, it is possible to constrain the latter by using numerical simulations. Based on the hypothesis for the formation of Aram Chaos, we conducted an extensive parameter study to determine the parameter settings leading to sufficient melting of the buried ice sheet. We find that the surface heat flux in the Aram Chaos region during the Hesperian was most likely between 20 and 45 mW m⁻² with a possible maximum value of up to 60 mW m⁻².     

Effect of atmospheric dust loading on martian albedo measurements

This work is devoted to the analysis of the variation of albedo measured by orbiting instruments with atmospheric opacity on Mars. The study has been conduced by analysing Mars Global Surveyor Thermal Emission Spectrometer (MGS-TES) data from martian regions with different surface albedo.In support of these data, synthetic spectra with different surface albedo and atmospheric opacities have been computed, so that a comparison has been performed. The synthetic spectra have been retrieved by using two different grain sizes for suspended dust (0.5 and 1.2 μm), allowing a comparison between the two models and the observations.Using the DCI, a parameter describing the quantity of dust deposited on the surface, the effectiveness of the single scattering approximation has been tested for low atmospheric opacity by analysing the quality of the linear fit up to different atmospheric opacity.For more opaque conditions two kinds of fits have been applied to the data, linear and second-order degree polynomial. In this case, we found that the polynomial fit better describes the observations.The analysis of these data made it possible to notice a peculiar trend, already reported by Christensen (1988), of the albedo over Syrtis Major after the occurrence of dust storms, but, differently from that work, now the study of DCI together with atmospheric opacity and albedo allowed us to robustly confirm the hypothesis made by Christensen.Finally, the comparison between observations and synthetic spectra computed with models with different particles grain sizes indicates that dust particles of 0.5 μm diameter are the most effective to change the aerosol atmospheric opacity on Mars.     

Assimilation of thermal emission spectrometer atmospheric data during the Mars Global Surveyor aerobraking period

The Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft has produced an extensive atmospheric data set, beginning during aerobraking and continuing throughout the extended scientific mapping phase. Temperature profiles for the atmosphere below about 40 km, surface temperatures and total dust and water ice opacities, can be retrieved from infrared spectra in nadir viewing mode. This paper describes assimilation of nadir retrievals from the spacecraft aerobraking period, L_S=190°–260°, northern hemisphere autumn to winter, into a Mars general circulation model. The assimilation scheme is able to combine information from temperature and dust optical depth retrievals, making use of a model forecast containing information from the assimilation of earlier observations, to obtain a global, time-dependent analysis. Given sufficient temperature retrievals, the assimilation procedure indicates errors in the a priori dust distribution assumptions even when lacking dust observations; in this case there are relatively cold regions above the poles compared to a model which assumes a horizontally-uniform dust distribution. One major reason for using assimilation techniques is in order to investigate the transient wave behavior on Mars. Whilst the data from the 2-h spacecraft mapping orbit phase is much more suitable for assimilation, even the longer (45–24 h) period aerobraking orbit data contain useful information about the three-dimensional synoptic-scale martian circulation which the assimilation procedure can reconstruct in a consistent way. Assimilations from the period of the Noachis regional dust storm demonstrate that the combined assimilation of temperature and dust retrievals has a beneficial impact on the atmospheric analysis.