Nature of the Martian uplands: Effect on Martian meteorite age distribution and secondary cratering

Abstract— Martian meteorites (MMs) have been launched from an estimated 5–9 sites on Mars within the last 20 Myr. Some 80–89% of these launch sites sampled igneous rock formations from only the last 29% of Martian time. We hypothesize that this imbalance arises not merely from poor statistics, but because the launch processes are dominated by two main phenomena: first, much of the older Martian surface is inefficient in launching rocks during impacts, and second, the volumetrically enormous reservoir of original cumulate crust enhances launch probability for 4.5 Gyr old rocks. There are four lines of evidence for the first point, not all of equal strength. First, impact theory implies that MM launch is favored by surface exposures of near‐surface coherent rock (≤10² m deep), whereas Noachian surfaces generally should have ≥10² m of loose or weakly cemented regolith with high ice content, reducing efficiency of rock launch. Second, similarly, both Mars Exploration Rovers found sedimentary strata, 1–2 orders of magnitude weaker than Martian igneous rocks, favoring low launch efficiency among some fluvial‐derived Hesperian and Noachian rocks. Even if launched, such rocks may be unrecognized as meteorites on Earth. Third, statistics of MM formation age versus cosmic‐ray exposure (CRE) age weakly suggest that older surfaces may need larger, deeper craters to launch rocks. Fourth, in direct confirmation, one of us (N. G. B.) has found that older surfaces need larger craters to produce secondary impact crater fields (cf. Barlow and Block 2004). In a survey of 200 craters, the smallest Noachian, Hesperian, and Amazonian craters with prominent fields of secondaries have diameters of ?45 km, ?19 km, and ?10 km, respectively. Because 40% of Mars is Noachian, and 74% is either Noachian or Hesperian, the subsurface geologic characteristics of the older areas probably affect statistics of recognized MMs and production rates of secondary crater populations, and the MM and secondary crater statistics may give us clues to those properties.     

Planetary cratering 1. The question of multiple impactor populations: Lunar evidence

Abstract— This paper addresses several current issues related to use of craters in interpreting planetary surface histories. The primary goal is to test the widely adopted hypothesis of multiple populations of impactors at different times or places in the Solar System. New data presented here revise a “lunar highland” crater diameter distribution that has been widely used as evidence of an early distinct population of impactors. This curve, which has a depression of the size distribution at mid-sizes, does not, in fact, represent the lunar highlands generally. I show that it is associated with regions of intercrater plains. The more extensive the obliteration by intercrater plains, the deeper the depression. Modeling indicates that the depression of the curve is caused by the obliteration process itself. The oldest, most cratered regions of lunar highlands do not show the depression. These findings call into question earlier interpretations of multiple populations of impactors in the Solar System and of a distinctive primordial population. The present work is consistent, instead, with (1) a relatively uniform size distribution of interplanetary impactors, of mixed origins, back to 4 Ga ago and throughout the sampled Solar System; (2) fragmentation as the process that produced that size distribution; (3) saturation equilibrium on the most heavily cratered surfaces; and (4) differences in structure in the size distribution caused not by distinct impactor populations but by episodes of endogenic obliteration. If accepted, these results would modify some studies of solar system evolution, including assertions of two to five distinct populations of impactors, assumptions of lack of saturation equilibrium, and identifications of specific heliocentric or planetocentric sources for impactors within outer planet satellite systems.     

The structural state of the former lower continental crust in Calabria (S. Italy)

In the Northern Serre (S. Italy), a continuous section through a former lower continental crust is exposed (Schenk 1984, 1985). These granulite facies rocks have not been penetratively deformed during their uplift to higher crustal levels. Therefore, they allow to study the structural state of the lower continental crust. Most rockforming minerals, among them ortho- and clinopyroxene and cordierite, show recrystallization structures. Quarz basal subgrain boundaries are developed in lower levels of the lower crust and disappear toward higher levels. In addition to dimensional orientations of the rockforming minerals, the preferred crystallographic orientations of plagioclases, cordierite, amphibole, pyroxenes and, rarely, quartz reflect the strain orientation of the prograde high-T deformation. During this deformation, the stretching direction was constant NNE-SSW throughout the exposed lower crustal section. Shear indicators are rare. This points to the fact that deformation has been homogeneous and coaxial flattening. There is only a small variation in orientation and strength of textures throughout the exposed lower crust. Therefore, neighbouring layers of different composition rather than different deformation might serve as seismic reflectors in the former lower crust.

