Numerical simulations of granular dynamics II: Particle dynamics in a shaken granular material

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties.We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration.We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.     

10Be content in clasts from fallout suevitic breccia in drill cores from the Bosumtwi impact crater,Ghana: Clues to preimpact target distribution

Rocks from drill cores LB‐07A (crater fill) and LB‐08A (central uplift) into the Bosumtwi impact crater, Ghana, were analyzed for the presence of the cosmogenic radionuclide ¹⁰Be. The aim of the study was to determine the extent to which target rocks of various depths were mixed during the formation of the crater‐filling breccia, and also to detect meteoric water infiltration within the impactite layer. ¹⁰Be abundances above background were found in two (out of 24) samples from the LB‐07A core, and in none of five samples from the LB‐08A core. After excluding other possible explanations for an elevated ¹⁰Be signal, we conclude that it is most probably due to a preimpact origin of those clasts from target rocks close to the surface. Our results suggest that in‐crater breccias were well mixed during the impact cratering process. In addition, the lack of a ¹⁰Be signal within the rocks located very close to the lake sediment–impactite boundary suggests that infiltration of meteoric water below the postimpact crater floor was limited. This may suggest that the infiltration of the meteoric water within the crater takes place not through the aerial pore‐space, but rather through a localized system of fractures.     

Clastic matrix in EH3 chondrites

Abstract— Patches of clastic matrix (15 to 730 μm in size) constitute 4.9 vol% of EH3 Yamato (Y‐) 691 and 11.7 vol% of EH3 Allan Hills (ALH) 81189. Individual patches in Y‐691 consist of 1) ?25 vol% relatively coarse opaque grain fragments and polycrystalline assemblages of kamacite, schreibersite, perryite, troilite (some grains with daubréelite exsolution lamellae), niningerite, oldhamite, and caswellsilverite; 2) ?30 vol% relatively coarse silicate grains including enstatite, albitic plagioclase, silica and diopside; and 3) an inferred fine nebular component (?45 vol%) comprised of submicrometer‐size grains. Clastic matrix patches in ALH 81189 contain relatively coarse grains of opaques (?20 vol%; kamacite, schreibersite, perryite and troilite) and silicates (?30 vol%; enstatite, silica and forsterite) as well as an inferred fine nebular component (?50 vol%). The O‐isotopic composition of clastic matrix in Y‐691 is indistinguishable from that of olivine and pyroxene grains in adjacent chondrules; both sets of objects lie on the terrestrial mass‐fractionation line on the standard three‐isotope graph. Some patches of fine‐grained matrix in Y‐691 have distinguishable bulk concentrations of Na and K, inferred to be inherited from the solar nebula. Some patches in ALH 81189 differ in their bulk concentrations of Ca, Cr, Mn, and Ni. The average compositions of matrix material in Y‐691 and ALH 81189 are similar but not identical‐matrix in ALH 81189 is much richer in Mn (0.23 ± 0.05 versus 0.07 ± 0.02 wt%) and appreciably richer in Ni (0.36 ± 0.10 versus 0.18 ± 0.05 wt%) than matrix in Y‐691. Each of the two whole‐rocks exhibits a petrofabric, probably produced by shock processes on their parent asteroid.     

Earth–Moon Weak Stability Boundaries in the restricted three and four body problem

This work deals with the structure of the lunar Weak Stability Boundaries (WSB) in the framework of the restricted three and four body problem. Geometry and properties of the escape trajectories have been studied by changing the spacecraft orbital parameters around the Moon. Results obtained using the algorithm definition of the WSB have been compared with an analytical approximation based on the value of the Jacobi constant. Planar and three-dimensional cases have been studied in both three and four body models and the effects on the WSB structure, due to the presence of the gravitational force of the Sun and the Moon orbital eccentricity, have been investigated. The study of the dynamical evolution of the spacecraft after lunar capture allowed us to find regions of the WSB corresponding to stable and safe orbits, that is orbits that will not impact onto lunar surface after capture. By using a bicircular four body model, then, it has been possible to study low-energy transfer trajectories and results are given in terms of eccentricity, pericenter altitude and inclination of the capture orbit. Equatorial and polar capture orbits have been compared and differences in terms of energy between these two kinds of orbits are shown. Finally, the knowledge of the WSB geometry permitted us to modify the design of the low-energy capture trajectories in order to reach stable capture, which allows orbit circularization using low-thrust propulsion systems.     

Lunar intrusive domes: Morphometric analysis and laccolith modelling

