East-west trending magnetic anomalies in the Southern Hemisphere of Mars: Modeling analysis and interpretation

Maps of the vector components of the martian crustal magnetic field over the strongly magnetized Terra Cimmeria/Sirenum region are constructed using Mars Global Surveyor magnetometer data. Although pronounced east-west trending anomalies are present on the radial and north field component maps at the mapping altitude (∼360-380 km), these trends are much less prominent at the lower aerobraking altitude (∼90-150 km). Comparisons with similar maps produced using artificial data at the aerobraking altitude indicate that elongated sources in this region may have maximum lengths along the martian surface of ∼500 km and maximum aspect ratios of ∼2. Iterative forward modeling of several relatively isolated anomalies in the mapped region yields paleomagnetic pole positions consistent with those estimated in previous studies of other anomalies using mapping phase and science phasing orbit data. On this basis, it is inferred that sources in the studied region are most probably magnetized primarily in northward or southward directions. Using this additional constraint, iterative forward modeling is then applied to determine a magnetization distribution that is consistent with data at both the aerobraking altitude and the mapping altitude. The model magnetization distribution, which includes 41 discrete sources, again indicates no highly elongated sources. An examination of surface geology in the region as well as a consideration of the global distribution of anomalies suggests that magmatic intrusions (e.g., subsurface dike swarms), cooling in the presence of water, are the most likely sources of the magnetic anomalies.     

Mechanical modeling of thrust faults in the Thaumasia region, Mars, and implications for the Noachian heat flux

Insight into the state of the early martian lithosphere is gained by modeling the topography above surface breaking thrust faults in the southern Thaumasia region. Crater counts of key surface units associated with the faulting indicate a scarp emplacement in the late Noachian-early Hesperian periods between 4.0 and 3.7 Gyr. The seismogenic layer thickness at the time of faulting is constrained to 27-35 km and 21-28 km for the two scarps investigated, implying paleo geothermal gradients of 12-18 and 15-23 K km⁻¹, corresponding to heat flows of 24-36 and 30-46 mW m⁻². The heat flow values obtained in this study are considerably lower than those derived from rift flank uplift at the close-by Coracis Fossae for a similar time period, indicating that surface heat flow is a strong function of regional setting. If viewed as representative for magmatically active and inactive regions, the thermal gradients at rifts and scarps span the range of admissible global mean values. This implies , with the true value probably being closer to the lower bound.     

An improved crustal magnetic field map of Mars from electron reflectometry: Highland volcano magmatic history and the end of the martian dynamo

We apply improved kinetic modeling of electron transport in the martian thermosphere to fit pitch angle distributions measured by the Mars Global Surveyor (MGS) Magnetometer/Electron Reflectometer (MAG/ER), together with appropriate filtering, binning, averaging and error correction techniques, to create the most reliable ER global map to date of crustal magnetic field magnitude at 185 km altitude, with twice the spatial resolution and considerably higher sensitivity to crustal fields than global maps of magnetic field components produced with MAG data alone. This map compares favorably to sparsely sampled dayside MAG data taken at similar altitudes, insofar as a direct comparison is meaningful. Using this map, we present two case studies. The first compares the magnetic signatures of two highland volcanoes, concluding that the comparatively greater thermal demagnetization at Syrtis Major compared with Tyrrhena Patera is likely due to a higher ratio of intruded to extruded magmas. The second uses the map along with topographic data to compare the magnetic signatures and crater retention ages of the demagnetized Hellas impact basin and magnetized Ladon impact basin. From this comparison, we determine that the martian global dynamo magnetic field went from substantial to very weak or nonexistent in the absolute model age time interval 4.15±0.05 to 4.07±0.05 Ga ago.     

Yearly comparisons of the martian polar caps: 1999-2003 Mars Orbiter Camera observations

The Mars Global Surveyor Mars Orbiter Camera wide-angle cameras were used to obtain images of the north and south seasonal and residual polar caps between 1999 and 2003. Wide-angle red camera images were used in assembling mosaics of the north and south polar recessions and regression rates were measured and compared. There are small variations in the north polar recession between 2000 and 2002, especially between LS=7° and LS=50°, however there is no evidence for the plateau in the recession curves that has been observed in some prior years. The south polar recession changes very little from year to year, and the 2001 dust storm had little if any effect on the average cap recession that year. Albedo values of the geographic north pole were measured using wide-angle red and blue camera images, and the residual south polar cap configuration was compared between the three years observed by MOC. The albedo of the geographic north pole generally varies between 0.5 and 0.6 as measured from MOC wide-angle red camera images. There were only minor variations near the edges of the residual south polar cap between the three years examined.     

