Crater clusters on Mars: Shedding light on martian ejecta launch conditions

We have identified two classes of crater clusters on Mars. One class is “small clusters” (crater diameter D∼ tens m, spread over few hundred m), fitting our earlier calculations for the breakup of weak stone meteoroids in the martian atmosphere [Popova, O.P., Nemtchinov, I.V., Hartmann, W.K., 2003. Meteorit. Planet. Sci. 38, 905-925]. The second class is “large clusters” (D∼ few hundred m, spread over 2 to 30 km), which do not fit any predictions for breakup of known meteoroid types. We consider a range of possible explanations. The best explanation relates to known, high-speed ejection of large, semi-coherent, fractured rock masses from the surface, as secondary debris from primary impacts. The clusters are probably due to breakup of partly fracture, few-hundred-meter scale weak blocks, especially during ascent (producing moderate lateral spreading velocities among the fragments during sub-orbital flight), and also during descent of the resulting swarm. These conclusions illuminate the launch conditions of martian meteorites, including fragmentation processes, although more work is needed on the lateral separation of fragments (during either atmosphere descent or ascent) due to the effects of volatiles in the projectiles. Martian meteorites probably come from smaller martian craters than the clusters' source craters. The latter probably have D?85 km, although we have not ruled out diameters as small as 15 km.     

Numerical experiments with rubble piles: equilibrium shapes and spins

We present numerical experiments investigating the shape and spin limits of self-gravitating “perfect” rubble piles that consist of identical, smooth, rigid, spherical particles with configurable normal coefficient of restitution and no sliding friction. Such constructs are currently employed in a variety of investigations, ranging from the formation of asteroid satellites to the dynamical properties of Saturn's densest rings. We find that, owing to cannonball stacking behavior, rubble piles can maintain non-spherical shapes without bulk spin, unlike a fluid, and can spin faster than a perfect fluid before shedding mass, consistent with the theory for the more general continuum rubble pile model (Holsapple, 2004, Icarus 172, 272-303). Rubble piles that reassemble following a catastrophic disruption reconfigure themselves to lie within stability limits predicted by the continuum theory. We also find that coarse configurations consisting of a small number of particles are more resistant to tidal disruption than fine configurations with many particles. Overall this study shows that idealized rubble piles behave qualitatively in a manner similar to certain granular materials, at least in the limit where global shape readjustments and/or mass shedding begins. The limits obtained here may provide constraints on the possible internal structure of some small Solar System bodies that have extreme shapes or are under high stress. Amalthea is presented as a case study.     

Two chondrule groups each with distinctive rims in Murchison recognized by cathodoluminescence

Abstract— Two groups of chondrules in the Murchison CM chondrite, which have previously been identified on the basis of FeO in the chondrule grains, are readily identified from cathodoluminescence (CL) and belong to those of the ordinary chondrite group A and B chondrules of Sears et al. (1992a). All chondrules are surrounded by fine-grained rims containing forsterite with bright red CL, but on group A chondrules an outer thin rim grades into a much thicker rim, with a lower density of forsterite grains, which in turn grades into the central chondrule. Group B chondrules have only the thin outer rim with a high density of small forsterite grains. This is the first time an unequivocal correlation has been observed between chondrule rim thickness and the composition of the object on which the rim is located. We suggest that while all objects in the meteorite (group B chondrules, refractory inclusions, mineral and chondrule fragments, clasts) acquired a very thin rim during processing in a wet regolith, the thick rims on group A chondrules were formed by aqueous alteration of precursor metal- and sulfide-rich rims which are a characteristic of group A chondrules in ordinary chondrites.     

The formation of asteroid satellites in large impacts: results from numerical simulations

We present results of 161 numerical simulations of impacts into 100-km diameter asteroids, examining debris trajectories to search for the formation of bound satellite systems. Our simulations utilize a 3-dimensional smooth-particle hydrodynamics (SPH) code to model the impact between the colliding asteroids. The outcomes of the SPH models are handed off as the initial conditions for N-body simulations, which follow the trajectories of the ejecta fragments to search for the formation of satellite systems. Our results show that catastrophic and large-scale cratering collisions create numerous fragments whose trajectories can be changed by particle-particle interactions and by the reaccretion of material onto the remaining target body. Some impact debris can enter into orbit around the remaining target body, which is a gravitationally reaccreted rubble pile, to form a SMAshed Target Satellite (SMATS). Numerous smaller fragments escaping the largest remnant may have similar trajectories such that many become bound to one another, forming Escaping Ejecta Binaries (EEBs). Our simulations so far seem to be able to produce satellite systems qualitatively similar to observed systems in the main asteroid belt. We find that impacts of 34-km diameter projectiles striking at 3 km s⁻¹ at impact angles of ∼30° appear to be particularly efficient at producing relatively large satellites around the largest remnant as well as large numbers of modest-size binaries among their escaping ejecta.     

