A 1 Gyr climate model for Mars: new orbital statistics and the importance of seasonally resolved polar processes

New results from a 1 Gyr integration of the martian orbit are presented along with a seasonally resolved energy balance climate model employed to illuminate the gross characteristics of the long-term atmospheric pressure evolution. We present a new analysis of the statistical variation of the martian obliquity and precession prior to and subsequent to the formation of the Tharsis uplift, and explore the long term effects on the martian climate. We find that seasonal polar cycles have a critical influence on the ability for the regolith to release CO₂ at high obliquities, and find that the atmospheric CO₂ actually decreases at high obliquities due to the cooling effect of polar deposits at latitudes where seasonal caps form. At low obliquity, the formation of massive, permanent polar caps depends critically on the values of the frost albedo, A_frost, and frost emissivity, ?_frost. Using our model with values of A_frost=0.67 and ?_frost=0.55, matched to the NASA Ames General Circulation Model (GCM) results (Haberle et al., 1993, J. Geophys. Res. 98, 3093-3123, and Haberle et al., 2003, Icarus 161, 66-89), we find that permanent caps only form at low obliquities (<13°), suggesting that any permanent deposits on the surface of Mars today may be residuals left over from a period of very low obliquity, or are the result of mechanisms not represented by this model. Thus, contrary to expectations, the martian atmospheric pressure is remarkable static over time, and decreases both at high and low obliquity. Also, from our one billion year orbital model, we present new results on the fraction of time Mars is expected to experience periods of low obliquity and high obliquity.     

Equilibrium rotational stability and figure of Mars

Studies extending over three decades have concluded that the current orientation of the martian rotation pole is unstable. Specifically, the gravitational figure of the planet, after correction for a hydrostatic form, has been interpreted to indicate that the rotation pole should move easily between the present position and a site on the current equator, 90° from the location of the massive Tharsis volcanic province. We demonstrate, using general physical arguments supported by a fluid Love number analysis, that the so-called non-hydrostatic theory is an inaccurate framework for analyzing the rotational stability of planets, such as Mars, that are characterized by long-term elastic strength within the lithosphere. In this case, the appropriate correction to the gravitational figure is the equilibrium rotating form achieved when the elastic lithospheric shell (of some thickness LT) is accounted for. Moreover, the current rotation vector of Mars is shown to be stable when the correct non-equilibrium theory is adopted using values consistent with recent, independent estimates of LT. Finally, we compare observational constraints on the figure of Mars with non-equilibrium predictions based on a large suite of possible Tharsis-driven true polar wander (TPW) scenarios. We conclude, in contrast to recent comparisons of this type based on a non-hydrostatic theory, that the reorientation of the pole associated with the development of Tharsis was likely less than 15° and that the thickness of the elastic lithosphere at the time of Tharsis formation was at least ∼50 km. Larger Tharsis-driven TPW is possible if the present-day gravitational form of the planet at degree 2 has significant contributions from non-Tharsis loads; in this case, the most plausible source would be internal heterogeneities linked to convection.     

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.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.     

Gully formation on Mars: Two recent phases of formation suggested by links between morphology, slope orientation and insolation history

The unusual 80 km diameter Noachian-aged Asimov crater in Noachis Terra (46°S, 5°E) is characterized by extensive Noachian-Hesperian crater fill and a younger superposed annulus of valleys encircling the margins of the crater floor. These valleys provide an opportunity to study the relationships of gully geomorphology as a function of changing slope orientation relative to solar insolation. We found that the level of development of gullies was highly correlated with slope orientation and solar insolation. The largest and most complex gully systems, with the most well-developed fluvial landforms, are restricted to pole-facing slopes. In contrast, gullies on equator-facing slopes are smaller, more poorly developed and integrated, more highly degraded, and contain more impact craters. We used a 1D version of the Laboratoire de Météorologie Dynamique GCM, and slope geometries (orientation and angle), driven by predicted spin-axis/orbital parameter history, to assess the distribution and history of surface temperatures in these valleys during recent geological history. Surface temperatures on pole-facing slopes preferential for water ice accumulation and subsequent melting are predicted to occur as recently as 0.5-2.1 Ma, which is consistent with age estimates of gully activity elsewhere on Mars. In contrast, the 1D model predicts that water ice cannot accumulate on equator-facing slopes until obliquities exceed 45°, suggesting they are unlikely to have been active over the last 5 Ma. The correlation of the temperature predictions and the geological evidence for age differences suggests that there were two phases of gully formation in the last few million years: an older phase in which top-down melting occurred on equator-facing slopes and a younger more robust phase on pole-facing slopes. The similarities of small-scale fluvial erosion features seen in the gullies on Mars and those observed in gullies cut by seasonal and perennial snowmelt in the Antarctic Dry Valleys supports a top-down melting origin for these gullies on Mars.     