Stability analysis of the martian obliquity during the Noachian era

We performed numerical simulations of the obliquity evolution of Mars during the Noachian era, at which time the giant planets were on drastically different orbits than today. For the preferred primordial configuration of the planets we find that there are two large zones where the martian obliquity is stable and oscillates with an amplitude lower than 20°. These zones occur at obliquities below 30° and above 60°; intermediate values show either resonant or chaotic behaviour depending on the primordial orbits of the terrestrial planets.     

Monitoring the perennial martian northern polar cap with MGS MOC

We have used the Mars Global Surveyor Mars Orbiter Camera Wide Angle (MGS MOC WA) dataset to study albedo trends on the martian northern residual cap. Six study regions were selected, the Chasma Boreale source region, three regions near the center of the cap (“fish hook” region, latitude = 87°; “bottle opener” region, latitude = 87°, “steep-shallow” region, latitude = 85°), and two lower latitude regions (crater, latitude = 77°, and polar outlier, latitude = 82°), and the albedos of these six regions were examined. These regions were chosen due to their good temporal coverage in the MOC dataset, as well as having been studied by other researchers (Bass et al., 2000, Icarus 144, 382-396; Calvin and Titus, 2004, Lunar Planet. Sci. XXXV, Abstract 1455). The picture which emerges is complex. Most areas experience a combination of darkening and brightening through the northern summer; only one area consistently brightens (the polar outlier region). A good deal of interannual repeatability in each region's albedo behavior is seen, however. Possible causes for the observed complex behaviors include dust deposition from late summer storms, sintering of frost grains over the course of the summer, and cold trapping of volatiles on bright, cold surfaces.     

Viscous relaxation of craters within the martian south polar layered deposits

The South Polar Layered Deposits (PLD) are of fundamental importance to martian climatology, as they may comprise the largest reservoir of near-surface water on Mars. The South PLD exhibit relatively young crater retention surface ages, which are widely attributed to recent resurfacing. However, we show that both constructional and destructional resurfacing mechanisms (such as dust deposition and water ice sublimation, respectively) are inconsistent with the size, depth, and spatial distributions of South PLD craters. We demonstrate that another process—viscous creep relaxation of dusty water ice—is more compatible with the observed cratering of South PLD surfaces. The results of our finite element relaxation simulations suggest that, despite their apparent youthfulness, the PLD have been stable for at least several hundred million and perhaps even over a billion years. Consequently, our modeling implies that the time scales for the formation and preservation of the layers characteristic of the South (and possibly North) PLD are much longer than generally assumed.     

Eolian controlled modification of the martian south polar layered deposits

Data from the Mars Orbiter Laser Altimeter (MOLA) and Mars Orbiter Camera (MOC) aboard the Mars Global Surveyor (MGS) mission and the Thermal Emission Imaging System (THEMIS) aboard the Mars Odyssey mission have revealed unique surface features in a particular region of the South Polar Layered Deposits (SPLD). The dominant morphology is large-scale quasi-parallel grooves that extend for hundreds of kilometers with only tens of meters of vertical relief, that we have termed here the “Wire Brush” terrain. The grooves are also transected by disjointed, yet roughly continuous, low-relief sinuous ridges that cross roughly perpendicular to the trend-direction of the large-scale grooves and show only tens of meters of relief. We interpret these ridges to be eroded remnants of folded layers. At the northern end of the large-scale grooves there are non-symmetric mounds. They are frequently preceded by a significant depression and/or trailing grooves that are parallel to the Wire Brush trend. We find that a two-stage process involving winds that intermittently remove a low-density crust exposing the underlying ice to ablation is the interpretation that best explains the multitude of features observed here. These features appear to be currently inactive indicating higher winds in previous epochs.     

The effects of martian orbital variations upon the sublimation and relaxation of north polar troughs and scarps

The North Polar Layered Deposits (PLD) of Mars are climatologically significant because they represent the largest actively-exchanging reservoir of martian water. The kilometer-scale topography of the North PLD is dominated by troughs and scarps, which exhibit highly-correlated surface slopes and total depths. The most widespread theories of PLD evolution presume that the asymmetrical nature of North PLD troughs (characterized by equatorward-facing slopes that are generally steeper than poleward-facing slopes) is the result of preferential H₂O sublimation from equatorward-facing trough walls. However, our orbitally-modulated simulations indicate that there is no long-term sublimation advantage of equatorward-facing trough walls, because of the effects of obliquity upon the slope dependence of sublimation rate. In contrast, we propose that viscous relaxation of subsurface water is consistent with the slope and depth distributions of North PLD troughs and scarps. The results of our finite element simulations suggest that a mere 2 K difference in the subsurface temperatures of opposing trough walls is sufficient to produce the observed slope disparity, due to the slower rate of uplift of colder poleward-facing trough walls. Both our sublimation and relaxation simulations indicate that present-day North PLD troughs have formed since 5 Ma and are not sites of long-term deposition; additionally, our results imply that glacial flow may govern the large-scale evolution of the North PLD, especially at high obliquity.     