Quantitative risk assessment of slope hazards along a section of railway in the Canadian Cordillera—a methodology considering the uncertainty in the results

Railway alignments through the Canadian Cordillera are constantly exposed to slope instabilities. Proactive mitigation strategies have been in place for a few decades now, and instability record keeping has been recognized as an important aspect of them. Such a proactive strategy has enhanced the industry’s capacity to manage slope risks, and some sections have been recognized as critical due to the frequency of instabilities. At these locations, quantification of the risks becomes necessary. Risk analysis requires knowledge of some variables for which statistical data are scarce or not available, and elicitation of subjective probabilities is needed. A limitation of such approaches lies in the uncertainty associated to those elicited probabilities. In this paper, a quantitative risk analysis is presented for a section of railway across the Canadian Cordillera. The analysis focused on the risk to life of the freight train crews working along this section. Upper and lower bounds were elicited to cope with the uncertainties associated with this approach. A Monte Carlo simulation technique was then applied to obtain the probability distribution of the estimated risks. The risk probability distribution suggests that the risk to life of the crews is below previously published evaluation criteria and within acceptable levels. The risk assessment approach proposed focuses on providing a measure of the uncertainty associated with the estimated risk and is capable of handling distributions that cover more than two orders of magnitude.     

Application of the Christensen Failure Criterion to Intact Rock

The Christensen criterion, originally introduced in materials science, has a simple mathematical form and uniaxial tensile and compressive strength as the only parameters, making it an attractive candidate for rock engineering purposes. In this study, the applicability of the criterion to rock materials is examined. Explicit equations for application of the criterion under biaxial, triaxial compression, triaxial extension, and polyaxial states of stresses are derived. A comprehensive strength data set including the results of tests on synthetic rock, chert dyke, Carrara marble and Westerly granite is utilized to examine the accuracy of the Christensen criterion to the failure of rock material. The two surprising findings about the Christensen criterion are the zero values of tensile strength and the very low slopes of the failure envelope obtained from fitting analyses for chert dyke and Westerly granite. It is shown that the two problems are interrelated and the values of tensile strength tend to zero to produce higher slopes. It is then mathematically proven that the maximum initial slope of the Christensen failure envelope is limited to 4 in triaxial compression and 2.5 in triaxial extension which is considerably lower than the slope of experimental data. The accuracy of the Christensen criterion was found to be significantly lower than the well-established Hoek–Brown criterion. The circular π-plane representations and brittle-to-ductile transition limits from the Christensen criterion are also inconsistent with the observed behavior of rocks.     

Magmatism at continental passive margins inferred from Ambient‐Noise Phase‐velocity in the Gulf of Aden

Non‐volcanic continental passive margins have traditionally been considered to be tectonically and magmatically inactive once continental breakup has occurred and seafloor spreading has commenced. We use ambient‐noise tomography to constrain Rayleigh‐wave phase‐velocity maps beneath the eastern Gulf of Aden (eastern Yemen and southern Oman). In the crust, we image low velocities beneath the Jiza‐Qamar (Yemen) and Ashawq‐Salalah (Oman) basins, likely caused by the presence of partial melt associated with magmatic plumbing systems beneath the rifted margin. Our results provide strong evidence that magma intrusion persists after breakup, modifying the composition and thermal structure of the continental margin. The coincidence between zones of crustal intrusion and steep gradients in lithospheric thinning, as well as with transform faults, suggests that magmatism post‐breakup may be driven by small‐scale convection and enhanced by edge‐driven flow at the juxtaposition of lithosphere of varying thickness and thermal age.     

Retreat dynamics of the eastern sector of the British–Irish Ice Sheet during the last glaciation

The findings of BRITICE-CHRONO Transect 2 through the North Sea Basin and eastern England are reported. We define ice-sheet marginal oscillation between ~31 and 16 ka, with seven distinctive former ice-sheet limits (L1–7) constrained by Bayesian statistical analysis. The southernmost limit of the North Sea Lobe is recorded by the Bolders Bank Formation (L1; 25.8–24.6 ka). L2 represents ice-sheet oscillation and early retreat to the northern edge of the Dogger Bank (23.5–22.2 ka), with the Garret Hill Moraine in north Norfolk recording a significant regional readvance to L3 at 21.5–20.8 ka. Ice-marginal oscillations at ~26–21 ka resulted in L1, L2 and L3 being partially to totally overprinted. Ice-dammed lakes related to L1–3, including Lake Humber, are dated at 24.1–22.3 ka. Ice-sheet oscillation and retreat from L4 to L5 occurred between 19.7 and 17.3 ka, with grounding zone wedges marking an important transition from terrestrial to marine tidewater conditions, triggered by the opening of the Dogger Lake spillway between 19.9 and 17.5 ka. L6 relates to ice retreat under glacimarine conditions and final ice retreat into the Firth of Forth by 15.8 ka. L7 (~15 ka) represents an ice retreat from Bosies Bank into the Moray Firth.