-
Ejection and Dynamics of Aggregates in the Coma of Comet 67P/Churyumov-Gerasimenko
Authors:
Pablo Lemos,
Jessica Agarwal,
Raphael Marschall,
Marius Pfeifer
Abstract:
The process of cometary activity continues to pose a challenging question in cometary science. The activity modeling of comet 67P/Churyumov-Gerasimenko, based on data from the Rosetta mission, has significantly enhanced our comprehension of cometary activity. But thermophysical models have difficulties in simultaneously explaining the production rates of various gas species and dust. It has been s…
▽ More
The process of cometary activity continues to pose a challenging question in cometary science. The activity modeling of comet 67P/Churyumov-Gerasimenko, based on data from the Rosetta mission, has significantly enhanced our comprehension of cometary activity. But thermophysical models have difficulties in simultaneously explaining the production rates of various gas species and dust. It has been suggested that different gas species might be responsible for the ejection of refractory material in distinct size ranges. This work focuses on investigating abundance and the ejection mechanisms of large ($\gtrsim$ 1 cm) aggregates from the comet nucleus. We aim to determine their properties and map the distribution of their source regions across the comet surface. This can place constraints on activity models for comets. We examined 189 images acquired at five epochs by the OSIRIS/NAC instrument. Our goal was to identify bright tracks produced by individual aggregates as they traversed the camera field of view. We generated synthetic images based on the output of dynamical simulations involving various types of aggregates. By comparing these synthetic images with the observations, we determine the characteristics of the simulated aggregates that most closely resembled the observations. We identified over 30000 tracks present in the OSIRIS images, derived constraints on the characteristics of the aggregates and mapped their origins on the nucleus surface. The aggregates have an average radius of $\simeq5$ cm, and a bulk density consistent with that of the comet's nucleus. Due to their size, gas drag exerts only a minor influence on their dynamical behavior, so an initial velocity is needed in order to bring them into the camera field of view. The source regions of these aggregates are predominantly located near the boundaries of distinct terrains on the surface.
△ Less
Submitted 13 May, 2024;
originally announced May 2024.
-
Dynamics and potential origins of decimeter-sized particles around comet 67P/Churyumov-Gerasimenko
Authors:
Marius Pfeifer,
Jessica Agarwal,
Raphael Marschall,
Björn Grieger,
Pablo Lemos
Abstract:
Methods. We algorithmically tracked thousands of individual particles through four OSIRIS/NAC image sequences of 67P's near-nucleus coma. We then traced concentrated particle groups back to the nucleus surface, and estimated their potential source regions, size distributions, and projected dynamical parameters. Finally, we compared the observed activity to dust coma simulations. Results. We traced…
▽ More
Methods. We algorithmically tracked thousands of individual particles through four OSIRIS/NAC image sequences of 67P's near-nucleus coma. We then traced concentrated particle groups back to the nucleus surface, and estimated their potential source regions, size distributions, and projected dynamical parameters. Finally, we compared the observed activity to dust coma simulations. Results. We traced back 409 decimeter-sized particles to four suspected source regions. The regions strongly overlap and are mostly confined to the Khonsu-Atum-Anubis area. The activity may be linked to rugged terrain, and the erosion of fine dust and the ejection of large boulders may be mutually exclusive. Power-law indices fitted to the particle size--frequency distributions range from $3.4 \pm 0.3$ to $3.8 \pm 0.4$. Gas drag fits to the radial particle accelerations provide an estimate for the local gas production rates ($Q_\text{g} = 3.6 \cdot 10^{-5}$ kg s$^{-1}$ m$^{-2}$), which is several times higher than our model predictions based on purely insolation-driven water ice sublimation. Our observational results and our modeling results both reveal that our particles were likely ejected with substantial nonzero initial velocities of around 0.5$-$0.6 m s$^{-1}$. Conclusions. Our findings strongly suggest that the observed ejection of decimeter-sized particles cannot be explained by water ice sublimation and favorable illumination conditions alone. Instead, the local structures and compositions of the source regions likely play a major role. In line with current ejection models of decimeter-sized particles, we deem an overabundance of CO$_2$ ice and its sublimation to be the most probable driver. In addition, because of the significant initial velocities, we suspect the ejection events to be considerably more energetic than gradual liftoffs.
△ Less
Submitted 28 February, 2024;
originally announced February 2024.
-
Origin and Evolution of Jupiter's Trojan Asteroids
Authors:
William F. Bottke,
Raphael Marschall,
David Nesvorný,
David Vokrouhlický
Abstract:
The origin of the Jupiter Trojan asteroids has long been a mystery. Dynamically, the population, which is considerably smaller than the main asteroid belt, librates around Jupiter's stable L4 and L5 Lagrange points, 60 deg ahead and behind Jupiter. It is thought that these bodies were captured into these orbits early in solar system history, but any capture mechanism must also explain why the Troj…
▽ More
The origin of the Jupiter Trojan asteroids has long been a mystery. Dynamically, the population, which is considerably smaller than the main asteroid belt, librates around Jupiter's stable L4 and L5 Lagrange points, 60 deg ahead and behind Jupiter. It is thought that these bodies were captured into these orbits early in solar system history, but any capture mechanism must also explain why the Trojans have an excited inclination distribution, with some objects reaching inclinations of 35 deg. The Trojans themselves, individually and in aggregate, also have spectral and physical properties that appear consistent with many small bodies found in the outer solar system (e.g., irregular satellites, Kuiper belt objects). In this review, we assemble what is known about the Trojans and discuss various models for their origin and collisional evolution. It can be argued that the Trojans are unlikely to be captured planetesimals from the giant planet zone, but instead were once denizens of the primordial Kuiper belt, trapped by the events taking place during a giant planet instability. The Lucy mission to the Trojans is therefore well positioned to not only answer questions about these objects, but also about their place in planet formation and solar system evolution studies.
△ Less
Submitted 5 December, 2023;
originally announced December 2023.
-
Dust Emission and Dynamics
Authors:
Jessica Agarwal,
Yoonyoung Kim,
Michael S. P. Kelley,
Raphael Marschall
Abstract:
When viewed from Earth, most of what we observe of a comet is dust. The influence of solar radiation pressure on the trajectories of dust particles depends on their cross-section to mass ratio. Hence solar radiation pressure acts like a mass spectrometer inside a cometary tail. The appearances of cometary dust tails have long been studied to obtain information on the dust properties, such as chara…
▽ More
When viewed from Earth, most of what we observe of a comet is dust. The influence of solar radiation pressure on the trajectories of dust particles depends on their cross-section to mass ratio. Hence solar radiation pressure acts like a mass spectrometer inside a cometary tail. The appearances of cometary dust tails have long been studied to obtain information on the dust properties, such as characteristic particle size and initial velocity when entering the tail. Over the past two decades, several spacecraft missions to comets have enabled us to study the dust activity of their targets at much greater resolution than is possible with a telescope on Earth or in near-Earth space, and added detail to the results obtained by the spacecraft visiting comet 1P/Halley in 1986. We now know that the dynamics of dust in the inner cometary coma is complex and includes a significant fraction of particles that will eventually fall back to the surface. The filamented structure of the near-surface coma is thought to result from a combination of topographic focussing of the gas flow, inhomogeneous distribution of activity across the surface, and projection effects. It is possible that some larger-than-centimetre debris contains ice when lifted from the surface, which can affect its motion. Open questions remain regarding the microphysics of the process that leads to the detachment and lifting of dust from the surface, the evolution of the dust while travelling away from the nucleus, and the extent to which information on the nucleus activity can be retrieved from remote observations of the outer coma and tail.
△ Less
Submitted 10 July, 2024; v1 submitted 22 September, 2023;
originally announced September 2023.
-
Neutral gas coma dynamics: modeling of flows and attempts to link inner coma structures to properties of the nucleus
Authors:
Raphael Marschall,
Björn J. R. Davidsson,
Martin Rubin,
Valeriy Tenishev
Abstract:
Deriving properties of cometary nuclei from coma data is of significant importance for our understanding of cometary activity and has implications beyond. Ground-based data represent the bulk of measurements available for comets. Yet, to date these observations only access a comet's gas and dust coma at rather large distances from the surface and do not directly observe its surface. In contrast, s…
▽ More
Deriving properties of cometary nuclei from coma data is of significant importance for our understanding of cometary activity and has implications beyond. Ground-based data represent the bulk of measurements available for comets. Yet, to date these observations only access a comet's gas and dust coma at rather large distances from the surface and do not directly observe its surface. In contrast, spacecraft missions are one of the only tools that gain direct access to surface measurements. However, such missions are limited to roughly one per decade. We can overcome these challenges by recognising that the coma contains information about the nucleus's properties. In particular, the near-surface gas environment is most representative of the nucleus. It can inform us about the composition, regionality of activity, and sources of coma features and how they link to the topography, morphology, or other surface properties. The inner coma data is a particularly good proxy because it has not yet, or only marginally, been contaminated by coma chemistry or secondary gas sources. Additionally, the simultaneous observation of the innermost coma with the surface provides the potential to make a direct link between coma measurements and the nucleus. If we hope to link outer coma measurements obtained by Earth-based telescopes to the surface, we must first understand how the inner coma measurements are linked to the surface. Numerical models that describe the flow from the surface into the immediate surroundings are needed to make this connection. This chapter focuses on the advances made to understand the flow of the neutral gas coma from the surface to distances up to a few tens of nuclei radii. We describe both simple/heuristic models and state-of-the-art physically consistent models. Model limitations and what they each are best suited for are discussed.
△ Less
Submitted 10 August, 2023;
originally announced August 2023.
-
An inflationary disk phase to explain extended protoplanetary dust disks
Authors:
Raphael Marschall,
Alessandro Morbidelli
Abstract:
Understanding planetesimal formation is an essential first step to understanding planet formation. The distribution of these first solid bodies will drive the locations where planetary embryos can grow. We seek to understand the parameter space of possible protoplanetary disk formation and evolution models of our Solar System. A good protoplanetary disk scenario for the Solar System must meet at l…
▽ More
Understanding planetesimal formation is an essential first step to understanding planet formation. The distribution of these first solid bodies will drive the locations where planetary embryos can grow. We seek to understand the parameter space of possible protoplanetary disk formation and evolution models of our Solar System. A good protoplanetary disk scenario for the Solar System must meet at least the following three criteria: 1) an extended dust disk (at least 45 au); 2) formation of planetesimals in at least two distinct locations; and 3) transport of high temperatures condensates (i.e., calcium-aluminium-rich inclusion, CAIs) to the outer disk. We explore a large parameter space to study the effect of the disk viscosity, the timescale of infall of material into the disk, the distance within which material is deposited into the disk, and the fragmentation threshold of dust particles.
We find that scenarios with a large initial disk viscosity ($α>0.05$), relatively short infall timescale ($T_{infall}<100-200$ kyr), and a small centrifugal radius ($R_C\sim0.4$~au; the distance within which material falls into the disk) result in disks that satisfy the criteria for a good protoplanetary disk of the Solar System. The large initial viscosity and short infall timescale result in a rapid initial expansion of the disk, which we dub the inflationary phase of the disk. Furthermore, a temperature-dependent fragmentation threshold, which mimics that cold icy particles break more easily, results in larger and more massive disks. This results in more "icy" than "rocky" planetesimals. Such scenarios are also better in line with our Solar System, which has small terrestrial planets and massive giant planet cores. Finally, we find that scenarios with large $R_C$ cannot transport CAIs to the outer disk and do not produce planetesimals at two locations within the disk.
△ Less
Submitted 3 July, 2023;
originally announced July 2023.
-
Activity distribution of comet 67P/Churyumov-Gerasimenko from combined measurements of non-gravitational forces and torques
Authors:
Nicholas Attree,
Laurent Jorda,
Olivier Groussin,
Jessica Agarwal,
Riccardo Lasagni Manghi,
Paolo Tortora,
Marco Zannoni,
Raphael Marschall
Abstract:
Aims. Understanding the activity is vital for deciphering the structure, formation, and evolution of comets. We investigate models of cometary activity by comparing them to the dynamics of 67P/Churyumov-Gerasimenko. Methods. We matched simple thermal models of water activity to the combined Rosetta datasets by fitting to the total outgassing rate and four components of the outgassing induced non-g…
▽ More
Aims. Understanding the activity is vital for deciphering the structure, formation, and evolution of comets. We investigate models of cometary activity by comparing them to the dynamics of 67P/Churyumov-Gerasimenko. Methods. We matched simple thermal models of water activity to the combined Rosetta datasets by fitting to the total outgassing rate and four components of the outgassing induced non-gravitational force and torque, with a final manual adjustment of the model parameters to additionally match the other two torque components. We parametrised the thermal model in terms of a distribution of relative activity over the surface of the comet, and attempted to link this to different terrain types. We also tested a more advanced thermal model based on a pebble structure. Results. We confirm a hemispherical dichotomy and non-linear water outgassing response to insolation. The southern hemisphere of the comet and consolidated terrain show enhanced activity relative to the northern hemisphere and dust-covered, unconsolidated terrain types, especially at perihelion. We further find that the non-gravitational torque is especially sensitive to the activity distribution, and to fit the pole-axis orientation in particular, activity must be concentrated (in excess of the already high activity in the southern hemisphere and consolidated terrain) around the south pole and on the body and neck of the comet over its head. This is the case for both the simple thermal model and the pebble-based model. Overall, our results show that water activity cannot be matched by a simple model of sublimating surface ice driven by the insolation alone, regardless of the surface distribution, and that both local spatial and temporal variations are needed to fit the data.
△ Less
Submitted 12 January, 2023;
originally announced January 2023.
-
Implications for the collisional strength of Jupiter Trojans from the Eurybates family
Authors:
Raphael Marschall,
David Nesvorný,
Rogerio Deienno,
Ian Wong,
Harold F. Levison,
William F. Bottke
Abstract:
In this work, we model the collisional evolution of the Jupiter Trojans and determined under which conditions the Eurybates-Queta system survives. We show that the collisional strength of the Jupiter Trojans and the age of the Eurybates family and by extension Queta are correlated.
The collisional grinding of the Jupiter Trojan population over 4.5 Gy results in a size-frequency distribution (SFD…
▽ More
In this work, we model the collisional evolution of the Jupiter Trojans and determined under which conditions the Eurybates-Queta system survives. We show that the collisional strength of the Jupiter Trojans and the age of the Eurybates family and by extension Queta are correlated.
The collisional grinding of the Jupiter Trojan population over 4.5 Gy results in a size-frequency distribution (SFD) that remains largely unaltered at large sizes (>10 km) but is depleted at small sizes (10 m to 1 km). This results in a turnover in the SFD, the location of which depends on the collisional strength of the material. It is to be expected that the Trojan SFD bends between 1 and 10 km.
Based on the SFD of the Eurybates family, we find that the family was likely the result of a catastrophic impact onto a 100 km rubble pile target. This corresponds to objects with a rather low collisional strength (10 times weaker than that of basaltic material studied in Benz & Asphaug1999).
Assuming this weak strength, and an initial cumulative slope of the size frequency distribution of 2.1 between diameters of 2 m and 100 km when the Trojans were captured, the existence of Queta, the satellite of Eurybates, implies an upper limit for the family age of 3.7 Gy. Alternatively, we demonstrate that an unconventional collisional strength with a minimum at 20 m is a plausible candidate to ensure the survival of Queta over the age of the Solar System.
Finally, we show how different collisional histories change the expected number of craters on the targets of the Lucy mission and that Lucy will be able to differentiate between them.
△ Less
Submitted 22 August, 2022;
originally announced August 2022.
-
Determining the dust environment of an unknown comet for a spacecraft fly-by: The case of ESA's Comet Interceptor mission
Authors:
Raphael Marschall,
Vladimir Zakharov,
Cecilia Tubiana,
Michael S. P. Kelley,
Carlos Corral van Damme,
Colin Snodgrass,
Geraint H. Jones,
Stavro L. Ivanovski,
Frank Postberg,
Vincenzo Della Corte,
Jean-Baptiste Vincent,
Olga Muñoz,
Fiorangela La Forgia,
Anny-Chantal Levasseur-Regourd,
the Comet Interceptor Team
Abstract:
We present a statistical approach to assess the dust environment for a yet unknown comet (or when its parameters are known only with large uncertainty). This is of particular importance for missions such as ESA's Comet Interceptor mission to a dynamically new comet.
We find that the lack of knowledge of any particular comet results in very large uncertainties (~3 orders of magnitude) for the dus…
▽ More
We present a statistical approach to assess the dust environment for a yet unknown comet (or when its parameters are known only with large uncertainty). This is of particular importance for missions such as ESA's Comet Interceptor mission to a dynamically new comet.
We find that the lack of knowledge of any particular comet results in very large uncertainties (~3 orders of magnitude) for the dust densities within the coma. The most sensitive parameters affecting the dust densities are the dust size distribution, the dust production rate and coma brightness, often quantified by Af$ρ$. Further, the conversion of a coma's brightness (Af$ρ$) to a dust production rate is poorly constrained. The dust production rate can only be estimated down to an uncertainty of ~0.5 orders of magnitude if the dust size distribution is known in addition to the Af$ρ$.
To accurately predict the dust environment of a poorly known comet, a statistical approach as we propose here needs to be taken to properly reflect the uncertainties. This can be done by calculating an ensemble of comae covering all possible combinations within parameter space as shown in this work.
△ Less
Submitted 9 August, 2022;
originally announced August 2022.
-
Binary Planetesimal Formation from Gravitationally Collapsing Pebble Clouds
Authors:
David Nesvorny,
Rixin Li,
Jacob B. Simon,
Andrew N. Youdin,
Derek C. Richardson,
Raphael Marschall,
William M. Grundy
Abstract:
Planetesimals are compact astrophysical objects roughly 1-1000 km in size, massive enough to be held together by gravity. They can grow by accreting material to become full-size planets. Planetesimals themselves are thought to form by complex physical processes from small grains in protoplanetary disks. The streaming instability (SI) model states that mm/cm-size particles (pebbles) are aerodynamic…
▽ More
Planetesimals are compact astrophysical objects roughly 1-1000 km in size, massive enough to be held together by gravity. They can grow by accreting material to become full-size planets. Planetesimals themselves are thought to form by complex physical processes from small grains in protoplanetary disks. The streaming instability (SI) model states that mm/cm-size particles (pebbles) are aerodynamically collected into self-gravitating clouds which then directly collapse into planetesimals. Here we analyze ATHENA simulations of the SI to characterize the initial properties (e.g., rotation) of pebble clouds. Their gravitational collapse is followed with the PKDGRAV N-body code, which has been modified to realistically account for pebble collisions. We find that pebble clouds rapidly collapse into short-lived disk structures from which planetesimals form. The planetesimal properties depend on the cloud's scaled angular momentum, l=L/(M R_H^2 Omega, where L and M are the angular momentum and mass, R_H is the Hill radius, and Omega is the orbital frequency. Low-l pebble clouds produce tight (or contact) binaries and single planetesimals. Compact high-l clouds give birth to binary planetesimals with attributes that closely resemble the equal-size binaries found in the Kuiper belt. Significantly, the SI-triggered gravitational collapse can explain the angular momentum distribution of known equal-size binaries -- a result pending verification from studies with improved resolution. About 10% of collapse simulations produce hierarchical systems with two or more large moons. These systems should be found in the Kuiper belt when observations reach the threshold sensitivity.
△ Less
Submitted 13 November, 2020;
originally announced November 2020.
-
The dust-to-gas ratio, size distribution, and dust fall-back fraction of comet 67P/Churyumov-Gerasimenko: Inferences from linking the optical and dynamical properties of the inner comae
Authors:
Raphael Marschall,
Johannes Markkanen,
Selina-Barbara Gerig,
Olga Pinzón-Rodríguez,
Nicolas Thomas,
Jong-Shinn Wu
Abstract:
In this work, we present results that simultaneously constrain the dust size distribution, dust-to-gas ratio, fraction of dust re-deposition, and total mass production rates for comet 67P/Churyumov-Gerasimenko. We use a 3D Direct Simulation Monte Carlo (DSMC) gas dynamics code to simulate the inner gas coma of the comet for the duration of the Rosetta mission. The gas model is constrained by ROSIN…
▽ More
In this work, we present results that simultaneously constrain the dust size distribution, dust-to-gas ratio, fraction of dust re-deposition, and total mass production rates for comet 67P/Churyumov-Gerasimenko. We use a 3D Direct Simulation Monte Carlo (DSMC) gas dynamics code to simulate the inner gas coma of the comet for the duration of the Rosetta mission. The gas model is constrained by ROSINA/COPS data. Further, we simulate for different epochs the inner dust coma using a 3D dust dynamics code including gas drag and the nucleus' gravity. Using advanced dust scattering properties these results are used to produce synthetic images that can be compared to the OSIRIS data set. These simulations allow us to constrain the properties of the dust coma and the total gas and dust production rates.
We determined a total volatile mass loss of $(6.1 \pm 1.5) \cdot 10^9$~kg during the 2015 apparition. Further, we found that power-laws with $q=3.7^{+0.57}_{-0.078}$ are consistent with the data. This results in a total of $5.1^{+6.0}_{-4.9}\cdot10^9$~kg of dust being ejected from the nucleus surface, of which $4.4^{+4.9}_{-4.2}\cdot10^9$~kg escape to space and $6.8^{+11}_{-6.8}\cdot10^8$~kg (or an equivalent of $14^{+22}_{-14}$~cm over the smooth regions) is re-deposited on the surface. This leads to a dust-to-gas ratio of $0.73^{+1.3}_{-0.70}$ for the escaping material and $0.84^{+1.6}_{-0.81}$ for the ejected material. We have further found that the smallest dust size must be strictly smaller than $\sim30μ$m and nominally even smaller than $\sim12μ$m.
△ Less
Submitted 27 May, 2020;
originally announced May 2020.
-
Stability of Jovian Trojans and their collisional families
Authors:
Timothy R. Holt,
David Nesvorny,
Jonathan Horner,
Rachel King,
Raphael Marschall,
Melissa Kamrowski,
Brad Carter,
Leigh Brookshaw,
Christopher Tylor
Abstract:
The Jovian Trojans are two swarms of objects located around the L$_4$ and L$_5$ Lagrange points. The population is thought to have been captured by Jupiter during the Solar system's youth. Within the swarms, six collisional families have been identified in previous work, with four in the L$_4$ swarm, and two in the L$_5$. Our aim is to investigate the stability of the two Trojan swarms, with a par…
▽ More
The Jovian Trojans are two swarms of objects located around the L$_4$ and L$_5$ Lagrange points. The population is thought to have been captured by Jupiter during the Solar system's youth. Within the swarms, six collisional families have been identified in previous work, with four in the L$_4$ swarm, and two in the L$_5$. Our aim is to investigate the stability of the two Trojan swarms, with a particular focus on these collisional families. We find that the members of Trojan swarms escape the population at a linear rate, with the primordial L$_4$ (23.35% escape) and L$_5$ (24.89% escape) population sizes likely 1.31 and 1.35 times larger than today. Given that the escape rates were approximately equal between the two Trojan swarms, our results do not explain the observed asymmetry between the two groups, suggesting that the numerical differences are primordial in nature, supporting previous studies. Upon leaving the Trojan population, the escaped objects move onto orbits that resemble those of the Centaur and short-period comet populations. Within the Trojan collisional families, the 1996 RJ and 2001 UV$_{209}$ families are found to be dynamically stable over the lifetime of the Solar system, whilst the Hektor, Arkesilos and Ennomos families exhibit various degrees of instability. The larger Eurybates family shows 18.81\% of simulated members escaping the Trojan population. Unlike the L4 swarm, the escape rate from the Eurybates family is found to increase as a function of time, allowing an age estimation of approximately $1.045\pm 0.364 \times 10^9$ years.
△ Less
Submitted 7 May, 2020;
originally announced May 2020.
-
Limitations in the determination of surface emission distributions on comets through modelling of observational data -- A case study based on Rosetta observations
Authors:
Raphael Marschall,
Ying Liao,
Nicolas Thomas,
Jong-Shinn Wu
Abstract:
The European Space Agency's (ESA) Rosetta mission has returned a vast data set of measurements of the inner gas coma of comet 67P/Churyumov-Gerasimenko. These measurements have been used by different groups to determine the distribution of the gas sources at the nucleus surface. The solutions that have been found differ from each other substantially and illustrate the degeneracy of this issue. It…
▽ More
The European Space Agency's (ESA) Rosetta mission has returned a vast data set of measurements of the inner gas coma of comet 67P/Churyumov-Gerasimenko. These measurements have been used by different groups to determine the distribution of the gas sources at the nucleus surface. The solutions that have been found differ from each other substantially and illustrate the degeneracy of this issue. It is the aim of this work to explore the limitations that current gas models have in linking the coma measurements to the surface. In particular, we discuss the sensitivity of Rosetta's ROSINA/COPS, VIRTIS, and MIRO instruments to differentiate between vastly different spatial distributions of the gas emission from the surface. We have applied a state of the art 3D DSMC gas dynamics code to simulate the inner gas coma of different models that vary in the fraction of the surface that contains ice and in different sizes of active patches. These different distributions result in jet interactions that differ in their dynamical behaviour. We have found that ROSINA/COPS measurements by themselves cannot detect the differences in our models. While ROSINA/COPS measurements are important to constrain the regional inhomogeneities of the gas emission, they can by themselves not determine the surface emission distribution of the gas sources to a spatial accuracy of better than a few hundred metres (400 m). Any solutions fitting the ROSINA/COPS measurements is hence fundamentally degenerate, be it through a forward or inverse model. Only other instruments with complementary measurements can potentially lift this degeneracy as we show here for VIRTIS and MIRO. Finally, as a by-product, we have explored the effect of our activity distributions on lateral flow at the surface that may be responsible for some of the observed aeolian features.
△ Less
Submitted 27 April, 2020;
originally announced April 2020.
-
Superparticle Method for Simulating Collisions
Authors:
David Nesvorny,
Andrew N. Youdin,
Raphael Marschall,
Derek C. Richardson
Abstract:
For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated to represent a collectively large collisional cross section of real particles. Here we develop a superparticle-based method that replicates the kinetic energy…
▽ More
For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated to represent a collectively large collisional cross section of real particles. Here we develop a superparticle-based method that replicates the kinetic energy loss during real-world collisions, implement it in an $N$-body code and test it. The tests provide interesting insights into dynamics of self gravitating collisional systems. They show how particle systems evolve over several free fall timescales to form central concentrations and equilibrated outer shells. The superparticle method can be extended to account for the accretional growth of objects during inelastic mergers.
△ Less
Submitted 14 April, 2020;
originally announced April 2020.
-
Migrating Scarps as a Significant Driver for Cometary Surface Evolution
Authors:
Samuel Birch,
Alexander Hayes,
Orkan Umurhan,
Yuhui Tang,
Jean-Baptiste Vincent,
Nilda Oklay,
Dennis Bodewits,
Bjorn Davidsson,
Raphael Marschall,
Jason Soderblom,
Jeff Moore,
Paul Corlies,
Steven Squyres
Abstract:
Rosetta observations of 67P/Churyumov-Gerasimenko (67P) reveal that most changes occur in the fallback-generated smooth terrains, vast deposits of granular material blanketing the comet's northern hemisphere. These changes express themselves both morphologically and spectrally across the nucleus, yet we lack a model that describes their formation and evolution. Here we present a self-consistent mo…
▽ More
Rosetta observations of 67P/Churyumov-Gerasimenko (67P) reveal that most changes occur in the fallback-generated smooth terrains, vast deposits of granular material blanketing the comet's northern hemisphere. These changes express themselves both morphologically and spectrally across the nucleus, yet we lack a model that describes their formation and evolution. Here we present a self-consistent model that thoroughly explains the activity and mass loss from Hapi's smooth terrains. Our model predicts the removal of dust via re-radiated solar insolation localized within depression scarps that are substantially more ice-rich than previously expected. We couple our model with numerous Rosetta observations to thoroughly capture the seasonal erosion of Hapi's smooth terrains, where local scarp retreat gradually removes the uppermost dusty mantle. As sublimation-regolith interactions occur on rocky planets, comets, icy moons and KBOs, our coupled model and observations provide a foundation for future understanding of the myriad of sublimation-carved worlds.
△ Less
Submitted 9 September, 2019;
originally announced September 2019.
-
Diurnal variation of dust and gas production in comet 67P/Churyumov-Gerasimenko at the inbound equinox as seen by OSIRIS and VIRTIS-M on board Rosetta
Authors:
C. Tubiana,
G. Rinaldi,
C. Güttler,
C. Snodgrass,
X. Shi,
X. Hu,
R. Marschall,
M. Fulle,
D. Bockelée-Morvan,
G. Naletto,
F. Capaccioni,
H. Sierks,
G. Arnold,
M. A. Barucci,
J. -L. Bertaux,
I. Bertini,
D. Bodewits,
M. T. Capria,
M. Ciarniello,
G. Cremonese,
J. Crovisier,
V. Da Deppo,
S. Debei,
M. De Cecco,
J. Deller
, et al. (31 additional authors not shown)
Abstract:
On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the…
▽ More
On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the one calculated using a thermo-physical model. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O) and 4200 nm (CO2). We measured the average signal in a circular annulus, to study spatial variation around the comet, and in a sector of the annulus, to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet centre. The spatial correlation between dust and water, both coming from the sun-lit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0deg sub-solar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production, with respect to the value calculated using a simple thermo-physical model, is observed when the head lobe and regions of the Southern hemisphere with strong seasonal variations are illuminated. A drastic decrease in dust production, when the water production (both measured and from the model) displays a maximum, happens when typical Northern consolidated regions are illuminated and the Southern hemisphere regions with strong seasonal variations are instead in shadow.
△ Less
Submitted 8 May, 2019;
originally announced May 2019.
-
Local manifestations of cometary activity
Authors:
Jean-Baptiste Vincent,
Tony Farnham,
Ekkehard Kührt,
Yuri Skorov,
Raphael Marschall,
Nilda Oklay,
Ramy El-Maarry,
Horst Uwe Keller
Abstract:
Comets are made of volatile and refractory material and naturally experience various degrees of sublimation as they orbit around the Sun. This gas release, accompanied by dust, represents what is traditionally described as activity. Although the basic principles are well established, most details remain elusive, especially regarding the mechanisms by which dust is detached from the surface and sub…
▽ More
Comets are made of volatile and refractory material and naturally experience various degrees of sublimation as they orbit around the Sun. This gas release, accompanied by dust, represents what is traditionally described as activity. Although the basic principles are well established, most details remain elusive, especially regarding the mechanisms by which dust is detached from the surface and subsequently accelerated by the gas flows surrounding the nucleus. During its 2 years rendez-vous with comet 67P/Churyumov-Gerasimenko, ESA's Rosetta has observed cometary activity with unprecedented details, in both the inbound and outbound legs of the comet's orbit. This trove of data provides a solid ground on which new models of activity can be built. In this chapter, we review how activity manifests at close distance from the surface, establish a nomenclature for the different types of observed features, discuss how activity is at the same time transforming and being shaped by the topography, and finally address several potential mechanisms.
△ Less
Submitted 11 April, 2019;
originally announced April 2019.
-
Constraining spatial pattern of early activity of comet 67P/C-G with 3D modeling of the MIRO observations
Authors:
Y. Zhao,
L. Rezac,
P. Hartogh,
J. Ji,
R. Marschall,
H. U. Keller
Abstract:
Our aim is to investigate early activity (July 2014) of 67P/CG with 3D coma and radiative transfer modeling of MIRO measurements, accounting for nucleus shape, illumination, and orientation of the comet. We investigate MIRO line shape information for spatial distribution of water activity on the nucleus during the onset of activity. During this period we show that MIRO line shape have enough infor…
▽ More
Our aim is to investigate early activity (July 2014) of 67P/CG with 3D coma and radiative transfer modeling of MIRO measurements, accounting for nucleus shape, illumination, and orientation of the comet. We investigate MIRO line shape information for spatial distribution of water activity on the nucleus during the onset of activity. During this period we show that MIRO line shape have enough information to clearly isolate contribution from Hapi and Inhotep independently, and compare it to the nominal case of activity from the entire illuminated surface. We also demonstrate that spectral line shapes differ from the 1D model for different viewing geometries and coma conditions relevant to this study. Specifically, line shapes are somewhat sensitive to the location of the terminator in the coma. At last, fitting the MIRO observations we show that the Imhotep region (possible extended source of H$_{2}$O due to CO$_{2}$ activity) contributes only a small fraction of the total number of water molecules into MIRO beam in the early activity. On the other hand, a strong enhancement of water activity from the Hapi region seems required to fit the MIRO line shapes. This is consistent with earlier Rosetta results. Nevertheless, within the assumption of our coma and surface boundary conditions, we cannot get a reasonable fit to all MIRO mapping observations in July 2014, which may illustrate that a more sophisticated coma model or more accurate temperature/velocity distribution is needed.
△ Less
Submitted 27 February, 2019;
originally announced February 2019.