---

Zusammenfassung In der nördlichen Serre (Süditalien) ist ein zusammenhängender Bereich einer ehemaligen kontinentalen Unterkruste aufgeschlossen (Schenk 1984, 1985). Diese granulitfaziellen Gesteine sind bei ihrem Aufstieg in höhere Krustenbereiche nicht durchgreifend überprägt worden. Deshalb läßt sich an ihnen noch der strukturelle Zustand der ehemaligen Unterkruste studieren. Außer in Feldspäten, Quarz und Amphibol tritt auch in Ortho- und Klinopyroxen und in Cordierit Rekristallisation auf. In tieferen Krustenniveaus sind im Quarz basisparallele Subkorngrenzen entwickelt, die in höheren Krustenbereichen nicht auftreten. Sowohl die Gestaltregelungen der gesteinsbildenden Minerale als auch die kristallographischen Plagioklas-, Cordierit-, Amphibol-, Pyroxen- und (seltener) Quarzregelungen spiegeln die Strain-orientierung der prograden hochtemperierten Deformation wider. Während dieser Deformation blieb die Streckung im gesamten Unterkrustenbereich konstant NNE-SSW orientiert. Scherindikatoren sind selten. Beides deutet darauf hin, daß die Deformation homogen und als koachsiale Plättung abgelaufen ist. Im gesamten aufgeschlossenen Unterkrustenbereich variieren weder Intensität noch Orientierung der kristallographischen und gestaltlichen Regelungen der gesteinsbildenden Minerale wesentlich. Deshalb sollten in der ehemaligen Unterkruste Lagen unterschiedlicher Deformation weniger als seismische Reflektoren in Frage kommen als Lagen unterschiedlicher Zusammensetzung.

---

