The authors and co-signers of the Terrestrial Planets Comparative Climatology (TPCC) mission concept white paper advocate that planetary science in the next decade would greatly benefit from comparatively studying the fundamental behavior of the atmospheres of Venus and Mars, contemporaneously and with the same instrumentation, to capture atmospheric response to the same solar forcing, and with a minimum of instrument-related variability. Thus, this white paper was created for the 2023-2032 Planetary Science Decadal Survey process. It describes the science rationale for such a mission, and a mission concept that could achieve such a mission.

Robust atmospheric and radiative transfer modeling will be required to properly interpret reflected light and thermal emission spectra of terrestrial exoplanets. This will help break observational degeneracies between the numerous atmospheric, planetary, and stellar factors that drive planetary climate. Here we simulate the climates of Earth-like worlds around the Sun with increasingly slow rotation periods, from Earth-like to fully Sun-synchronous, using the ROCKE-3D general circulation model. We then provide these results as input to the Spectral Planet Model (SPM), which employs the SMART radiative transfer model to simulate the spectra of a planet as it would be observed from a future space-based telescope. We find that the primary observable effects of slowing planetary rotation rate are the altered cloud distributions, altitudes, and opacities which subsequently drive many changes to the spectra by altering the absorption band depths of biologically-relevant gas species (e.g., H2O, O2, and O3). We also identify a potentially diagnostic feature of synchronously rotating worlds in mid-infrared H2O absorption/emission lines.

This paper constrains the scattering phase function of water ice clouds (WICs) found within Mars' Aphelion Cloud Belt (ACB), determined from orbit by processing publicly available raw Mars Color Imager (MARCI) data spanning solar longitudes (Ls) 42-170 during Mars Years (MYs) 28 and 29. MARCI visible wavelength data were calibrated and then pipeline-processed to select the pixels most likely to possess clouds. Mean phase function curves for the MARCI blue filter data were derived, and for all seasons investigated, modeled aggregates, plates, solid and hollow columns, bullet rosettes, and droxtals were all found to be plausible habits. Spheres were found to be the least plausible, but still possible. Additionally, this work probed the opposition surge to examine the slope of the linear relationship between column ice water content and cloud opacity on Mars, and found a significant dependence on particle radius. The half-width-half-maxima (HWHM) of the visible 180 degree peak of five MARCI images were found to agree better with modeled HWHM for WICs than with modeled HWHM for dust.

Northern line-of-sight extinction within Gale Crater during the 2018 global dust storm was monitored daily using MSL's Navcam. Additional observations with Mastcam (north) and Navcam (all directions) were obtained at a lower cadence. Using feature identification and geo-referencing, extinction was estimated in all possible directions. Peak extinction of $>1.1$ km$^{-1}$ was measured between sols 2086 and 2090, an order of magnitude higher than previous maxima. Northern and western directions show an initial decrease, followed by a secondary peak in extinction, not seen in column opacity measurements. Due to foreground topography, eastern direction results are provided only as limits, and southern results were indeterminable. Mastcam red and green filter results agree well, but blue filter results show higher extinctions, likely due to low signal-to-noise. Morning results are systematically higher than afternoon results, potentially indicative of atmospheric mixing.

While recently discovered exotic new planet-types have both challenged our imaginations and broadened our knowledge of planetary system workings, perhaps the most compelling objective of exoplanet science is to detect and characterize habitable and possibly inhabited worlds orbiting in other star systems. For the foreseeable future, characterizations of extrasolar planets will be made via remote sensing of planetary spectroscopic and temporal signals, along with careful fitting of this data to advanced models of planets and their atmospheres. Terrestrial planets are small and significantly more challenging to observe compared to their larger gaseous brethren; however observatories coming on-line in the coming decade will begin to allow their characterization. Still, it is not enough to invest only in observational endeavors. Comprehensive modeling of planetary atmospheres is required in order to fully understand what it is that our grand telescopes see in the night-sky. In our quest to characterize habitable, and possibly inhabited worlds, 3D general circulation models (GCMs) should be used to evaluate potential climate states and their associated temporal and spatial dependent observable signals. 3D models allow for coupled, self-consistent, multi-dimensional simulations, which can realistically simulate the climates of terrestrial extrasolar planets. A complete theoretical understanding of terrestrial exoplanetary atmospheres, gained through comprehensive 3D modeling, is critical for interpreting spectra of exoplanets taken from current and planned instruments, and is critical for designing future missions that aim to measure spectra of potentially habitable exoplanets as one of their key science goals. We recommend continued institutional support for 3D GCM modeling teams that focus on planetary and exoplanetary applications.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

