Search

We identify a new mechanism that can lead to the destruction of small, close-in planetary satellites. If a small moon close to the planet has a sizable eccentricity and inclination, its ejecta that escape to planetocentric orbit would often re-impact with much higher velocity due to the satellite's and the fragment's orbits precessing out of alignment. If the impacts of returning ejecta result in net erosion, a runaway process can occur which may end in disruption of the satellite, and we term this process ``sesquinary catastrophe''. We expect the moon to re-accrete, but on an orbit with significantly lower eccentricity and inclination. We find that the large majority of small close-in moons in the Solar System, have orbits that are immune to sesquinary catastrophe. The exceptions include a number of resonant moonlets of Saturn for which resonances may affect the velocities of re-impact of their own debris. Additionally, we find that Neptune's moon Naiad (and to a lesser degree, Jupiter's Thebe) must have substantial internal strength, in line with prior estimates based on Roche limit stability. We also find that sesquinary instability puts important constraints on the plausible past orbits of Phobos and Deimos or their progenitors., Comment: Accepted for ApJ

Pathogenic fungi are a leading cause of crop disease and primarily spread through microscopic, durable spores adapted differentially for both persistence and dispersal. Computational Earth System Models and air pollution models have been used to simulate atmospheric spore transport for aerial-dispersal-adapted (airborne) rust diseases, but the importance of atmospheric spore transport for soil-dispersal-adapted (soilborne) diseases remains unknown. This study adapts the Community Atmosphere Model, the atmospheric component of the Community Earth System Model, to simulate the global transport of the plant pathogenic soilborne fungus Fusarium oxysporum, F. oxy. Our sensitivity study assesses the model's accuracy in long-distance aerosol transport and the impact of deposition rate on long-distance spore transport in Summer 2020 during a major dust transport event from Northern Sub-Saharan Africa to the Caribbean and southeastern U.S. We find that decreasing wet and dry deposition rates by an order of magnitude improves representation of long distance, trans-Atlantic dust transport. Simulations also suggest that a small number of viable spores can survive trans-Atlantic transport to be deposited in agricultural zones. This number is dependent on source spore parameterization, which we improved through a literature search to yield a global map of F. oxy spore distribution in source agricultural soils. Using this map and aerosol transport modeling, we show how viable spore numbers in the atmosphere decrease with distance traveled and offer a novel danger index for viable spore deposition in agricultural zones.

NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $\beta$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $\beta$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous., Comment: 23 pages, 13 figures, published in PSJ

We report on the energy dependence of galactic cosmic rays (GCRs) in the very local interstellar medium (VLISM) as measured by the Low Energy Charged Particle (LECP) instrument on the Voyager 1 (V1) spacecraft. The LECP instrument includes a dual-ended solid state detector particle telescope mechanically scanning through 360 deg across eight equally-spaced angular sectors. As reported previously, LECP measurements showed a dramatic increase in GCR intensities for all sectors of the >=211 MeV count rate (CH31) at the V1 heliopause (HP) crossing in 2012, however, since then the count rate data have demonstrated systematic episodes of intensity decrease for particles around 90{\deg} pitch angle. To shed light on the energy dependence of these GCR anisotropies over a wide range of energies, we use V1 LECP count rate and pulse height analyzer (PHA) data from >=211 MeV channel together with lower energy LECP channels. Our analysis shows that while GCR anisotropies are present over a wide range of energies, there is a decreasing trend in the amplitude of second-order anisotropy with increasing energy during anisotropy episodes. A stronger pitch-angle scattering at the higher velocities is argued as a potential cause for this energy dependence. A possible cause for this velocity dependence arising from weak rigidity dependence of the scattering mean free path and resulting velocity-dominated scattering rate is discussed. This interpretation is consistent with a recently reported lack of corresponding GCR electron anisotropies.