Long term wind erosion on Mars

Wind erosion and deposition are powerful agents of surface change on Mars. Erosion is sensitive to the atmospheric density, so feedback between orbit variations and atmospheric density can enhance the sensitivity of erosion rates to orbital parameters. We use statistics derived from a 1 Gyr integration of the spin axis of Mars, coupled with runs of the NASA Ames 3-D general circulation model (GCM) at a variety of orbital conditions and pressures, to explore this feedback. We find that wind erosion in the GCM is associated with two factors: baroclinic zone winds, and strong cross-equator solstice flows. Both of these factors are influenced by topography, producing an asymmetry in the erosion pattern between the north and the south. The erosion model, averaged over 1 Gyr, produces potential erosion rates of in the north and in the south, which increase by an order of magnitude in an early 28 mbar atmosphere. The stability of these potential erosion patterns over geological time indicates that the lowland regions of Mars are continuously eroded, and that wind is capable of revealing ancient sedimentary deposits, or eroding substantial deposits that may have otherwise been preserved.     

MOC observations of the 2001 Mars planet-encircling dust storm

From 15 September 1997 through 21 January 2006, only a single planet-encircling martian dust storm was observed by MGS-MOC. The onset of the storm occurred on 26 June 2001 (Ls=184.7°), earliest recorded to date. It was initiated in the southern mid-to-low latitudes by a series of local dust storm pulses that developed along the seasonal cap edge in Malea and in Hellas basin (Ls=176.2°-184.4°). The initial expansion of the storm, though asymmetric, was very rapid in all directions (3-32 m s⁻¹). The main direction of propagation, however, was to the east, with the storm becoming planet encircling in the southern hemisphere on Ls=192.3°. Several distinct centers of active dust lifting were associated with the storm, with the longest persisting for 86 sols (Syria-Claritas). These regional storms helped generate and sustain a dust cloud (“haze”), which reached an altitude of about 60 km and a peak opacity of τdust∼5.0. By Ls=197.0°, the cloud had encircled the entire planet between 59.0° S and 60.0° N, obscuring all but the largest volcanoes. The decay phase began around Ls∼200.4° with atmospheric dust concentrations returning to nominal seasonal low-levels at Ls∼304.0°. Exponential decay time constants ranged from 30-117 sols. The storm caused substantial regional albedo changes (darkening and brightening) as a result of the redistribution (removal and deposition) of a thin veneer of surface dust at least 0.1-11.1 μm thick. It also caused changes in meteorological phenomena (i.e., dust storms, dust devils, clouds, recession of the polar caps, and possibly surface temperatures) that persisted for just a few weeks to more than a single Mars year. The redistribution of dust by large annual regional storms might help explain the long period (∼30 years) between the largest planet-encircling dust storms events.     

The Mars 1:2 resonant population

An excess of around 400 asteroids in the distribution of the semimajor axes of the asteroids is identified by means of numerical integrations as generated by a population of approximately 1000 asteroids evolving inside the exterior resonance 1:2 with Mars. Approximately 200 asteroids are librating around the asymmetric libration centers and their evolution in a time-scale of 1 million years appears stable but with a strong influence of Mars' eccentricity. The biggest Mars 1:2 resonant asteroid is (142) Polana.     

Mars high resolution gravity fields from MRO, Mars seasonal gravity, and other dynamical parameters

