Modelling the behaviour of granular material on the surface of asteroids and under different gravity conditions (e.g., Mars, the Moon)

This thesis, at the interface between the scientific disciplines of planetary science and granular physics, has two key components, both of which intend to increase our understanding of granular dynamics in varying gravitational conditions. The dynamics of granular materials are involved in the evolution of solid planets and small bodies in our Solar System, whose surfaces are generally covered with regolith. Understanding granular dynamics is also critical for the design and/or operations of landers, sampling devices and rovers to be included in space missions. The first component of this thesis is the validation of the hard-sphere discrete element method implementation in the N-body code pkdgrav to model the dynamics of granular material. By direct comparison with results from laboratory experiments, it is demonstrated that the hard-sphere discrete element method implementation in pkdgrav is valid for modelling granular material in dilute regimes and is capable of reproducing the complex dynamical behaviour of a specific dense system as well. The second component is focussed on the AstEx parabolic flight experiment. This experiment, with the aim of characterising the response of granular material to rotational shear forces in a microgravity environment, was designed, constructed, flown and the data were analysed as part of this thesis. It was found that the effect of constant shearing on a granular material in a direction perpendicular to the gravity field is not strongly influenced by gravity. The AstEx experiment has demonstrated, for the first time, that the efficiency of granular convection may decrease in the presence of a weak gravitational field, similar to that on the surface of small bodies. The first measurements of transient weakening of granular material after shear reversal in microgravity are also reported. Results suggest that the force contact network may be weaker in microgravity, although the influence of any change in the contact network is felt by the granular material over much larger distances. This may have important implications for our interpretation of asteroid surfaces. Continued advancement of our understanding of granular materials in varying gravitational conditions requires futher experiments and the development of the soft-sphere discrete element method implementation in pkdgrav in order to model the granular regimes that are inaccessbile to the hard-sphere implementation

Micro-meteoroid seismic uplift and regolith concentration on kilometric scale asteroids

Seismic shaking is an attractive mechanism to explain the destabilisation of regolith slopes and the regolith migration found on the surfaces of asteroids (Richardson et al. 2004; Miyamoto et al. 2007). Here, we use a continuum mechanics method to simulate the seismic wave propagation in an asteroid. Assuming that asteroids can be described by a cohesive core surrounded by a thin non-cohesive regolith layer, our numerical simulations of vibrations induced by micro-meteoroids suggest that the surface peak ground accelerations induced by micro-meteoroid impacts may have been previously under-estimated. Our lower bound estimate of vertical accelerations induced by seismic waves is about 50 times larger than previous estimates. It suggests that impact events triggering seismic activity are more frequent than previously assumed for asteroids in the kilometric and sub-kilometric size range. The regolith lofting is also estimated by a first order ballistic approximation. Vertical displacements are small, but lofting times are long compared to the duration of the seismic signals. The regolith movement has a non-linear dependence on the distance to the impact source which is induced by the type of seismic wave generating the first movement. The implications of regolith concentration in lows of surface acceleration potential are also discussed. We suggest that the resulting surface thermal inertia variations of small fast rotators may induce an increased sensitivity of these objects to the Yarkovsky effect.Comment: Accepted for publication in Icaru

Astex microgravity experiment: simulated asteroid regoliths

AstEx is a microgravity experiment selected to fly on ESA's 51st Microgravity Research Campaign in November 2009. The experiment will investigate the dynamics of regolith on asteroid surfaces. Despite their very low surface gravities, asteroids exhibit a number of different geological processes involving granular matter. Understanding the mechanical response of this granular material subject to external forces in microgravity conditions is vital to the design of a successful asteroid sub-surface sampling mechanism, and in the interpretation of the fascinating geology on an asteroid. The AstEx experiment uses a microgravity modified Taylor-Couette shear cell to investigate granular flow caused by shear forces under the conditions of parabolic flight microgravity. It is intended to determine how a steady state granular flow is achieved in microgravity conditions, and what effect prior shear history has on the timescales involved in initiating a steady state flow in a granular material. Presented are the technical details of the AstEx experimental design with particular emphasis on how the team have designed the equipment specifically for the parabolic flight microgravity environment

