Atacama Large Aperture Submillimeter Telescope (AtLAST) science: Planetary and cometary atmospheres

The study of planets and small bodies within our Solar System is fundamental for understanding the formation and evolution of the Earth and other planets. Compositional and meteorological studies of the giant planets provide a foundation for understanding the nature of the most commonly observed exoplanets, while spectroscopic observations of the atmospheres of terrestrial planets, moons, and comets provide insights into the past and present-day habitability of planetary environments, and the availability of the chemical ingredients for life. While prior and existing (sub)millimeter observations have led to major advances in these areas, progress is hindered by limitations in the dynamic range, spatial and temporal coverage, as well as sensitivity of existing telescopes and interferometers. Here, we summarize some of the key planetary science use cases that factor into the design of the Atacama Large Aperture Submillimeter Telescope (AtLAST), a proposed 50-m class single dish facility: (1) to more fully characterize planetary wind fields and atmospheric thermal structures, (2) to measure the compositions of icy moon atmospheres and plumes, (3) to obtain detections of new, astrobiologically relevant gases and perform isotopic surveys of comets, and (4) to perform synergistic, temporally-resolved measurements in support of dedicated interplanetary space missions. The improved spatial coverage (several arcminutes), resolution (~ 1.2′′ − 12′′), bandwidth (several tens of GHz), dynamic range (~ 10 5) and sensitivity (~ 1 mK km s −1) required by these science cases would enable new insights into the chemistry and physics of planetary environments, the origins of prebiotic molecules and the habitability of planetary systems in general.


Introduction
Since the latter half of the 20th century, (sub)millimeter (∼ 0.3 − 3 mm; 100-950 GHz) telescope facilities have been increasingly employed for the study of planetary surfaces, atmospheres, and ring/moon systems.Brightness temperature measurements can be readily derived for large solar system objects using traditional single-dish radio facilities, while the high spectral resolution combined with increasing bandwidth and sensitivity of modern (sub)millimeter wave telescopes and interferometers has enabled pioneering spectroscopic surveys and spatialspectral mapping of planetary bodies and comets.The measurement of thermal emission in the shallow subsurfaces of natural satellites, small bodies, and terrestrial planets fills the gap between radar sounding of the deep (>1 km) surface and the near surface temperatures derived from infrared observations.The (sub)millimeter wavelength regime also provides access to rotational transitions from gas-phase molecules in planetary atmospheres, the observation of which leads to robust, spectroscopic molecular detections as well as abundance, temperature, and wind measurements with well-defined, and often high, degrees of quantitative accuracy.The available molecular species tend to be dominated by the more abundant reactive elements in our solar system (C, H, N, O, P and S), and therefore include hydrocarbons, oxides, hydrides, nitriles and other members of the class of organic molecules associated with terrestrial biology.The Atacama Large Aperture Submillimeter Telescope (AtLAST) is proposed as a more powerful, ground-based, single-dish submillimeter facility to begin operating in the 2030s [83,120,110,111].Combining improved sensitivity, angular resolution and field-of-view, AtLAST would open new frontiers in our understanding of planetary at-mospheres and surfaces, as well as providing more detailed characterization of the gases produced by comets.More sensitive access to the long-wavelength region of the spectrum is also crucial for providing supporting observations and large-scale context for the more focused measurements performed by interplanetary spacecraft missions.In particular, the availability of improved mapping capabilities, wider instantaneous bandwidth, and higher spatial resolution made possible by AtLAST would result in improvements in our ability to characterize and monitor dynamic and complex physical phenomena, such as winds and storms, as well as chemical processes (revealed by rotational spectroscopy), occurring in the atmospheres and on the surfaces of planetary bodies.The study of 'local' planets in our Solar System also provides a fundamental baseline for interpreting the properties of the vast number of exoplanets that are now known to exist around other stars throughout the galaxy [144].
The paradigm within which we understand the environmental conditions, climate and possible habitability of exoplanets exists by virtue of the active development of an extensive knowledge-base on atmospheric (and bulk) chemical inventories and physical processes occurring in the planets, moons, and other minor bodies of our Solar System.In this era of rapid discovery, it is thus more important than ever to develop and build new telescope facilities such as AtLAST, to continue to advance our understanding of the compositions, spatial and temporal features, and evolutionary processes occurring on the bodies within our Solar System.

Prior/existing instruments and results
Single-dish radio/submillimeter facilities provide a wealth of unique and complementary observations to spacecraft and facilities at other wavelengths.The Owens Valley Radio Observatory (OVRO), Institut de Radioastronomie Millimétrique (IRAM) 30-m telescope, Green Bank Observatory (GBO), Arecibo radar facility, among others, were previously used to study the brightness temperatures of the Giant planets, Mars, Venus, and large satellites.They enabled study of the deep atmospheric composition of Gas and Ice Giants, characterization of the near subsurface properties of the Galilean Satellites, as well as compositional and dynamical studies of terrestrial planet and satellite atmospheres.These were complemented by early interferometric observations by the Very Large Array (VLA), Berkeley-Illinois-Maryland Association (BIMA) array, Submillimeter Array (SMA), and Atacama Large Millimeter/submillimeter Array (ALMA), enabling higher sensitivity and angular resolution studies [55,56,112,2,86,142,99,73,104,59,57,53,54].Cometary science -the study of our Solar System's oldest, yet most pristine materials -has benefited hugely from advances in single-dish millimeter-wave instrumentation.The first detections of key organics HCN and CH 3 OH in cometary comae were obtained through observations using the National Radio Astronomy Observatory (NRAO) 11-m and IRAM 30-m telescopes in the 3 mm and 2 mm bands, respectively [80,16].Millimeter-wave spectroscopy continues to provide the primary method for remotely detecting and characterizing new molecules (including complex organic molecules) in cometary comae [13], while also providing crucial insights into the origins of our solar system's primitive materials via detailed studies of molecular isotopic ratios [75,14].ALMA is currently revolutionizing the study of comets at (sub)millimeter wavelengths, revealing the coma and nucleus outgassing sources in unprecedented detail from the ground [46,49,52], yet snapshot maps of the largest coma scales remain out of reach with present facilities.A new, highly-sensitive single-dish (sub)millimeter telescope with wide-field imaging capabilities will help advance coma mapping studies and help break new ground in our understanding of the chemical composition of comets.Observations at (sub)millimeter wavelengths have enabled the detection and mapping of new molecular species in the atmospheres of various planets in our solar system.CO and HCN were first detected on Neptune in the early 1990s using the James Clerk Maxwell Telescope (JCMT) and CalTech Submillimeter Observatory (CSO) telescopes [98], with abundances ∼ 1000 times higher than predicted from thermochemical models.Similarly, CO, HCN, and HNC were detected on Pluto with recent ALMA observations [88], and a variety of trace gases have been discovered on Titan using millimeterwave spectroscopy [11,47,42,118,115,133].Sub-millimeter observations can also be used to probe gas and ice giant interiors.For example, Neptune's very high observed CO abundance in the stratosphere and troposphere, combined with thermochemical models, has been used to argue for an ice-dominated interior [91,84,23].However, a comet impact in the last few hundred years, as suggested by (sub)millimeter observations of stratospheric CO [84] and CS [102], would imply significant external flux into Neptune's atmosphere.Therefore, a rock dominated interior is also plausible [126] and more consistent with D/H ratio measurements if the interior is well mixed [63].The Herschel Space Observatory provided numerous demonstrations of the power of (sub)millimeter singledish heterodyne spectroscopy for characterizing planetary atmospheres, producing the first detection of water in the extended gaseous torus around Enceladus [74].The Herschel Heterodyne Instrument for the Far-Infrared (HIFI) instrument also obtained the first detection of hydrogen isocyanide (HNC) in Titan's atmosphere [105].Stratospheric H 2 O measurements with HIFI have been used to demonstrate the cometary origin of water in Jupiter's stratosphere [28] and to constrain the background flux of interplanetary dust particles into Uranus' and Neptune's atmospheres [127].Fletcher et al. [65] used the Spectral and Photometric Imaging Receiver (SPIRE) instrument on Herschel to constrain the vertical distributions of phosphine, ammonia and methane on Saturn, as well as providing measurements of stratospheric water and precise upper limits for a range of exotic compounds including halides.In general, the improvements in (sub)millimeter wave sensitivity and mapping capabilities ushered in by ALMA have led to major strides in our understanding of the compositions and dynamical states of planetary atmospheres, as well as satellite and small body surfaces [44,48,41,43,40,131,132,130,133,129,128,102,82,59,53,54,134,135,31,32,8,9,19,20].However, the lack of simultaneous, high-sensitivity total power observations hinders ALMA's ability to characterize localised/dynamical phenomena such as storms and plumes on rotating planetary surfaces, due to temporal smearing.A new, single-dish telescope with improved total power sensitivity and widefield mapping capabilities is therefore required.Instantaneous (snapshot) mapping over the entire area of a body is critical for detailed studies of the rapidly rotating and evolving atmospheres of solar system objects.Such studies have historically been the realm of interferometry, though the lack of short baselines in a typical interferometric array leads to spatial filtering of the resulting images, which can preclude reliable measurements of extended emission from larger bodies such as the giant planets and extended cometary comae.Further, the repositioning of antennas to larger configurations (such as the extended baselines of ALMA) may impede mapping of the giant planets during transient events, such as impacts and storms: impact events tell us both about the population, bulk properties, and volatile composition of potential impactors in the outer Solar System, while storms provide us access to the chemical composition deep below the topmost clouds, and reveal how dramatic meteorology evolves in non-terrestrial and cold environments.Combined with their large angular scales (Figure 1) and total flux, the planets present challenging targets for observation of small brightness temperature variations or weak spectral lines against a strong background continuum at (sub)millimeter wavelengths, and therefore place stringent demands on the stability, calibration, and dynamic range of heterodyne telescope instruments.The availability of a new, cutting edge single-dish (sub)millimeter facility at a dry, high-altitude site will help to fill this gap, allowing us to address important outstanding questions in planetary and cometary science as outlined in this article.
The Need for a Larger Aperture, Multi-Beam, Single Dish (Sub)Millimeter Facility