With 2 years of tracking data collection from the MRO spacecraft, there is noticeable improvement in the high frequency portion of the spherical harmonic Mars gravity field. The new JPL Mars gravity fields, MRO110B and MRO110B2, show resolution near degree 90. Additional years of MGS and Mars Odyssey tracking data result in improvement for the seasonal gravity changes which compares well to global circulation models and Odyssey neutron data and Mars rotation and precession (). Once atmospheric dust is accounted for in the spacecraft solar pressure model, solutions for Mars solar tide are consistent between data sets and show slightly larger values (k₂ = 0.164 ± 0.009, after correction for atmospheric tide) compared to previous results, further constraining core models. An additional 4 years of Mars range data improves the Mars ephemeris, determines 21 asteroid masses and bounds solar mass loss (dGM_Sun/dt < 1.6 × 10⁻¹³ GM_Sun year⁻¹).     

Preservation of Late Amazonian Mars ice and water-related deposits in a unique crater environment in Noachis Terra: Age relationships between lobate debris tongues and gullies

The Amazonian period of Mars has been described as static, cold, and dry. Recent analysis of high-resolution imagery of equatorial and mid-latitude regions has revealed an array of young landforms produced in association with ice and liquid water; because near-surface ice in these regions is currently unstable, these ice-and-water-related landforms suggest one or more episodes of martian climate change during the Amazonian. Here we report on the origin and evolution of valley systems within a degraded crater in Noachis Terra, Asimov Crater. The valleys have produced a unique environment in which to study the geomorphic signals of Amazonian climate change. New high-resolution images reveal Hesperian-aged layered basalt with distinctive columnar jointing capping interior crater fill and providing debris, via mass wasting, for the surrounding annular valleys. The occurrence of steep slopes (>20°), relatively narrow (sheltered) valleys, and a source of debris have provided favorable conditions for the preservation of shallow-ice deposits. Detailed mapping reveals morphological evidence for viscous ice flow, in the form of several lobate debris tongues (LDT). Superimposed on LDT are a series of fresh-appearing gullies, with typical alcove, channel, and fan morphologies. The shift from ice-rich viscous-flow formation to gully erosion is best explained as a shift in martian climate, from one compatible with excess snowfall and flow of ice-rich deposits, to one consistent with minor snow and gully formation. Available dating suggests that the climate transition occurred >8 Ma, prior to the formation of other small-scale ice-rich flow features identified elsewhere on Mars that have been interpreted to have formed during the most recent phases of high obliquity. Taken together, these older deposits suggest that multiple climatic shifts have occurred over the last tens of millions of years of martian history.     

Origin and sustainability of the population of asteroids captured in the exterior resonance 1:2 with Mars

At present, approximately 1500 asteroids are known to evolve inside or sticked to the exterior 1:2 resonance with Mars at a ? 2.418 AU, being (142) Polana the largest member of this group. The effect of the forced secular modes superposed to the resonance gives rise to a complex dynamical evolution. Chaotic diffusion, collisions, close encounters with massive asteroids and mainly orbital migration due to the Yarkovsky effect generate continuous captures to and losses from the resonance, with a fraction of asteroids remaining captured over long time scales and generating a concentration in the semimajor axis distribution that exceeds by 20% the population of background asteroids. The Yarkovsky effect induces different dynamics according to the asteroid size, producing an excess of small asteroids inside the resonance. The evolution in the resonance generates a signature on the orbits, mainly in eccentricity, that depends on the time the asteroid remains captured inside the resonance and on the magnitude of the Yarkovsky effect. The greater the asteroids, the larger the time they remain captured in the resonance, allowing greater diffusion in eccentricity and inclination. The resonance generates a discontinuity and mixing in the space of proper elements producing misidentification of dynamical family members, mainly for Vesta and Nysa-Polana families. The half-life of resonant asteroids large enough for not being affected by the Yarkovsky effect is about 1 Gyr. From the point of view of taxonomic classes, the resonant population does not differ from the background population and the excess of small asteroids is confirmed.     

Polar wander on Mars: Evidence in the geoid

