Abstract

Crater counting on the Moon and other bodies is crucial to constrain the dynamical history of the Solar System. This has traditionally been done by visual inspection of images, thus limiting the scope, efficiency, and/or accuracy of retrieval. In this paper we demonstrate the viability of using convolutional neural networks (CNNs) to determine the positions and sizes of craters from a Lunar digital elevation map (DEM). We recover 92% of craters from the human-generated test set and almost double the total number of crater detections. Of these new craters, 15% are smaller in diameter than the minimum crater size in the ground-truth dataset. Our median fractional longitude, latitude and radius errors are 11% or less, representing good agreement with the human-generated datasets. From a manual inspection of 361 new craters we estimate the false positive rate of new craters to be 11%. Moreover, our Moon-trained CNN performs well when tested on DEM images of Mercury, detecting a large fraction of craters in each map. Our results suggest that deep learning will be a useful tool for rapidly and automatically extracting craters on various Solar System bodies. Code is publicly available via the links above

Abstract

Anorthosites that comprise the bulk of the lunar crust are believed to have formed during solidification of a lunar magma ocean (LMO) in which these rocks would have floated to the surface. This early flotation crust would have formed a thermal blanket over the remaining LMO, prolonging solidification. Geochronology of lunar anorthosites indicates a long timescale of LMO cooling, or remelting and recrystallization in one or more late events. To better interpret this geochronology, we model LMO solidification in a scenario where the Moon is being continuously bombarded by returning projectiles released from the Moon-forming giant impact. More than one lunar mass of material escaped the Earth-Moon system onto heliocentric orbits following the giant impact, much of it to come back on returning orbits for a period of 100 Myr. If large enough, these projectiles would have punctured holes in the nascent floatation crust of the Moon, exposing the LMO to space and causing more rapid cooling. We model these scenarios using a thermal evolution model of the Moon that allows for production (by cratering) and evolution (solidification and infill) of holes in the flotation crust that insulates the LMO. For effective hole production, solidification of the magma ocean can be significantly expedited, decreasing the cooling time by more than a factor of 5. If hole production is inefficient, but shock conversion of projectile kinetic energy to thermal energy is efficient, then LMO solidification can be somewhat prolonged, lengthening the cooling time by 50% or more.

Abstract