Improved Sensitivity and Dynamic Range for Molecular Detections
The spectral and imaging dynamic range limitations of currently-leading facilities such as ALMA hinder our ability to detect and map new chemical species in planetary atmospheres.ALMA's nominal spectroscopic dynamic range (the ratio of strongest to weakest signals in the spectrum -typically expressed as the line-to-continuum intensity ratio), is ∼ 1000 at 1.3 mm, but falls to 150-400 at shorter wavelengths (∼ 0.7 mm) using standard calibration procedures.Such limitations make it difficult to reliably measure weak spectral lines and detect new molecules in the (sub)millimeter range.This is a particular problem for planetary targets that have a relatively bright (sub-)mm continuum.The maximum achievable dynamic range is dictated by the ability to accurately calibrate (or flatten) the spectroscopic baseline, which, for interferometers like ALMA, requires the observation of a point-like calibration source with a strong, line-free continuum.Improved calibration accuracy of a single-dish ('total power') facility such as AtLAST could be achieved through observations of brighter, more extended, spectrally featureless solar system objects (such as the Moon or Mercury).This would enable more precise calibration of the intrinsic receiver bandpass.
Advances in instrument design and optics since the construction of present-generation (sub)millimeter facilities such as the IRAM 30-m and JCMT will further facilitate the acquisition of spectra at high dynamic range, allowing the detection of new, spatially distributed trace gases in planetary atmospheres, such as complex organic molecules and other (perhaps biologically relevant) species.For optimal science return on AtLAST, a spectroscopic dynamic range ∼ 10 5 across large bandwidths should be the goal (see Science Case ii).To realize this objective, close attention should be paid to achieving the flattest possible spectral bandpass, through careful choice of optical design, receiver components, and optimal calibration strategies.To this end, it will also be important to mitigate internal reflections and standing waves between the telescope components as much as possible.
Cometary comae are highly diffuse and extended objects.Their molecular excitation is mostly governed by a combination of thermal and fluorescent processes, leading to rotational emission lines with peak intensities in the (sub)millimeter band.Nevertheless, molecular column densities remain low, so that detection of new molecules and mapping of known species can be challenging in all but the brightest comets.Studies of cometary chemical compositions therefore benefit in proportion to the available spectral line sensitivity; ideally, AtLAST would provide a factor of at least a few improvement compared with existing facilities.This can be achieved through a combination of ∼ 2× smaller beam size for a given frequency (focusing in on the denser, innercoma), combined with lower receiver temperatures, and reduced sky opacities compared with the IRAM 30-m, which is the current state-of-the art instrument for total power mm-wave spectroscopy of comets.Our goal sensitivity for new spectroscopic detections of cometary molecules with AtLAST should be ∼ 1 mK (per km s −1 of bandwidth) in ≲ 8 h observing time (see Science Case iv).

High Spatial Resolution
The relatively low spatial resolution of current singledish (sub)millimeter facilities has limited the potential for mapping of planetary atmospheres and surfaces.The angular diameters of the major planets (and Titan) are shown in Figure 1, demonstrating the spread in their apparent sizes over time due to orbital motions.For comparison, the diffraction-limited angular resolution (primary beam FWHM) for a 50 m diameter antenna at a variety of frequencies is also shown.The maximum number of resolution elements within the area of the major Solar System bodies as a function of these frequencies (and wavelengths) are given in Table 1.These data demonstrate that a 50-m class single-dish (sub)millimeter telescope can readily resolve Venus, Mars, Jupiter and Saturn -often, throughout their orbits, in a number of frequency ranges -while the Ice Giants (Uranus and Neptune) can be marginally resolved at terahertz frequencies.Titan remains smaller than 1 ′′ throughout its year, but is still an important target for total power spectroscopy with At-LAST, since a 50 m single dish would allow high-frequency observations without severe beam dilution (particularly above ∼ 800 GHz).Large planets such as Jupiter and Saturn are a particular challenge for interferometers like ALMA, NOEMA and SMA as they possess large angular structures that are resolved out at long baselines.As shown in Figure 1 and Table 1, a large-aperture, single-dish telescope would enable the observation of planetary objects and surface/atmospheric features that span a wider range of spatial scales.Providing planetary observers with a consistent range of angular resolutions at all times, would facilitate studies of temporal phenomena, contemporaneous observations with other facilities, and support for spacecraft measurements.This latter task can prove difficult for interferometers, which tend to have variable configurations (and thus, variable angular sensitivities) throughout the observing schedule.Further, a single-dish instrument with focal-plane receiver array at least an arcminute in diameter would permit mapping of both large and small-scale planetary features.This is particularly valuable for moderately well mixed gases that may only vary on global scales.A small beam size of ≲ 5 ′′ is required to resolve the majority of planets in our Solar System, which is achievable at frequencies >250 GHz for a 50 m facility.Additionally, a small beam would help minimize beam dilution for the small (∼ 0.1 ′′ -2 ′′ diameter) bodies such as the icy moons of Jupiter, and will be required to ensure Saturn's moon Enceladus can be reliably separated from Saturn's bright continuum emission (see Science Case iii).

Instantaneous Multi-beam Spectral Mapping
AtLAST has the potential to capitalize on recent advances in heterodyne multiplexing and multi-beam receiver technology to produce unprecedented spectro-spatial sub-mm maps of the larger planets, as well as cometary comae, the latter of which often extend over hundreds of arcseconds on the sky.The ability to map in two spatial dimensions the chemical distributions and dynamical motions of Solar System objects at a given instant in time is crucial for properly characterizing short-term phenomena such as planetary impacts and storms, and cometary jets and outbursts.Spatial resolution of a few arcseconds over a field of view up to a few arcminutes, with ∼ 10-30 beams across a single axis (for example, in a square or hexagonal, ∼ 25 × 25 pixel array), at a spectral resolution ∼ 0.1 km s −1 , would allow detailed imaging and Doppler studies of molecular species in planetary and cometary atmospheres to study the combined effects of chemical and dynamical processes in these bodies.