Electron reflectometry in the martian atmosphere

The technique of electron reflectometry, a method for remote estimation of planetary magnetic fields, is expanded from its original use of mapping crustal magnetic fields at the Moon to achieving the same purpose at Mars, where the presence of a substantial atmosphere complicates matters considerably. The motion of solar wind electrons, incident on the martian atmosphere, is considered in detail, taking account of the following effects: the electrons' helical paths around the magnetic field lines to which they are bound, the magnetic mirror force they experience due to converging field lines in the vicinity of crustal magnetic anomalies, their acceleration/deceleration by electrostatic potentials, their interactions with thermal plasma, their drifts due to magnetic field line curvature and perpendicular electric fields and their scattering off, and loss of energy through a number of different processes to, atmospheric neutrals. A theoretical framework is thus developed for modeling electron pitch angle distributions expected when a spacecraft is on a magnetic field line which is connected to both the martian crust and the interplanetary magnetic field. This framework, along with measured pitch angle distributions from the Mars Global Surveyor (MGS) Magnetometer/Electron Reflectometer (MAG/ER) experiment, can be used to remotely measure crustal magnetic field magnitudes and atmospheric neutral densities at ∼180 km above the martian datum, as well as estimate average parallel electric fields between 200 and 400 km altitude. Detailed analysis and full results, concerning the crustal magnetic field and upper thermospheric density of Mars, are left to two companion papers.     

Predicted and observed magnetic signatures of martian (de)magnetized impact craters

The current morphology of the martian lithospheric magnetic field results from magnetization and demagnetization processes, both of which shaped the planet. The largest martian impact craters, Hellas, Argyre, Isidis and Utopia, are not associated with intense magnetic fields at spacecraft altitude. This is usually interpreted as locally non- or de-magnetized areas, as large impactors may have reset the magnetization of the pre-impact material. We study the effects of impacts on the magnetic field. First, a careful analysis is performed to compute the impact demagnetization effects. We assume that the pre-impact lithosphere acquired its magnetization while cooling in the presence of a global, centered and mainly dipolar magnetic field, and that the subsequent demagnetization is restricted to the excavation area created by large craters, between 50- and 500-km diameter. Depth-to-diameter ratio of the transient craters is set to 0.1, consistent with observed telluric bodies. Associated magnetic field is computed between 100- and 500-km altitude. For a single-impact event, the maximum magnetic field anomaly associated with a crater located over the magnetic pole is maximum above the crater. A 200-km diameter crater presents a close-to-1-nT magnetic field anomaly at 400-km altitude, while a 100-km diameter crater has a similar signature at 200-km altitude. Second, we statistically study the 400-km altitude Mars Global Surveyor magnetic measurements modelled locally over the visible impact craters. This approach offers a local estimate of the confidence to which the magnetic field can be computed from real measurements. We conclude that currently craters down to a diameter of 200 km can be characterized. There is a slight anti-correlation of −0.23 between magnetic field intensity and impact crater diameters, although we show that this result may be fortuitous. A complete low-altitude magnetic field mapping is needed. New data will allow predicted weak anomalies above craters to be better characterized, and will bring new constraints on the timing of the martian dynamo and on Mars’ evolution.     

A top-down origin for martian mantle plumes

The two main volcanic centers on Mars, Tharsis and Elysium, are often interpreted in terms of mantle plume hotspots, even though there are several problems with the plume hypothesis for Mars. We present results of 2D cylindrical shell numerical mantle convection experiments in which we try to ascertain whether flushing of the hot lower mantle could provide a mechanism for the generation of a small number of plume-like features, i.e., localized upwelling of hot material. In this scenario the formation of hot upwellings is driven from the top by cold downwellings rather than from a hot thermal boundary layer at the CMB. First we construct a range of Mars interior structure models consistent with observations in order to demonstrate that the presence of a thin lower mantle in the martian interior is a viable scenario. Then we use a series of numerical convection experiments to investigate the effects of solid-state phase transitions, different stratified and temperature-dependent viscosity models, and the presence of a thick southern hemisphere crust on the operation of such a mechanism. Our results show that it is possible to generate hot strong localized upwellings from top-down dynamics if the lithosphere is thin or actively involved in the convective pattern. The presence of a thick, immobile, insulating southern hemisphere crust reduces the number of upwellings, and the perovskite phase transition causes a focusing of the upwellings. Further experiments demonstrate that an initial 500 Myr phase of mobile lid is sufficient to start this process create an upwelling which is stable for billions of years.     

