No Arabic abstract
Gravity inversion allows us to constrain the interior mass distribution of a planetary body using the observed shape, rotation, and gravity. Traditionally, techniques developed for gravity inversion can be divided into Monte Carlo methods, matrix inversion methods, and spectral methods. Here we employ both matrix inversion and Monte Carlo in order to explore the space of exact solutions, in a method which is particularly suited for arbitrary shape bodies. We expand the mass density function using orthogonal polynomials, and map the contribution of each term to the global gravitational field generated. This map is linear in the density terms, and can be pseudo-inverted in the under-determined regime using QR decomposition, to obtain a basis of the affine space of exact interior structure solutions. As the interior structure solutions are degenerate, assumptions have to be made in order to control their properties, and these assumptions can be transformed into scalar functions and used to explore the solutions space using Monte Carlo techniques. Sample applications show that the range of solutions tend to converge towards the nominal one as long as the generic assumptions made are correct, even in the presence of moderate noise. We present the underlying mathematical formalism and an analysis of how to impose specific features on the global solution, including uniform solutions, gradients, and layered models. Analytical formulas for the computation of the relevant quantities when the shape is represented using several common methods are included in the Appendix.
Impacts between planetary-sized bodies can explain the origin of satellites orbiting large ($R>500$~km) trans-Neptunian objects. Their water rich composition, along with the complex phase diagram of water, make it important to accurately model the wide range of thermodynamic conditions material experiences during an impact event and in the debris disk. Since differences in the thermodynamics may influence the system dynamics, we seek to evaluate how the choice of an equation of state (EOS) alters the systems evolution. Specifically, we compare two EOSs that are constructed by different approaches: either by a simplified analytic description (Tillotson), or by interpolation of tabulated data (Sesame). Approximately $50$ pairs of Smoothed Particle Hydrodynamics impact simulations were performed, with similar initial conditions but different EOSs, in the parameter space in which the Pluto-Charon binary is thought to form (slow impacts between Pluto-size, water rich bodies). Generally, we show that impact outcomes (e.g., circumplanetary debris disk) are consistent between EOSs. Some differences arise, importantly in the production of satellitesimals (large intact clumps) that form in the post-impact debris disk. When utilizing an analytic EOS, the emergence of satellitesimals is highly certain, while when using the tabulated EOS it is less common. This is because for the typical densities and energies experienced in these impacts, the analytic EOS predicts very low pressure values, leading to particles artificially aggregating by a tensile instability.
During the formation of terrestrial planets, volatile loss may occur through nebular processing, planetesimal differentiation, and planetary accretion. We investigate iron meteorites as an archive of volatile loss during planetesimal processing. The carbon contents of the parent bodies of magmatic iron meteorites are reconstructed by thermodynamic modelling. Calculated solid/molten alloy partitioning of C increases greatly with liquid S concentration and inferred parent body C concentrations range from 0.0004 to 0.11 wt.%. Parent bodies fall into 2 compositional clusters characterized by cores with medium, and low C/S. Both of these require significant planetesimal degassing, as metamorphic devolatilization on chondrite-like precursors is insufficient to account for their C depletions. Planetesimal core formation models, ranging from closed system extraction to degassing of a wholly molten body, show that significant open system silicate melting and volatile loss is required to match medium and low C/S parent body core compositions. Greater depletion in C relative to S is the hallmark of silicate degassing, indicating that parent body core compositions record processes that affect composite silicate/iron planetesimals. Degassing of bare cores stripped of their silicate mantles would deplete S with negligible C loss, and could not account for inferred parent body core compositions. Devolatilization during small-body differentiation is thus a key process in shaping the volatile inventory of terrestrial planets derived from planetesimals and planetary embryos.
It is known that the so-called problem of solar power pacemaker related to possible existence of some hidden but key mechanism of energy influence of the Sun on fundamental geophysical processes is one of the principal and puzzling problems of modern climatology. The tracks of this mechanism have been shown up in different problems of solar-terrestrial physics for a long time and, in particular, in climatology, where the solar-climate variability is stably observed. However, the mechanisms by which small changes in the Suns energy (solar irradiance or insolation) output during the solar cycle can cause change in the weather and climate are still unknown. We analyze possible causes of the solar-climate variability concentrating ones attention on the physical substantiation of strong correlation between the temporal variations of magnetic flux of the solar tachocline zone and the Earth magnetic field (Y-component). We propose an effective mechanism of solar dynamo-geodynamo connection which plays the role of the solar power pacemaker of the Earth global climate.
The porosity of an asteroid is important when studying the evolution of our solar system through small bodies and for planning mitigation strategies to avoid disasters due to asteroid impacts. Our knowledge of asteroid porosity largely relies on meteorites sampled on Earth. However, chondrites sampled on Earth are suggested to be sorted by strength. In this study, we obtained an estimate of the most porous structure of primordial granular chondrite parent bodies based on measurements of the compaction behavior of chondrite component analogs. We measured compaction curves of dust and dust-bead mixture samples. The dust sample consisted of various spherical and irregular particles with diameters on the order of 10^0-10^1 $mu$m. The mixture sample consisted of dust and beads with different dust volume fractions (~0.2-1). We used 1.5 and 4.8 $mu$m particles as dust as a first step, although the typical size of materials in matrix may be much smaller. We approximated the compaction curve of each sample with a power-law form and calculated the porosity structure of the primordial chondrite parent bodies using the experimental results. Our results show that the primordial parent bodies are likely to have higher porosity than the chondrites. Moreover, the relatively higher volume fraction of the matrix may be one of the reasons why most meteorites with high porosity are carbonaceous chondrites.
Collisions are the core agent of planet formation. In this work, we derive an analytic description of the dynamical outcome for any collision between gravity-dominated bodies. We conduct high-resolution simulations of collisions between planetesimals; the results are used to isolate the effects of different impact parameters on collision outcome. During growth from planetesimals to planets, collision outcomes span multiple regimes: cratering, merging, disruption, super-catastrophic disruption, and hit-and-run events. We derive equations (scaling laws) to demarcate the transition between collision regimes and to describe the size and velocity distributions of the post-collision bodies. The scaling laws are used to calculate maps of collision outcomes as a function of mass ratio, impact angle, and impact velocity, and we discuss the implications of the probability of each collision regime during planet formation. The analytic collision model presented in this work will significantly improve the physics of collisions in numerical simulations of planet formation and collisional evolution. (abstract abridged)