We present spatially resolved millimeter maps of Neptune between 95 and 242 GHz taken with the Atacama Large Millimeter/submillimeter Array (ALMA) in $2016-2017$. The millimeter weighting functions peak between 1 and 10 bar on Neptune, lying in between the altitudes probed at visible/infrared and centimeter wavelengths. Thus, these observations provide important constraints on the atmospheric structure and dynamics of Neptune. We identify seven well-resolved latitudinal bands of discrete brightness temperature variations, on the order of $0.5-3$K in all three observed ALMA spectral bands. We model Neptune's brightness temperature using the radiative transfer code Radio-BEAR and compare how various H$_2$S, CH$_4$, and \textit{ortho/para} H$_2$ abundance profiles can fit the observed temperature variations across the disk. We find that observed variations in brightness temperature with latitude can be explained by variations in the H$_2$S profile that range from sub- to super-saturations at altitudes above the 10-bar pressure level, while variations in CH$_4$ improve the quality of fit near the equator. At the south polar cap, our best fit model has a depleted deep atmospheric abundance of H$_2$S from 30 to only 1.5 times the protosolar value, while simultaneously depleting the CH$_4$ abundance. This pattern of enhancement and depletion of condensable species is consistent with a global circulation structure where enriched air rises at the mid-latitudes ($32^{\circ}-12^{\circ}$S) and north of the equator ($2^{\circ}-20^{\circ}$N), and dry air descends at the poles ($90^{\circ}-66^{\circ}$S) and just south of the equator ($12^{\circ}$S$-2^{\circ}$N). Our analysis finds more complex structure near the equator than accounted for in previous circulation models.