Continuous monitoring of nightside upper thermospheric mass densities in the martian southern hemisphere over 4 martian years using electron reflectometry

Details are presented of an improved technique to use atmospheric absorption of magnetically reflecting solar wind electrons to constrain neutral mass densities in the nightside martian upper thermosphere. The helical motion of electrons on converging magnetic field lines, through an extended neutral atmosphere, is modeled to enable prediction of loss cone pitch angle distributions measured by the Magnetometer/Electron Reflectometer (MAG/ER) experiment on Mars Global Surveyor at 400 km altitude. Over the small fraction of Mars' southern hemisphere (∼2.5%) where the permanent crustal magnetic fields are both open to the solar wind and sufficiently strong as to dominate the variable induced martian magnetotail field, spherical harmonic expansions of the crustal fields are used to prescribe the magnetic field along the electron's path, allowing least-squares fitting of measured loss cones, in order to solve for parameters describing the vertical neutral atmospheric mass density profile from 160 to 230 km. Results are presented of mass densities in the southern hemisphere at 2 a.m. LST at the mean altitude of greatest sensitivity, 180 km, continuously over four martian years. Seasonal variability in densities is largely explained by orbital and latitudinal changes in dayside insolation that impacts the nightside through the resulting thermospheric circulation. However, the physical processes behind repeatable rapid, late autumnal cooling at mid-latitudes and near-aphelion warming at equatorial latitudes is not fully clear. Southern winter polar warming is generally weak or nonexistent over several Mars years, in basic agreement with MGS and MRO accelerometer observations. The puzzling response of mid-latitude densities from 160° to 200° E to the 2001 global dust storm suggests unanticipated localized nightside upper thermospheric lateral and vertical circulation patterns may accompany such storms. The downturn of the 11-year cycle of solar EUV flux is likely responsible for lower aphelion densities in 2004 and 2006 (Mars years 27 and 28).     

The evolution of the martian elastic lithosphere and implications for crustal and mantle rheology

Estimates of the martian elastic lithosphere thickness Te imply that Te increased from around 20 km in the Noachian to about 70 km in the Amazonian period. A phase of rapid lithospheric growth is observed during the Hesperian and we propose that this elastic thickness history is a consequence of the martian crustal rheology and its thermal evolution. A wet crustal rheology is found to generate a mechanically incompetent layer in the lower crust during the early evolution and the rapid growth of Te during the Hesperian results from the disappearance of this layer due to planetary cooling. The incompetent layer and the related rapid lithospheric growth are absent for a dry basaltic crustal rheology, which is therefore incompatible with the observations. Furthermore, we find that the observed elastic thickness evolution is best compatible with a wet mantle rheology, although a dry mantle cannot be ruled out. It therefore seems likely that rheologically significant amounts of water were retained in the Martian crust and mantle after planetary accretion.     

Radiative habitable zones in martian polar environments

The biologically damaging solar ultraviolet (UV) radiation (quantified by the DNA-weighted dose) reaches the martian surface in extremely high levels. Searching for potentially habitable UV-protected environments on Mars, we considered the polar ice caps that consist of a seasonally varying CO2 ice cover and a permanent H2O ice layer. It was found that, though the CO2 ice is insufficient by itself to screen the UV radiation, at approximately 1 m depth within the perennial H2O ice the DNA-weighted dose is reduced to terrestrial levels. This depth depends strongly on the optical properties of the H2O ice layers (for instance snow-like layers). The Earth-like DNA-weighted dose and Photosynthetically Active Radiation (PAR) requirements were used to define the upper and lower limits of the northern and southern polar Radiative Habitable Zone (RHZ) for which a temporal and spatial mapping was performed. Based on these studies we conclude that photosynthetic life might be possible within the ice layers of the polar regions. The thickness varies along each martian polar spring and summer between approximately 1.5 and 2.4 m for H2O ice-like layers, and a few centimeters for snow-like covers. These martian Earth-like radiative habitable environments may be primary targets for future martian astrobiological missions. Special attention should be paid to planetary protection, since the polar RHZ may also be subject to terrestrial contamination by probes.     

Accretion, loss, and fractionation of martian water