Science Case i: Giant Planet Atmospheres
The complex, extended atmospheres of the Gas and Ice giant planets provide interesting, yet challenging targets for remote observations due to (1) small-scale vertical (radial) variations in temperature and composition, (2) the brightness of their continuum emission (particularly in contrast to trace atmospheric species), and (3) their high atmospheric pressures and fast rotation speeds, which can significantly broaden the spectral lines.Millimeter/submillimeter facilities provide powerful probes of the composition, structure, and dynamics of the tropospheres and stratospheres of giant planets through the combination of rotational emission/absorption line measurements, studies of highly pressure-broadened (i.e.>1 GHz wide) pseudo-continuum absorption features, and spatial mapping of brightness temperature variations.The continuum emission from giant planets is dominated by the temperature structure and deep abundances of H 2 , He, NH 3 , H 2 S, PH 3 , CH 4 , and other gases, which are difficult properties to observe directly through remote sensing and often require in situ measurements through interplanetary probes [143,108,109].With the advancement of larger aperture and interferometric facilities, the spatial variability of the giant planet continuum emission can be mapped to better understand the connection between deeper atmospheric composition and dynamical activity, and will help to contextualize spacecraft observations including those made by Juno and the recently launched JUpiter Icy Moons Explorer (JUICE) mission [68].The development of a larger-aperture, single-dish telescope (At-LAST) will further help to improve our understanding of gas and ice giant planet atmospheres -and by extension, exoplanets with similar sizes and compositions -in several key areas: • The spatial and temporal variation of temperature structure throughout the troposphere and stratosphere, which relates to the global and local circulation, dynamics and composition of the atmosphere.
• Variability of atmospheric winds and vertical wind shear on short-term and seasonal timescales, and how that relates to climatological and other forcing mechanisms.• The composition and distribution of trace species throughout the stratosphere and troposphere.• The changes in chemical abundances over short and long timescales, and their connection to seasonal and transient events such as storms, aurorae and infall from space (cometary impacts, micrometeorites and dust).

