We present the Strategic Scientific Plan (SSP) for the direction and activities of the Gemini Observatory in the 2020s. The overarching goal is to ensure that Gemini best uses the available resources to serve the needs of its international user community throughout the coming decade. The actionable items fall into three general categories: (1) preserving Gemini's current facilities and strengths; (2) developing instrumentation and software systems, including data pipelines, to enable new scientific capabilities that build on those strengths; (3) strategizing how visiting instruments can deliver additional valuable capabilities. We provide a high-level timeline (schematically illustrated in one figure) for the main developments discussed in this SSP. The schedule is ambitious, but in light of the recent Gemini in the Era of Multi-Messenger Astronomy (GEMMA) award from the NSF, the plan becomes achievable. Lists of milestones are given for gauging progress. As these milestones are reached and new instruments become available, some current instruments will need to be retired; we make recommendations in this regard. The final section concludes by reemphasizing the importance of a strong partnership committed to the needs of all members.

As part of the National Science Foundation funded "Gemini in the Era of MultiMessenger Astronomy" (GEMMA) program, Gemini Observatory is developing GNAO, a widefield adaptive optics (AO) facility for Gemini-North on Maunakea, the only 8m-class open-access telescope available to the US astronomers in the northern hemisphere. GNAO will provide the user community with a queue-operated Multi-Conjugate AO (MCAO) system, enabling a wide range of innovative solar system, Galactic, and extragalactic science with a particular focus on synergies with JWST in the area of time-domain astronomy. The GNAO effort builds on institutional investment and experience with the more limited block-scheduled Gemini Multi-Conjugate System (GeMS), commissioned at Gemini South in 2013. The project involves close partnerships with the community through the recently established Gemini AO Working Group and the GNAO Science Team, as well as external instrument teams. The modular design of GNAO will enable a planned upgrade to a Ground Layer AO (GLAO) mode when combined with an Adaptive Secondary Mirror (ASM). By enhancing the natural seeing by an expected factor of two, GLAO will vastly improve Gemini North's observing efficiency for seeing-limited instruments and strengthen its survey capabilities for multi-messenger astronomy.

A combination band due to a mechanism whereby a photon excites two or more vibrational modes ({\it e.g.} a bend and a stretch) of an individual molecule is commonly seen in laboratory and astronomical spectroscopy. Here, we present evidence of a much less commonly seen combination band $-$ one where a photon simultaneously excites two adjacent molecules in an ice. In particular, we present near-infrared spectra of laboratory CO/N$_2$ ice samples where we identify a band at 4467.5 cm$^{-1}$ (2.239 $\mu$m) that results from single photons exciting adjacent pairs of CO and N$_2$ molecules. We also present a near-infrared spectrum of Neptune's largest satellite Triton taken with the Gemini-South 8.1 meter telescope and the Immersion Grating Infrared Spectrograph (IGRINS) that shows this 4467.5 cm$^{-1}$ (2.239 $\mu$m) CO-N$_2$ combination band. The existence of the band in a spectrum of Triton indicates that CO and N$_2$ molecules are intimately mixed in the ice rather than existing as separate regions of pure CO and pure N$_2$ deposits. Our finding is important because CO and N$_2$ are the most volatile species on Triton and so dominate seasonal volatile transport across its surface. Our result will place constraints on the interaction between the surface and atmosphere of Triton. 1

Two groundbreaking new facilities will commence operations early in the 2020s and thereafter define much of the broad landscape of US optical-infrared astronomy in the remaining decade. The Large Synoptic Survey Telescope (LSST), perched atop Cerro Pachon in the Chilean Andes, will revolutionize the young field of Time Domain Astronomy through its wide-field, multi-band optical imaging survey. At the same time, the James Webb Space Telescope (JWST), orbiting at the Sun-Earth L2 Lagrange point, will provide stunningly high-resolution views of selected targets from the red end of the optical spectrum to the mid-infrared. However, the spatial resolution of the LSST observations will be limited by atmospheric seeing, while JWST will be limited in its time-domain capabilities. This paper highlights the scientific opportunities lying between these two landmark missions, i.e., science enabled by systems capable of astronomical observations with both high cadence in the time domain and high resolution in the spatial domain. The opportunities range from constraining the late phases of stellar evolution in nearby resolved populations to constraining dark matter distributions and cosmology using lensed transient sources. We describe a system that can deliver the required capabilities.