The geoid of Mars is dominated by its equilibrium figure and by the effect of the Tharsis rise. To investigate the rotational stability of Mars prior to the rise of Tharsis, we produced a residual non-hydrostatic geoid without Tharsis. First the hydrostatic component of the present-day flattening was removed. This procedure was performed using a 6% non-hydrostatic component of flattening, a value set by the spin axis precession rate of Mars. Then zonal spherical harmonics up to degree 6 centered on Tharsis were removed. Finally, the resultant residual geoid was evaluated for rotational stability by comparing polar and equatorial moments at 4050 trial pole positions. If the spin axis of ancient Mars was secularly stable, our analysis indicates that substantial polar wander has occurred with the rise of Tharsis. Stable spin axis positions on the non-hydrostatic residual figure of Mars are 15° to 90° from the present-day poles. This result is consistent with previously proposed paleopoles based on magnetic anomalies, geomorphology, and grazing impacts.     

Evolution of the obliquities of the giant planets in encounters during migration

Tsiganis et al. [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461] have proposed that the current orbital architecture of the outer Solar System could have been established if it was initially compact and Jupiter and Saturn crossed the 2:1 orbital resonance by divergent migration. The crossing led to close encounters among the giant planets, but the orbital eccentricities and inclinations were damped to their current values by interactions with planetesimals. Brunini [Brunini, A., 2006. Nature 440, 1163-1165] has presented widely publicized numerical results showing that the close encounters led to the current obliquities of the giant planets. We present a simple analytic argument which shows that the change in the spin direction of a planet relative to an inertial frame during an encounter between the planets is very small and that the change in the obliquity (which is measured from the orbit normal) is due to the change in the orbital inclination. Since the inclinations are damped by planetesimal interactions on timescales much shorter than the timescales on which the spins precess due to the torques from the Sun, especially for Uranus and Neptune, the obliquities should return to small values if they are small before the encounters. We have performed simulations using the symplectic integrator SyMBA, modified to include spin evolution due to the torques from the Sun and mutual planetary interactions. Our numerical results are consistent with the analytic argument for no significant remnant obliquities.     

Reorientation of planets with lithospheres: The effect of elastic energy

It is commonly assumed that internal energy dissipation will ultimately drive planets to principal axis rotation, i.e., where the rotation vector is aligned with the maximum principle axis, since this situation corresponds to the minimum rotational energy state. This assumption simplifies long-term true polar wander (TPW) studies since the rotation pole can then be found by diagonalizing the appropriate (non-equilibrium) inertia tensor. We show that for planets with elastic lithospheres the minimum energy state does not correspond to principal axis rotation. As the planet undergoes reorientation elastic energy is stored in the deforming lithosphere, and the state of minimum total energy is achieved before principal axis rotation. We find solutions for the TPW of planets that include this effect by calculating the elastic stresses associated with deformation, and then minimizing the total (rotational and elastic) energy. These expressions indicate that the stored elastic energy acts to reduce the effective size of the driving load (relative to predictions which do not include this energy term). Our derivation also yields expressions for the TPW-induced stress field that generalizes several earlier results. As an illustration of the new theory, we consider TPW driven by the development of the Tharsis volcanic province on Mars. Once the size of the Tharsis load and the Mars model is specified, the extended theory yields a more limited range on the possible TPW.     

The cause for the north-south orientation of the crustal dichotomy and the equatorial location of Tharsis on Mars

The crustal dichotomy and the Tharsis rise are the most prominent topographic features on Mars. The dichotomy is largely an expression of different crustal thicknesses in the northern and southern hemispheres, while Tharsis is centered near the equator at the dichotomy boundary. However, the cause for the orientation of the dichotomy and the equatorial location of Tharsis remains poorly understood. Here we show that the crustal thickness variations associated with the dichotomy may have driven true polar wander, establishing the north-south orientation of the dichotomy very early in martian history. Such a reorientation that placed the dichotomy boundary near the equator would also have constrained the Tharsis region on the dichotomy boundary to have originated near the equator. We present a scenario for the early generation and subsequent reorientation of the hemispheric dichotomy, although the reorientation is independent of the formation mechanism. Our results also have implications for the sharply different remanent magnetizations between the two hemispheres.     