In single star systems like our own Solar system, comets dominate the mass budget of bodies that are ejected into interstellar space, since they form further away and are less tightly bound. However 1I/`Oumuamua, the first interstellar object detected, appears asteroidal in its spectra and in its lack of detectable activity. We argue that the galactic budget of interstellar objects like 1I/`Oumuamua should be dominated by planetesimal material ejected during planet formation in circumbinary systems, rather than in single star systems or widely separated binaries. We further show that in circumbinary systems, rocky bodies should be ejected in comparable numbers to icy ones. This suggests that a substantial fraction of additional interstellar objects discovered in the future should display an active coma. We find that the rocky population, of which 1I/`Oumuamua seems to be a member, should be predominantly sourced from A-type and late B-star binaries.

Abstract

We provide a fast method for computing constraints on impactor pre-impact orbits, applying this to the late giant impacts in the Solar system. These constraints can be used to make quick, broad comparisons of different collision scenarios, identifying some immediately as low-probability events, and narrowing the parameter space in which to target follow-up studies with expensive N-body simulations. We benchmark our parameter space predictions, finding good agreement with existing N-body studies for the Moon. We suggest that high-velocity impact scenarios in the inner Solar system, including all currently proposed single impact scenarios for the formation of Mercury, should be disfavoured. This leaves a multiple hit-and-run scenario as the most probable currently proposed for the formation of Mercury.

Abstract

As discoveries of multiple planets in the habitable zone of their parent star mount, developing analytical techniques to quantify extrasolar intra-system panspermia will become increasingly important. Here, we provide user-friendly prescriptions that describe the asteroid impact characteristics which would be necessary to transport life both inwards and outwards within these systems within a single framework. Our focus is on projectile generation and delivery and our expressions are algebraic, eliminating the need for the solution of differential equations. We derive a probability distribution function for life-bearing debris to reach a planetary orbit, and describe the survival of micro-organisms during planetary ejection, their journey through interplanetary space, and atmospheric entry.

Abstract

We report 885 μm ALMA continuum flux densities for 24 Taurus members spanning the stellar/substellar boundary with spectral types from M4 to M7.75. Of the 24 systems, 22 are detected at levels ranging from 1.0 to 55.7 mJy. The two nondetections are transition disks, though other transition disks in the sample are detected. Converting ALMA continuum measurements to masses using standard scaling laws and radiative transfer modeling yields dust mass estimates ranging from ∼0.3 to 20 M_⊕. The dust mass shows a declining trend with central object mass when combined with results from submillimeter surveys of more massive Taurus members. The substellar disks appear as part of a continuous sequence and not a distinct population. Compared to older Upper Sco members with similar masses across the substellar limit, the Taurus disks are brighter and more massive. Both Taurus and Upper Sco populations are consistent with an approximately linear relationship in M_dust to M_star, although derived power-law slopes depend strongly upon choices of stellar evolutionary model and dust temperature relation. The median disk around early-M stars in Taurus contains a comparable amount of mass in small solids as the average amount of heavy elements in Kepler planetary systems on short-period orbits around M-dwarf stars, with an order of magnitude spread in disk dust mass about the median value. Assuming a gas-to-dust ratio of 100:1, only a small number of low-mass stars and brown dwarfs have a total disk mass amenable to giant planet formation, consistent with the low frequency of giant planets orbiting M dwarfs.

Abstract

This paper considers the dynamics of the scattering of planetesimals or planetary embryos by a planet on a circumstellar orbit. We classify six regions in the planet's mass versus semimajor axis parameter space according to the dominant outcome for scattered objects: ejected, accreted, remaining, escaping, Oort Cloud, and depleted Oort Cloud. We use these outcomes to consider which planetary system architectures maximize the observability of specific signatures, given that signatures should be detected first around systems with optimal architectures (if such systems exist in nature). Giant impact debris is most readily detectable for 0.1-10 M_⊕ planets at 1-5 au, depending on the detection method and spectral type. While A stars have putative giant impact debris at 4-6 au consistent with this sweet spot, that of FGK stars is typically ≪1 au contrary to expectations; an absence of 1-3 au giant impact debris could indicate a low frequency of terrestrial planets there. Three principles maximize the cometary influx from exo-Kuiper belts: a chain of closely separated planets interior to the belt, none of which is a Jupiter-like ejector; planet masses not increasing strongly with distance (for a net inward torque on comets); and ongoing replenishment of comets, possibly by embedded low-mass planets. A high Oort Cloud comet influx requires no ejectors and architectures that maximize the Oort Cloud population. Cold debris discs are usually considered classical Kuiper belt analogues. Here we consider the possibility of detecting scattered disc analogues, which could be betrayed by a broad radial profile and lack of small grains, as well as spherical 100-1000 au mini-Oort Clouds. Some implications for escaping planets around young stars, detached planets akin to Sedna, and the formation of super-Earths are also discussed.

Abstract

Discs of dusty debris around main-sequence stars indicate fragmentation of orbiting planetesimals, and for a few A-type stars, a gas component is also seen that may come from collisionally-released volatiles. Here we find the sixth example of a CO-hosting disc, around the ∼30 Myr-old A0-star HD 32997. Two more of these CO-hosting stars, HD 21997 and 49 Cet, have also been imaged in dust with SCUBA-2 within the SONS project. A census of 27 A-type debris hosts within 125 pc now shows 7/16 detections of carbon-bearing gas within the 5-50 Myr epoch, with no detections in 11 older systems. Such a prolonged period of high fragmentation rates corresponds quite well to the epoch when most of the Earth was assembled from planetesimal collisions. Recent models propose that collisional products can be spatially asymmetric if they originate at one location in the disc, with CO particularly exhibiting this behaviour as it can photodissociate in less than an orbital period. Of the six CO-hosting systems, only β Pic is in clear support of this hypothesis. However, radiative transfer modelling with the ProDiMo code shows that the CO is also hard to explain in a proto-planetary disc context.

Abstract

Many asteroids are likely rubble-piles that are a collection of smaller objects held together by gravity and possibly cohesion. These asteroids are seismically shaken by impacts, which leads to excitation of their constituent particles. As a result it has been suggested that their surfaces and sub-surface interiors may be governed by a size sorting mechanism known as the Brazil Nut Effect. We study the behavior of a model asteroid that is a spherical, self-gravitating aggregate with a binary size-distribution of particles under the action of applied seismic shaking. We find that above a seismic threshold, larger particles rise to the surface when friction is present, in agreement with previous studies that focussed on cylindrical and rectangular box configurations. Unlike previous works we also find that size sorting takes place even with zero friction, though the presence of friction does aid the sorting process above the seismic threshold. Additionally we find that while strong size sorting can take place near the surface, the innermost regions remain unsorted under even the most vigorous shaking.

Abstract

The Oort cloud is usually thought of as a collection of icy comets inhabiting the outer reaches of the Solar system, but this picture is incomplete. We use simulations of the formation of the Oort cloud to show that ~4 per cent of the small bodies in the Oort cloud should have formed within 2.5 au of the Sun, and hence be ice-free rock-iron bodies. If we assume that these Oort cloud asteroids have the same size distribution as their cometary counterparts, the Large Synoptic Survey Telescope should find roughly a dozen Oort cloud asteroids during 10 years of operations. Measurement of the asteroid fraction within the Oort cloud can serve as an excellent test of the Solar system's formation and dynamical history. Oort cloud asteroids could be of particular concern as impact hazards as their high mass density, high impact velocity, and low visibility make them both hard to detect and hard to divert or destroy. However, they should be a rare class of object, and we estimate globally catastrophic collisions should only occur about once per billion years.

Abstract

Many stars are surrounded by disks of dusty debris formed in the collisions of asteroids, comets, and dwarf planets, but is gas also released in such events? Observations at submillimeter wavelengths of the archetypal debris disk around β Pictoris show that 0.3% of a Moon mass of carbon monoxide orbits in its debris belt. The gas distribution is highly asymmetric, with 30% found in a single clump 85 astronomical units from the star, in a plane closely aligned with the orbit of the inner planet, β Pictoris b. This gas clump delineates a region of enhanced collisions, either from a mean motion resonance with an unseen giant planet or from the remnants of a collision of Mars-mass planets.

Abstract

We consider the observational signatures of giant impacts between planetary embryos. While the debris released in the impact remains in a clump for only a single orbit, there is a much longer lasting asymmetry caused by the fact that all debris must pass through the collision-point. The resulting asymmetry is stationary, it does not orbit the star. The debris is concentrated in a clump at the collision-point, with a more diffuse structure on the opposite side. The asymmetry lasts for typically around 1000 orbital periods of the progenitor, which can be several Myr at distances of ~50 au. We describe how the appearance of the asymmetric disc depends on the mass and eccentricity of the progenitor, as well as viewing orientation. The wavelength of observation, which determines the grain sizes probed, is also important. Notably, the increased collision rate of the debris at the collision-point makes this the dominant production site for any secondary dust and gas created. For dust small enough to be removed by radiation pressure, and gas with a short lifetime, this causes their distribution to resemble a jet emanating from the (stationary) collision-point. We suggest that the asymmetries seen at large separations in some debris discs, like Beta Pictoris, could be the result of giant impacts. If so, this would indicate that planetary embryos are present and continuing to grow at several tens of au at ages of up to tens of Myr.

Abstract

We study the evolution of debris created in the giant impacts expected during the final stages of terrestrial planet formation. The starting point is the debris created in a simulation of the Moon-forming impact. The dynamical evolution is followed for 10 Myr including the effects of Earth, Venus, Mars and Jupiter. The spatial distribution evolves from a clump in the first few months to an asymmetric ring for the first 10 kyr and finally becoming an axisymmetric ring by about 1 Myr after the impact. By 10 Myr after the impact 20 per cent of the particles have been accreted on to Earth and 17 per cent on to Venus, with 8 per cent ejected by Jupiter and other bodies playing minor roles. However, the fate of the debris also depends strongly on how fast it is collisionally depleted, which depends on the poorly constrained size distribution of the impact debris. Assuming that the debris is made up of 30 per cent by mass mm-cm-sized vapour condensates and 70 per cent boulders up to 500 km, we find that the condensates deplete rapidly on ~1000 yr time-scales, whereas the boulders deplete predominantly dynamically. By considering the luminosity of dust produced in collisions within the boulder-debris distribution we find that the Moon-forming impact would have been readily detectable around other stars in Spitzer 24 μm surveys for around 25 Myr after the impact, with levels of emission comparable to many known hot dust systems. The vapour condensates meanwhile produce a short-lived, optically thick, spike of emission. We use these surveys to make an estimate of the fraction of stars that form terrestrial planets, F_TPF. Since current terrestrial planet formation models invoke multiple giant impacts, the low fraction of 10-100 Myr stars found to have warm (>~150 K) dust implies that F_TPF≲10 per cent. For this number to be higher, it would require that either terrestrial planets are largely fully formed when the protoplanetary disc disperses, or that impact generated debris consists purely of sub-km objects such that its signature is short-lived.

Abstract

We consider the evaporation of close-in planets by the star's intrinsic extreme-ultraviolet (EUV) and X-ray radiation. We calculate evaporation rates by solving the hydrodynamical problem for planetary evaporation including heating from both X-ray and EUV radiation. We show that most close-in planets (a < 0.1 au) are evaporating hydrodynamically, with the evaporation occurring in two distinct regimes: X-ray driven, in which the X-ray heated flow contains a sonic point, and EUV driven, in which the X-ray region is entirely sub-sonic. The mass-loss rates scale as L_X/a² for X-ray driven evaporation, and as Φ_*^1/2/a for EUV driven evaporation at early times, with mass-loss rates of the order of 10¹⁰-10¹⁴ g s^-1. No exact scaling exists for the mass-loss rate with planet mass and planet radius; however, in general evaporation proceeds more rapidly for planets with lower densities and higher masses. Furthermore, we find that in general the transition from X-ray driven to EUV driven evaporation occurs at lower X-ray luminosities for planets closer to their parent stars and for planets with lower densities.

Coupling our evaporation models to the evolution of the high-energy radiation - which falls with time - we are able to follow the evolution of evaporating planets. We find that most planets start off evaporating in the X-ray driven regime, but switch to EUV driven once the X-ray luminosity falls below a critical value. The evolution models suggest that while 'hot Jupiters' are evaporating, they are not evaporating at a rate sufficient to remove the entire gaseous envelope on Gyr time-scales. However, we do find that close in Neptune mass planets are more susceptible to complete evaporation of their envelopes. Thus we conclude that planetary evaporation is more important for lower mass planets, particularly those in the 'hot Neptune'/'super Earth' regime.

Abstract

We study the relationship between coronal X-ray emission and stellar age for late-type stars, and the variation of this relationship with spectral type. We select 717 stars from 13 open clusters and find that the ratio of X-ray to bolometric luminosity during the saturated phase of coronal emission decreases from 10^-3.1 for late K dwarfs to 10^-4.3 for early-F-type stars [across the range 0.29 ≤ (B-V)₀ < 1.41]. Our determined saturation time-scales vary between 10^7.8 and 10^8.3 yr, though with no clear trend across the whole FGK range. We apply our X-ray emission-age relations to the investigation of the evaporation history of 121 known transiting exoplanets using a simple energy-limited model of evaporation and taking into consideration Roche lobe effects and different heating/evaporation efficiencies. We confirm that a linear cut-off of the planet distribution in the M2/R3 versus a-2 plane is an expected result of population modification by evaporation and show that the known transiting exoplanets display such a cut-off. We find that for an evaporation efficiency of 25 per cent we expect around one in ten of the known transiting exoplanets to have lost ≥5 per cent of their mass since formation. In addition we provide estimates of the minimum formation mass for which a planet could be expected to survive for 4 Gyr for a range of stellar and planetary parameters. We emphasize the importance of the earliest periods of a planet's life for its evaporation history with 75 per cent expected to occur within the first Gyr. This raises the possibility of using evaporation histories to distinguish between different migration mechanisms. For planets with spin-orbit angles available from measurements of the Rossiter-McLaughlin effect, no difference is found between the distributions of planets with misaligned orbits and those with aligned orbits. This suggests that dynamical effects accounting for misalignment occur early in the life of the planetary system, although additional data are required to test this.