In this work we investigated the longitudinal (zonal) variability of H$_2$O and ammonia (NH$_3$) hydrate ices on the surface of Charon through analysis of the 1.65 $\mu$m and 2.21 $\mu$m absorption features, respectively. Near-infrared spectra presented here were obtained between 2015-07-14 and 2015-08-30 UT with the OSIRIS integral field spectrograph on Keck I. Spectra centered on six different sub-observer longitudes were obtained through the Hbb (1.473-1.803 $\mu$m) and Kbb (1.965-2.381 $\mu$m) filters. Gaussian functions were fit to the aforementioned bands to obtain information on band center, band depth, full width at half maximum, and band area. The shift in the band center of the temperature-dependent 1.65 $\mu$m feature was used to calculate the H$_2$O ice temperature. The mean temperature of the ice on the observable portion of Charon's surface is 45$\pm$14 K and we report no statistically significant variations in temperature across the surface. We hypothesize that the crystalline and amorphous phases of water ice reached equilibrium over 3.5 Gyr ago, with thermal recrystallization balancing the effects of irradiation amorphization. We do not believe that cryovolcanism is necessary to explain the presence of crystalline water ice on the surface of Charon. Absorption from ammonia species is detected between 12$^{\circ}$ and 290$^{\circ}$, in agreement with results from New Horizons. Ongoing diffusion of ammonia through the rocky mantle and upper layer of water ice is one possible mechanism for maintaining its presence in Charon's surface ice. Reduced Charon spectra corrected for telluric and solar absorption are available as supplementary online material.

From observations by the Hubble Space Telescope, Keck II Telescope, and Gemini North Telescope, we have determined the mutual orbit of the large transneptunian object (174567) Varda and its satellite Ilmar\"e. These two objects orbit one another in a highly inclined, circular or near-circular orbit with a period of 5.75 days and a semimajor axis of 4810 km. This orbit reveals the system mass to be (2.664 +/- 0.064) x 10^20 kg, slightly greater than the mass of the second most massive main-belt asteroid (4) Vesta. The dynamical mass can in turn be combined with estimates of the surface area of the system from Herschel Space Telescope thermal observations to estimate a bulk density of 1.24 +0.50 -0.35 g cm^-3. Varda and Ilmar\"e both have colors similar to the combined colors of the system, B-V = 0.886 +/- 0.025 and V-I = 1.156 +/- 0.029.

We present three improved and five new mutual orbits of transneptunian binary systems (58534) Logos-Zoe, (66652) Borasisi-Pabu, (88611) Teharonhiawako-Sawiskera, (123509) 2000 WK183, (149780) Altjira, 2001 QY297, 2003 QW111, and 2003 QY90 based on Hubble Space Telescope and Keck 2 laser guide star adaptive optics observations. Combining the five new orbit solutions with 17 previously known orbits yields a sample of 22 mutual orbits for which the period P, semimajor axis a, and eccentricity e have been determined. These orbits have mutual periods ranging from 5 to over 800 days, semimajor axes ranging from 1,600 to 37,000 km, eccentricities ranging from 0 to 0.8, and system masses ranging from 2 x 10^17 to 2 x 10^22 kg. Based on the relative brightnesses of primaries and secondaries, most of these systems consist of near equal-sized pairs, although a few of the most massive systems are more lopsided. The observed distribution of orbital properties suggests that the most loosely-bound transneptunian binary systems are only found on dynamically cold heliocentric orbits. Of the 22 known binary mutual or-bits, orientation ambiguities are now resolved for 9, of which 7 are prograde and 2 are retro-grade, consistent with a random distribution of orbital orientations, but not with models predicting a strong preference for retrograde orbits. To the extent that other perturbations are not dominant, the binary systems undergo Kozai oscillations of their eccentricities and inclinations with periods of the order of tens of thousands to millions of years, some with strikingly high amplitudes.