A giant impact origin for the Moon is generally accepted, but many aspects of lunar formation remain poorly understood and debated. \'Cuk et al. (2016) proposed that an impact that left the Earth-Moon system with high obliquity and angular momentum could explain the Moon's orbital inclination and isotopic similarity to Earth. In this scenario, instability during the Laplace Plane transition, when the Moon's orbit transitions from the gravitational influence of Earth's figure to that of the Sun, would both lower the system's angular momentum to its present-day value and generate the Moon's orbital inclination. Recently, Tian and Wisdom (2020) discovered new dynamical constraints on the Laplace Plane transition and concluded that the Earth-Moon system could not have evolved from an initial state with high obliquity. Here we demonstrate that the Earth-Moon system with an initially high obliquity can evolve into the present state, and we identify a spin-orbit secular resonance as a key dynamical mechanism in the later stages of the Laplace Plane transition. Some of the simulations by Tian and Wisdom (2020) did not encounter this late secular resonance, as their model suppressed obliquity tides and the resulting inclination damping. Our results demonstrate that a giant impact that left Earth with high angular momentum and high obliquity ($\theta > 61^{\circ}$) is a promising scenario for explaining many properties of the Earth-Moon system, including its angular momentum and obliquity, the geochemistry of Earth and the Moon, and the lunar inclination., Comment: Accepted for the Planetary Science Journal

We investigate what can be learned about a population of distant KBOs by studying the statistical properties of their light curves. Whereas others have successfully inferred the properties of individual, highly variable KBOs, we show that the fraction of KBOs with low amplitudes also provides fundamental information about a population. Each light curve is primarily the result of two factors: shape and orientation. We consider contact binaries and ellipsoidal shapes, with and without flattening. After developing the mathematical framework, we apply it to the existing body of KBO light curve data. Principal conclusions are as follows. (1) When using absolute magnitude H as a proxy for size, it is more accurate to use the maximum of the light curve rather than the mean. (2) Previous investigators have noted that smaller KBOs have higher-amplitude light curves, and have interpreted this as evidence that they are systematically more irregular in shape than larger KBOs; we show that a population of flattened bodies with uniform proportions could also explain this result. (3) Our analysis indicates that prior assessments of the fraction of contact binaries in the Kuiper Belt may be artificially low. (4) The pole orientations of some KBOs can be inferred from observed changes in their light curves; however, these KBOs constitute a biased sample, whose pole orientations are not representative of the population overall. (5) Although surface topography, albedo patterns, limb darkening, and other surface properties can affect individual light curves, they do not have a strong influence on the statistics overall. (6) Photometry from the OSSOS survey is incompatible with previous results and its statistical properties defy easy interpretation. We also discuss the promise of this approach for the analysis of future, much larger data sets such as the one anticipated from the Rubin Observatory.

The D68 ringlet is the innermost feature in Saturn's rings. Four clumps that appeared in D68 around 2014 remained evenly spaced about 30 degrees apart and moved very slowly relative to each other from 2014 up until the last measurements were taken in 2017. D68's narrowness and the distribution of clumps could either indicate that we have a collection of source bodies in a co-orbital configuration or imply that an outside force confines the observed dust and any source bodies. In this paper we explore the possibility that these four clumps arose from four source bodies in a co-orbital configuration. We find that there are no solutions with four masses that produce the observed spacings. We therefore consider whether an unseen fifth co-orbital object could account for the discrepancies in the angular separations and approach a stable stationary configuration. We find a range of solutions for five co-orbital objects where their mass ratios depend on the assumed location of the fifth mass. Numerical simulations of five co-orbitals are highly sensitive to initial conditions, especially for the range of masses we would expect the D68 clumps to have. The fragility of our D68 co-orbital system model implies that there is probably some outside force confining the material in this ringlet., Comment: 7 pages, 7 figures, accepted for publication by the Planetary Science Journal

In its 16 years of scientific measurements, the Spitzer Space Telescope performed a number of ground-breaking infrared measurements of Solar System objects. In this second of two papers, we describe results from Spitzer observations of asteroids, dust rings, and planets that provide new insight into the formation and evolution of our Solar System. The key Spitzer results presented here can be grouped into three broad classes: characterizing the physical properties of asteroids, notably including a large survey of Near Earth Objects; detection and characterization of several dust/debris disks in the Solar System; and comprehensive characterization of ice giant (Uranus, Neptune) atmospheres. Many of these observations provide critical foundations for future infrared space-based observations., Comment: Published in Nature Astronomy, Volume 4, page 940. This arXiv version pre-dates the proofs corrections; the final published version is available at https://rdcu.be/b8fgs [Updated on October 28 to get the right figures in this version]