Résumé Le Massif de la Serre Septentrionale en Calabre (Sud de l'Italie) représente un tronçon d'une ancienne croûte continentale profonde (Schenk 1984, 1985). Les roches du faciès granulite qui le constituent n'ont pas subi de déformation penetrative lors de leur ascension dans les niveaux plus superficiels de la croûte. Leur étude permet donc de reconstruire les caractéristiques structurales de l'ancienne croûte profonde. Non seulement les feldspaths, le quartz et l'amphibole, mais encore l'orthopyroxène, le clinopyroxène et la cordiérite ont recristallisé. Dans les niveaux les plus profonds de la croûte on observe dans les cristaux de quartz des bordures de sous-grains parallèles à la base. Celles-ci n'apparaissent plus dans les niveaux plus superficiels. L'orientation des cristaux des minéraux constitutifs des roches, de même que l'orientation des réseaux cristallins du plagioclase, de la cordiérite, de l'amphibole, du pyroxène et (mais plus rarement) du quartz réflètent la direction de contrainte de la phase de déformation prograde de haute température. Pendant cette phase de déformation, la direction d'extension est demeurée NNE-SSW. Les témoins de cisaillement sont rares, ce qui montre que la déformation s'est effectuée de façon homogène, en relation avec un aplatissement coaxial. Les textures varient très peu en orientation comme en intensité dans toute la coupe de croûte profonde exposée. En conséquence, l'existence de réflecteurs sismiques devrait refléter des différences dans les compositions minéralogiques des roches plutôt que dans leurs caractéristiques structurales.

---

Serre (Schenk, 1984, 1985). . . , , . , , . , , , , () . NNE-SSW. . , . , , . , , .

     

Remote comets and related bodies: VJHK colorimetry and surface materials

Spectrophotometric data show that major compositional groups among outer solar system (OSS) surfaces include bright ices and at least two distinct classes of blackish carbonaceous-like materials, called C-type and RD-type. VJHK colorimetry of asteroids, satellites, and laboratory samples shows that these three classes can be distinguished by VJHK colors. We define an “α index” that denotes the position of objects in VJHK color - color diagrams; it empirically increases with albedo and ice/dirt ratio. We use the above data to define color fields that may be useful in interpreting our observations of eight comets (1980–1981). All eight comets have colors generally resembling RD asteroids and are inconsistent with reflection off clean ice surfaces. The observations suggest that these comets' halos contain RD dirt or dirty ice grains colored by RD dirt, supporting J. Gradie and J. Veverka's [Nature283, 840–842 (1980)] prediction of RD, rather than C, material in comets. Remote Comet P/Schwassmann-Wachmann 1 was observed both during outburst and quiescence and had the highest α index of any observed comet. Comet α indices appear to be correlated with solar distance. Further work will be needed to clarify possible coloring effects due to particle size, dispersal, and composition. We suggest a number of physical interpretations based on a single two-component mixing model, which assumes that all OSS planetesimals formed primarily from bright ices and dark carboneceous-like dirt, consistent with condensation theory. We discuss differentiation processes that concentrated one component or the other at the surface. All measured OSS interplanetary bodies have surfaces of dark dirt or dark dirty ice colored by the dirt component. Comets, consistent with the Whipple dirty iceberg model, are such objects close enough to the Sun for volatilization to throw dirty ice grains into the coma. In remote comets, the ice component of the grains remains stable, and we see dirty ice grains; in near comets, the ice component vaporizes, and we see dirt grains. A volatile-depleted dusty regolith on P/Schwassmann-Wachmann 1 and other remote comets could explain their eruptive behavior by means of gas pressure buildup in the porous, weakly bonded dust.     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.     

Roll Vortices and Boundary-Layer Development during a Cold Air Outbreak

The development of the boundary layer during a cold air outbreak in the FramStrait is documented by aircraft measurements. The convection was organisedinto roll vortices with aspect ratios increasing from 2.9 near the ice edgeto more than 6 at 100 km further downstream. This increase coincides with anincrease of the latent heat release in the cloud layer. The stability parameter-zi/L varies from about zero at the ice edge to 30 at a distance of 200 kmdownstream over open water where the satellite picture still shows cloudstreets. The increase is mainly due to the deepening of the boundary layer.The turbulent vertical sensible and latent heat fluxes near the surface amountto 400 W m^-2 within a 300 km off-ice zone. 25% of the upward heat fluxin the subcloud layer is carried out by organised roll motions. Experimentswith a 2-dimensional non-hydrostatic model show a similar roll aspect ratio inthe first 50 km, but further downstream where condensational heating is moreimportant the modelled roll wavelengths are distinctly smaller than the observedones.     

Reconstruction of atmospheric CO2 from ice-core data and the deep-sea record of ontong Java plateau: the Milankovitch chron

We provide a reconstruction of atmospheric CO₂ from deep-sea sediments, for the past 625000 years (Milankovitch chron). Our database consists of a Milankovitch template of sea-level variation in combination with a unique data set for the deep-sea record for Ontong Java plateau in the western equatorial Pacific. We redate the Vostok ice-core data of Barnola et al. (1987). To make the reconstructions we employ multiple regression between deep-sea data, on one hand, and ice-core CO₂ data in Antarctica, on the other. The patterns of correlation suggest that the main factors controlling atmospheric CO₂ can be described as a combination of sea-level state and sea-level change. For best results squared values of state and change are used. The square-of-sea-level rule agrees with the concept that shelf processes are important modulators of atmospheric CO₂ (e.g., budgets of shelf organic carbon and shelf carbonate, nitrate reduction). The square-of-change rule implies that, on short timescales, any major disturbance of the system results in a temporary rise in atmospheric CO₂.     

Gravity reduction with three-dimensional atmospheric pressure data for precise ground gravity measurements

The redistribution of air masses induces gravity variations (atmospheric pressure effect) up to about 20 μgal. These variations are disturbing signals in gravity records and they must be removed very carefully for detecting weak gravity signals. In the past, different methods have been developed for modelling of the atmospheric pressure effect. These methods use local or two-dimensional (2D) surface atmospheric pressure data and a standard height-dependent air density distribution. The atmospheric pressure effect is consisting of the elastic deformation and attraction term. The deformation term can be well modelled with 2D surface atmospheric pressure data, for instance with the Green's function method. For modelling of the attraction term, three-dimensional (3D) data are required. Results with 2D data are insufficient.From European Centre for Medium-Range Weather Forecasts (ECMWF) 3D atmospheric pressure data are now available. The ECMWF data used here are characterised by a spacing of Δ and Δλ = 0.5°, 60 pressure levels up to a height of 60 km and an interval of 6 h. These data are used for modelling of the atmospheric attraction term. Two attraction models have been developed based on the point mass attraction of air segments and the gravity potential of the air masses. The modelling shows a surface pressure-independent part of gravity variations induced by mass redistribution of the atmosphere in the order of some μgal. This part can only be determined by using 3D atmospheric pressure data. It has been calculated for the Vienna Superconducting Gravimeter site.From this follows that the gravity reduction can be improved by applying the 3D atmospheric attraction model for analysing long-periodic tidal waves including the polar tide. The same improvement is expected for reduction of long-term absolute gravity measurements or comparison of gravity measurements at different seasonal times. By using 3D atmospheric pressure data, the gravity correction can be improved up to some μgal.     

A parameterisation for the threshold shear velocity to initiate deflation of dry and wet sediment

A new parameterisation for the threshold shear velocity to initiate deflation of dry and wet particles is presented. It is based on the balance of moments acting on particles at the instant of particle motion. The model hence includes a term for the aerodynamic forces, including the drag force, the lift force and the aerodynamic-moment force, and a term for the interparticle forces. The effect of gravitation is incorporated in both terms. Rather than using an implicit function for the effect of the aerodynamic forces as reported earlier in literature, a constant aerodynamic coefficient was introduced. From consideration of the van der Waals force between two particles, it was further shown that the effect of the interparticle cohesion force between two dry particles on the deflation threshold should be inversely proportional to the particle diameter squared. The interparticle force was further extended to include wet bonding forces. The latter were considered as the sum of capillary forces and adhesive forces. A model that expresses the capillary force as a function of particle diameter squared and the inverse of capillary potential was deduced from consideration of the well-known model of Fisher and the Young–Laplace equation. The adhesive force was assumed to be equal to tensile strength, and a function which is proportional to particle diameter squared and the inverse of the potential due to adhesive forces was derived. By combining the capillary-force model and the adhesive force model, the interparticle force due to wet bonding was simplified and written as a function of particle diameter squared and the inverse of matric potential. The latter was loglinearly related to the gravimetric moisture content, a relationship that is valid in the low-moisture content range that is important in the light of deflation of sediment by wind. By introducing a correction to force the relationship to converge to zero moisture content at oven dryness, the matric potential–moisture content relationship contained only one unknown model parameter, viz. moisture content at −1.5 MPa. Working out the model led to a rather simple parameterisation containing only three coefficients. Two parameters were incorporated in the term that applies to dry sediment and were determined by using experimental data as reported by Iversen and White [Sedimentology 29 (1982) 111]. The third parameter for the wet-sediment part of the model was determined from wind-tunnel experiments on prewetted sand and sandy loam aggregates. The model was validated using data from wind-tunnel experiments on the same but dry sediment, and on data obtained from simulations with the model of Chepil [Soil Sci. Soc. Am. Proc. 20 (1956) 288]. The experiments showed that soil aggregates should be treated as individual particles with a density equal to their bulk density. Furthermore, it was shown that the surface had to dry to a moisture content of about 75% of the moisture content at −1.5 MPa before deflation became sustained. The threshold shear velocities simulated with our model were found to be in good agreement with own observations and with simulations using Chepil's model.