Images from instruments on Cassini as well as from telescopes on the ground reveal the presence of sporadic small-scale cloud activity in the cold late-winter north polar of Saturn's large moon Titan. These clouds lie underneath the previously discovered uniform polar cloud attributed to a quiescent ethane cloud at ~40 km and appear confined to the same latitudes as those of the largest known hydrocarbon lakes at the north pole of Titan. The physical properties of these clouds suggest that they are due to methane convection and condensation. Such convection has not been predicted for the cold winter pole, but can be caused by a process in many ways analogous to terrestrial lake-effect clouds. The lakes on Titan are a key connection between the surface and the meteorological cycle.

We report a multi-week sequence of B-band photometric measurements of the dwarf planet Eris using the {\it Swift} satellite. The use of an observatory in low-Earth orbit provides better temporal sampling than is available with a ground-based telescope. We find no compelling evidence for an unusually slow rotation period of multiple days, as has been suggested previously. A $\sim$1.08 day rotation period is marginally detected at a modest level of statistical confidence ($\sim$97%). Analysis of the combination of the $Swift$ data with the ground-based B-band measurements of \citet{2007AJ....133...26R} returns the same period ($\sim$1.08 day) at a slightly higher statistical confidence ($\sim$99%).

We report observations of the ro-vibronic transition of SO at 1.707 microns on Io. These data were taken while Io was eclipsed by Jupiter, on four nights between July 2000 and March 2003. We analyze these results in conjunction with a previously published night to investigate the temporal behavior of these emissions. The observations were all conducted using the near-infrared spectrometer NIRSPEC on the W.M. Keck II telescope. The integrated emitted intensity for this band varies from 0.8 x 10^27 to 2.4 x 10^27 photons/sec, with a possible link to variations in Loki's infrared brightness. The band-shapes imply rotational temperatures of 550-1000K for the emitting gas, lending further evidence to a volcanic origin for sulfur monoxide. An attempt to detect the ro-vibronic transition of SO at 0.97 microns was unsuccessful; simultaneous detection with the 1.707 micron band would permit determination of the SO column abundance.

We have used the Hubble Space Telescope to directly measure the angular size of the large Kuiper belt object 2003 UB313. By carefully calibrating the point spread function of a nearby field star, we measure the size of 2003 UB313 to be 34.3$\pm$1.4 milliarcseconds, corresponding to a diameter of 2400$\pm$100 km or a size $\sim5$% larger than Pluto. The V band geometric albedo of 2003 UB313 is $86\pm7$%. The extremely high albedo is consistent with the frosty methane spectrum, the lack of red coloring, and the lack of observed photometric variation on the surface of 2003 UB313. Methane photolysis should quickly darken the surface of 2003 UB313, but continuous evaporation and redeposition of surface ices appears capable of maintaining the extreme alebdo of this body.

We present measurements at optical wavelengths of the spectral reflectance, rotational light curve, and solar phase curve of 2003 EL61. With apparent visual magnitude 17.5 at 51 AU from the sun, this newly discovered member of the classical Kuiper Belt is now the third brightest KBO after Pluto and 2005 FY9. Our observations reveal an unambiguous, double-peaked rotational light curve with period 3.9154 +/- 0.0002 hours and peak to peak amplitude 0.28 +/- 0.04 mag. This is the fastest rotation period reliably determined for any body in the solar system larger than 100 km. Assuming the body has relaxed over time to the shape taken by a homogenous fluid body, our observations tightly constrain the shape and density. Given the mass we recently determined for 2003 EL61 from the orbit of a small satellite, we also constrain the size and albedo. We find a total length of 1960 to 2500 km, a mean density of 2600 to 3340 kg m-3, and a visual albedo greater than 0.6. We also measure a neutral reflectance at visible wavelengths and a linear phase curve with slope varying from 0.09 mag deg-1 in the B band to 0.13 mag deg-1 in the I band. The absolute V-band magnitude is 0.444+/-0.021.

We present the first observations of propane (C$_3$H$_8$) on Titan that unambiguously resolve propane features from other numerous stratospheric emissions. This is accomplished using a $R=\lambda/\delta\lambda\approx10^5$ spectrometer (TEXES) to observe propane's $\nu_{26}$ rotation-vibration band near 748 cm$^{-1}$. We find a best-fit fractional abundance of propane in Titan's stratosphere of $(6.2\pm1.2)\times10^{-7}$ in the altitude range to which we are sensitive (90-250 km or 13-0.24 mbar).

