Charge transfer processes from low to intermediate energies

The theoretical treatment of charge transfer processes is developed using ab‐initio molecular calculations. Semi‐classical and quantal dynamical approaches are presented for the determination of cross sections and rate constants which are important data for space chemistry models. Accurate cross section values may be determined, taking account of rotational couplings with regard to the collision energy. Such theoretical approaches provide besides an insight into the mechanism of these processes with consideration of anisotropic and vibrational effects for collisions with diatomic molecular targets (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

High-resolution simulations of the impacts of asteroids into the venusian atmosphere III: further 3D models

We report on high-resolution three-dimensional calculations of oblique impacts into planetary atmospheres, specifically the atmosphere of Venus, extending the results of Korycansky et al. (2000, Icarus 146, 387-403; 2002, Icarus 157, 1-23). We have made calculations for impacts at 0°, 45°, and 60° from the vertical, different impactor velocities (10, 20, and 40 km s⁻¹), and different impactor masses and orientations. We present results for porous impactors using a simple model of porosity. We have investigated the sensitivity to initial conditions of the calculations [as a follow-up to the results found in Korycansky et al. (2002)] and resolution effects. For use in cratering calculations, we fit simple functions to the numerical results for mass and momentum that penetrate to a given altitude (column mass) and investigate the behavior of the fit coefficients as functions of impactor parameters such as mass, velocity, and impact angle. Generally speaking, the mass and momentum (and hence resulting crater diameters) depend primarily on impactor mass and mass of atmosphere encountered and weakly or not at all on other parameters such as impactor velocity, impact angle, or porosity. The column mass to which the last portion of the impactor penetrates is approximately equal to the mass of impactor at the top of the atmosphere before the impact takes place. Finally, we present the beginnings of a simplified but physically based model for the impactor and its fragments to reproduce the mass and momentum fluxes as a function of height during the impact.     

High-Resolution Simulations of the Impacts of Asteroids into the Venusian Atmosphere II: 3D Models

We compare high-resolution 2D and 3D numerical hydrocode simulations of asteroids striking the atmosphere of Venus. Our focus is on aerobraking and its effect on the size of impact craters. We consider impacts both by spheres and by the real asteroid 4769 Castalia, a severely nonspherical body in a Venus-crossing orbit. We compute mass and momentum fluxes as functions of altitude as global measures of the asteroid's progress. We find that, on average, the 2D and 3D simulations are in broad agreement over how quickly an asteroid slows down, but that the scatter about the average is much larger for the 2D models than for the 3D models. The 2D models appear to be rather strongly susceptible to the “butterfly effect,” in which tiny changes in initial conditions (e.g., 0.05% change in the impact velocity) produce quite different chaotic evolutions. By contrast, the global properties of the 3D models appear more reproducible despite seemingly large differences in initial conditions. We argue that this difference between 2D and 3D models has its root in the greater geometrical constraints present in any 2D model, and in particular in the global conservation of enstrophy in 2D that forces energy to pool in large-scale structures. It is the interaction of these artificial large-scale structures that causes slightly different 2D models to diverge so greatly. These constraints do not apply in 3D and large scale structures are not observed to form. A one-parameter modified pancake model reproduces the expected crater diameters of the 3D Castalias reasonably well.     

Measurement of Impact Ejecta from Regolith Targets in Oblique Impacts

This paper reports the results of experiments on projectile impact into regolith targets at various impact angles. Copper projectiles of 240 mg are accelerated to 197 to 272 m s⁻¹ using an electromagnetic gun. The ejecta are detected by thin Al foil targets as secondary targets, and the resulting holes on the foil are measured to derive the spatial distribution of the ejecta. The ejecta that penetrated the foil are concentrated toward the downrange azimuths of impacting projectiles in oblique impacts. In order to investigate the ejecta velocity distribution, the nondimensional volume of ejecta with velocities higher than a given value is calculated from the spatial distribution. In the case of the vertical impact of the projectile, most ejecta have velocities lower than 24% of the projectile speed (∼50 m s⁻¹), and there are only several ejecta with velocities higher than 72 m s⁻¹. This result confirms the existence of an upper limit to the ejection velocity in the ejecta velocity distribution (Hartmann cutoff velocity) (W. K. Hartmann, 1985, Icarus63, 69-98). On the other hand, it is found that, in the oblique impacts, there are a large number of ejecta with velocities higher than the Hartmann cutoff velocity. The relative quantity of ejecta above the Hartmann cutoff velocity increases as the projectile impact angle decreases. Taking these results with the results of S. Yamamoto and A. M. Nakamura (1997, Icarus128, 160-170) from impact experiments using an impact angle of 30°, it can be concluded that the ejecta from these regolith targets exhibit a bimodal velocity distribution. Below a few tens of m s⁻¹, we see the expected velocity distribution of ejecta, but above this velocity we see a separate group of high-velocity ejecta.     

The formation of floor-fractured craters in Xanthe Terra

Floor-fractured craters (FFC) are a peculiar form of degradation of impact craters defined by the presence of crevice networks and mesas affecting crater floors. They are preferentially distributed near chaotic terrains and outflow channels. The scope of this paper is to present a detailed systematic analysis of FFC at Xanthe Terra. FFC morphologies in this region are classified into five types making a picture of different stages of the same degradation process. FFC are geographically intermixed with un-fractured normal craters (non-FFC). Young craters are less prone to show this type of degradation, as suggested by fresh ejecta layer with preserved crater floor. Size distributions of FFC and non-FFC indicate that larger craters are preferentially fractured. Careful examinations of the crater floor elevations reveal that the crevices often extend deeper than the original crater cavity. Furthermore, an onset depth for the formation of FFC is evidenced from the difference of spatial distributions between FFC and non-FFC. Roof-collapsed depressions observed in the same region have been also documented and their characteristics suggest the removal of subsurface material at depth from about 1200 to 4000 m. These observations taken together suggest a subsurface zone of volume deficit at depth from 1 to 2 km down to several kilometers responsible for FFC formation. Then a scenario of FFC formations is presented in the context of groundwater discharge events at the late Hesperian. This scenario involves two key processes, Earth fissuring and piping erosion, known to occur with rapid groundwater migrations on Earth.     

Origin of the structure and planform of small impact craters in fractured targets: Endurance Crater at Meridiani Planum, Mars

We present observations and models that together explain many hallmarks of the structure and growth of small impact craters forming in targets with aligned fractures. Endurance Crater at Meridiani Planum on Mars (diameter ≈ 150 m) formed in horizontally-layered aeolian sandstones with a prominent set of wide, orthogonal joints. A structural model of Endurance Crater is assembled and used to estimate the transient crater planform. The model is based on observations from the Mars Exploration Rover Opportunity: (a) bedding plane orientations and layer thicknesses measured from stereo image pairs; (b) a digital elevation model of the whole crater at 0.3 m resolution; and (c) color image panoramas of the upper crater walls. This model implies that the crater’s current shape was mostly determined by highly asymmetric excavation rather than long-term wind-mediated erosion. We show that modal azimuths of conjugate fractures in the surrounding rocks are aligned with the square component of the present-day crater planform, suggesting excavation was carried farther in the direction of fracture alignments. This was previously observed at Barringer Crater in Arizona and we show the same relationship also holds for Tswaing Crater in South Africa. We present models of crater growth in which excavation creates a “stellate” transient cavity that is concave-cuspate in planform. These models reproduce the “lenticular-crescentic” layering pattern in the walls of some polygonal impact craters such as Endurance and Barringer Craters, and suggest a common origin for tear faults and some crater rays. We also demonstrate a method for detailed error analysis of stereogrammetric measurements of bedding plane orientations.     

Effective collision strengths for forbidden transitions among the 3s23p3 fine-structure levels of Cl IIIIII

Effective collision strengths for the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^o, ²D^o and ²P^o levels in the 3s²3p³ configuration of Cl  iii are presented. The calculation employs the multichannel R-matrix method to compute the electron-impact excitation collision strengths in a close-coupling expansion, which incorporates the lowest 23 LS target eigenstates of Cl  iii . These states are formed from the 3s²3p³, 3s3p⁴, 3s²3p²3d and 3s²3p²4s configurations. The Maxwellian-averaged effective collision strengths are presented graphically for all 10 fine-structure transitions over a wide range of electron temperatures appropriate for astrophysical applications [log  T (K)=3.3−log  T (K)=5.9]. Comparisons are made with the earlier seven-state close-coupling calculation of Butler & Zeippen, and in general excellent agreement is found in the low-temperature region where a comparison is possible [log  T (K)=3.3−log  T (K)=4.7]. However, discrepancies of up to 30 per cent are found to occur for the forbidden transitions which involve the ⁴S^o ground state level, particularly for the lowest temperatures considered. At the higher temperatures, the present data are the only reliable results currently available.     

Scaling of oblique impacts in frictional targets: Implications for crater size and formation mechanisms

Almost every meteorite impact occurs at an oblique angle of incidence, yet the effect of impact angle on crater size or formation mechanism is only poorly understood. This is, in large part, due to the difficulty of inferring impactor properties, such as size, velocity and trajectory, from observations of natural craters, and the expense and complexity of simulating oblique impacts using numerical models. Laboratory oblique impact experiments and previous numerical models have shown that the portion of the projectile’s kinetic energy that is involved in crater excavation decreases significantly with impact angle. However, a thorough quantification of planetary-scale oblique impact cratering does not exist and the effect of impact angle on crater size is not considered by current scaling laws. To address this gap in understanding, we developed iSALE-3D, a three-dimensional multi-rheology hydrocode, which is efficient enough to perform a large number of well-resolved oblique impact simulations within a reasonable time. Here we present the results of a comprehensive numerical study containing more than 200 three-dimensional hydrocode-simulations covering a broad range of projectile sizes, impact angles and friction coefficients. We show that existing scaling laws in principle describe oblique planetary-scale impact events at angles greater than 30° measured from horizontal. The displaced mass of a crater decreases with impact angle in a sinusoidal manner. However, our results indicate that the assumption that crater size scales with the vertical component of the impact velocity does not hold for materials with a friction coefficient significantly lower than 0.7 (sand). We found that increasing coefficients of friction result in smaller craters and a formation process more controlled by impactor momentum than by energy.     

Electron impact excitation of the ground-state 3P fine-structure levels in atomic oxygen

Fine-structure collision strengths are calculated for transitions between the ground-state levels of atomic oxygen. The R -matrix method is used in which the (2p⁴) ³P, ¹D and ¹S terms are included as well as three pseudo-states chosen to represent the dipole polarizability of the ³P ground state. Very sophisticated configuration–interaction wavefunctions are used for the target states and a recoupling transformation to pair coupling is applied to the inner-region Hamiltonian matrix, with the ³P fine-structure energy levels being adjusted to match the observed splitting prior to diagonalization. Effective collision strengths are obtained by integrating over a Maxwellian distribution and the rate coefficient of the cooling of electron gas is determined. The cooling rates are significantly lower than those currently available, in confirmation of the result deduced from measurements of the Earth's ionosphere electron gas.     

Time-resolved studies of hypervelocity vertical impacts into porous particulate targets: Effects of projectile density on early-time coupling and crater growth

The depth and duration of energy and momentum coupling in an impact shapes the formation of the crater. The earliest stages of crater growth (when the projectile transfers its energy and momentum to the target) are unrecoverable when the event is described by late stage parameters, which collapse the initial conditions of the impact into a singular point in time and space. During the coupling phase, the details of the impact are mapped into the ejecta flow field. In this experimental study, we present new experimental and computational measurements of the ejecta distribution and crater growth extending from early times into main-stage ballistic flow for hypervelocity impacts over a range of projectile densities. Specifically, we assess the effect of projectile density on coupling depth and location in porous particulate (sand) targets. A non-invasive high-speed imaging technique is employed to capture the velocity of individual ejecta particles very early in the cratering event as a function of both time and launch position. These data reveal that the effects of early-stage coupling, such as non-constant ejection angles, manifest not only in early-time behavior but also extend to main-stage crater growth. Time-resolved comparisons with hydrocode calculations provide both benchmarking and insight into the parameters controlling the ejection process. Measurements of the launch position and metrics for the transient diameter to depth ratio as a function of time demonstrate non-proportional crater growth throughout much of excavation. Low-density projectiles couple closer to the surface, thereby leading to lower ejection angles and larger effective diameter to depth ratios. These results have implications for the ballistic emplacement of ejecta on planetary surfaces, and are essential to interpreting temporally resolved data from impact missions.     

Non-intrusive measurements of crater growth

An experimental technique to measure crater growth is presented whereby a high speed video captures profiles of a crater forming after impact obtained using a vertical laser sheet centered on the impact point. Unlike previous so called “quarter-space experiments,” where projectiles were launched along a transparent Plexiglas sheet so that growth of half a crater could be viewed, the use of the laser sheet permits viewing changes in crater shape without any physical interference to the cratering process. This technique indicates that for low velocity impacts (<300 m/s) into 220 μm glass beads that are without cohesion and where the projectile is not disrupted, craters initially grow somewhat proportionally, but that later their depths remain essentially constant while their diameters continue to expand. In addition, these experiments indicate that as the impact velocity increases, the rate of growth and the transient depth to diameter ratio at the end of ejecta excavation decreases. These last two observations are probably due to the large time of penetration of the projectile, which becomes a significant fraction of the time of crater formation. This is contrary to the expectations for the scaling rules, which assumes a point source. Very high curtain angles (>45°) are also seen, and could be due to the low friction angle of the target. Significant crater modification, which is rarely seen in “quarter-space experiments,” is also observed and appears to be controlled by the dynamic angle of repose of the target. These latter observations indicate that differences in target friction angles may need to be considered when determining near rim ejecta-mass distributions and large-scale crater modification processes on the planets.     

Early-stage ejecta velocity distribution for vertical hypervelocity impacts into sand

The ejecta dynamics during main-stage excavation flow in a cratering event have previously been well characterized, particularly for vertical impacts. In this experimental study, we present new results addressing the early-time, low-angle, high-speed component of the ejecta velocity distribution as a function of time for hypervelocity vertical impacts into sand. Although this regime represents a very small portion of total ejected mass in laboratory experiments, it comprises a greater percentage of growth for larger craters.     

Observations of martian gullies and constraints on potential formation mechanisms

The discovery of presumably geologically recent gully features on Mars (Malin and Edgett, 2000, Science 288, 2330-2335) has spawned a wide variety of proposed theories of their origin including hypotheses of the type of erosive material. To test the validity of gully formation mechanisms, data from the Mars Global Surveyor spacecraft has been analyzed to uncover trends in the dimensional and physical properties of the gullies and their surrounding terrain. We located 106 Mars Orbiter Camera (MOC) images that contain clear evidence of gully landforms, distributed in the southern mid and high latitudes, and analyzed these images in combination with Mars Orbiter Laser Altimeter (MOLA) and Thermal Emission Spectrometer (TES) data to provide quantitative measurements of numerous gully characteristics. Parameters we measured include apparent source depth and distribution, vertical and horizontal dimensions, slopes, orientations, and present-day characteristics that affect local ground temperatures. We find that the number of gully systems normalized to the number of MOC images steadily declines as one moves poleward of 30° S, reaches a minimum value between 60°-63° S, and then again rises poleward of 63° S. All gully alcove heads occur within the upper one-third of the slope encompassing the gully and the alcove bases occur within the upper two-thirds of the slope. Also, the gully alcove heads occur typically within the first 200 meters of the overlying ridge with the exception of gullies equatorward of 40° S where some alcove heads reach a maximum depth of 1000 meters. While gullies exhibit complex slope orientation trends, gullies are found on all slope orientations at all the latitudes studied. Assuming thermal conductivities derived from TES measurements as well as modeled surface temperatures, we find that 79% of the gully alcove bases lie at depths where subsurface temperatures are greater than 273 K and 21% of the alcove bases lie within the solid water regime. Most of the gully alcoves lie outside the temperature-pressure phase stability of liquid CO₂. Based on a comparison of measured gully features with predictions from the various models of gully formation, we find that models involving carbon dioxide, melting ground ice in the upper few meters of the soil, dry landslide, and surface snowmelt are the least likely to describe the formation of the martian gullies. Although some discrepancies still exist between prediction and observation, the shallow and deep aquifer models remain as the most plausible theories. Interior processes involving subsurface fluid sources are generally favored over exogenic processes such as wind and snowfall for explaining the origin of the martian gullies.     

HiRISE views enigmatic deposits in the Sirenum Fossae region of Mars

HiRISE images together with other recent orbital data from Mars define new characteristics of enigmatic Hesperian-aged deposits in Sirenum Fossae that are mostly 100-200 m thick, drape kilometers of relief, and often display generally low relief surfaces. New characteristics of the deposits, previously mapped as the “Electris deposits,” include local detection of meter-scale beds that show truncating relationships, a generally light-toned nature, and a variably blocky, weakly indurated appearance. Boulders shed by erosion of the deposits are readily broken down and contribute little to talus. Thermal inertia values for the deposits are ∼200 J m⁻² K⁻¹ s^−1/2 and they may incorporate hydrated minerals derived from weathering of basalt. The deposits do not contain anomalous amounts of water or water ice. Deflation may dominate degradation of the deposits over time and points to an inventory of fine-grained sediment. Together with constraints imposed by the regional setting on formation processes, these newly resolved characteristics are most consistent with an eolian origin as a loess-like deposit comprised of redistributed and somewhat altered volcanic ash. Constituent sediments may be derived from airfall ash deposits in the Tharsis region. An origin directly related to airfall ash or similar volcanic materials is less probable and emplacement by alluvial/fluvial, impact, lacustrine, or relict polar processes is even less likely.     

Secondary chaotic terrain formation in the higher outflow channels of southern circum-Chryse, Mars

Higher outflow channel dissection in the martian region of southern circum-Chryse appears to have extended from the Late Hesperian to the Middle Amazonian Epoch. These outflow channels were excavated within the upper 1 km of the cryolithosphere, where no liquid water is expected to have existed during these geologic epochs. In accordance with previous work, our examination of outflow channel floor morphologies suggests the upper crust excavated by the studied outflow channels consisted of a thin (a few tens of meters) layer of dry geologic materials overlying an indurated zone that extends to the bases of the investigated outflow channels (1 km in depth). We find that the floors of these outflow channels contain widespread secondary chaotic terrains (i.e., chaotic terrains produced by the destruction of channel-floor materials). These chaotic terrains occur within the full range of outflow channel dissection and tend to form clusters. Our examination of the geology of these chaotic terrains suggests that their formation did not result in the generation of floods. Nevertheless, despite their much smaller dimensions, these chaotic terrains are comprised of the same basic morphologic elements (e.g., mesas, knobs, and smooth deposits within scarp-bound depressions) as those located in the initiation zones of the outflow channels, which suggests that their formation must have involved the release of ground volatiles. We propose that these chaotic terrains developed not catastrophically but gradually and during multiple episodes of nested surface collapse. In order to explain the formation of secondary chaotic terrains within zones of outflow channel dissection, we propose that the regional Martian cryolithosphere contained widespread lenses of volatiles in liquid form. In this model, channel floor collapse and secondary chaotic terrain formation would have taken place as a consequence of instabilities arising during their exhumation by outflow channel dissection. Within relatively warm upper crustal materials in volcanic settings, or within highly saline crustal materials where cryopegs developed, lenses of volatiles in liquid form within the cryolithosphere could have formed, and/or remained stable.In addition, our numerical simulations suggest that low thermal conductivity, dry fine-grained porous geologic materials just a few tens of meters in thickness (e.g., dunes, sand sheets, some types of regolith materials), could have produced high thermal anomalies resulting in subsurface melting. The existence of a global layer of dry geologic materials overlying the cryolithosphere would suggest that widespread lenses of fluids existed (and may still exist) at shallow depths wherever these materials are fine-grained and porous. The surface ages of the investigated outflow channels and chaotic terrains span a full 500 to 700 Myr. Chaotic terrains similar in dimensions and morphology to secondary chaotic terrains are not observed conspicuously throughout the surface of Mars, suggesting that intra-cryolithospheric fluid lenses may form relatively stable systems. The existence of widespread groundwater lenses at shallow depths of burial has tremendous implications for exobiological studies and future human exploration. We find that the clear geomorphologic anomaly that the chaotic terrains and outflow channels of southern Chryse form within the Martian landscape could have been a consequence of large-scale resurfacing resulting from anomalously extensive subsurface melt in this region of the planet produced by high concentrations of salts within the regional upper crust. Crater count statistics reveal that secondary chaotic terrains and the outflow channels within which they occur have overlapping ages, suggesting that the instabilities leading to their formation rapidly dissipated, perhaps as the thickness of the cryolithosphere was reset following the disruption of the upper crustal thermal structure produced during outflow channel excavation.     

Evidence of gully formation by regional groundwater flow in the Gorgonum-Newton region (Mars)

The discovery of gullies on Mars suggests liquid water activity near the surface of the planet in recent times. Since liquid water is unstable under the present-day P-T martian conditions, the formation mechanisms of the gullies, and the source of the putative water, have been a matter of debate for the last years. To provide new insights into these matters, we have approached the problem studying the gullies in relation to their regional setting. A major point in our study relates to the geographic orientation of gullies, an aspect that has been previously regarded as a crucial matter in different models, and has profound implications regarding their origin. We present a comprehensive and detailed survey of the Gorgonum-Newton region, and a study of the Dao and Nirgal Vallis regions. The survey was carried out with the aid of 965 high-resolution MOC images (752 for Gorgonum-Newton, 102 for Nirgal Vallis and 111 for Dao Vallis regions), and MOLA-derived DEMs. We found that gullies display a clear regional pattern, geographically and topographically consistent with a decreasing regional slope. We interpret the results in terms of the existence of several groundwater flow systems operating at different scales, which ultimately may have led to gully formation by seepage at the slopes of craters and canyons. We suggest that aquifers discharging at gully systems may have recharged from the surface, in response to the melting of young partially eroded ice-rich deposits.     

Segmented lineaments on Europa: Implications for the formation of ridge complexes and bright bands

We describe several segmented lineaments on Europa’s surface. These lineaments are extensive, stretching for 100s-1000s of km, and have ridge complex or bright band morphologies. The geometries of the segmented portions of these features are diagnostic of the remote normal and shear stress environment in which they formed and, therefore, constrain ridge complex and bright band formation mechanisms. Analysis of four ridge complexes indicates that they formed in a remote normal stress environment that was tensile and isotropic (or nearly so) and that these lineaments may have formed in a manner more analogous to bands on Europa than to ridges. The stress environment associated with these ridge complexes may also explain the anastomosing nature of their interior morphology. Analysis of two bright bands indicate that one formed in a remote normal stress environment that was tensile and the other was reactivated under a combination of remote tensile normal stress and remote sinistral shear stress. Aspects of the morphologies of these features also indicate that bright bands likely have complex deformation histories that can include multiple episodes of reactivation.     

Physical properties of Titan's surface at the Huygens landing site from the Surface Science Package Acoustic Properties sensor (API-S)

We present the results from the first sonar to be deployed outside of Earth, and the first active acoustic instrument on Titan, onboard the Huygens probe, and the implications of its data for the geomorphology and characteristics of the Huygens landing site. Signals were recorded from 90 m downwards until impact, with a maximum sensor footprint diameter at the ground of 39.2 m. Probe impact speed was measured to be 4.67 m s⁻¹. Derivation of terrain topography in a transect beneath the probe may indicate a ridge-trough terrain with an amplitude of about 1 m and a wavelength of about 10 m, although a flat surface is also consistent with the results. Modelling of the returned signal indicates that the surface acoustic properties at the landing site must be specular in nature, which may have two possible (not incompatible) causes—the surface may consist of sorted interlocking grains, smooth on the centimetre scale, which would imply either fluvial sorting or the infill of small particles interstitial to the larger particles (similar to a terrestrial playa). Alternatively, specularity may indicate the presence of methane as an interstitial liquid or as very small pools. Due to mission constraints, tens of metres around the landing site were not well-imaged by Huygens' cameras except for the narrow azimuth observed after impact (the camera did not look straight down, and was not in imaging mode during the last few hundred metres of descent). Thus the data presented are among the few direct observations of the landing site surroundings.     

Observations of martian gullies and constraints on potential formation mechanisms: II. The northern hemisphere

The formation process(es) responsible for creating the observed geologically recent gully features on Mars has remained the subject of intense debate since their discovery. We present new data and analysis of northern hemisphere gullies from Mars Global Surveyor data which is used to test the various proposed mechanisms of gully formation. We located 137 Mars Orbiter Camera (MOC) images in the northern hemisphere that contain clear evidence of gully landforms and analyzed these images in combination with Mars Orbiter Laser Altimeter (MOLA) and Thermal Emission Spectrometer (TES) data to provide quantitative measurements of numerous gully characteristics. Parameters we measured include apparent source depth and distribution, vertical and horizontal dimensions, slopes, orientations, and present-day characteristics that affect local ground temperatures. Northern hemisphere gullies are clustered in Arcadia Planitia, Tempe Terra, Acidalia Planitia, and Utopia Planitia. These gullies form in craters (84%), knobby terrain (4%), valleys (3%), other/unknown terrains (9%) and are found on all slope orientations although the majority of gullies are equator-facing. Most gullies (63%) are associated with competent rock strata, 26% are not associated with strata, and 11% are ambiguous. Assuming thermal conductivities derived from TES measurements as well as modeled surface temperatures, we find that 95% of the gully alcove bases with adequate data coverage lie at depths where subsurface temperatures are greater than 273 K and 5% of the alcove bases lie within the solid water regime. The average alcove length is 470 m and the average channel length is 690 m. Based on a comparison of measured gully features with predictions from the various models of gully formation, we find that models involving carbon dioxide, melting ground ice in the upper few meters of the soil, dry landslide, and surface snowmelt are the least likely to describe the formation of the martian gullies. Although some discrepancies still exist between prediction and observation, the shallow and deep aquifer models remain as the most plausible theories. Interior processes involving subsurface fluid sources are generally favored over exogenic processes such as wind and snowfall for explaining the origin of the martian gullies. These findings gleaned from the northern hemisphere data are in general agreement with analyses of gullies in the southern hemisphere [Heldmann, J.L., Mellon, M.T., 2004. Icarus 168, 285-304].     

Possible ancient giant basin and related water enrichment in the Arabia Terra province, Mars

A circular albedo feature in the Arabia Terra province was first hypothesized as an ancient impact basin using Viking-era information. To test this unpublished hypothesis, we have analyzed the Viking era-information together with layers of new data derived from the Mars Global Surveyor (MGS) and Mars Odyssey (MO) missions. Our analysis indicates that Arabia Terra is an ancient geologic province of Mars with many distinct characteristics, including predominantly Noachian materials, a unique part of the highland-lowland boundary, a prominent paleotectonic history, the largest region of fretted terrain on the planet, outflow channels with no obvious origins, extensive exposures of eroded layered sedimentary deposits, and notable structural, albedo, thermal inertia, gravity, magnetic, and elemental signatures. The province also is marked by special impact crater morphologies, which suggest a persistent volatile-rich substrate. No one characteristic provides definitive answers to the dominant event(s) that shaped this unique province. Collectively the characteristics reported here support the following hypothesized sequence of events in Arabia Terra: (1) an enormous basin, possibly of impact origin, formed early in martian history when the magnetic dynamo was active and the lithosphere was relatively thin, (2) sediments and other materials were deposited in the basin during high erosion rates while maintaining isostatic equilibrium, (3) sediments became water enriched during the Noachian Period, and (4) basin materials were uplifted in response to the growth of the Tharsis Bulge, resulting in differential erosion exposing ancient stratigraphic sequences. Parts of the ancient basin remain water-enriched to the present day.