Neutral exospheres of the Moon, Mercury, and several other solar system bodies consist of particles on ballistic trajectories. Here, the vertical density profile of a surface-bounded exosphere is calculated using thermodynamic averages of an ensemble of ballistic trajectories. When the initial velocities follow the so-called "Maxwell-Boltzmann Flux" (MBF) distribution, the classical density profile results. For many other probability distributions, the density approaches infinity near the surface and has a more ground-hugging vertical profile than the classical solution. Even MBF on a rough surface results in a ground-hugging component. Observed vertical density profiles that were interpreted as a superposition of a hot and a cold (ground hugging) population may in fact be consistent with a population at a single temperature.

Recent dynamical analyses suggest that some Jupiter family comets (JFCs) may originate in the main asteroid belt instead of the outer solar system. This possibility is particularly interesting given evidence that icy main-belt objects are known to be present in the Themis asteroid family. We report results from dynamical analyses specifically investigating the possibility that icy Themis family members could contribute to the observed population of JFCs. Numerical integrations show that such dynamical evolution is indeed possible via a combination of eccentricity excitation apparently driven by the nearby 2:1 mean-motion resonance with Jupiter, gravitational interactions with planets other than Jupiter, and the Yarkovsky effect. We estimate that, at any given time, there may be tens of objects from the Themis family on JFC-like orbits with the potential to mimic active JFCs from the outer solar system, although not all, or even any, may necessarily be observably active. We find that dynamically evolved Themis family objects on JFC-like orbits have semimajor axes between 3.15 au and 3.40 au for the vast majority of their time on such orbits, consistent with the strong role that the 2:1 mean-motion resonance with Jupiter likely plays in their dynamical evolution. We conclude that a contribution from the Themis family to the active JFC population is plausible, although further work is needed to better characterize this contribution.

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.

We report a comprehensive analysis of the global spectrophotometric properties of Ceres using Dawn Framing Camera images collected from April to June 2015 during the RC3 and Survey mission phases. The single-scattering albedo of Ceres at 555 nm is 0.14$\pm$0.04, the geometric albedo is 0.096$\pm$0.006, and the Bond albedo is 0.037$\pm$0.002. The asymmetry factors calculated from the best-fit two-term Henyey-Greenstein (HG) single-particle phase function (SPPF) show a wavelength dependence, suggesting that the phase reddening of Ceres is dominated by single-particle scattering rather than multiple scattering or small-scale surface roughness. The Hapke roughness parameter of Ceres is derived to be 20$^\circ\pm$6$^\circ$ with no wavelength dependence. The phase function of Ceres shows appreciably strong scattering around 90$^\circ$ phase angle that cannot be fitted with a single-term HG SPPF, suggesting possible stronger forward scattering than other asteroids previously analyzed with spacecraft data. We speculate that such a scattering characteristic of Ceres might be related to its unique surface composition. We grouped the reflectance data into a 1$^\circ$ latitude-longitude grid and fitted each grid independently to study the spatial variations of photometric properties. The albedo and color maps are consistent with previous studies. The SPPF over the surface of Ceres shows stronger backscattering associated with lower albedo and vice versa, consistent with the general trend among asteroids. The Hapke roughness parameter does not vary much across the surface of Ceres, except for the ancient Vendimia Planitia region that has a slightly higher roughness. Based on the wavelength dependence of the SPPF of Ceres, we hypothesize that its regolith grains either contain a considerable fraction of $\lessapprox\mu$m-sized particles, or are strongly affected by internal scatterers of this size.

We analyze 760475 observations of 333026 main-belt objects obtained by the Pan-STARRS1 (PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discovered two main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 comet detection procedures currently consist of the comparison of the point spread functions (PSFs) of moving objects to those of reference stars, and the flagging of objects that show anomalously large radial PSF widths. Based on the number of missed discovery opportunities among comets discovered by other observers, we estimate an upper limit comet discovery efficiency rate of ~70% for PS1. Additional analyses that could improve comet discovery yields in future surveys include linear PSF analysis, modeling of trailed stellar PSFs for comparison to trailed moving object PSFs, searches for azimuthally localized activity, comparison of point-source-optimized photometry to extended-source-optimized photometry, searches for photometric excesses in objects with known absolute magnitudes, and crowd-sourcing. Analysis of PS1 survey statistics indicates an expected fraction of 59 MBCs per 10^6 outer main-belt asteroids, and a 95% confidence upper limit of 96 MBCs per 10^6 outer main-belt asteroids. We note that more sensitive future surveys could detect many more MBCs than estimated here. We find an excess of high eccentricities (0.1 < e < 0.3) among all known MBCs relative to the background asteroid population. Theoretical calculations show that, given these eccentricities, the sublimation rate for a typical MBC is orders of magnitude larger at perihelion than at aphelion, providing a plausible physical explanation for the observed behavior of MBCs peaking in observed activity strength near perihelion.