Dynamics
The dynamical state of planetary atmospheres -often described by their temperature structure and wind fieldsallows for the characterization of their atmospheric circulation, energy budget, and seasonal variability.Previous studies at long wavelengths have measured brightness temperature and compositional variations in the giant  2 shows the comparison of ALMA observations of Jupiter at 1.3 and 3.0 mm to a visible wavelength composite from the Hubble Space Telescope (HST), allowing for the comparison of brightness temperatures at multiple pressure levels from various zonal features.At frequencies >700 GHz, a 50 m single-dish facility with large focal plane receiver array will be able to provide comparable latitudinal coverage to ALMA Band 3 observations (Figure 2, right panel) in a single integration.Sensitivity ≲ 1 K over short integration times (minutes) would allow for the precise determination of brightness temperature variations compared to the bright continuum emission while avoiding longitudinal smearing of localized features due to planetary rotation and winds; higher sensitivity and dynamic range considerations pertain mostly to spectral line surveys, as discussed below.Wide bandwidth (>a few GHz) settings enable the pressure-broadened rotational transitions (e.g., CO at ∼ 115, 231, 346, 461, 576, 691, 807, and 922 GHz) and pseudo-continuum features (e.g., H 2 S at ∼ 169, 408, and 736 GHz; PH 3 at ∼ 267 GHz; and NH 3 at ∼ 572 GHz) to be modeled so the deep temperature and composition can be characterized, which complements observations in the IR down to ∼ 1 bar (e.g., 67).Observing the giant planets multiple times throughout their extensive orbital periods will improve our understanding of the influence of seasonal variability on the circulation and dynamical state of the atmosphere and allow for comparison to seasonal photochemical models [107,106,78,79].
In contrast to the meridional circulation of the atmosphere, which is too slow to readily detect using spectral line shifts, the direct measurement of horizontal zonal atmospheric winds is possible from the Doppler shifts of molecular rotational transitions at high (∼ 0.1 km s −1 ) spectral resolution.Wind speeds of the giant planets have typically been inferred directly from cloud tracking techniques (e.g. through Voyager imaging; 90), which rely on correlation analyses of small-scale features observed over multiple rotations.Zonal wind speeds can also be inferred through infrared measurements of the thermal wind shear [66].Direct observation of spectral line (e.g., HCN at ∼ 89, 265, 355, and 709 GHz) Doppler shifts are required in order to obtain precise wind speed measurements at altitudes above the clouds, avoiding the problems associated with changes in the shapes of resolved atmospheric features over multiple rotations.Observations of strong stratospheric rotational transitions of CO and HCN using ALMA has resulted in wind measurements for all giant planets apart from Uranus [31,9,20].
Figure 3 shows examples of winds derived from ALMA observations of Jupiter, Saturn, and Neptune compared to cloud tracking and auroral measurements.Measuring planetary wind speeds as a function of latitude and at multiple altitudes requires high spectral and spatial resolution combined with a large field of view, which is enabled by heterodyne and (sub)millimeter instrumentation.While the variability in spatial resolution of the giant planets (Figure 1, Table 1) precludes complete latitudinal coverage for all planets at all frequencies, even hemispheric comparisons of Doppler wind measurements over time provide interesting results -see, for example, the abrupt change in Titan's stratospheric jet [89,40].As a single-dish facility would not be confined to time-variable configuration (resolution) constraints, episodic changes in wind speeds and wave propagation can be investigated frequently, and over irregular timescales as needed.
Jupiter and Saturn exhibit large-scale weather events where dramatic changes in the planet's visible appearance are accompanied by energetic changes in cloud color and morphology.On Jupiter, such events periodically take place at low latitudes [64], in the typically white band know as the Equatorial Zone (±7 • latitude), and in the adjacent, darker North and South Equatorial Belts between ±7 − 18 • latitude.Each region undergoes changes in its normal cloud color and morphology, exhibiting waves, eddies, and wakes, initially beginning over a small lo-(e.g.HCN; [1]).However, the exceptional spatial and spectral resolution of ALM e utilized in a search for additional, undetected molecules produced in Jupiter's atm ted through cometary impacts.In order to continue tracing the species generated act and better deduce to what extent Jupiter's atmospheric composition can be us f impact history, we must investigate the full extent of minor atmospheric products pr high energy interactions or dredged up from Jupiter's deeper atmosphere.
temperature of the 1-3 mm ALMA data with best-fit models to the 2013-2014 VLA data (from dP19).We stress here that no new models were produced; the existing models were merely extended into the millimeter-wavelength range.Hence, as in dP19, we ignored opacity by clouds.The latter authors justified this assumption based upon disk-averaged spectra at millimeter to centimeter wavelengths.They argued that if cloud aps.The scans were created by median averaging over 60°of longitude, centered on d.Because a limb-darkened disk had been subtracted from the data, the background re offset for clarity by 10 K each, while each set is separated by 20 K.The spatial A maps, which lowers the feature contrast.The vertical dashed lines (at, e.g., the EZ, e top is the (eastward) wind profile as measured from the HST data; the scale is given ure 4. (B) Longitude-smeared ALMA maps of Jupiter's thermal emission at 1.3 and map from dP19, after subtraction of a limb-darkened disk.
ver the solar values, ctor of 3.2, and the an adiabat (typically e 165 K at the 1 bar ation profile (Lindal profile follows that CIRS) observations observed brightness variations in opacity perature.They show kely than changes in 9, we attribute all latter authors also temperature due to temperature of the 1-3 mm ALMA data with the 2013-2014 VLA data (from dP19).We new models were produced; the existing m extended into the millimeter-wavelength in dP19, we ignored opacity by clouds.justified this assumption based upon disk-a millimeter to centimeter wavelengths.They a opacity were important, the brightness tem meter wavelengths should be affected much centimeter range, because the mass absorp inversely proportional to wavelength for part compared to the wavelength (Gibson et al. Figure 7 shows the zonal-mean brigh spectra of the ALMA data together with 2013-2014 VLA data, superposed on the m meared ALMA and VLA maps.The scans were created by median averaging over 60°o on a longitude/latitude grid.Because a limb-darkened disk had been subtracted from th 24 GHz scan.The scans are offset for clarity by 10 K each, while each set is separate than at 1.3 mm and the VLA maps, which lowers the feature contrast.The vertical dashe atures.The green line at the top is the (eastward) wind profile as measured from the HST h the HST image from Figure 4. (B) Longitude-smeared ALMA maps of Jupiter's therm ectively), and a VLA 3 cm map from dP19, after subtraction of a limb-darkened disk.).Finally, despite the northern aurora being located on the western side, mostly behind the terminator, we see a broad signal on the eastern limb at polar latitudes with an average LOS velocity of about +100 m s 1 , for which we lack a clear explanation.A more favorable observation geometry of the northern polar region is thus required to improve our understanding of the stratospheric circulation in this region.

Discussion
The branch of the northern auroral jet we tentatively detect lies below the electrojet discovered at p < 1 µbar from infrared observations of H +  where u is the zonal component of the wind speed at the limb.at the equator in the latitude range from 3 S to 3 N to the average winds of the large equatorial jet that sp 20 S and 25 N, this narrow jet is relatively weak i vations at 350-500 mbar pressure, but it is very in observations at 60-250 mbar pressure (Fig. 9).Our tions between 0.1 mbar and 20 mbar do not reveal a narrow intense jet at the equator.This is quite surp expected to detect this peak which could be related t That we do not see any evidence for such a pe posed over the broad eastward jet in the CO data result from the large vertical extent of the CO wi tion function (Fig. 5), which may cancel any contrib SSAO by encompassing opposite phases of the SSA ingly, the HCN wind contribution function is more the CO one (around the peak level), and looking car HCN wind profile (Figs. 8 and 9), we detect a narro imum in the velocities between 5 S and 1 N, wit amplitude of 50 ± 20 m s 1 with respect to the aver 10 S and 10 N shown in Fig. 10.This is consistent w of magnitude of the SSAO peaks, according to G (2018).

Northern hemisphere mid-latitudes
Between 25 N and 60 N, we find tentative evide first time of a global westward wind with an avera -50 ± 30 m s 1 in both limbs from the CO data (Fi over, our HCN and CO wind measurements sh tropospheric eastward jet seen at 42 N (Fig. 9) has vanished in the stratosphere.We also find westwa with HCN between 25 N and 50 N, but only mar noteworthy that some of the only features that coul in Saturn's stratosphere were the hot vortices, nic beacons, that were formed in the stratosphere follow Great White Spot of 2010-2011.cal region and expanding to encompass the whole region for some time before returning back to its nominal state.These events on Jupiter happen on a timescale of between 2-18 years, with little correlation to adjacent regional activity, season, or orbital position [5,6].On Saturn, however, the northern hemisphere experiences a large convective storm event approximately every 30 years.The previous 'Great Saturn Storm' in 1990 occurred prior to the advent of the ALMA telescope, spanning ∼ 15 • latitude and eventually encircling the whole planet.This event led to significant stratospheric effects that could have been readily characterized in the (sub)millimeter [65].Additionally, 'Great White Spots' have previously been observed with Herschel [26].The sizes of these types of spectacular atmospheric events (∼ 2 − 5") would make them observable by AtLAST.As extensive meteorological events offer a window into the planets' deeper atmospheres, contemporaneous (sub)millimeter observations (which can probe different altitudes by highly resolving line profiles and measuring multiple rotational transitions), provide complementary information to other facilities (e.g., HST, JWST, the Juno Microwave Radiometer instrument).Since a large, single-dish facility would provide adequate spatial resolution to discern between the ∼ 1 − 5" zonal regions on the giant planets (see Table 1, Figure 2) over their orbital period, strategic and regular observations by AtLAST would serve as an excellent resource to investigate these spectacular episodic events and provide global context for other ground-based and spacecraft measurements.

Chemistry
The detection of "non-equilibrium" C, N, O and S-bearing molecular species -such as HCN, CO, and CS -in the upper atmospheres of the giant planets has prompted questions pertaining to their origin and chemical evolution.Chemistry and circulation may result in the production and distribution of trace species in the upper atmosphere, as may the infall of cometary, interplanetary dust, or satellite material, or some combination thereof [24,25,27,28,29,30,87,81,127].While previous observations at millimeter wavelengths have allowed for the investigation of these processes (Figure 4), higher sensitivity measurements are required to more fully characterize the governing sources and sinks of nonequilibrium species in the atmospheres of the giant planets.In addition, long-term observing campaigns allow for the measurement of temporal variability in these species, as well as potential sporadic events (e.g.storm outbreaks; cometary impacts) that may alter the upper atmospheric composition and can be measured by AtLAST.To perform these investigations requires wideband (≳ 16 GHz), high-resolution (≲ 1 MHz) spectral imaging, as the spectral transitions from minor atmospheric constituents can manifest as significantly extended (several GHz wide), pressure broadened pseudo-continuum absorption wings around a narrow (∼ MHz-wide) line core (see e.g., CO on Neptune; [93,94]).Finally, the detection of new, trace chemical species in the atmospheres of the giant planets will inform our understanding of their atmospheric photochemistry, and elucidate connections between the external sources from the surrounding system (e.g.dust, satellite infall) or beyond (e.g.comets, galactic cosmic rays).While ALMA has greatly facilitated the detection of heretofore unknown molecules in the atmosphere of Titan, the background continuum emission from the giant planets presents a challenge to the detection of all but the strongest emission lines (CO, HCN, CS, etc.) in their atmospheres.To further expand the known chemical inventory of molecular species, both predicted and serendipitous, the nominal dynamic range of AtLAST needs to be sufficiently high so that emission lines over 1000 times weaker than the continuum level can be reliably detected against large (∼ 20 − 30"), bright (several hundred Kelvin) continuum sources.Chemical species of interest include HC 3 N (∼ 136, 264, 318, and 672 GHz), CH 3 CN (bands at ∼ 129, 221, 312, and 698 GHz), CH 3 CCH (bands at ∼ 154, 222, 325, 632 GHz), C 3 H 8 , and other nitriles and hydrocarbons.AtLAST therefore holds the promise of detecting new chemical species on giant planets through advancements in sensitivity and dynamic range, while a large focal plane array (≳ 25 × 25 pixels in size) is required to map their spatial distributions without rotational smearing.Such observations would provide major insights into the origins of non-equilibrium species, through spatial associations with storms, aurorae or infall phenomena, and would lead to a new understanding of the extent of photochemical processing of gases in giant planetary atmospheres.In addition, studies of molecular isotopic ratios provide profound insights into the formation and evolution of the giant planets ( [116] and references therein).The measurement of C-, N-, and O-bearing isotopes on the giant planets enables the determination of the main reservoirs for these species, providing crucial inputs for planet formation models.With a spectroscopic dynamic range ≳ 10 4 , measurement of more accurate isotopic ratios for 13 C, 17 O, 18 O, and 15 N in CO, HCN and other gases will be possible, resulting in improved understanding of how planetary systems form and evolve from the isotopically diverse ice, gas and dust reservoirs in the protosolar accretion disk.

Science Case ii: Terrestrial-Type Atmospheres
Terrestrial-type planets with thick atmospheres are well suited to molecular spectroscopic studies at (sub)millimeter wavelengths.
With atmospheric molecular column densities in the range 2.4 × 10 23 -1.5 × 10 27 cm −2 , and angular diameters ∼ 1 − 60 ′′ , Venus, Mars and Titan present compelling targets for spectroscopic studies using a sensitive (sub)millimeter facility.As shown in Figure 1 and Table 1, Mars and Venus can be readily resolved using a 50 m (sub)millimeter telescope, allowing spectral mapping of trace atmospheric gases to improve our understanding of their atmospheric physical and chemical processes.For each source, the models that best fit the emission core are displayed: an internal source yielding a mole fraction of 1.9 ˆ10 ´8 in the upper troposphere (red line), a steady external flux (due to IDP or a local source) of 2.2 ˆ10 5 cm ´2 s ´1 (blue line), and a comet with a diameter of 640 m depositing 3.4 ˆ10 13 g of CO above the 0.1 mbar level "300 years ago (green line).These models were computed with the thermal profile of Feuchtgruber et al. (2013).The internal source model overestimates the line core width and produces a broad absorption that is not observed in the data.The external source models can barely be differentiated.Right: vertical profiles associated with the spectra.
time down to 1 mbar in Uranus in our model, but other combinations of deposition time and level are possible.To infer the mass and diameter of the comet, we assumed the comet density was 0.5 g cm ´3 (Weissman et al. 2004;Davidsson et al. 2007) and that the comet yielded 50% CO at impact (Lellouch et al. 1997).
The vertical profiles and resulting spectra corresponding to the three sources, as obtained with the thermal profile of Feuchtgruber et al. (2013), are displayed in Fig. 2. The best fits to the spectrum are obtained for external source models.Despite resulting in different vertical profiles, a steady flux φ CO " 2.2 ˆ10 5 cm ´2 s ´1 and a 640 m diameter comet depositing 3.5 ˆ10 13 g of CO (y 0 " 3.1 ˆ10 ´7) result in lines that are indistinguishable from the standpoint of our observations.Such impact at Uranus occurs every "500 years with a factor of 6 uncertainty (Zahnle et al. 2003).Such timescales are fully compatible with our assumption on ∆t.With the thermal profile of Orton et al. (2013a), we obtain slightly higher values because of lower stratospheric temperatures: φ CO " 2.7 ˆ10 5 cm ´2 s ´1 and y 0 " 3.9 ˆ10 ´7 (i.e., a 700 m diameter comet).All fit parameters are listed in Table 2.These values remain to be confirmed by more rigorous (photochemical) modeling and higher S/N data.
The amplitude of the CO emission peak is reproduced with an internal source model in which y CO " 1.9 ˆ10 ´8 (see Fig. 2).With the thermal profile of Orton et al. (2013a), y CO " 2.7 ˆ10 ´8.We note that "three times more tropospheric CO is needed in this model, compared to the "uniform" empirical model value derived in Sect.3.2.This is due to the fact that the observed emission line probes the mbar level, i.e., where the CO vertical profile sharply decreases because of the low homopause in the atmosphere of Uranus.As a result, a stronger internal source is required to reach a sufficient level of abundance of CO around the mbar level.The main outcome of this model is that it now overestimates the line core width and results in additional broad absorption because CO is much more abundant in the lower stratosphere than in the external source models (by as much as two orders of magnitude at 10 mbar).
The absence of such a broad CO absorption in the data cannot be caused by our sinusoidal ripple removal procedure because we have removed sine waves of much shorter period than the total width of such broad absorption wings.We can rule out the internal source model because the width of the line core is not fitted, there is no broad absorption in the spectrum, and the derived y CO values are an order of magnitude larger than the upper limit set by Herschel-SPIRE observations of this region of the atmosphere (Teanby & Irwin 2013).Thus, as long as there is no significant photochemical source of CO in the stratosphere, the HIFI line is caused by external CO.The upper tropospheric mole fraction of CO is fixed at the level where the thermochemical equilibrium is quenched by vertical diffusion.
The upper limit of Teanby & Irwin (2013) on the internal source (y CO " 2.1 ˆ10 ´9) can be further used to try and constrain the deep atmospheric O/H ratio in Uranus.Their observations probe between 10 and 2000 mbar, i.e., well below the homopause level (see Fig. 2 right).As a consequence, this upper limit is valid even if the authors have not accounted for the low homopause of Uranus.
We have adapted the thermochemical model developed by Venot et al.

Venus
Mapping of spectral lines from Venus with strong transitions (such as SO 2 , SO and CO) can provide insights into the chemistry and dynamics of our sister planet's atmosphere, much of which remains to be fully understood.The high spectral resolving power (∼ 10 6 -10 7 ) offered by heterodyne spectroscopy results in measurements mostly sensitive to altitudes ≳ 60 km -above the thick Venusian cloud layers -allowing detailed studies of the mesosphere and above.Previous space missions to Venus including Pioneer Venus, Venera, Venus Express, and Akatsuki have allowed us to better understand the sulphur and water cycles and global circulation of the atmosphere [97,125].As described by [62], ground based studies provide a valuable complement to these datasets, thanks to their ability to monitor the behavior of gases as function of time, both in the short term (∼hours) and the long term (∼years or decades).With respect to space-based data, ground-based observations have the advantage of instantaneous snap-shot imaging over the whole disk of Venus, enabling global atmospheric studies as a function of latitude, longitude, and local time.The Venus Express and Akatsuki orbiters detected strong, long-term (multi-year) variations in the planet-scale zonal wind speeds and atmospheric SO 2 and CO abundances that remain to be fully understood.In the wake of these space missions, continued ground-based (sub)millimeter monitoring is therefore important to better understand these phenomena, which will help us understand the evolution of Venus's global climate.Spectral-spatial mapping of the entire (∼ 10 − 60 ′′ diameter) Venusian disk at a spectral resolution ∼ 0.1 km s −1 and spatial resolution of ∼ a few arcseconds is required to fulfill this science objective.The complex behaviour of sulphur-bearing molecules is a particularly important topic for further study due to their central role in Venus's climate, their close relation-ship with other volcanic gases (such as the ubiquitous clouds of sulfuric acid that blanket Venus), and their ability to trace photochemical and dynamical phenomena [97].The strong temporal and spatial variability of SO and SO 2 abundances at high altitudes remains largely mysterious, so further mapping studies are required to elucidate the sources and and sinks of these molecules in the upper atmosphere, and relate them to the lower altitude (volcanic) gases and aerosols.Such temporal studies are difficult using ALMA due to its variable array configuration that resolves out most of Venus's disk for a large part of the year, so rapid, multi-beam total-power mapping is required.
CO plays a key role in the carbon dioxide cycle, which maintains the chemical stability of the primary constituent of Venus' atmosphere.Previous (sub)millimeter studies [34,36,35,61] disagree on the extent of spatial variations in the CO distribution, so additional, temporally resolved mapping studies of this molecule are needed.In particular, future, temporally and spatially resolved, ground-based measurements of the atmospheric composition of Venus will be required to contextualize in situ data from upcoming spacecraft missions such as DAVINCI+ [69], VERITAS [21] and EnVision [70].
The unparalleled sensitivity of AtLAST would enable searches for new molecules (including sulfur, phosphorous and chlorine-bearing species, as well as organic molecules) that would help improve our understanding of the chemical processes occurring in the Venusian atmosphere, and their link to the volcanic gases that are thought to be so critical in governing Venus's hot global climate.Due to the bright (several hundred Kelvin) Venus continuum, and the comparative weakness of the mm-wave absorption lines from any yet-to-be-detected molecules (∼ 10 mK), searches for new molecules will place strong demands on the spectroscopic dynamic range (DR) of the instrument (which would ideally reach DR ∼ 10 5 to allow in-depth studies of Venus's trace gas inventory).Close attention will therefore need to be paid to the bandpass calibration accuracy and spectroscopic baseline stability, in order to achieve very flat baselines, to allow us to realise the full potential of such a powerful telescope for planetary science.A precise measurement or upper limit of the phosphine (PH 3 ) abundance would be just one of the ways in which AtLAST would help revolutionize Venus atmospheric science [72,137,7,51].

Mars
The divergent histories of Earth and Mars contain valuable lessons in our quest to determine the conditions required for the origin and persistence of life on planetary surfaces.Understanding the present-day climate and atmospheric composition of Mars, and how it evolved over the history of the Solar System, is therefore an important objective in the fields of planetary science and astrobiology.Mapping and detection of trace atmospheric gases such as CO and H 2 O allows us to constrain photochemical networks and climate models, while measurements of molecular isotopic ratios provide unique insights into the planet's atmospheric loss history [136,116].Quantifying the Martian atmospheric circulation/global wind field is another essential objective towards revealing the past climate of Mars.Detailed studies of the composition and wind field of Mars's upper atmosphere were recently performed using data from the Mars Atmosphere and Volatile Evolution (MAVEN) spacecraft/orbiter [10].These measurements were obtained in the thermospheric altitude range (∼ 140-240 km), using data acquired over several Martian seasons and spanning a large number of spacecraft orbits.Discrepancies between the observed wind speeds and those predicted by state-of-the-art climate models indicate that the global air circulation on Mars remains to be fully understood [121]; further studies of these phenomena, including how the thermospheric wind speeds are influenced by those at lower altitudes, are therefore needed.Spatially complete measurements of Mars's 4"-26" disk at high (∼ 0.1 km s −1 ) spectral resolution are required in order to derive the zonal wind speeds as a function of latitude, similar to those described for the giant planets (Science Case i).The circulation in the middle atmosphere (40-80 km) is also affected by atmospheric transport and climatic processes in the lower atmosphere (<40 km).Although the lower atmospheric dynamics of Mars have been characterized by previous spacecraft missions (e.g.[77,71]), and are reasonably well reproduced by general circulation models (GCMs), the middle atmosphere has been less well studied, so models are less constrained in this region [103].(Sub)millimeter heterodyne observations are uniquely well suited to probing the Martian middle atmosphere, yielding temperature and wind measurements from spectrally and spatially resolved molecular line profiles [85,22].As shown by [103], comparison of groundbased Mars CO mapping observations with GCM predictions (see Fig. 5) indicates significant discrepancies in the wind speeds and their spatial structure, and the GCMs often underestimate the temperatures in the middle (20-50 km) atmosphere.No Mars space mission to-date has included highresolution heterodyne instrumentation.Mars science would thus benefit from an improved-sensitivity groundbased (sub)millimeter telescope.With an angular resolution of a few arcseconds in the (sub)millimeter range (to cover e.g. the CO J = 3-2, 6-5, 7-6 or 8-7 lines between 300-1000 GHz), combined with a large field of view, At-LAST would perform global monitoring of the middle atmosphere of Mars (over timescales from hours to years), enabling critical benchmarking and subsequent improvement of climate models.

Titan
Since the end of the Cassini mission, ALMA has risen the forefront of ground-based Titan studies, thanks to its ability to perform spatial/spectral imaging of a large number of atmospheric gases, measure their (3D) distributions, and detect new species, including molecules of possible astrobiological relevance [45,47,42,118,115,132,133].This has demonstrated the power of long wavelength astronomy to provide new insights into the chemistry, dynamics, climate, and potential habitability of small icy moons in our Solar System.Although a 50-m class facility such as AtLAST will not spatially resolve Titan (which has an angular diameter ≈ 1 ′′ , including its extended atmosphere; Table 1), with sufficient collecting area, bandwidth and atmospheric transmission, AtLAST would be able to compete with large (sub)millimeter interferometers such as ALMA and NOEMA in terms of flux sensitivity per beam.Consequently, broadband molecular line surveys with AtLAST would be able to reach unprecedented sensitivity towards objects such as Titan.Broad bandwidths (10's of GHz) have been shown to be crucial for the detection of increasingly complex organic molecules in extraterrestrial environments, due to their large number of individual rotational transitions spread across the cm/(sub)millimeter range [13,100,92,39].Combining the fluxes from multiple emission lines in a broadband survey (e.g. between 200-300 GHz at sub-mK RMS and 1 km s −1 spectral resolution), would therefore provide sufficient sensitivity to detect new molecules on Titan, including molecules of possible (pre-)biotic relevance, and help improve the accuracy of abundance ratio derivations.Similarly, AtLAST could enable the detection and characterisation of new isotopologues of known molecules (such as 13 C, 15 N and deuterium-substituted forms of C, N and H-bearing organics), which have the power to elucidate the long-term physico-chemical evolution of Titan's mysterious atmosphere [96,95,116].

Science Case iii: Plumes from Icy Moons
Advances in our knowledge regarding the interior structure of icy moons have demonstrated the likely presence of subsurface oceans throughout the Solar System [114,123].Enceladus, Europa, Ganymede, Callisto and Triton are known or suspected ocean worlds.Cryovolcanic activity can launch plumes high into their atmospheres, and their molecular emission could be detected by a sensitive (sub)millimeter instrument.Plumes linking the subsurface with the atmosphere provide a potential means for probing the compositions of subsurface liquid reservoirs, and open up the possibility of remote habitability studies.Sensitive (sub)millimeter single-dish spectroscopy can probe the chemical compositions of the tenuous atmospheres of the icy Jovian satellites Europa, Ganymede and Callisto, the Saturnian satellite Enceladus, and the Neptunian satellite Triton to help determine the surface and/or interior compositions, from which insights into the chemistry and habitability of the subsurface oceans of these moons can be obtained.Subsurface oceans are primary targets in the search for extraterrestrial life [76], but in-situ survey missions are extremely costly and difficult.One of the most outstanding results of the Cassini mission to the Saturnian system was the discovery of plumes emanating from the south polar region of Enceladus, believed to originate from a warm, subsurface ocean [119,124].However, Cassini mass spectrometry of the plume was relatively low resolution leading to ambiguity in the detection of molecules with equivalent integer atomic mass (e.g.N 2 /C 2 H 4 and HCN/C 2 H 3 ; [139]).The possibility of chemical reactions and molecular dissociation within the Cassini mass spectrometer chamber is also difficult to quantify or rule out, as is the potential contamination of mass spectra due to Cassini fly-throughs of the Titan atmosphere.Unambiguous spectroscopic followup of the Enceladus plume composition is therefore essential.Following the end of the Cassini spacecraft mission, ground-based spectroscopic studies present a feasible alternative for studying the composition of the Enceladus plume(s), and complement observations from JWST [138].However, present observatories suffer from a lack of sufficient sensitivity and/or angular resolution to readily detect the faint, ∼ 1 ′′ long plumes.Contamination from the bright, nearby (< 25 ′′ away) Saturn is also a major issue.Even large interferometers such as ALMA struggle due to the relatively large primary beam size provided by a 12 m sized antenna, which makes it difficult to exclude emission from Saturn and its rings from the field of view, leading to troublesome image artifacts.To avoid Saturn and exclude it from the 2nd and 3rd sidelobes of the primary beam, AtLAST would require a HPBW < 8 ′′ at 1.1 mm, and as such requires an aperture diameter >40 m.
As an example of what may be detectable, assuming a plume H 2 O production rate of ∼ 10 28 s −1 , constant outflow velocity of 0.5 km s −1 , and rotational temperature of 25 K [140,138], with a jet opening angle of 30 • and HCN abundance of 10 −4 (10× less than observed in typical comets [15]), the peak antenna temperature for the HCN J = 4 − 3 transition is calculated to be ∼ 25 mK (after convolution with a Gaussian beam of FWHM 3.5 ′′ ).Combining a small beam size with milli-Kelvin sensitivity, AtLAST would therefore be able to open up the possibility of direct ground-based studies of the Enceladus ocean chemistry.Specific target lines for this science case include HCN (354 GHz), H 2 CO (351 GHz), H 2 O (183 GHz), CS (243 GHz), SO 2 (252 GHz), H 2 S (369 GHz) and transitions of various salts and complex organic molecules including alcohols, carbon chains, aromatics and amino acids in the range 125-373 GHz.Detections (or useful upper limits) on the abundances of these species in the Enceladus and Europa plumes could thus be obtained.Searches for new molecules in the extended Enceladus torus [75] would also be possible.This would not only provide a crucial followup and validation of the Cassini mass spectrometry studies, but also allow changes in the plume composition as a function of time to be observed.

Science Case iv: Comets
Comets accreted in the Solar System at around the same time as the planets, and are believed to contain pristine (largely unprocessed) material from the protosolar accretion disk and prior interstellar cloud [113].Studies of cometary ices therefore provide unique information on the physical and chemical conditions prevalent during the earliest history of the Solar System.Due to their chemical compositions, rich in water and organic molecules, cometary impacts could also have been important for delivering the ingredients of life to otherwise barren planetary surfaces throughout the Solar System, and may therefore have played a role in initiating prebiotic chemistry [33].Examining the molecular and isotopic content of comets provide unique insights into the relationship between interstellar and planetary material.Through broadband molecular surveys, we gain insights into the diversity of cometary compositions, improving our knowledge of the chemistry that occurred during planetary system formation.Based on prior experience with the CSO 10-m and IRAM 30-m telescopes [17,12,14], detections of new molecular species (previously undetected in comets), including complex organic molecules relevant to the origin of life, can be expected using a next-generation, single-dish (sub)millimeter facility.Observations of deuterated and 15 N/ 13 C-substituted isotopologues of known molecules in a statistical sample of comets will provide fundamental new insights into the chemical origins of cometary (and thus, planetary) material.Spatial-spectral coma studies (see Figure 6) will provide fundamental information on the physics of cometary outgassing, as well as the complex thermal and photochemical processes occurring in the coma.Using currentlyavailable telescopes, these crucial measurements can only be obtained in the very brightest comets (of which only a few appear every 10 years), so a more powerful (sub)millimeter facility (optimized for fast mapping) will significantly expand our knowledge and generate reliable statistics for the (chemically diverse) comet population.Specific target lines include HCN (354 GHz), DCN (362 GHz), HC 15 N (344 GHz), H 2 CO (351 GHz), HDCO (335 GHz), D 2 CO (342 GHz), H 2 O (183 GHz), HDO (465 GHz, 894 GHz), D 2 O (607 GHz, 898 GHz) and transitions of various complex organic molecules including alcohols, carbon chains, aromatics and amino acids in the range 125-373 GHz (see Figure 7 for a selection of the available lines from known cometary species in the millimeter waveband).Isotopic studies provide a powerful means for understanding the evolution, origins and inter-relationships between different Solar System bodies.In particular, materials throughout the Solar System exhibit differing degrees of deuterium (heavy hydrogen) enrichment, the measurement of which can reveal their chemical and thermal histories.Comets are thought to contain a significant amount of pristine, deuterium-enriched interstellar material, predating the birth of the Solar System [4].However, D-enrichment in water ice may also be possible due to fractionation reactions in the protosolar accretion disk, before the gas dissipates [116].The fraction of pristine (interstellar) vs. reprocessed (disk) ices incorporated into planetary materials is an important open question that can be addressed via improved measurements of cometary D/H ratios.In particular, comparison of cometary D/H ratios with protoplanetary disk models can help unravel the details of where, and from which materials cometary ices were accreted [141,1,38].As shown in Figure 8, the strongest cometary HDO line accessible using ground-based sub-mm telescopes the 1 1,1 − 0 0,0 line at 894 GHz.The next strongest lines at 509 and 600 GHz are on the flanks of the high-opacity 557 GHz H 2 O line, making them undesirable targets, so the 465 GHz line presents the next best alternative after the 894 GHz line.A simulated HDO 894 GHz flux map (and spectral line profile) for a typical comet (assuming HPBW = 1.4 ′′ ) is shown in the right panel of Figure 8. Access to this line would provide an unprecedented opportunity to survey HDO in a large number of comets (several per year), spanning different levels (and modes) of outgassing activity, which would provide new insights into the origin of water on Earth, and elsewhere in the Solar System [75,38,3,116].Due to their spatially extended nature (up to several hundred arcseconds), and ∼ 1/r brightness profile, optimal sensitivity to weak cometary lines is achieved using a large, single-aperture antenna rather than an interferometer, since the latter is blind to coma structures larger than the angular scale probed by the shortest baselines.To expand the statistics for cometary abundance measurements in a meaningful way (observing 2 moderately bright comets per year over a 10-year timespan), with sensitivity to all the species in Figure 6, will require milli-Kelvin sensitivity across ∼ 1 km s −1 spectral line width.Since comets are dynamic, time-variable objects (on timescales of days or even hours), simultaneous observation of the species of interest is desirable for reliable molecular abundance comparisons.An instantaneous bandwidth of several 10's of GHz is therefore required, at a spectral resolution ∼ 0.1 km s −1 in order to resolve the complex kinematical structure of the coma (e.g.[122]).Simultaneous, multi-chroic observations in several frequency bands are required to probe the complete range of species of interest.The ability to rapidly respond to transient cometary phenomena such as outbursts or surprise apparitions will necessitate flexible telescope scheduling, for example, with a turnaround time of a few days between identification of the target phenomenon and acquisition of the observations.It will also be important to accommodate time-constrained observing requests (ideally, to an accuracy of up to a few minutes), to perform synergistic observations with ground or spacebased observatories at other wavelengths, or to capture specific, temporally isolated events such as rapid perihelion/perigee passages or spacecraft flybys.In terms of frequency coverage, sensitive access to the (sub)millimeter range is required, ideally up to and including the HDO 1 1,1 − 0 0,0 line at 894 GHz line.

Conclusion
We have highlighted some the most important planetary science use cases for a new, large-aperture, singledish submillimeter telescope.In particular, the improvements in spectral line sensitivity and dynamic range offered by such a cutting-edge facility will permit the detection of new molecules and isotopologues that will directly enhance our understanding of chemical complexity and abundances in the upper atmospheres of the terrestrial and giant planets, as well as in comets.In turn, these observations will improve our understanding of exoplanetary atmospheric dynamics and compositions, and allow for the improved assessment of planets and moons within our solar system to harbor past or present habitable environments.Instantaneous (snapshot), high resolution, broadband spectral-spatial molecular mapping will provide new insights into transient phenomena such as planetary winds and storms, as well as jets, outbursts and ambient molecular distributions in cometary comae and icy moon exospheres.The relevant instrument specifications required to achieve this science are summarized below.

Summary of Instrument Requirements
Sensitivity RMS sensitivities of ∼ 1 mK per 1 km s −1 spectral resolution element in a few hours on-source will be required for spectroscopic characterization of Solar System bodies (which are often time-variable).This will allow spectral mapping of known trace species, isotopic measurements, and new molecular detections in planetary atmospheres and comets.

Spectroscopic Dynamic Range
Spectroscopic dynamic range (ratio of continuum to line intensity) ≳ 10 5 is required for detection of weak atmospheric spectral lines from trace gases, against the bright continuum of Venus, Mars, Jupiter, Saturn, Uranus, Neptune and Titan.

Spectral Coverage and Resolution
(Sub)millimeter spectroscopic observations of planetary and cometary atmospheres require a combination of broad bandwidth (several GHz or more to properly characterize the profiles of strongly pressure-broadened rotational lines) and high spectral resolution (∼ 100 kHz to fully resolve thermally-broadened line profiles ∼ 1 km s −1 wide).Broad (10's of GHz) of bandwidth in the (sub)millimeter range -preferably, simultaneously across multiple receiver bands -is also necessary to observe multiple molecules and rotational transitions simultaneously, and would enable spectral line stacking analyses for more complex molecules whose lines are often distributed over a large frequency range.This will lead to orders-of-magnitude improvements in sensitivity for the more complex molecules, which can have hundreds of lines within the spectral range of interest.Instantaneous bandwidth of 32-64 GHz would be desirable to avoid having to change receiver tunings during time-critical observations, facilitating relative abundance and isotopic ratio analyses.

Spatial Resolution
A small beam size of ∼ 5 ′′ at 300 GHz is required for the majority of planets in our Solar System to be spatially resolved (see Table 1), paving the way toward more detailed atmospheric studies than have previously been possible with single-dish facilities.Such a small beam is also required to minimize beam dilution for the small (∼ 0.1 ′′ -2 ′′ diameter) bodies such as icy moons and their plumes/exospheres.A small primary beam (≲ 8 ′′ ) is also strictly required for Enceladus observations, to avoid contamination due to its close proximity (< 25 ′′ ) from the bright emission sources of Saturn and its rings.

Mapping Requirements
Spatial resolution of a few arcseconds over a field of view up to a few arcminutes, with ∼ 20-50 beams across a single axis (for example, in a square or hexagonal array, ∼ 25 pixels in across) is required for snapshot mapping of terrestrial and gas-giant planetary atmospheres.A 50-m single-dish telescope with large focal plane array could achieve this without resolving out the larger planets, which is a problem for interferometers.

Time-domain Considerations
Observations of Solar System objects are often timecritical, in the case of observing transient atmospheric phenomena, cometary outbursts etc., and to account for moons being in and out of eclipse/occultation, or to avoid observing when a satellite is too close to its parent planet.Target-of-opportunity observations, coordination with other telescopes/space missions, and accommodation for regular monitoring of science targets with a predefined cadence (over a timescale of days, weeks, months or years), will all be crucial for maximum science return on time-variable Solar System objects.Scheduling the most critical observations to within an accuracy of a few minutes will therefore be desirable.Responding to unexpected transient phenomena such as cometary outbursts or surprise apparitions will necessitate flexible telescope scheduling, for example, with a turnaround time of a few days.

Figure 1 .
Figure 1.Angular diameters of major Solar System bodies and Saturn's largest moon, Titan.The vertical extents of the coloured bars for each body represent their range of angular sizes due to differing geocentric distances throughout the year.The second y-axis shows the angular resolution a 50 m single dish facility would achieve for a range of representative frequencies.

HST:
Jupiter as observed by the Hubble Space Telescope in April, 2017 (left), and ALM 17 at 1.3 mm (center) and 3 mm (right).Image credits: (left) NASA, ESA, A. right) [2].dition, these observations will help to provide global context to the recent aft, whose mission has been extended until 2021.The Jovian Infrared Aurora AM) and Microwave Radiometer (MWR) instruments onboard Juno provide high 1

Figure 2 .
Figure 2. HST (left) and longitude-smeared ALMA (center, right) observations of Jupiter from de Pater et al. [59] (used with permission).The ALMA observations show data after the subtraction of a limb-darkened disk model, enabling the high contrast (≲ 10 K) differences in brightness temperature between Jupiter's zonal structure to be easily observed.

Fig. 3 .
Fig. 3. Jupiter's UV aurora and stratospheric HCN winds.This composite image shows the LOS wind velocities (in m s 1 ) derived from the ALMA observations and the statistical emission of the aurorae (Clarke et al. 2009) in the configuration of the ALMA observations.The northern and southern aurora regions are best seen in the dedicated zoomed-in quadrants.The M = 30 footprints of the magnetic field model from Connerney et al. (2018) are good markers of the positions of the main ovals as seen by Juno-UVS (Gladstone et al. 2017) and are plotted in orange.The white ellipses indicate the spatial resolution of the ALMA observations.The directions of the strongest winds in the equatorial and auroral regions are indicated with the red and ⌦ symbols.
was expected to be close to tangential to the limb on its poleward edge (see Appendix F and Fig. F.1).It is thus no surprise that we find no clear evidence of the jet on the northern edge of the oval.Within the framework of our simplified model, assuming a 300 m s 1 counterrotation wind inside the northern oval nonetheless provides a good fit to the measured wind speeds poleward of 55 N on the western limb where the northern oval was rising (see Fig. E.1 3 emission by Rego et al. (1999) and further constrained by Stallard et al. (2001) and Johnson et al. (2017).This 1 ern oval.Models by Majeed et al. (2016) and Yates et al. (2020) predict that neutrals have higher velocities below the southern oval than below the northern one.Although we find relatively similar velocities underneath the two ovals, our detection in the northern oval remains tentative such that we cannot conclude on the relative magnitude between the two auroral jets.This particular point thus needs to be confirmed with new observations.Majeed et al. (2016) and Yates et al. (2020) also predict that the southern jets are expected to disappear around the µbar level.On the contrary, our data demonstrate that the neutrals are still flowing with a substantial counterrotation velocity at the sub-millibar level below the southern oval (and probably also below the northern one), that is, ⇠900 km below the corresponding ionospheric winds of Rego et al. (1999) and 100 500 km below the tentative H2 flow of Chaufray et al. (2011).Despite the strong signal-tonoise limitations of our CO observations at 3 mbar, we find that the southern auroral jets are at least twice slower in the millibar range than at sub-millibar levels, possibly disappearing between the sub-millibar and the millibar levels.The detection of these auroral vortices down to the submillibar level may bear crucial implications for Jovian atmo-B.Benmahi et al.: First absolute wind measurements in Saturn's stratosphere from ALMA observations

Fig. 9 .
Fig. 9. Zonally averaged eastward winds from CO and HCN observations compared with the García-Melendo et al. (2011) measurements.The eastward winds for HCN and CO are obtained by averaging the winds between the two limbs of the planet, i.e. u east u west 2 Fletcher et al. (2 found that the post-merger beacon had a westwar 1.6 ± 0.2 per day (i.e.⇠-15 m s 1 at ⇠35 N).These are thus consistent with the average wind obtained in range from our data.In Fig. 8, the eastward and westward peaks seen winds at 61 N, 55 N, 50 N, and 45 N at the weste at 55 N, 59 N, and 67 N at the eastern limb, an amplitudes exceeding 100 m s 1 do not correspond to culation because they do not have a symmetrical co the other limb.In this latitudinal range in the trop dynamics are often perturbed by the presence of other instabilities due to meridional shear of the spheric jets at mid-latitudes (see Trammell et al. 2 this pressure level, in the stratosphere, and at the the question is whether a circulation similar to th by Trammell et al. (2014) in the upper troposphere whether or not a 150 m s 1 velocity is realistic for thetical eddies.This could explain our results of non Carrión-Gonzále

Fig. 3 .
Fig. 3. Retrieval results.Left: retrieved best-fitting wind profi polynomial.Red and blue shadowed regions contain the ensem grey rectangle indicates the unobservable northern latitudes.T (Sromovsky et al. 1993).Right: wind variations with altitude at in the literature (see Sect. 5 for details).Cloud-tracking winds f as they have uncertainties of about 1000 m s 1 .Winds from Flet thermal wind equation applied to a reanalysis of the 1989 IRIS/V Our equatorial wind is about 200 +100 80 m s 1 less intens the Voyager reference.Assuming nominal probed lev ⇠1 mbar and ⇠1 bar, respectively, this indicates a +70 +30 20 wind shear per pressure decade (or @u/@z = +30 ± 10 per scale height), where the positive sign is related to the grade wind direction.At 70 S, our winds are about 70 +180 170 less intense than Voyager's.We find a much smaller wind at 70 S, although the uncertainties are larger than for e rial winds: 20 ± 60 m s 1 per pressure decade (or @u 9 ± 25 m s 1 per scale height).Our results compare we the estimates from French et al. (1998, their Fig.11b) studied the wind-shear between Voyager's cloud-tracking (which they assumed at 100 mbar) and their occultation d 0.38 mbar.French et al. (1998) determined a wind shear of +30 m s 1 per scale height at the equator and 15 m s 1 pe height at 70 S.In contrast, cloud-tracking measurements Tollefson & Pater (2018) appear somewhat at odds our estimated wind shear, as their H-band measurem assumed by the authors to probe deeper levels -indica

Figure 3 .
Figure 3. (Left) Jupiter zonal wind velocities derived from ALMA observations of HCN at ∼ 1 ′′ resolution (color map); ultraviolet auroral emission at the south pole is also shown [37, 31].(Center) Zonal wind speeds as a function of latitude in Saturn's stratosphere from Benmahi et al. [9], compared to winds from Cassini imaging.(Right) Comparison of ALMA wind speed measurements of Neptune from CO and HCN emission lines compared to Voyager cloud tracking measurements [20].Images used with permission from their respective copyright holders.

Fig. 2 .
Fig.2.Left: Herschel-HIFI observation of the CO(8´7) line in Uranus on June 15, 2012, expressed in terms of line-to-continuum ratio (l{c, black line).For each source, the models that best fit the emission core are displayed: an internal source yielding a mole fraction of 1.9 ˆ10 ´8 in the upper troposphere (red line), a steady external flux (due to IDP or a local source) of 2.2 ˆ10 5 cm ´2 s ´1 (blue line), and a comet with a diameter of 640 m depositing 3.4 ˆ10 13 g of CO above the 0.1 mbar level "300 years ago (green line).These models were computed with the thermal profile ofFeuchtgruber et al. (2013).The internal source model overestimates the line core width and produces a broad absorption that is not observed in the data.The external source models can barely be differentiated.Right: vertical profiles associated with the spectra.

3. 4 .
An upper limit on the deep O/H ratio in Uranus Thermochemistry in the deep interior of Uranus links the CO abundance to H 2 O abundance and thus to the internal O/H ratio (Fegley & Prinn 1988; Lodders & Fegley 1994) with the following net thermochemical equilibrium reaction, H 2 O `CH 4 " CO `3H 2 .
(2012) to Uranus to constrain the O/H ratio.This model accounts for C, N, and O species.We extended our thermal profile to high pressures following the dry adiabat.The profiles of Feuchtgruber et al. (2013) and Orton et al. (2013a) are similar in the upper troposphere and thus give similar deep tropospheric profiles.We constrained the O/H and C/H ratios by fitting the following upper tropospheric mole fractions with errors A33, page 4 of 6

Figure 4 .
Figure 4. Figure from Cavalié et al. [29] with permission) showing Herschel/HIFI data of the CO (J=8-7) spectrum on Uranus compared to internal, cometary, and steady external infall source models (left), and the corresponding vertical profiles (right).The comparison of Herschel data to the CO source models allows for the inference of an external source of Uranus's CO.

Figure 5 .
Figure 5. Spectral line map of CO J = 1 − 0 absorption towards Mars, obtained using the IRAM Plateau de Bure Interferometer with a beam size ∼ 7 ′′ (adapted from [103]; overlaid on an optical map of the Martian surface from the NASA Solar System Simulator; https://space.jpl.nasa.gov/).Black lines: Observations.Red lines: Line profile expected from GCM predictions.Green lines: Fit of observations with retrieved thermal profiles.

Table 1 .
Maximum Number of Beams Per Planet Area for 50 m Diffraction-Limited Aperture [59,58,101]34,135]plement studies in the optical and infrared.The temperature and circulation of Neptune's troposphere and stratosphere have been studied using the the Combined Array for Research in Millimeter-wave Astronomy (CARMA), VLA, and ALMA[93,94,82,134,135]through observations of CO and continuum features.Similarly, the NH 3 and brightness temperature distribution on Jupiter from VLA and ALMA observations has been used to contextualize storm outbreaks and observations with the Juno spacecraft[59,58,101].Figure