Many adaptive optics systems operate by measuring the distortion of the wavefront in one wavelength range and performing the scientific observations in a second, different wavelength range. One common technique is to measure wavefront distortions at wavelengths <~1 micron while operating the science instrument at wavelengths >~1 micron. The index of refraction of air decreases sharply from shorter visible wavelengths to near-infrared wavelengths. Therefore, because the adaptive optics system is measuring the wavefront distortion in one wavelength range and the science observations are performed at a different wavelength range, residual image motion occurs and the maximum exposure time before smearing of the image can be significantly limited. We demonstrate the importance of atmospheric differential refraction, present calculations to predict the effect of atmospheric differential refraction, and finally discuss the implications of atmospheric differential refraction for several current and proposed observatories.

Using adaptive optics on the W.M. Keck II telescope we imaged Titan several times during 1999 to 2001 in narrowband near-infrared filters selected to probe Titan's stratosphere and upper troposphere. We observed a bright feature around the south pole, possibly a collar of haze or clouds. Further, we find that solar phase angle explains most of the observed east-west brightness asymmetry of Titan's atmosphere, although the data do not preclude the presence of a `morning fog' effect at small solar phase angle.

We present 2.03-2.30 micron near-infrared spectroscopy of Neptune taken 1999 June 2 (UT) with the W.M. Keck Observatory's near-infrared spectrometer (NIRSPEC) during the commissioning of the instrument. The spectrum is dominated by a bright cloud feature, possibly a storm or upwelling, in the southern hemisphere at approximately 50 degrees S latitude. The spectrum also includes light from a dimmer northern feature at approximately 30 degrees N latitude. We compare our spectra (R ~ 2000) of these two features with a simple model of Neptune's atmosphere. Given our model assumption that the clouds are flat reflecting layers, we find that the top of the bright southern cloud feature sat at a pressure level of 0.14 (+0.05, -0.03) bar, and thus this cloud did not extend into the stratosphere (P < 0.1 bar). A similar analysis of the dimmer northern feature gives a cloud-top pressure of 0.084 +/- 0.026 bar. This suggests that the features we observed efficiently transport methane to the base of the stratosphere, but do not directly transport methane to the upper stratosphere (P < 0.01-0.001 bar) where photolysis occurs. Our observations do not constrain how far these clouds penetrate down into the troposphere. We find that our model fits to the data restrict the fraction of molecular hydrogen in ortho/para thermodynamic equilibrium to greater than 0.8.

New abundances for neutron-capture (n-capture) elements in a large sample of metal-poor giants from the Bond survey are presented. The spectra were acquired with the KPNO 4-m echelle and coude feed spectrographs, and have been analyzed using LTE fine-analysis techniques with both line analysis and spectral synthesis. Abundances of eight n-capture elements (Sr, Y, Zr, Ba, La, Nd, Eu, Dy) in 43 stars have been derived from blue (lambda = 4070--4710, R~20,000, S/N ratio~100-200) echelle spectra and red (lambda = 6100--6180, R~22,000, S/N ratio~100-200) coude spectra, and the abundance of Ba only has been derived from the red spectra for an additional 27 stars. Overall, the abundances show clear evidence for a large star-to-star dispersion in the heavy element-to-iron ratios. The new data also confirm that at metallicities [Fe/H] <~ --2.4, the abundance pattern of the heavy (Z >= 56) n-capture elements in most giants is well-matched to a scaled Solar System r-process nucleosynthesis pattern. The onset of the main r-process can be seen at [Fe/H] ~ --2.9. Contributions from the s-process can first be seen in some stars with metallicities as low as [Fe/H] ~ --2.75, and are present in most stars with metallicities [Fe/H] > --2.3. The lighter n-capture elements (Sr-Y-Zr) are enhanced relative to the heavier r-process element abundances. Their production cannot be attributed solely to any combination of the Solar System r- and main s-processes, but requires a mixture of material from the r-process and from an additional n-capture process which can operate at early Galactic time.