Water ice is thought to be trapped in large permanently shadowed regions (PSRs) in the Moon's polar regions, due to their extremely low temperatures. Here, we show that many unmapped cold traps exist on small spatial scales, substantially augmenting the areas where ice may accumulate. Using theoretical models and data from the Lunar Reconnaissance Orbiter, we estimate the contribution of shadows on scales from 1 km down to 1 cm, the smallest distance over which we find cold-trapping to be effective for water ice. Approximately 10-20\% of the permanent cold trap area for water is found to be contained in these "micro cold traps," which are the most numerous cold traps on the Moon. Consideration of all spatial scales therefore substantially increases the number of cold traps over previous estimates, for a total area of ~40,000 km^2. A majority of cold traps for water ice is found at latitudes >80{\deg} because permanent shadows equatorward of 80{\deg} are typically too warm to support ice accumulation. Our results show that water trapped at the lunar poles may be more accessible as a resource for future missions than previously thought.

Asteroid pairs had a single progenitor that split due to rotational-fission of a weak, rubble-pile structured body. By constructing shape models of asteroid pairs from multiple-apparition observations and using a lightcurve inversion technique, we mapped the gravitational and rotational accelerations on the surfaces of these asteroids. This allows us to construct a map of local slopes on the asteroids' surfaces. In order to test for frictional failure, we determine the maximum rotation rate at which an area larger than half the surface area of the secondary member (assumed to be the ejected component) has a slope value greater than 40 degrees, the angle of friction of lunar regolith, where loose material will begin sliding. We use this criterion to constrain the failure stress operating on the body, just before disruption at the commonly observed spin barrier of 2.2 h. Our current sample includes shape models of eleven primary members of asteroid pairs, observed from the Wise Observatory in the last decade. In the studied parameter space we find that the shape models only reach the spin barrier when their bulk density is larger than the ~2 gr cm-3 measured for the rubble pile structured 25143 Itokawa, suggesting that km-sized asteroid pairs are dense compared to sub-km bodies. Assuming ejection of secondary components that are larger than those observed (up to the maximal size allowing separation), can also increase the spin barrier of the asteroids, thus supporting the previously suggested scenario of continuous disruption of the secondary. In addition, cohesion levels of hundreds of Pascals are also required to prevent these shape models from disrupting at spin rates slower than the usual spin barrier.

Current lunar origin scenarios suggest that Earth's Moon may have resulted from the merger of two (or more) smaller moonlets. Dynamical studies of multiple moons find that these satellite systems are not stable, resulting in moonlet collision or loss of one or more of the moonlets. We perform Smoothed Particle Hydrodynamic (SPH) impact simulations of two orbiting moonlets inside the planetary gravitational potential and find that the classical outcome of two bodies impacting in free space is altered as erosive mass loss is more significant with decreasing distance to the planet. Depending on the conditions of accretion, each moonlet could have a distinct isotopic signature, therefore, we assess the initial mixing during their merger, in order to estimate whether future measurements of surface variations could distinguish between lunar origin scenarios (single vs. multiple moonlets). We find that for comparable-size impacting bodies in the accretionary regime, surface mixing is efficient, but in the hit-and-run regime, only a small amount of material is transferred between the bodies. However, sequences of hit-and-run impacts are expected, which will enhance the surface mixing. Overall, our results show that large scale heterogeneities can arise only from the merger of drastically different component masses. Surfaces of moons resulting from the merger of comparable-sized components have little material heterogeneities, and such impacts are preferred, as the relatively massive impactor generates more melt, extending the lunar magma ocean phase.

The hypothesis of lunar origin by a single giant impact can explain some aspects of the Earth-Moon system. However, it is difficult to reconcile giant impact models with the compositional similarity of the Earth and Moon without violating angular momentum constraints. Furthermore, successful giant impact scenarios require very specific conditions such that they have a low probability of occurring. Here we present numerical simulations suggesting that the Moon could instead be the product of a succession of a variety of smaller collisions. In this scenario, each collision forms a debris disk around the proto-Earth that then accretes to form a moonlet. The moonlets tidally advance outward, and may coalesce to form the Moon. We find that sub-lunar moonlets are a common result of impacts expected onto the proto-Earth in the early solar system and find that the planetary rotation is limited by impact angular momentum drain. We conclude that, assuming efficient merger of moonlets, a multiple impact scenario can account for the formation of the Earth-Moon system with its present properties.