We present a novel explanation of the well known $M_{\bullet}-\sigma$ relation. In a triaxial potential binaries with chaotic orbits within a sphere that encompass $\sim100$ times the mass of the super-massive black-hole (SMBH) have a finite probability to be tidally disrupted by the SMBH. As a result one component loses energy and might itself break apart tidally and accreted onto the SMBH. More significantly, the other component, which gains energy, returns to the bulge and equilibrates its excess energy with the environment thereby changing the kinetic temperature, hence the velocity dispersion. We develop a mathematical model and find that its results are in agreement with the observed relation., Comment: Comments are welcomed

In this paper, we investigate whether Uranus's 98$^{\circ}$ obliquity was a by-product of a secular spin-orbit resonance assuming that the planet originated closer to the Sun. In this position, Uranus's spin precession frequency is fast enough to resonate with another planet located beyond Saturn. Using numerical integration, we show that resonance capture is possible in a variety of past solar system configurations, but that the timescale required to tilt the planet to 90$^{\circ}$ is of the order $\sim\!10^{8}$ yr -- a timespan that is uncomfortably long. A resonance kick could tilt the planet to a significant 40$^{\circ}$ in $\sim\!10^{7}$ yr only if conditions were ideal. We also revisit the collisional hypothesis for the origin of Uranus's large obliquity. We consider multiple impacts with a new collisional code that builds up a planet by summing the angular momentum imparted from impactors. Because gas accretion imparts an unknown but likely large part of the planet's spin angular momentum, we compare different collisional models for tilted, untilted, spinning, and nonspinning planets. We find that a 1 $M_{\oplus}$ strike is sufficient to explain the planet's current spin state, but that two $0.5\,M_{\oplus}$ collisions produce better likelihoods. Finally, we investigate hybrid models and show that resonances must produce a tilt of at least $\sim\!40^{\circ}$ for any noticeable improvements to the collision model. Because it is difficult for spin-orbit resonances to drive Uranus's obliquity to 98$^{\circ}$ even under these ideal conditions, giant impacts seem inescapable., Comment: 33 pages, 17 figures, accepted to PSJ

The Moon is known to have a small liquid core, and it is thought that in the distant past the core may have produced strong magnetic fields recorded in lunar samples. Here we implement a numerical model of lunar orbital and rotational dynamics that includes the effects of a liquid core. In agreement with previous work, we find that the lunar core is dynamically decoupled from the lunar mantle, and that this decoupling happened very early in lunar history. Our model predicts that the lunar core rotates sub-synchronously, and the difference between the core and the mantle rotational rates was significant when the Moon had a high forced obliquity during and after the Cassini State transition. We find that the presence of the lunar liquid core further destabilizes synchronous rotation of the mantle for a wide range of semimajor axes centered around the Cassini State transition. CMB torques make it even more likely that the Moon experienced large-scale inclination damping during the Cassini State transition. We present estimates for the mutual core-mantle obliquity as a function of Earth-Moon distance, and we discuss plausible absolute time-lines for this evolution. We conclude that our results are consistent with the hypothesis of a precession-driven early lunar dynamo and may explain the variability of the inferred orientation of the past lunar dynamo., Comment: Accepted for JGR:Planets

Giant collisions can account for Uranus's and Neptune's large obliquities, yet generating two planets with widely different tilts and strikingly similar spin rates is a low-probability event. Trapping into a secular spin-orbit resonance, a coupling between spin and orbit precession frequencies, is a promising alternative, as it can tilt the planet without altering its spin period. We show with numerical integrations that if Uranus harbored a massive circumplanetary disk at least three times the mass of its satellite system while it was accreting its gaseous atmosphere, then its spin precession rate would increase enough to resonate with its own orbit, potentially driving the planet's obliquity to 70${^\circ}$. We find that the presence of a massive disk moves the Laplace radius significantly outward from its classical value, resulting in more of the disk contributing to the planet's pole precession. Although we can generate tilts greater than 70${^\circ}$ only rarely and cannot drive tilts beyond 90${^\circ}$, a subsequent collision with an object about $0.5\,M_{\oplus}$ could tilt Uranus from 70${^\circ}$ to 98${^\circ}$. Minimizing the masses and number of giant impactors from two or more to just one increases the likelihood of producing Uranus's spin states by about an order of magnitude. Neptune, by contrast, needs a less massive disk to explain its 30${^\circ}$ tilt, eliminating the need for giant collisions altogether., Comment: 16 pages, 8 figures, accepted for publication in ApJ, minor edits