This study examines a set of lunar domes with very low flank slopes which differ in several respects from the frequently occurring lunar effusive domes. Some of these domes are exceptionally large, and most of them are associated with faults or linear rilles of presumably tensional origin. Accordingly, they might be interpreted as surface manifestations of laccolithic intrusions formed by flexure-induced vertical uplift of the lunar crust (or, alternatively, as low effusive edifices due to lava mantling of highland terrain, or kipukas, or structural features). All of them are situated near the borders of mare regions or in regions characterised by extensive effusive volcanic activity. Clementine multispectral UVVIS imagery indicates that they do not preferentially occur in specific types of mare basalt. Our determination of their morphometric properties, involving a combined photoclinometry and shape from shading technique applied to telescopic CCD images acquired at oblique illumination, reveals large dome diameters between 10 and more than 30 km, flank slopes below 0.9°, and volumes ranging from 0.5 to 50 km³. We establish three morphometric classes. The first class, In1, comprises large domes with diameters above 25 km and flank slopes of 0.2°-0.6°, class In2 is made up by smaller and slightly steeper domes with diameters of 10-15 km and flank slopes between 0.4° and 0.9°, and domes of class In3 have diameters of 13-20 km and flank slopes below 0.3°. While the morphometric properties of several candidate intrusive domes overlap with those of some classes of effusive domes, we show that a possible distinction criterion are the characteristic elongated outlines of the candidate intrusive domes. We examine how they differ from typical effusive domes of classes 5 and 6 defined by Head and Gifford [Head, J.W., Gifford, A., 1980. Lunar mare domes: classification and modes of origin. Moon Planets 22, 235-257], and show that they are likely no highland kipukas due to the absence of spectral contrast to their surrounding. These considerations serve as a motivation for an analysis of the candidate intrusive domes in terms of the laccolith model by Kerr and Pollard [Kerr, A.D., Pollard, D.D., 1998. Toward more realistic formulations for the analysis of laccoliths. J. Struct. Geol. 20(12), 1783-1793], to estimate the geophysical parameters, especially the intrusion depth and the magma pressure, which would result from the observed morphometric properties. Accordingly, domes of class In1 are characterised by intrusion depths of 2.3-3.5 km and magma pressures between 18 and 29 MPa. For the smaller and steeper domes of class In2 the magma intruded to shallow depths between 0.4 and 1.0 km while the inferred magma pressures range from 3 to 8 MPa. Class In3 domes are similar to those of class In1 with intrusion depths of 1.8-2.7 km and magma pressures of 15-23 MPa. As an extraordinary feature, we describe in some detail the concentric crater Archytas G associated with the intrusive dome Ar1 and discuss possible modes of origin. In comparison to the candidate intrusive domes, terrestrial laccoliths tend to be smaller, but it remains unclear if this observation is merely a selection effect due to the limited resolution of our telescopic CCD images. An elongated outline is common to many terrestrial laccoliths and the putative lunar laccoliths, while the thickness values measured for terrestrial laccoliths are typically higher than those inferred for lunar laccoliths, but the typical intrusion depths are comparable.     

Planar deformation features in quartz from impact‐produced polymict breccia of the Xiuyan crater,China

Abstract– The 1.8 km‐diameter Xiuyan crater is an impact structure in northeastern China, exposed in a Proterozoic metamorphic rock complex. The major rocks of the crater are composed of granulite, hornblendite, gneiss, tremolite marble, and marble. The bottom at the center of the crater covers about 100 m thick lacustrine sediments underlain by 188 m thick crater‐fill breccia. A layer of polymict breccia composed of clasts of granulite, gneiss, hornblendite, and fragments of glass as well as clastic matrix, occurs near the base, in the depth interval from 260 to 295 m. An investigation in quartz from the polymict breccia in the crater‐fill units reveals abundant planar deformation features (PDFs). Quartz with multiple sets of PDFs is found in clasts of granulite that consist of mainly quartz and feldspar, and in fine‐grained matrix of the impact‐produced polymict breccia. A universal stage was used to measure the orientation of PDFs in 70 grains of quartz from five thin sections made from the clasts of granulite of polymict breccia recovered at the depth of 290 m. Forty‐four percent of the quartz grains contain three sets of PDFs, and another 40% contain two sets of PDFs. The most abundant PDFs are rhombohedron forms of , , and with frequency of 33.5, 22.3, and 9.6%, respectively. A predominant PDF form of in quartz suggests a shock pressure >20 GPa. The occurrence of PDFs in quartz from the polymict breccia provides crucial evidence for shock metamorphism of target rocks and confirms the impact origin of this crater, which thus appears to be the first confirmed impact crater in China.     