The layered polar caps of Mars have long been thought to be related to variations in orbit and axial tilt. We dynamically link Mars's past climate variations with the stratigraphy and isotopic composition of its ice by modeling the exchange of H2O and HDO among three reservoirs. The model shows that the interplay among equatorial, mid-latitude, and north-polar layered deposits (NPLD) induces significant isotopic changes in the cap. The diffusive properties of the sublimation lags and dust content in our model result in a cap size consistent with current Mars. The layer thicknesses are mostly controlled by obliquity variations, but the precession period of 50 kyr dominates the variations in the isotopic composition during epochs of relatively low and nearly constant obliquity such as at present. Isotopic sampling of the top 100 meters may reveal climate oscillations unseen in the layer thicknesses and would thus probe recent precession-driven climate cycles.

The Earth-Moon system is suggested to have formed through a single giant collision, in which the Moon accreted from the impact-generated debris disk. However, such giant impacts are rare, and during its evolution the Earth experienced many more smaller impacts, producing smaller satellites that potentially coevolved. In the multiple-impact hypothesis of lunar formation, the current Moon was produced from the mergers of several smaller satellites (moonlets), each formed from debris disks produced by successive large impacts. In the Myrs between impacts, a pre-existing moonlet tidally evolves outward until a subsequent impact forms a new moonlet, at which point both moonlets will tidally evolve until a merger or system disruption. In this work, we examine the likelihood that pre-existing moonlets survive subsequent impact events, and explore the dynamics of Earth-moonlet systems that contain two moonlets generated Myrs apart. We demonstrate that pre-existing moonlets can tidally migrate outward, remain stable during subsequent impacts, and later merge with newly created moonlets (or re-collide with the Earth). Formation of the Moon from the mergers of several moonlets could therefore be a natural byproduct of the Earth's growth through multiple impacts. More generally, we examine the likelihood and consequences of Earth having prior moons, and find that the stability of moonlets against disruption by subsequent impacts implies that several large impacts could post-date Moon formation.

Unraveling the stratigraphic record is the key to understanding ancient climate and past climate changes on Mars. River deposits when placed in stratigraphic order could constrain the number, magnitudes, and durations of the wettest climates in Mars history. We establish the stratigraphic context of river deposits in Aeolis Dorsa sedimentary basin, 10E of Gale crater. Here, wind has exhumed a stratigraphic section of >=4 unconformity-bounded sedimentary rock packages, recording >=3 distinct episodes of surface runoff. Early deposits (>700m thick) are embayed by river deposits (>400m), which are in turn unconformably draped by fan-shaped deposits (<100m) which we interpret as alluvial fans. Yardang-forming deposits (>900 m) unconformably drape all previous deposits. River deposits embay a dissected sedimentary-rock landscape, and comprise >=2 distinguishable units. The total interval spanned by river deposits is >(1x10^6-2x10^7) yr; more if we include alluvial-fan deposits. Alluvial-fan deposits unconformably postdate thrust faults which crosscut river deposits. We infer a relatively dry interval of >4x10^7 yr after river deposits formed and before fan-shaped deposits formed. The time gap between the end of river deposition and the onset of yardang-forming deposits is constrained to >10^8 yr by the density of impact craters embedded at the unconformity. We correlate yardang-forming deposits to the upper layers of Gale crater's mound (Mt. Sharp/Aeolis Mons), and fan-shaped deposits to Peace Vallis fan. Alternations between periods of low vs. high mean obliquity may have modulated erosion-deposition cycling in Aeolis. This is consistent with results from an ensemble of simulations of Solar System orbital evolution and the resulting history of Mars obliquity. Almost all simulations yield intervals of continuously low mean Mars obliquity that are long enough to match our unconformity data.