Seismometer Detection of Dust Devil Vortices by Ground Tilt

We report seismic signals on a desert playa caused by convective vortices and dust devils. The long-period (10-100s) signatures, with tilts of ~10$^{-7}$ radians, are correlated with the presence of vortices, detected with nearby sensors as sharp temporary pressure drops (0.2-1 mbar) and solar obscuration by dust. We show that the shape and amplitude of the signals, manifesting primarily as horizontal accelerations, can be modeled approximately with a simple quasi-static point-load model of the negative pressure field associated with the vortices acting on the ground as an elastic half space. We suggest the load imposed by a dust devil of diameter D and core pressure {\Delta}Po is ~({\pi}/2){\Delta}PoD$^2$, or for a typical terrestrial devil of 5 m diameter and 2 mbar, about the weight of a small car. The tilt depends on the inverse square of distance, and on the elastic properties of the ground, and the large signals we observe are in part due to the relatively soft playa sediment and the shallow installation of the instrument. Ground tilt may be a particularly sensitive means of detecting dust devils. The simple point-load model fails for large dust devils at short ranges, but more elaborate models incorporating the work of Sorrells (1971) may explain some of the more complex features in such cases, taking the vortex winds and ground velocity into account. We discuss some implications for the InSight mission to Mars.Comment: Contributed Article for Bulletin of the Seismological Society of America, Accepted 29th August 201

Estimations of the Seismic Pressure Noise on Mars Determined from Large Eddy Simulations and Demonstration of Pressure Decorrelation Techniques for the Insight Mission

The atmospheric pressure fluctuations on Mars induce an elastic response in the ground that creates a ground tilt, detectable as a seismic signal on the InSight seismometer SEIS. The seismic pressure noise is modeled using Large Eddy Simulations (LES) of the wind and surface pressure at the InSight landing site and a Green’s function ground deformation approach that is subsequently validated via a detailed comparison with two other methods: a spectral approach, and an approach based on Sorrells’ theory (Sorrells,Geophys. J. Int. 26:71–82, 1971; Sorrells et al., Nat. Phys. Sci. 229:14–16, 1971). The horizontal accelerations as a result of the ground tilt due to the LES turbulence-induced pressure fluctuations are found to be typically ∼ 2–40 nm/s2 in amplitude, whereas the direct horizontal acceleration is two orders of magnitude smaller and is thus negligible in comparison. The vertical accelerations are found to be ∼ 0.1–6 nm/s2 in amplitude. These are expected to be worst-case estimates for the seismic noise as we use a half-space approximation; the presence at some (shallow) depth of a harder layer would significantly reduce quasi-static displacement and tilt effects. We show that under calm conditions, a single-pressure measurement is representative of the large-scale pressure field (to a distance of several kilometers), particularly in the prevailing wind direction. However, during windy conditions, small-scale turbulence results in a reduced correlation between the pressure signals, and the single-pressure measurement becomes less representative of the pressure field. The correlation between the seismic signal and the pressure signal is found to be higher for the windiest period because the seismic pressure noise reflects the atmospheric structure close to the seismometer. In the same way that we reduce the atmospheric seismic signal by making use of a pressure sensor that is part of the InSight Auxiliary Payload Sensor Suite, we also the use the synthetic noise data obtained from the LES pressure field to demonstrate a decorrelation strategy. We show that our decorrelation approach is efficient, resulting in a reduction by a factor of ∼ 5 in the observed horizontal tilt noise (in the wind direction) and the vertical noise. This technique can, therefore, be used to remove the pressure signal from the seismic data obtained on Mars during the InSight mission

An experimental study of low-velocity impacts into granular material in reduced gravity