Unexpected dynamic phenomena have surprised solar system observers in the past and have led to important discoveries about solar system workings. Observations at the initial stages of these events provide crucial information on the physical processes at work. We advocate for long-term/permanent programs on ground-based and space-based telescopes of all sizes - including Extremely Large Telescopes (ELTs) - to conduct observations of high-priority dynamic phenomena, based on a predefined set of triggering conditions. These programs will ensure that the best initial dataset of the triggering event are taken; separate additional observing programs will be required to study the temporal evolution of these phenomena. While not a comprehensive list, the following are notional examples of phenomena that are rare, that cannot be anticipated, and that provide high-impact advances to our understandings of planetary processes. Examples include: new cryovolcanic eruptions or plumes on ocean worlds; impacts on Jupiter, Saturn, Uranus, or Neptune; extreme eruptions on Io; convective superstorms on Saturn, Uranus, or Neptune; collisions within the asteroid belt or other small-body populations; discovery of an interstellar object passing through our solar system (e.g. 'Oumuamua); and responses of planetary atmospheres to major solar flares or coronal mass ejections.

Imaging and resolved spectroscopy reveal varying environmental conditions in our dynamic solar system. Many key advances have focused on how these conditions change over time. Observatory-level commitments to conduct annual observations of solar system bodies would establish a long-term legacy chronicling the evolution of dynamic planetary atmospheres, surfaces, and rings. Science investigations will use these temporal datasets to address potential biosignatures, circulation and evolution of atmospheres from the edge of the habitable zone to the ice giants, orbital dynamics and planetary seismology with ring systems, exchange between components in the planetary system, and the migration and processing of volatiles on icy bodies, including Ocean Worlds. The common factor among these diverse investigations is the need for a very long campaign duration, and temporal sampling at an annual cadence.

We report the discovery of a large ($\sim$8500 km diameter) infrared-bright storm at Neptune's equator in June 2017. We tracked the storm over a period of 7 months with high-cadence infrared snapshot imaging, carried out on 14 nights at the 10 meter Keck II telescope and 17 nights at the Shane 120 inch reflector at Lick Observatory. The cloud feature was larger and more persistent than any equatorial clouds seen before on Neptune, remaining intermittently active from at least 10 June to 31 December 2017. Our Keck and Lick observations were augmented by very high-cadence images from the amateur community, which permitted the determination of accurate drift rates for the cloud feature. Its zonal drift speed was variable from 10 June to at least 25 July, but remained a constant $237.4 \pm 0.2$ m s$^{-1}$ from 30 September until at least 15 November. The pressure of the cloud top was determined from radiative transfer calculations to be 0.3-0.6 bar; this value remained constant over the course of the observations. Multiple cloud break-up events, in which a bright cloud band wrapped around Neptune's equator, were observed over the course of our observations. No "dark spot" vortices were seen near the equator in HST imaging on 6 and 7 October. The size and pressure of the storm are consistent with moist convection or a planetary-scale wave as the energy source of convective upwelling, but more modeling is required to determine the driver of this equatorial disturbance as well as the triggers for and dynamics of the observed cloud break-up events.

Since 2013, observations of Neptune with small telescopes have resulted in several detections of long-lived bright atmospheric features that have also been observed by large telescopes such as Keck II or Hubble. The combination of both types of images allows the study of the long term evolution of major cloud systems in the planet. In 2013 and 2014 two bright features were present on the planet at southern mid latitudes. These may have merged in late 2014, possibly leading to the formation of a single bright feature observed during 2015 at the same latitude. This cloud system was first observed in January 2015 and nearly continuously from July to December 2015 in observations with telescopes in the 2 to 10 meter class and in images from amateur astronomers. These images show the bright spot as a compact feature at 40.1 deg South planetographic latitude well resolved from a nearby bright zonal band that extended from 42 deg South to 20 deg South. Tracking its motion from July to November 2015 suggests a longitudinal oscillation of 16 deg in amplitude with a 90 day period, typical of dark spots on Neptune and similar to the Great Red Spot oscillation in Jupiter. The limited time covered by high-resolution observations only covers one full oscillation and other interpretations of the changing motions could be possible. HST images in September 2015 show the presence of a dark spot at short wavelengths in the southern flank of the bright cloud observed throughout 2015.

We present and analyze three-dimensional data cubes of Neptune from the OSIRIS integral-field spectrograph on the 10-m Keck telescope, from July 2009. These data have a spatial resolution of 0.035"/pixel and spectral resolution of R~3800 in the H and K broad bands. We focus our analysis on regions of Neptune's atmosphere that are near-infrared dark- that is, free of discrete bright cloud features. We use a forward model coupled to a Markov chain Monte Carlo algorithm to retrieve properties of Neptune's aerosol structure and methane profile above ~4 bar in these near-infrared dark regions. Using a set of high signal-to-noise spectra in a cloud-free band from 2-12N, we find that Neptune's cloud opacity is dominated by a compact, optically thick cloud layer with a base near 3 bar and composed of low albedo, forward scattering particles, with an assumed characteristic size of ~1$\mu$m. Above this cloud, we require a vertically extended haze of smaller (~0.1 $\mu$m) particles, which reaches from the upper troposphere (~0.6 bar) into the stratosphere. The particles in this haze are brighter and more isotropically scattering than those in the deep cloud. When we extend our analysis to 18 cloud-free locations from 20N to 87S, we observe that the optical depth in aerosols above 0.5 bar decreases by a factor of 2-3 or more at mid- and high-southern latitudes relative to low latitudes. We also consider Neptune's methane (CH$_4$) profile, and find that our retrievals indicate a strong preference for a low methane relative humidity at pressures where methane is expected to condense. Our preferred solution at most locations is for a methane relative humidity below 10% near the tropopause in addition to methane depletion down to 2.0-2.5 bar. We tentatively identify a trend of lower CH$_4$ columns above 2.5 bar at mid- and high-southern latitudes over low latitudes.

We present the first near infra-red spectrum of the faint white dwarf companion HD 114174 B, obtained with Project 1640. Our spectrum, covering the Y, J and H bands, combined with previous TRENDS photometry measurements, allows us to place further constraints on this companion. We suggest two possible scenarios; either this object is an old, low mass, cool H atmosphere white dwarf with Teff ~ 3800 K or a high mass white dwarf with Teff > 6000 K, potentially with an associated cool (Teff ~ 700 K) brown dwarf or debris disk resulting in an infra-red excess in the L band. We also provide an additional astrometry point for 2014 June 12 and use the modelled companion mass combined with the RV and direct imaging data to place constraints on the orbital parameters for this companion.

The nearby Sun-like star GJ 758 hosts a cold substellar companion, GJ 758 B, at a projected separation of $\lesssim$30 AU, previously detected in high-contrast multi-band photometric observations. In order to better constrain the companion's physical characteristics, we acquired the first low-resolution ($R \sim 50$) near-infrared spectrum of it using the high-contrast hyperspectral imaging instrument Project 1640 on Palomar Observatory's 5-m Hale telescope. We obtained simultaneous images in 32 wavelength channels covering the $Y$, $J$, and $H$ bands ($\sim$952-1770 nm), and used data processing techniques based on principal component analysis to efficiently subtract chromatic background speckle-noise. GJ 758 B was detected in four epochs during 2013 and 2014. Basic astrometric measurements confirm its apparent northwest trajectory relative to the primary star, with no clear signs of orbital curvature. Spectra of SpeX/IRTF observed T dwarfs were compared to the combined spectrum of GJ 758 B, with ${\chi}^2$ minimization suggesting a best fit for spectral type T7.0$\pm$1.0, but with a shallow minimum over T5-T8. Fitting of synthetic spectra from the BT-Settl13 model atmospheres gives an effective temperature $T_{\text{eff}}=741 \pm 25$ K and surface gravity $\log g = 4.3 \pm 0.5$ dex (cgs). Our derived best-fit spectral type and effective temperature from modeling of the low-resolution spectrum suggest a slightly earlier and hotter companion than previous findings from photometric data, but do not rule out current results, and confirm GJ 758 B as one of the coolest sub-stellar companions to a Sun-like star to date.

$\mu$ Her is a nearby quadruple system with a G-subgiant primary and several low mass companions arranged in a 2+2 architecture. While the BC components have been well characterized, the Ab component has been detected astrometrically and with direct imaging but there has been some confusion over its nature, in particular whether the companion is stellar or substellar. Using near-infrared spectroscopy we are able to estimate the spectral type of the companion as a M4$\pm$1V star. In addition, we have measured the astrometry of the system for over a decade. We combined the astrometry with archival radial velocity measurements to compute an orbit of the system. From the combined orbit, we are able to compute the mass sum of the system. Using the estimated mass of the primary, we estimate the mass of the secondary as 0.32 M_sun, which agrees with the estimated spectral type. Our computed orbit is preliminary due to the incomplete orbital phase coverage, but it should be sufficient to predict ephemerides over the next decade.

HD 177830 is an evolved K0IV star with two known exoplanets. In addition to the planetary companions it has a late-type stellar companion discovered with adaptive optics imagery. We observed the binary star system with the PHARO near-IR camera and the Project 1640 coronagraph. Using the Project 1640 coronagraph and integral field spectrograph we extracted a spectrum of the stellar companion. This allowed us to determine that the spectral type of the stellar companion is a M4$\pm$1V. We used both instruments to measure the astrometry of the binary system. Combining these data with published data, we determined that the binary star has a likely period of approximately 800 years with a semi-major axis of 100-200 AU. This implies that the stellar companion has had little or no impact on the dynamics of the exoplanets. The astrometry of the system should continue to be monitored, but due to the slow nature of the system, observations can be made once every 5-10 years.

We present an analysis of the orbital motion of the four sub-stellar objects orbiting HR8799. Our study relies on the published astrometric history of this system augmented with an epoch obtained with the Project 1640 coronagraph + Integral Field Spectrograph (IFS) installed at the Palomar Hale telescope. We first focus on the intricacies associated with astrometric estimation using the combination of an Extreme Adaptive Optics system (PALM-3000), a coronagraph and an IFS. We introduce two new algorithms. The first one retrieves the stellar focal plane position when the star is occulted by a coronagraphic stop. The second one yields precise astrometric and spectro-photometric estimates of faint point sources even when they are initially buried in the speckle noise. The second part of our paper is devoted to studying orbital motion in this system. In order to complement the orbital architectures discussed in the literature, we determine an ensemble of likely Keplerian orbits for HR8799bcde, using a Bayesian analysis with maximally vague priors regarding the overall configuration of the system. While the astrometric history is currently too scarce to formally rule out coplanarity, HR8799d appears to be misaligned with respect to the most likely planes of HR8799bce orbits. This misalignment is sufficient to question the strictly coplanar assumption made by various authors when identifying a Laplace resonance as a potential architecture. Finally, we establish a high likelihood that HR8799de have dynamical masses below 13 M_Jup using a loose dynamical survival argument based on geometric close encounters. We illustrate how future dynamical analyses will further constrain dynamical masses in the entire system.

We report observations of Neptune made in H-(1.4-1.8 {\mu}m) and K'-(2.0-2.4 {\mu}m) bands on 14 and 16 July 2009 from the 10-m W.M. Keck II Telescope using the near-infrared camera NIRC2 coupled to the Adaptive Optics (AO) system. We track the positions of 54 bright atmospheric features over a few hours to derive their zonal and latitudinal velocities, and perform radiative transfer modeling to measure the cloud-top pressures of 50 features seen simultaneously in both bands. We observe one South Polar Feature (SPF) on 14 July and three SPFs on 16 July at ~65 deg S. The SPFs observed on both nights are different features, consistent with the high variability of Neptune's storms. There is significant dispersion in Neptune's zonal wind velocities about the smooth Voyager wind profile fit of Sromovsky et al., Icarus 105, 140 (1993), much greater than the upper limit we expect from vertical wind shear, with the largest dispersion seen at equatorial and southern mid-latitudes. Comparison of feature pressures vs. residuals in zonal velocity from the smooth Voyager wind profile also directly reveals the dominance of mechanisms over vertical wind shear in causing dispersion in the zonal winds. Vertical wind shear is not the primary cause of the difference in dispersion and deviation in zonal velocities between features tracked in H-band on 14 July and those tracked in K'-band on 16 July. Dispersion in the zonal velocities of features tracked over these short time periods is dominated by one or more mechanisms, other than vertical wind shear, that can cause changes in the dispersion and deviation in the zonal velocities on timescales of hours to days.

We present maps of Neptune in and near the CO (2-1) rotation line at 230.538 GHz. These data, taken with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) represent the first published spatially-resolved maps in the millimeter. At large (~5 GHz) offsets from the CO line center, the majority of the emission originates from depths of 1.1-4.7 bar. We observe a latitudinal gradient in the brightness temperature at these frequencies, increasing by 2-3 K from 40 degrees N to the south pole. This corresponds to a decrease in the gas opacity of about 30% near the south pole at altitudes below 1 bar, or a decrease of order a factor of 50 in the gas opacity at pressures greater than 4 bar. We look at three potential causes of the observed gradient: variations in the tropospheric methane abundance, variations in the H2S abundance, and deviations from equilibrium in the ortho/para ratio of hydrogen. At smaller offsets (0-213 MHz) from the center of the CO line, lower atmospheric pressures are probed, with contributions from mbar levels down to several bars. We find evidence of latitudinal variations at the 2-3% level in the CO line, which are consistent with the variations in zonal-mean temperature near the tropopause found by Conrath et al. (1998) and Orton et al. (2007).