The decay of the martian atmosphere - which is dominated by carbon dioxide - is a component of the long-term environmental change on Mars from a climate that once allowed rivers to flow to the cold and dry conditions of today. The minimum size of craters serves as a proxy for palaeopressure of planetary atmospheres, because thinner atmospheres permit smaller objects to reach the surface at high velocities and form craters. The Aeolis Dorsa region near Gale crater on Mars contains a high density of preserved ancient craters interbedded with river deposits and thus can provide constraints on atmospheric density around the time of fluvial activity. Here we use high-resolution orthophotos and digital terrain models from the Mars Reconnaissance Orbiter to identify ancient craters in Aeolis Dorsa that date to about 3.6 Gyr ago and compare their size distribution with models of atmospheric filtering of impactors. We obtain an upper limit of 0.9$\pm$0.1 bar, rising to 1.9$\pm$0.2 bar if rimmed circular mesas - interpreted to be erosionally-resistant fills of floors of impact craters - are excluded. We assume target properties appropriate for desert alluvium: if sediment had rock-like rock-mass strength similar to bedrock at the time of impact, the upper limit increases by a factor of up to two. If Mars did not have a stable multibar atmosphere at the time that the rivers were flowing - as suggested by our results - then the warm and wet CO2/H2O greenhouse is ruled out, and long-term average temperatures were most likely below freezing.

We used infrared data from the Lunar Reconnaissance Orbiter (LRO) Diviner Lunar Radiometer Experiment to globally map thermophysical properties of the Moon's regolith fines layer. Thermal conductivity varies from 7.4$\times$10$^{-4}$ W m$^{-1}$ K$^{-1}$ at the surface, to 3.4$\times$10$^{-3}$ W m$^{-1}$ K$^{-1}$ at depths of $\sim$1 m, given density values of 1100 kg m$^{-3}$ at the surface, to 1800 kg m$^{-3}$ at 1-m depth. On average, the scale height of these profiles is $\sim$7 cm, corresponding to a thermal inertia of 55 $\pm$2 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$ at 273 K, relevant to the diurnally active near-surface layer, $\sim$4-7 cm. The temperature-dependence of thermal conductivity and heat capacity leads to a $\sim$2$\times$ diurnal variation in thermal inertia at the equator. On global scales, the regolith fines are remarkably uniform, implying rapid homogenization by impact gardening of this layer on timescales $<$ 1 Gyr. Regional and local scale variations show prominent impact features $<$ 1 Gyr old, including higher thermal inertia ($>$ 100 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) in the interiors and ejecta of Copernican-aged impact craters, and lower thermal inertia ($<$ 50 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) within the lunar cold spots identified by Bandfield et al. (2014). Observed trends in ejecta thermal inertia provide a potential tool for age-dating craters of previously unknown age, complementary to the approach suggested by Ghent et al. (2014). Several anomalous regions are identified in the global 128 pixels-per-degree maps presented here, including a high-thermal inertia deposit near the antipode of Tycho crater.

The high planetary multiplicity revealed by Kepler implies that Transit Time Variations (TTVs) are intrinsically common. The usual procedure for detecting these TTVs is biased to long-period, deep transit planets whereas most transiting planets have short periods and shallow transits. Here we introduce the Spectral Approach to TTVs technique that allows expanding the TTVs catalog towards lower TTV amplitude, shorter orbital period, and shallower transit depth. In the Spectral Approach we assume that a sinusoidal TTV exists in the data and then calculate the improvement to $\chi^2$ this model allows over that of linear ephemeris model. This enables detection of TTVs even in cases where the transits are too shallow so individual transits cannot be timed. The Spectral Approach is more sensitive due to the reduced number of free parameters in its model. Using the Spectral Approach, we: (a) detect 131 new periodic TTVs in Kepler data (an increase of ~2/3 over a previous TTV catalog); (b) Constrain the TTV periods of 34 long-period TTVs and reduce amplitude errors of known TTVs; (c) Identify cases of multi-periodic TTVs, for which absolute planetary mass determination may be possible. We further extend our analysis by using perturbation theory assuming small TTV amplitude at the detection stage, which greatly speeds up our detection (to a level of few seconds per star). Our extended TTVs sample shows no deficit of short period or low amplitude transits, in contrast to previous surveys in which the detection schemes were significantly biased against such systems.