Fourier Power Spectra of Solar Noise Storms

We analyzed three noise storms recorded on 200?–?400 MHz Trieste Callisto radio spectra on 2 July 2012, 8 July 2012, and 16 July 2012 by the Fourier method. We divided intervals of the noise storms into five-minute intervals, and in these intervals we computed the mean Fourier spectra as a function of the wave numbers in the frequency and height-scale spaces. We found that these Fourier spectra, where the spectrum from the quiet-activity interval was subtracted, are power-law spectra. The mean power-law index of these spectra in the range \(\ln(k_{z}) = [1.8, 2.9]\) (where \(k_{z}\) is the wave number in the height-scale space) is \(-1.7\pm0.14\), \(-1.6\pm0.14\), and \(-1.5 \pm0.12\) for the 2 July 2012, the 8 July 2012, and the 16 July 2012 noise storms, respectively. It appears that as the number of Type-I bursts in the studied interval increases, the power-law index becomes closer to \(-5/3\); this is known as the Kolmogorov spectral index. The power-law index of the noise storms is very similar to that of the narrowband dm-spikes found in our previous studies. Furthermore, we found a break in the power spectra at \(\ln(k_{z}) \approx2.9\), and the mean power-law index values above this break are \(-2.9\pm0.46\), \(-3.1\pm0.65\), and \(-3.4\pm0.98\), respectively.     

Chemical composition and iron oxidation state of amorphous matrix silicates in the carbonaceous chondrite Acfer 094

Nanoscale amorphous silicates are a major component in primitive carbonaceous chondrite matrices and anhydrous interplanetary dust particles. Owing to their metastability and sensitive response to reactions with water, this material is of particular interest in understanding nebular and parent body processes in the early solar system. Here we investigated the amorphous silicate matrix (ASM) in the ungrouped carbonaceous chondrite Acfer 094 regarding its texture, chemical composition, and Fe oxidation state. We applied transmission electron microscopy techniques on six, focused ion beam technique-prepared, electron-transparent lamellae of Acfer 094 to determine the textures of this material. Furthermore, we used energy-dispersive X-ray analysis and electron energy loss spectroscopy to quantify the Fe content and the Fe oxidation state of the ASM. Textural investigations reveal differences in sulfide content, porosity, and distribution of the ASM among the samples, as well as evidence for rare recrystallization of phyllosilicate fibers. The chemical composition reveals mobilization of Fe. Furthermore, the determined Fe³⁺/ΣFe ratios of the ASM in the six samples display a homogeneously high oxidation state (0.66–0.73). This high and homogeneous Fe oxidation state in the ASM of Acfer 094 disagrees with its formation as a primary phase in a reduced solar gas and must have been induced in a later stage process. Most likely, this process was aqueous alteration on the Acfer 094 parent body, which led to hydration and oxidation of the ASM, which is supported by textural and chemical evidence of aqueous alteration.     

Experimental impact cratering: A summary of the major results of the MEMIN research unit

This paper reviews major findings of the Multidisciplinary Experimental and Modeling Impact Crater Research Network (MEMIN). MEMIN is a consortium, funded from 2009 till 2017 by the German Research Foundation, and is aimed at investigating impact cratering processes by experimental and modeling approaches. The vision of this network has been to comprehensively quantify impact processes by conducting a strictly controlled experimental campaign at the laboratory scale, together with a multidisciplinary analytical approach. Central to MEMIN has been the use of powerful two-stage light-gas accelerators capable of producing impact craters in the decimeter size range in solid rocks that allowed detailed spatial analyses of petrophysical, structural, and geochemical changes in target rocks and ejecta. In addition, explosive setups, membrane-driven diamond anvil cells, as well as laser irradiation and split Hopkinson pressure bar technologies have been used to study the response of minerals and rocks to shock and dynamic loading as well as high-temperature conditions. We used Seeberger sandstone, Taunus quartzite, Carrara marble, and Weibern tuff as major target rock types. In concert with the experiments we conducted mesoscale numerical simulations of shock wave propagation in heterogeneous rocks resolving the complex response of grains and pores to compressive, shear, and tensile loading and macroscale modeling of crater formation and fracturing. Major results comprise (1) projectile–target interaction, (2) various aspects of shock metamorphism with special focus on low shock pressures and effects of target porosity and water saturation, (3) crater morphologies and cratering efficiencies in various nonporous and porous lithologies, (4) in situ target damage, (5) ejecta dynamics, and (6) geophysical survey of experimental craters.