We present observations of Neptune's 1- and 3-mm spectrum from the Combined Array for Research in Millimeter-wave Astronomy (CARMA). Radiative transfer analysis of the CO (2-1) and (1-0) rotation lines was performed to constrain the CO vertical abundance profile. We find that the data are well matched by a CO mole fraction of 0.1^+0.2_-0.1 parts per million (ppm) in the troposphere, and 1.1^+0.2_-0.3 ppm in the stratosphere. A flux of 0.5-20 times 10^8 CO molecules cm-2 s-1 to the upper stratosphere is implied. Using the Zahnle et al. (2003) estimate for cometary impact rates at Neptune, we calculate the CO flux that could be formed from (sub)kilometer-sized comets; we find that if the diffusion rate near the tropopause is small (200 cm2 s-1), these impacts could produce a flux as high as 0.5^+0.8_-0.4 times 10^8 CO molecules cm-2 s-1. We also revisit the calculation of Neptune's internal CO contribution using revised calculations for the CO ->CH4 conversion timescale in the deep atmosphere (Visscher et al. 2011). We find that an upwelled CO mole fraction of 0.1 ppm implies a global O/H enrichment of at least 400, and likely more than 650, times the protosolar value.

Keck near-infrared images of Neptune from UT 26 July 2007 show that the cloud feature typically observed within a few degrees of Neptune's south pole had split into a pair of bright spots. A careful determination of disk center places the cloud centers at -89.07 +/- 0 .06 and -87.84 +/- 0.06 degrees planetocentric latitude. If modeled as optically thick, perfectly reflecting layers, we find the pair of features to be constrained to the troposphere, at pressures greater than 0.4 bar. By UT 28 July 2007, images with comparable resolution reveal only a single feature near the south pole. The changing morphology of these circumpolar clouds suggests they may form in a region of strong convection surrounding a Neptunian south polar vortex.