This study constrains the lower bound of the scattering phase function of Martian water ice clouds (WICs) through the implementation of a new observation aboard the Mars Science Laboratory (MSL). The Phase Function Sky Survey (PFSS) was a multiple pointing all-sky observation taken with the navigation cameras (Navcam) aboard MSL. The PFSS was executed 35 times during the Aphelion Cloud Belt (ACB) season of Mars Year 34 over a solar longitude range of L_s=61.4{\deg}-156.5{\deg}. Twenty observations occurred in the morning hours between 06:00 and 09:30 LTST, and 15 runs occurred in the evening hours between 14:30 and 18:00 LTST, with an operationally required 2.5 hour gap on either side of local noon due the sun being located near zenith. The resultant WIC phase function was derived over an observed scattering angle range of 18.3{\deg} to 152.61{\deg}, normalized, and compared with 9 modeled phase functions: seven ice crystal habits and two Martian WIC phase functions currently being implemented in models. Through statistical chi-squared probability tests, the five most probable ice crystal geometries observed in the ACB WICs were aggregates, hexagonal solid columns, hollow columns, plates, and bullet rosettes with p-values greater than or equal to 0.60, 0.57,0.56,0.56, and 0.55, respectively. Droxtals and spheres had p-values of 0.35, and 0.2, making them less probable components of Martian WICs, but still statistically possible ones. Having a better understanding of the ice crystal habit and phase function of Martian water ice clouds directly benefits Martian climate models which currently assume spherical and cylindrical particles.

The potential habitability of known exoplanets is often categorized by a nominal equilibrium temperature assuming a Bond albedo of either 0.3, similar to Earth, or 0. As an indicator of habitability, this leaves much to be desired, because albedo on other planets can be very different, and because surface temperature exceeds equilibrium temperature due to the atmospheric greenhouse effect. We use an ensemble of 3-dimensional general circulation model simulations to show that for a range of habitable planets, much of the variability of Bond albedo, equilibrium temperature, and even surface temperature can be predicted with useful accuracy from incident stellar flux and stellar temperature, two known external parameters for every confirmed exoplanet.

Highly volcanic exoplanets, which can be variously characterized as 'lava worlds', 'magma ocean worlds', or 'super-Ios' are high priority targets for investigation. The term 'lava world' may refer to any planet with extensive surface lava lakes, while the term 'magma ocean world' refers to planets with global or hemispherical magma oceans at their surface. 'Highly volcanic planets', including super-Ios, may simply have large, or large numbers of, active explosive or extrusive volcanoes of any form. They are plausibly highly diverse, with magmatic processes across a wide range of compositions, temperatures, activity rates, volcanic eruption styles, and background gravitational force magnitudes. Worlds in all these classes are likely to be the most characterizable rocky exoplanets in the near future due to observational advantages that stem from their preferential occurrence in short orbital periods and their bright day-side flux in the infrared. Transit techniques should enable a level of characterization of these worlds analogous to hot Jupiters. Understanding processes on highly volcanic worlds is critical to interpret imminent observations. The physical states of these worlds are likely to inform not just geodynamic processes, but also planet formation, and phenomena crucial to habitability. Volcanic and magmatic activity uniquely allows chemical investigation of otherwise spectroscopically inaccessible interior compositions. These worlds will be vital to assess the degree to which planetary interior element abundances compare to their stellar hosts, and may also offer pathways to study both the very young Earth, and the very early form of many silicate planets where magma oceans and surface lava lakes are expected to be more prevalent. We suggest that highly volcanic worlds may become second only to habitable worlds in terms of both scientific and public long-term interest.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