In the giant impact hypothesis for lunar origin, the Moon accreted from an equatorial circum-terrestrial disk; however the current lunar orbital inclination of 5 degrees requires a subsequent dynamical process that is still debated. In addition, the giant impact theory has been challenged by the Moon's unexpectedly Earth-like isotopic composition. Here, we show that tidal dissipation due to lunar obliquity was an important effect during the Moon's tidal evolution, and the past lunar inclination must have been very large, defying theoretical explanations. We present a new tidal evolution model starting with the Moon in an equatorial orbit around an initially fast-spinning, high-obliquity Earth, which is a probable outcome of giant impacts. Using numerical modeling, we show that the solar perturbations on the Moon's orbit naturally induce a large lunar inclination and remove angular momentum from the Earth-Moon system. Our tidal evolution model supports recent high-angular momentum giant impact scenarios to explain the Moon's isotopic composition and provides a new pathway to reach Earth's climatically favorable low obliquity., Comment: Preprint version of Nature vol. 539, pp-402-406 (2016)

We searched for dust or debris rings in the Pluto-Charon system before, during, and after the New Horizons encounter. Methodologies included searching for back-scattered light during the approach to Pluto (phase $\sim15^\circ$), in situ detection of impacting particles, a search for stellar occultations near the time of closest approach, and by forward-scattered light during departure (phase $\sim165^\circ$). A search using HST prior to the encounter also contributed to the results. No rings, debris, or dust features were observed, but our detection limits provide an improved picture of the environment throughout the Pluto-Charon system. Searches for rings in back-scattered light covered 35,000-250,000 km from the system barycenter, a zone that starts interior to the orbit of Styx, and extends to four times the orbital radius of Hydra. We obtained our firmest limits using the NH LORRI camera in the inner half of this region. Our limits on the normal $I/F$ of an unseen ring depends on the radial scale of the rings: $2\times10^{-8}$ ($3\sigma$) for 1500 km wide rings, $1\times10^{-8}$ for 6000 km rings, and $7\times10^{-9}$ for 12,000 km rings. Beyond $\sim100,000$ km from Pluto, HST observations limit normal $I/F$ to $\sim8\times10^{-8}$. Searches for dust from forward-scattered light extended from the surface of Pluto to the Pluto-Charon Hill sphere ($r_{\rm Hill}=6.4\times10^6$ km). No evidence for rings or dust was detected to normal $I/F$ limits of $\sim8.9\times10^{-7}$ on $\sim10^4$ km scales. Four occulation observations also probed the space interior to Hydra, but again no dust or debris was detected. Elsewhere in the solar system, small moons commonly share their orbits with faint dust rings. Our results suggest that small grains are quickly lost from the system due to solar radiation pressure, whereas larger particles are unstable due to perturbations by the known moons., Comment: Submitted to Icarus, 38 pages, 24 figures

The WASP-18 system, with its massive and extremely close-in planet, WASP-18b (M$_{p}$ = 10.3M$_{J}$, a = 0.02 AU, P = 22.6 hours), is one of the best known exoplanet laboratories to directly measure Q', the modified tidal quality factor and proxy for efficiency of tidal dissipation, of the host star. Previous analysis predicted a rapid orbital decay of the planet toward its host star that should be measurable on the time scale of a few years, if the star is as dissipative as is inferred from the circularization of close-in solar-type binary stars. We have compiled published transit and secondary eclipse timing (as observed by WASP, TRAPPIST, and Spitzer) with more recent unpublished light curves (as observed by TRAPPIST and HST) with coverage spanning nine years. We find no signature of a rapid decay. We conclude that the absence of rapid orbital decay most likely derives from Q' being larger than was inferred from solar-type stars, and find that Q' $\geq $1$\times10^{6}$, at 95\,\% confidence; this supports previous work suggesting that F-stars, with their convective cores and thin convective envelopes, are significantly less tidally dissipative than solar-type stars, with radiative cores and large convective envelopes., Comment: Accepted to ApJL