Seven of the nine known Mars Trojan asteroids belong to an orbital cluster named after its largest member 5261 Eureka. Eureka is likely the progenitor of the whole cluster, which formed at least 1 Gyr ago. It was suggested that the thermal YORP effect spun-up Eureka resulting with fragments being ejected by the rotational-fission mechanism. Eureka's spectrum exhibits a broad and deep absorption band around 1 {\mu}m, indicating an olivine-rich composition. Here we show evidence that the Trojan Eureka cluster progenitor could have originated as impact debris excavated from the Martian mantle. We present new near-infrared observations of two Trojans (311999 2007 NS2 and 385250 2001 DH47) and find that both exhibit an olivine-rich reflectance spectrum similar to Eureka's. These measurements confirm that the progenitor of the cluster has an achondritic composition. Olivine-rich reflectance spectra are rare amongst asteroids but are seen around the largest basins on Mars. They are also consistent with some Martian meteorites (e.g. Chassigny), and with the material comprising much of the Martian mantle. Using numerical simulations, we show that the Mars Trojans are more likely to be impact ejecta from Mars than captured olivine-rich asteroids transported from the main belt. This result links directly specific asteroids to debris from the forming planets.

At a mean diameter of ~650 m, the near-Earth asteroid (455213) 2001 OE84 (OE84 for short) has a rapid rotation period of 0.486542+-0.000002 hours, which is uncommon for asteroids larger than ~200 m. We revisited OE84 14 years after it was first, and last, observed by Pravec et al. (2002) in order to measure again its spin rate and to search for changes. We have confirmed the rapid rotation and, by fitting the photometric data from 2001 and 2016 using the lightcurve inversion technique, we determined a retrograde sense of rotation, with the spin axis close to the ecliptic south pole; an oblate shape model of a/b=1.32+-0.04 and b/c=1.8+-0.2; and no change in spin rate between 2001 and 2016. Using these parameters we constrained the body's internal strength, and found that current estimations of asteroid cohesion (up to ~80 Pa) are insufficient to maintain an intact rubble pile at the measured spin rate of OE84. Therefore, we argue that a monolithic asteroid, that can rotate at the rate of OE84 without shedding mass and without slowing down its spin rate, is the most plausible for OE84, and we give constraints on its age, since the time it was liberated from its parent body, between 2-10 million years.

The rubble pile spin barrier is an upper limit on the rotation rate of asteroids larger than ~200-300 m. Among thousands of asteroids with diameters larger than ~300 m, only a handful of asteroids are known to rotate faster than 2.0 h, all are in the sub-km range (<=0.6 km). Here we present photometric measurements suggesting that (60716) 2000 GD65, an S-complex, inner-main belt asteroid with a relatively large diameter of 2.3 +0.6-0.7 km, completes one rotation in 1.9529+-0.0002 h. Its unique diameter and rotation period allow us to examine scenarios about asteroid internal structure and evolution: a rubble pile bound only by gravity; a rubble-pile with strong cohesion; a monolithic structure; an asteroid experiencing mass shedding; an asteroid experiencing YORP spin-up/down; and an asteroid with a unique octahedron shape results with a four-peak lightcurve and a 3.9 h period. We find that the most likely scenario includes a lunar-like cohesion that can prevent (60716) 2000 GD65 from disrupting without requiring a monolithic structure or a unique shape. Due to the uniqueness of (60716) 2000 GD65, we suggest that most asteroids typically have smaller cohesion than that of lunar regolith.

We model the primary crater production of small (D < 100 m) primary craters on Mars and the Moon using the observed annual flux of terrestrial fireballs. From the size-frequency distribution (SFD) of meteor diameters, with appropriate velocity distributions for Mars and the Moon, we are able to reproduce martian and lunar crater-count chronometry systems (isochrons) in both slope and magnitude. We include an atmospheric model for Mars that accounts for the deceleration, ablation, and fragmentation of meteors. We find that the details of the atmosphere or the fragmentation of the meteors do not strongly influence our results. The downturn in the crater SFD from atmospheric filtering is predicted to occur at D ~ 10-20 cm, well below the downturn observed in the distribution of fresh craters detected by the Mars Global Surveyor (MGS) Mars Orbiter Camera (MOC) or the Mars Reconnaissance Orbiter (MRO) Context Camera (CTX). Crater counts conducted on the ejecta blanket of Zunil crater on Mars and North Ray crater on the Moon yielded crater SFDs with similar slopes and ages (~1 Ma, and ~58 Ma, respectively) to our model, indicating that the average cratering rate has been constant on these bodies over these time periods. Since our Monte Carlo simulations demonstrate that the existing crater chronology systems can be applied to date young surfaces using small craters on the Moon and Mars, we conclude that the signal from secondary craters in the isochrons must be relatively small, as our Monte Carlo model only generates primary craters.