A recently published model for accretion of water by Mars from asteroids and comets beyond 2.5 AU calls for Mars to have accreted 6-27% of Earth's ocean (equivalent to 600-2700 meters depth on Mars' surface) with a D/H ratio 1.2-1.6 times Standard Mean Ocean Water. This model is shown to provide enough outgassed water to account for the subsequent evolution of D/H because of loss of hydrogen and the D/H ratio in water of hydration found in martian meteoric minerals, but enhanced rate loss of hydrogen from early Mars is necessary to permit adequate fractionation to occur during the lifetime of the planet.     

Sorted stone circles in Elysium Planitia, Mars: Implications for recent martian climate

We have found sorted stone circles and polygons near the equator of Mars, using new 25 cm/pixel NASA HiRISE (High Resolution Imaging Science Experiment) images. The sorted circles occur in geologically recent catastrophic flood deposits in the equatorial Elysium Planitia region, and are diagnostic of periglacial processes: sorted polygons do not form from volcanic activity, as has been suggested for non-sorted polygons in this region. These landforms indicate that (i) a long-lived, geologically recent, active cryoturbation layer of ground ice was present in the regolith, (ii) there was some degree of freeze-thaw, and thus (iii) there were sustained period(s), likely within the last 10 Ma, in which the martian climate was 40 to 60 K warmer than current models predict.     

Interpretation of Mars southern highlands high amplitude magnetic field with total gradient and fractal source modeling: New insights into the magnetic mystery of Mars

Using model studies, the total gradient (TG) of the Z-component magnetic field is shown to be a useful quantity for delineating sources of satellite-altitude magnetic anomalies; this field is used to constrain the location and lateral boundaries of sources of high amplitude magnetic anomalies of southern highlands of Mars. The TG field suggests two parallel linear and oppositely magnetized sources of 1000 and 1800 km length separated by 1000 km of region of intervening non-parallel sources. The simplest interpretation of the long, linear features is that they are zones of multitudinous crustal scale dikes formed in separate episodes of rifting, and not features associated with the mechanism of seafloor spreading. Forward modeling with uniformly magnetized sources suggests that magnetizations of the order of 10-50 A/m (40 km thickness) over ∼100 km width in the case of the southern source and of 12.5-27.5 A/m (40 km thickness) and ∼200 km width for the northern source are necessary to explain the Z-component amplitudes and features of the TG field. If the crustal magnetization on Mars were to be distributed fractally as on Earth, magnetizations matching the largest amplitude features on Mars may be spatially correlated from a 50-100 km distance range (β ∼ 3) to approaching nearly uniform magnetization (β ∼ 5) values. To keep magnetization intensity as small as possible, the higher end of β values are preferred, whereas, small amplitude anomaly features could be generated from sources with β ∼ 3. Many of the Mars anomaly features could be coalescence effects similar to the coalescence of anomaly features observed on http://icarus.cornell.edu/information/keywords.html.Earth.     

Valley network-fed, open-basin lakes on Mars: Distribution and implications for Noachian surface and subsurface hydrology

A new catalog of 210 open-basin lakes (lakes with outlet valleys) fed by valley networks shows that they are widely distributed in the Noachian uplands of Mars. In order for an outlet valley to form, water must have ponded in the basin to at least the level of the outlet. We use this relationship and the present topography to directly estimate the minimum amount of water necessary to flood these basins in the past. The volumes derived for the largest lakes (∼3×¹⁰⁴ to ∼2×¹⁰⁵ km³) are comparable to the largest lakes and small seas on modern Earth, such as the Caspian Sea, Black Sea, and Lake Baikal. We determine a variety of other morphometric properties of these lakes and their catchments (lake area, mean depth, volume, shoreline development, outlet elevation, and watershed area). Most candidate lakes have volumes proportional to and commensurate with their watershed area, consistent with precipitation as their primary source. However, other lakes have volumes that are anomalously large relative to their watershed areas, implying that groundwater may have been important in their filling. Candidate groundwater-sourced lakes are generally concentrated in the Arabia Terra region but also include the Eridania basin [Irwin, R.P., Howard, A.D., Maxwell, T.A., 2004a. J. Geophys. Res. 109, doi: 10.1029/2004JE002287. E12009; Irwin, R.P., Watters, T.R., Howard, A.D. Zimbelman, J.R., 2004b. J. Geophys. Res. 109, doi: 10.1029/2004JE002248. E09011] and several lakes near the dichotomy boundary. This areal distribution is broadly consistent with where groundwater should have reached the surface as predicted by current models. Both surface runoff and groundwater flow appear to have been important sources for lakes and lake chains, suggesting a vertically integrated hydrological system, the absence of a global cryosphere, and direct communication between the surface and subsurface hydrosphere of early Mars.     

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.