Penitentes are ablative features observed in snow and ice that, on Earth, are characterized by regular cm to tens of cm spaced bowl-shaped depressions whose edges grade into tall spires up to several meters tall (Nichols et al. 1939, Lliboutry et al. 1954, Claudin et al. 2015). While penitentes have been suggested as an explanation for anomalous radar data on Europa (Hobley et al. 2013), hitherto no penitentes have been identified conclusively on other planetary bodies. Regular ridges with spacing of 3000 m to 5000 m and a depth of $\sim500$ m with morphologies that resemble penitentes (Fig. 1) have been observed by the New Horizons spacecraft (Moore et al. 2016, Stern et al. 2015, Gladstone et al. 2016, Moore et al. 2017) in the Tartarus Dorsa (TD) region of Pluto (approximately 220-250$^\circ$ E, 0-20$^\circ$ N). Here we report simulations, based upon a recent model (Claudinet al. 2015) adapted to conditions on Pluto (Gladstone et al. 2016, Toigo et al. 2015), that reproduce both the tri-modal orientation and the spacing of these features by deepening penitentes. These penitentes deepen by of order 1 cm per orbital cycle in the present era and grow only during periods of relatively high atmospheric pressure, suggesting a formation timescale of several tens of millions of years, consistent with cratering ages and the current atmospheric loss rate of methane. This time scale, in turn, implies that the penitentes formed from initial topographic variations of no more than a few 10s of meters, consistent with Pluto's youngest terrains.

Those that do not sow care little about such mundane things as equinoxes or planting seasons, or even crop rotation for that matter. Wherever and whenever the reavers reave, the mood is always foul and the nights are never warm or pleasant. For the rest of the good folks of Westeros, however, a decent grasp of the long-term weather forecast is a necessity. Many a maester have tried to play the Game of Weather Patterns and foretell when to plant those last turnip seeds, hoping for a few more years of balmy respite. Tried and failed. For other than the somewhat vague (if not outright meaningless) omens of "Winter is Coming", their meteorological efforts have been worse than useless. To right that appalling wrong, here we attempt to explain the apparently erratic seasonal changes in the world of G.R.R.M. A natural explanation for such phenomena is the unique behavior of a circumbinary planet. Thus, by speculating that the planet under scrutiny is orbiting a pair of stars, we utilize the power of numerical three-body dynamics to predict that, unfortunately, it is not possible to predict either the length, or the severity of any coming winter. We conclude that, alas, the Maesters were right -- one can only throw their hands in the air in frustration and, defeated by non-analytic solutions, mumble "Coming winter? May be long and nasty (~850 days, T<268K) or may be short and sweet (~600 days, T~273K). Who knows..."

The 2013 Planetary Science Decadal Survey identified a detailed investigation of the Trojan asteroids occupying Jupiter's L4 and L5 Lagrange points as a priority for future NASA missions. Observing these asteroids and measuring their physical characteristics and composition would aid in identification of their source and provide answers about their likely impact history and evolution, thus yielding information about the makeup and dynamics of the early Solar System. We present a conceptual design for a mission to the Jovian Trojan asteroids: the Trojan ASteroid Tour, Exploration, and Rendezvous (TASTER) mission, that is consistent with the NASA New Frontiers candidate mission recommended by the Decadal Survey and the final result of the 2011 NASA-JPL Planetary Science Summer School. Our proposed mission includes visits to two Trojans in the L4 population: a 500 km altitude fly-by of 1999 XS143, followed by a rendezvous with and detailed observations of 911 Agamemnon at orbital altitudes of 1000 - 100 km over a 12 month nominal science data capture period. Our proposed instrument payload - wide- and narrow-angle cameras, a visual and infrared mapping spectrometer, and a neutron/gamma ray spectrometer - would provide unprecedented high-resolution, regional-to-global datasets for the target bodies, yielding fundamental information about the early history and evolution of the Solar System. Although our mission design was completed as part of an academic exercise, this study serves as a useful starting point for future Trojan mission design studies. In particular, we identify and discuss key issues that can make large differences in the complex trade-offs required when designing a mission to the Trojan asteroids.