Cometary activity in main-belt asteroids probes the ice content of these objects and provides clues to the history of volatiles in the inner solar system. We search the Palomar Transient Factory (PTF) survey to derive upper limits on the population size of active main-belt comets (MBCs). From data collected March 2009 through July 2012, we extracted 2 million observations of 220 thousand known main-belt objects (40% of the known population, down to 1-km diameter) and discovered 626 new objects in multi-night linked detections. We formally quantify the extendedness of a small-body observation, account for systematic variation in this metric (e.g., due to on-sky motion) and evaluate this method's robustness in identifying cometary activity using observations of 115 comets, including two known candidate MBCs and six newly-discovered non-main-belt comets (two of which were originally designated as asteroids by other surveys). We demonstrate a 66% detection efficiency with respect to the extendedness distribution of the 115 sampled comets, and a 100% detection efficiency with respect to extendedness levels greater than or equal to those we observed in the known candidate MBCs P/2010 R2 (La Sagra) and P/2006 VW139. Using a log-constant prior, we infer 95% confidence upper limits of 33 and 22 active MBCs (per million main-belt asteroids down to 1-km diameter), for detection efficiencies of 66% and 100%, respectively. In a follow-up to this morphological search, we will perform a photometric (disk-integrated brightening) search for MBCs.

The time-variable electromagnetic sky has been well-explored at a wide range of wavelengths. Numerous high-energy space missions take advantage of the dark Gamma-ray and X-ray sky and utilize very wide field detectors to provide almost continuous monitoring of the entire celestial sphere. In visible light, new wide-field ground-based surveys cover wide patches of sky with ever decreasing cadence, progressing from monthly-weekly time scale surveys to sub-night sampling. In the radio, new powerful instrumentation offers unprecedented sensitivity over wide fields of view, with pathfinder experiments for even more ambitious programs underway. In contrast, the ultra-violet (UV) variable sky is relatively poorly explored, even though it offers exciting scientific prospects. Here, we review the potential scientific impact of a wide-field UV survey on the study of explosive and other transient events, as well as known classes of variable objects, such as active galactic nuclei and variable stars. We quantify our predictions using a fiducial set of observational parameters which are similar to those envisaged for the proposed ULTRASAT mission. We show that such a mission would be able to revolutionize our knowledge about massive star explosions by measuring the early UV emission from hundreds of events, revealing key physical parameters of the exploding progenitor stars. Such a mission would also detect the UV emission from many tens of tidal-disruption events of stars by super massive black holes at galactic nuclei and enable a measurement of the rate of such events. The overlap of such a wide-field UV mission with existing and planned gravitational-wave and high-energy neutrino telescopes makes it especially timely.

The Palomar Transient Factory (PTF) is a synoptic survey designed to explore the transient and variable sky in a wide variety of cadences. We use PTF observations of fields that were observed multiple times (>=10) per night, for several nights, to find asteroids, construct their lightcurves and measure their rotation periods. Here we describe the pipeline we use to achieve these goals and present the results from the first four (overlapping) PTF fields analyzed as part of this program. These fields, which cover an area of 21 deg^2, were observed on four nights with a cadence of ~20 min. Our pipeline was able to detect 624 asteroids, of which 145 (~20%) were previously unknown. We present high quality rotation periods for 88 main-belt asteroids and possible period or lower limit on the period for an additional 85 asteroids. For the remaining 451 asteroids, we present lower limits on their photometric amplitudes. Three of the asteroids have lightcurves that are characteristic of binary asteroids. We estimate that implementing our search for all existing high-cadence PTF data will provide rotation periods for about 10,000 asteroids mainly in the magnitude range ~14 to ~20.