Evolution of Mercury's obliquity

Mercury has a near-zero obliquity, i.e. its spin axis is nearly perpendicular to its orbital plane. The value of the obliquity must be known precisely in order to constrain the size of the planet's core within the framework suggested by Peale [Peale, S.J., 1976. Nature 262, 765-766]. Rambaux and Bois [Rambaux, N., Bois, E., 2004. Astron. Astrophys. 413, 381-393] have suggested that Mercury's obliquity varies on thousand-year timescales due to planetary perturbations, potentially ruining the feasibility of Peale's experiment. We use a Hamiltonian approach (free of energy dissipation) to study the spin-orbit evolution of Mercury subject to secular planetary perturbations. We can reproduce an obliquity evolution similar to that of Rambaux and Bois [Rambaux, N., Bois, E., 2004. Astron. Astrophys. 413, 381-393] if we integrate the system with a set of initial conditions that differs from the Cassini state. However the thousand-year oscillations in the obliquity disappear if we use initial conditions corresponding to the equilibrium position of the Cassini state. This result indicates that planetary perturbations do not force short-period, large amplitude oscillations in the obliquity of Mercury. In the absence of excitation processes on short timescales, Mercury's obliquity will remain quasi-constant, suggesting that one of the important conditions for the success of Peale's experiment is realized. We show that interpretation of data obtained in support of this experiment will require a precise knowledge of the spin-orbit configuration, and we provide estimates for two of the critical parameters, the instantaneous Laplace plane orientation and the orbital precession rate from numerical fits to ephemeris data. Finally we provide geometrical relationships and a scheme for identifying the correct initial conditions required in numerical integrations involving a Cassini state configuration subject to planetary perturbations.     

Did tidal deformation power the core dynamo of Mars?

We first show that 7 out of the 20 giant impact basins of Mars recently reported by Frey [Frey, H., 2008. Geophys. Res. Lett. 35. L13203] trace a great circle on Mars. The other five basins trace another great circle and still the other three basins trace yet another great circle. The latter great circle is in good agreement with the pre-Tharsis equator of Mars that is estimated from modeling crustal magnetic anomalies [Arkani-Hamed, J., 2001. Geophys. Res. Lett. 28, 3409-3412] and diagonalizing the moment of inertia of Mars after removing the loading effects of Tharsis bulge [Sprenke, K.F., Baker, L.L., Williams, A.F., 2005. Icarus 174, 486-489]. It is shown in this paper that the three great circles were likely the equatorial plane of Mars at certain times and Mars experienced appreciable polar wander. The great circles also indicate that the asteroids that created the basins were satellites of Mars whose orbits decayed in time through spin-orbit coupling with tidally deforming Mars, and eventually impacted on the planet creating the giant basins at around 4 Ga. The orbital dynamics of four largest asteroids show that they could have orbited Mars for several hundred million years if they were retrograde satellites. Continual elliptical straining of otherwise circular fluid streamlines of the liquid core of Mars by tidal deformation could have exerted a strong strain that was large enough to overcome dissipation and excite the elliptical instability inside the core. We investigate the physical properties of the martian core that are required to allow the tidal deformation to power the core dynamo, i.e., the growth time of the elliptical instability to become shorter than the dissipation time. The tidal energy dissipation rate inside Mars caused by even only one of the 4 largest asteroids is found to be over two orders of magnitude greater than the magnetic energy dissipation rate in the core, indicating that if only one of the 4 largest asteroids were orbiting in retrograde sense, it would have likely powered the core dynamo of Mars for several hundred million years.     

Pole wandering of Mars: Evidence from paleomagnetic poles