In order to understand the rate of merger of stellar-mass black hole binaries (BHBs) by gravitational wave (GW) emission it is important to determine the major pathways to merger. We use numerical simulations to explore the evolution of BHBs inside the radius of influence of supermassive black holes (SMBHs) in galactic centers. In this region the evolution of binaries is dominated by perturbations from the central SMBH. In particular, as first pointed out by Antonini and Perets, the Kozai-Lidov (KL) mechanism trades relative inclination of the BHB to the SMBH for eccentricity of the BHB, and for some orientations can bring the BHB to an eccentricity near unity. At very high eccentricities, GW emission from the BHB can become efficient, causing the members of the BHB to coalesce. We use a novel combination of two N-body codes to follow this evolution. We are forced to simulate small systems to follow the behavior accurately. We have completed 400 simulations that range from $\sim$ 300 stars around a $10^{3}$ M$_{\odot}$ black hole to $\sim$ 4500 stars around a $10^{4}$ M$_{\odot}$ black hole. These simulations are the first to follow the internal orbit of a binary near a SMBH while also following the changes to its external orbit self-consistently. We find that this mechanism could produce mergers at a maximum rate per volume of $\sim 100$ Gpc$^{-3}$ yr$^{-1}$ or considerably less if the inclination oscillations of the binary remain constant as the BHB inclination to the SMBH changes, or if the binary black hole fraction is small., Comment: 15 pages, 8 figures, submitted to ApJ

We present observations at optical wavelengths with the Cassini Spacecraft's Imaging Science System of the Phoebe ring, a vast debris disk around Saturn that seems to be collisionally generated by its irregular satellites. The analysis reveals a radial profile from 80-260 Saturn radii ($R_S$) that changes behavior interior to $\approx 110 R_S$, which we attribute to either the moon Iapetus sweeping up small particles, or to orbital instabilities that cause the ring to flare up vertically. Our study yields an integrated I/F at 0.635 $\mu$m along Saturn's shadow in the Phoebe ring's midplane from 80-250 $R_S$ of $2.7^{+0.9}_{-0.3} \times 10^{-9}$. We develop an analytical model for the size-dependent secular dynamics of retrograde Phoebe ring grains, and compare this model to the observations. This analysis implies that 1) the "Phoebe" ring is partially sourced by debris from irregular satellites beyond Phoebe's orbit and 2) the scattered light signal is dominated by small grains ($\lesssim 20\mu$m in size). If we assume that the Phoebe ring is generated through steady-state micrometeoroid bombardment, this implies a power-law size distribution with index $> 4$, which is unusually steep among solar system rings. This suggests either a steep size distribution of ejecta when material is initially released, or a subsequent process that preferentially breaks up large grains., Comment: 17 pages, 14 figures. Accepted for publication in Icarus

Out of nine known stable Mars Trojans, seven appear to be members of an orbital grouping including the largest Trojan, Eureka. In order to test if this could be a genetic family, we simulated the long term evolution of a tight orbital cluster centered on Eureka. We explored two cases: cluster dispersal through planetary gravity alone over 1 Gyr, and a 1 Gyr evolution due to both gravity and the Yarkovsky effect. We find that the dispersal of the cluster in eccentricity is primarily due to dynamical chaos, while the inclinations and libration amplitudes are primarily changed by the Yarkovsky effect. Current distribution of the cluster members orbits are indicative of an initially tight orbital grouping that was affected by a negative acceleration (i.e. one against the orbital motion) consistent with the thermal Yarkovsky effect. We conclude that the cluster is a genetic family formed either in a collision or through multiple rotational fissions. The cluster's age is on the order of 1 Gyr, and its long-term orbital evolution is likely dominated by the seasonal, rather than diurnal, Yarkovsky effect. If confirmed, Gyr-scale dominance of the seasonal Yarkovsky effect may indicate suppression of the diurnal Yarkovsky drift by the related YORP effect. Further study of Mars Trojans is essential for understanding the long-term orbital and rotational dynamics of small bodies in the absence of frequent collisions., Comment: Submitted to Icarus