In order to improve our understanding of landing on small bodies and of asteroid evolution, we use our novel drop tower facility \citep{sunday2016} to perform low-velocity (2 - 40 cm/s), shallow impact experiments of a 10 cm diameter aluminum sphere into quartz sand in low effective gravities (~0.2 - 1 m/s^2). Using in-situ accelerometers we measure the acceleration profile during the impacts and determine the peak accelerations, collision durations and maximum penetration depth. We find that the penetration depth scales linearly with the collision velocity but is independent of the effective gravity for the experimental range tested, and that the collision duration is independent of both the effective gravity and the collision velocity. No rebounds are observed in any of the experiments. Our low-gravity experimental results indicate that the transition from the quasi-static regime to the inertial regime occurs for impact energies two orders of magnitude smaller than in similar impact experiments under terrestrial gravity. The lower energy regime change may be due to the increased hydrodynamic drag of the surface material in our experiments, but may also support the notion that the quasi-static regime reduces as the effective gravity becomes lower

The Noise Model of the SEIS Seismometer of the InSight Mission to Mars

The SEIS (Seismic Experiment for Interior Structures) instrument on board the InSight mission to Mars is the critical instrument for determining the interior structure of Mars, the current level of tectonic activity and the meteorite flux. Meeting the performance requirements of the SEIS instrument is vital to successfully achieve these mission objectives. The InSight noise model is a key tool for the InSight mission and SEIS instrument requirement setup. It will also be used for future operation planning. This paper presents the analyses made to build a model of the Martian seismic noise as measured by the SEIS seismometer, around the seismic bandwidth of the instrument (from 0.01 Hz to 1 Hz). It includes the instrument self-noise, but also the environment parameters that impact the measurements. We present the general approach for the model determination, the environment assumptions, and we analyze the major and minor contributors to the noise model

Evaluating the Wind-Induced Mechanical Noise on the InSight Seismometers

The SEIS (Seismic Experiment for Interior Structures) instrument onboard the InSight mission to Mars is the critical instrument for determining the interior structure of Mars, the current level of tectonic activity and the meteorite flux. Meeting the performance requirements of the SEIS instrument is vital to successfully achieve these mission objectives. Here we analyse in-situ wind measurements from previous Mars space missions to understand the wind environment that we are likely to encounter on Mars, and then we use an elastic ground deformation model to evaluate the mechanical noise contributions on the SEIS instrument due to the interaction between the Martian winds and the InSight lander. Lander mechanical noise maps that will be used to select the best deployment site for SEIS once the InSight lander arrives on Mars are also presented. We find the lander mechanical noise may be a detectable signal on the InSight seismometers. However, for the baseline SEIS deployment position, the noise is expected to be below the total noise requirement > 97 % of the time and is, therefore, not expected to endanger the InSight mission objectives

Seismic Coupling of Short-Period Wind Noise Through Mars’ Regolith for NASA’s InSight Lander

NASA’s InSight lander will deploy a tripod-mounted seismometer package onto the surface of Mars in late 2018. Mars is expected to have lower seismic activity than the Earth, so minimisation of environmental seismic noise will be critical for maximising observations of seismicity and scientific return from the mission. Therefore, the seismometers will be protected by a Wind and Thermal Shield (WTS), also mounted on a tripod. Nevertheless, wind impinging on the WTS will cause vibration noise, which will be transmitted to the seismometers through the regolith (soil). Here we use a 1:1-scale model of the seismometer and WTS, combined with field testing at two analogue sites in Iceland, to determine the transfer coefficient between the two tripods and quantify the proportion of WTS vibration noise transmitted through the regolith to the seismometers. The analogue sites had median grain sizes in the range 0.3–1.0 mm, surface densities of 1.3-1.8 gcm−3, and an effective regolith Young’s modulus of 2.5+1.9−1.4 MPa. At a seismic frequency of 5 Hz the measured transfer coefficients had values of 0.02–0.04 for the vertical component and 0.01–0.02 for the horizontal component. These values are 3–6 times lower than predicted by elastic theory and imply that at short periods the regolith displays significant anelastic behaviour. This will result in reduced short-period wind noise and increased signal-to-noise. We predict the noise induced by turbulent aerodynamic lift on the WTS at 5 Hz to be ∼2×10−10 ms−2 Hz−1/2 with a factor of 10 uncertainty. This is at least an order of magnitude lower than the InSight short-period seismometer noise floor of 10−8 ms−2 Hz−1/2