We use the mapping-phase high-altitude magnetic measurements provided by Mars Global Surveyor (MGS) between March 1999 and April 2003 to model nine relatively isolated magnetic anomalies of Mars. Each anomaly is modeled with an elliptical prism. Each component of the observed magnetic field is modeled independently using an elliptical prism in order to assess the reliability of the results and suppress non-crustal and nearby crustal source contaminations. The paleomagnetic pole positions are obtained from the magnetization vectors of the model source bodies. We clean the data by removing the bad tracks and then divide the entire data into two sets that are measured at different times. Applying covariance analysis in the Fourier domain to two maps of the same magnetic component that are derived from the two sets provides a means to extract the most common features of the maps. The quality of a model is evaluated and only good models are used in the final geophysical interpretation. Most poles that come from good models cluster in the Tharsis region, suggesting that Mars experienced polar motion since the magnetic source bodies were magnetized.     

The formation and evolution of youthful gullies on Mars: Gullies as the late-stage phase of Mars’ most recent ice age

Gullies are extremely young erosional/depositional systems on Mars that have been carved by an agent that was likely to have been comprised in part by liquid water [Malin, M.C., Edgett, K.S., 2000. Evidence for recent groundwater seepage and surface runoff on Mars. Science 288, 2330-2335; McEwen, A.S. et al., 2007. A closer look at water-related geologic activity on Mars. Science 317, 1706-1709]. The strong latitude and orientation dependencies that have been documented for gullies require (1) a volatile near the surface, and (2) that insolation is an important factor for forming gullies. These constraints have led to two categories of interpretations for the source of the volatiles: (1) liquid water at depth beneath the melting isotherm that erupts suddenly (“groundwater”), and (2) ice at the surface or within the uppermost layer of soil that melts during optimal insolation conditions (“surface/near-surface melting”). In this contribution we synthesize global, hemispheric, regional and local studies of gullies across Mars and outline the criteria that must be met by any successful explanation for the formation of gullies. We further document trends in both hemispheres that emphasize the importance of top-down melting of recent ice-rich deposits and the cold-trapping of atmospherically-derived H₂O frost/snow as important components in the formation of gullies. This provides context for the incorporation of high-resolution multi-spectral and hyper-spectral data from the Mars Reconnaissance Orbiter that show that (1) cold-trapping of seasonal H₂O frost occurs at the alcove/channel-level on contemporary Mars; (2) gullies are episodically active systems; (3) gullies preferentially form in the presence of deposits plausibly interpreted as remnants of the Late Amazonian emplacement of ice-rich material; and (4) gully channels frequently emanate from the crest of alcoves instead of the base, showing that alcove generation is not necessarily a product of undermining and collapse at these locations, a prediction of the groundwater model. We interpret these various lines of evidence to mean that the majority of gullies on Mars are explained by the episodic melting of atmospherically emplaced snow/ice under spin-axis/orbital conditions characteristic of the last several Myr.     

MOC observations of four Mars year variations in the south polar residual cap of Mars

The residual south polar cap of Mars (RSPC) is distinct from the residual north polar cap both in composition and in morphology. CO₂ frost in the RSPC is stabilized by its high albedo during southern spring and summer despite the relatively large insolation during that period. The morphology of the RSPC in summer displays a bewildering variety of depressions that are formed in relatively thin layers of CO₂. The increase of the size of these depressions between each of the first three years of Mars Global Surveyor (MGS) observations may possibly signal some sort of climate change on the planet. For example, the erosion of the bright plateaus might reduce the RSPC albedo and affect the energy balance. The Mars Orbiter Cameras (MOC) on MGS observed Mars for four consecutive martian years before contact with the spacecraft was lost in late 2006. During this period coverage of the polar regions was particularly dense because MGS flew over them on every orbit. In this paper we report on the four-year behavior of the morphological features in the RSPC and on the large-scale variability in RSPC albedo over the period. The changes in the size of the surface features in the RSPC due to backwasting that were first observed between Mars years (MY) 24 and 25 and subsequently between MY25 and M26 was observed to continue at the same rate through MY 27. The results indicate that on average thicker layers in the RSPC retreat faster than thinner ones, roughly in proportion to their thickness. We argue that a simple difference in porosity between the A and B layers can explain this difference although other factors could be involved. The large-scale albedo of the RSPC decreases as the depressions are uncovered by sublimation of seasonal CO₂. However, any interannual differences in albedo due to the backwasting process are masked by interannual differences in the summer dust opacity in the RSPC region.