Between 2004 and 2017, spectral observations have been gathered by the Visual and Infrared Mapping Spectrometer (VIMS) on-board Cassini (Brown et al., 2004) during 23 Enceladus close encounters, in addition to more distant surveys. The objective of the present study is to produce a global hyperspectral mosaic of the complete VIMS data set of Enceladus in order to highlight spectral variations among the different geological units. This requires the selection of the best observations in terms of spatial resolution and illumination conditions. We have carried out a detailed investigation of the photometric behavior at several key wavelengths (1.35, 1.5, 1.65, 1.8, 2.0, 2.25, 2.55 and 3.6 ${\mu}$m), characteristics of the infrared spectra of water ice. We propose a new photometric function, based on the model of Shkuratov et al. (2011). When combined, corrected mosaics at different wavelengths reveal heterogeneous areas, in particular in the terrains surrounding the Tiger Stripes on the South Pole and in the northern hemisphere around 30{\deg}N, 90{\deg}W. Those areas appear mainly correlated to tectonized units, indicating an endogenous origin, potentially driven by seafloor hotspots.

Studying the atmospheric Planetary Boundary Layer (PBL) is crucial to understand the climate of a planet. The meteorological measurements by the instruments onboard InSight make a uniquely rich dataset to study the active turbulent dynamics of the daytime PBL on Mars. Here we use the high-sensitivity continuous pressure, wind, temperature measurements in the first 400 sols of InSight operations to analyze convective gusts, cells, and vortices in Mars' daytime PBL. We compare the InSight measurements to turbulence-resolving Large-Eddy Simulations (LES). The daytime PBL turbulence at the InSight landing site is found to be very active, with clearly identified signatures of convective cells and a vast population of almost 10,000 recorded vortex encounters, adequately represented by a power-law with a 3.5 exponent in agreement with LES results. While the daily variability of vortex encounters at the InSight landing site can be explained by the statistical nature of turbulence, the seasonal variability is strongly correlated with ambient wind speed, which is supported by LES. However, wind gustiness is more correlated to surface temperature than ambient wind speed, confirming the radiative control of the daytime martian PBL; and fewer convective vortices are forming in LES when the background wind is doubled. Thus, the long-term seasonal variability of vortex encounters at the InSight landing site is mainly controlled by the advection of convective vortices by the ambient wind speed. Typical tracks followed by vortices forming in the LES show a similar distribution in direction and length as orbital imagery of the InSight region.

We use the EAGLE simulations to study the effects of the intra-cluster medium (ICM) on the spatially resolved star-formation activity in galaxies. We study three cases of galaxy asymmetry dividing each galaxy in two halves using the plane (i) perpendicular to the $\rm \texttt{velocity}$ direction, differentiating the galaxy part approaching to the cluster center, hereafter dubbed as the "leading half", and the opposite one "trailing half", (ii) perpendicular to the $\rm \texttt{radial}$ position of the satellite to the centre of the cluster, (iii) that maximizes the star-formation rate ($\rm SFR$) difference between the two halves. For (i), we find an enhancement of the $\rm SFR$, star formation efficiency ($\rm SFE$), and interstellar medium pressure in the leading half with respect to the trailing one and normal star-forming galaxies in the EAGLE simulation, and a clear overabundance of gas particles in their trailing. These results suggest that ram pressure (RP) is boosting the star formation by gas compression in the leading half, and transporting the gas to the trailing half. This effect is more pronounced in satellites of intermediate stellar masses $\rm 10^{9.5-10.5} M_{\odot}$, with gas masses above $\rm 10^{9} M_{\odot}$, and located within one virial radius or in the most massive clusters. In (iii) we find an alignment between the velocity and the vector perpendicular to the plane that maximizes the $\rm SFR$ difference between the two halves. It suggests that finding this plane in real galaxies can provide insights into the velocity direction.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will survey galaxy and star formation history over cosmological time scales. Rather than identifying individual objects, EXCLAIM will be a pathfinder to demonstrate an intensity mapping approach, which measures the cumulative redshifted line emission. EXCLAIM will operate at 420-540 GHz with a spectral resolution R=512 to measure the integrated CO and [CII] in redshift windows spanning 0 < z < 3.5. CO and [CII] line emissions are key tracers of the gas phases in the interstellar medium involved in star-formation processes. EXCLAIM will shed light on questions such as why the star formation rate declines at z < 2, despite continued clustering of the dark matter. The instrument will employ an array of six superconducting integrated grating-analog spectrometers (micro-spec) coupled to microwave kinetic inductance detectors (MKIDs). Here we present an overview of the EXCLAIM instrument design and status.

The texture, composition, and morphology of dunes observed in the equatorial regions of Titan may reflect present and/or past climatic conditions. Determining the physio-chemical properties and the morphodynamics of Titan's dunes is therefore essential to understanding of the climatic and geological history of the largest moon of Saturn. We quantitatively derived average surface properties of dune and interdune areas (texture, composition) from modeling of the microwave backscattered signal and Monte-Carlo inversion of the despeckled Cassini/SAR data over Titan's three largest sand seas: Belet, Shangri-La and Fensal. We present the first analysis of the backscatter functions extracted from despeckled SAR images that cover such a large range in incidence angles, including data from the beginning of the Cassini mission up to its Grand Finale. We show that dunes and interdunes have significantly different physical properties. Dunes are found to be more microwave absorbent than interdunes. Additionally, potential secondary bedforms, such as ripples and avalanches, may have been detected, providing potential evidence for currently active dunes and sediment transport. Our modelling shows that the interdunes have multi-scale roughnesses with higher dielectric constants than the dunes which have a low dielectric constant consistent with organic sand. The radar brightness of the interdunes can be explained by the presence of a shallow layer of significantly larger organic grains, possibly non-mobilized by the winds. {Together, our} findings suggest that Titan's sand seas evolve under the current multi-directional wind regimes with dunes that elongate with their crests aligned in the residual drift direction.

The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne telescope mission to search for inflationary gravitational waves from the early universe. PIPER employs two 32x40 arrays of superconducting transition-edge sensors, which operate at 100 mK. An open bucket dewar of liquid helium maintains the receiver and telescope optics at 1.7 K. We describe the thermal design of the receiver and sub-kelvin cooling with a continuous adiabatic demagnetization refrigerator (CADR). The CADR operates between 70-130 mK and provides ~10 uW cooling power at 100 mK, nearly five times the loading of the two detector assemblies. We describe electronics and software to robustly control the CADR, overall CADR performance in flight-like integrated receiver testing, and practical considerations for implementation in the balloon float environment.

Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of numerical simulations for merging black holes. The catalog contains 2018 distinct configurations (a factor of 11 increase compared to the 2013 SXS catalog), including 1426 spin-precessing configurations, with mass ratios between 1 and 10, and spin magnitudes up to 0.998. The median length of a waveform in the catalog is 39 cycles of the dominant $\ell=m=2$ gravitational-wave mode, with the shortest waveform containing 7.0 cycles and the longest 351.3 cycles. We discuss improvements such as correcting for moving centers of mass and extended coverage of the parameter space. We also present a thorough analysis of numerical errors, finding typical truncation errors corresponding to a waveform mismatch of $\sim 10^{-4}$. The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% ($90^{\text{th}}$ percentile), about an order of magnitude better than analytical models for remnant properties. The full catalog is publicly available at https://www.black-holes.org/waveforms .

In this white paper, we present a cross-section of important scientific questions that remain partially or completely unanswered, ranging from Titan exosphere to the deep interior, and we detail which instrumentation and mission scenarios should be used to answer them. Our intention is to formulate the science goals for the next generation of planetary missions to Titan in order to prepare the future exploration of the moon. The ESA L-class mission concept that we propose is composed of a Titan orbiter and at least an in situ element (lake lander and/or drone(s)).

We have analyzed the complete Visual and Infrared Mapping Spectrometer (VIMS) data archive of Titan. Our objective is to build global surface cartographic products, by combining all the data gathered during the 127 targeted flybys of Titan into synthetic global maps interpolated on a grid at 32 pixels per degree (~1.4 km/pixel at the equator), in seven infrared spectral atmospheric windows. Multispectral summary images have been computed for each single VIMS cube in order to rapidly identify their scientific content and assess their quality. These summary images are made available to the community on a public website (vims.univ-nantes.fr). The global mapping work faced several challenges due to the strong absorbing and scattering effects of the atmosphere coupled to the changing observing conditions linked to the orbital tour of the Cassini mission. We determined a surface photometric function which accounts for variations in incidence, emergence and phase angles, and which is able to mitigate brightness variations linked to the viewing geometry of the flybys. The atmospheric contribution has been reduced using the subtraction of the methane absorption band wings, considered as proxies for atmospheric haze scattering. We present a new global three color composite map of band ratios (red: 1.59/1.27 {\mu}m; green: 2.03/1.27 {\mu}m; blue: 1.27/1.08 {\mu}m), which has also been empirically corrected from an airmass (the solar photon path length through the atmosphere) dependence. This map provides a detailed global color view of Titan's surface partially corrected from the atmosphere and gives a global insight of the spectral variability, with the equatorial dunes fields appearing in brownish tones, and several occurrences of bluish tones localized in areas such as Sinlap, Menvra and Selk craters. (...)

We have analyzed the complete archive of the Visual and Infrared Mapping Spectrometer (VIMS) data in order to monitor and analyze the evolution of the clouds and haze coverage at both poles of Titan during the entire Cassini mission. Our objective is to give a cartographic synopsis from a VIMS perspective, to provide a global view of the seasonal evolution of Titan's atmosphere over the poles. We leave the detailed comparison with the Imaging Science Subsystem (ISS) and the Composite Infrared Spectrometer (CIRS) data sets to further studies. We have computed global hyperspectral mosaics for each of the 127 targeted flybys of Titan to produce synthetic color maps emphasizing the main atmospheric features. The north pole appears fully covered by a huge cloud as soon as the first observations in 2004 and up to the equinox in 2009 (Le Mou\'elic et al. 2012). The northern skies then became progressively clearer, after the circulation turnover in 2009, revealing the underlying lakes and seas to the optical instruments up to 2017. The reverse situation is observed over the south pole, which was mostly clear of such a high obscuring cloud during the first years of the mission, but started to develop a polar cloud in 2012. This feature grew up month after month until the end of the mission in 2017, with a poleward latitudinal extent of 75$^\circ$S in 2013 up to 58$^\circ$S in April 2017. Thanks to the spectral capabilities of VIMS, we have detected HCN spectral signatures over the north pole in almost all flybys between 2004 and 2008. These HCN signatures started then to show up over the south pole in almost all flybys between 2012 and 2017, so perfectly matching the timing and spatial extent of the northern and southern polar atmospheric features.

The Primordial Inflation Polarization ExploreR (PIPER) is a balloon-borne telescope designed to measure the polarization of the Cosmic Microwave Background on large angular scales. PIPER will map 85% of the sky at 200, 270, 350, and 600 GHz over a series of 8 conventional balloon flights from the northern and southern hemispheres. The first science flight will use two 32x40 arrays of backshort-under-grid transition edge sensors, multiplexed in the time domain, and maintained at 100 mK by a Continuous Adiabatic Demagnetization Refrigerator. Front-end cryogenic Variable-delay Polarization Modulators provide systematic control by rotating linear to circular polarization at 3 Hz. Twin telescopes allow PIPER to measure Stokes I, Q, U, and V simultaneously. The telescope is maintained at 1.5 K in an LHe bucket dewar. Cold optics and the lack of a warm window permit sensitivity at the sky-background limit. The ultimate science target is a limit on the tensor-to-scalar ratio of r ~ 0.007, from the reionization bump to l ~ 300. PIPER's first flight will be from the Northern hemisphere, and overlap with the CLASS survey at lower frequencies. We describe the current status of the PIPER instrument.

The James Webb Space Telescope (JWST), scheduled for launch in 2018, is the successor to the Hubble Space Telescope (HST) but with a significantly larger aperture (6.5 m) and advanced instrumentation focusing on infrared science (0.6-28.0 $\mu$m ). In this paper we examine the potential for scientific investigation of Titan using JWST, primarily with three of the four instruments: NIRSpec, NIRCam and MIRI, noting that science with NIRISS will be complementary. Five core scientific themes are identified: (i) surface (ii) tropospheric clouds (iii) tropospheric gases (iv) stratospheric composition and (v) stratospheric hazes. We discuss each theme in depth, including the scientific purpose, capabilities and limitations of the instrument suite, and suggested observing schemes. We pay particular attention to saturation, which is a problem for all three instruments, but may be alleviated for NIRCam through use of selecting small sub-arrays of the detectors - sufficient to encompass Titan, but with significantly faster read-out times. We find that JWST has very significant potential for advancing Titan science, with a spectral resolution exceeding the Cassini instrument suite at near-infrared wavelengths, and a spatial resolution exceeding HST at the same wavelengths. In particular, JWST will be valuable for time-domain monitoring of Titan, given a five to ten year expected lifetime for the observatory, for example monitoring the seasonal appearance of clouds. JWST observations in the post-Cassini period will complement those of other large facilities such as HST, ALMA, SOFIA and next-generation ground-based telescopes (TMT, GMT, EELT).

Titan's equatorial regions are covered by eastward propagating linear dunes. This direction is opposite to mean surface winds simulated by Global Climate Models (GCMs), which are oriented westward at these latitudes, similar to trade winds on Earth. Different hypotheses have been proposed to address this apparent contradiction, involving Saturn's gravitational tides, large scale topography or wind statistics, but none of them can explain a global eastward dune propagation in the equatorial band. Here we analyse the impact of equinoctial tropical methane storms developing in the superrotating atmosphere (i.e. the eastward winds at high altitude) on Titan's dune orientation. Using mesoscale simulations of convective methane clouds with a GCM wind profile featuring superrotation, we show that Titan's storms should produce fast eastward gust fronts above the surface. Such gusts dominate the aeolian transport, allowing dunes to extend eastward. This analysis therefore suggests a coupling between superrotation, tropical methane storms and dune formation on Titan. Furthermore, together with GCM predictions and analogies to some terrestrial dune fields, this work provides a general framework explaining several major features of Titan's dunes: linear shape, eastward propagation and poleward divergence, and implies an equatorial origin of Titan's dune sand.

Water-ice-poor, 5-$\mu$m-bright material on Saturn's moon Titan has previously been geomorphologically identified as evaporitic. Here we present a global distribution of the occurrences of the 5-$\mu$m-bright spectral unit, identified with Cassini's Visual Infrared Mapping Spectrometer (VIMS) and examined with RADAR when possible. We explore the possibility that each of these occurrences are evaporite deposits. The 5-$\mu$m-bright material covers 1\% of Titan's surface and is not limited to the poles (the only regions with extensive, long-lived surface liquid). We find the greatest areal concentration to be in the equatorial basins Tui Regio and Hotei Regio. Our interpretations, based on the correlation between 5-$\mu$m-bright material and lakebeds, imply that there was enough liquid present at some time to create the observed 5-$\mu$m-bright material. We address the climate implications surrounding a lack of evaporitic material at the south polar basins: if the south pole basins were filled at some point in the past, then where is the evaporite?

The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne cosmic microwave background (CMB) polarimeter designed to search for evidence of inflation by measuring the large-angular scale CMB polarization signal. BICEP2 recently reported a detection of B-mode power corresponding to the tensor-to-scalar ratio r = 0.2 on ~2 degree scales. If the BICEP2 signal is caused by inflationary gravitational waves (IGWs), then there should be a corresponding increase in B-mode power on angular scales larger than 18 degrees. PIPER is currently the only suborbital instrument capable of fully testing and extending the BICEP2 results by measuring the B-mode power spectrum on angular scales $\theta$ = ~0.6 deg to 90 deg, covering both the reionization bump and recombination peak, with sensitivity to measure the tensor-to-scalar ratio down to r = 0.007, and four frequency bands to distinguish foregrounds. PIPER will accomplish this by mapping 85% of the sky in four frequency bands (200, 270, 350, 600 GHz) over a series of 8 conventional balloon flights from the northern and southern hemispheres. The instrument has background-limited sensitivity provided by fully cryogenic (1.5 K) optics focusing the sky signal onto four 32x40-pixel arrays of time-domain multiplexed Transition-Edge Sensor (TES) bolometers held at 140 mK. Polarization sensitivity and systematic control are provided by front-end Variable-delay Polarization Modulators (VPMs), which rapidly modulate only the polarized sky signal at 3 Hz and allow PIPER to instantaneously measure the full Stokes vector (I, Q, U, V) for each pointing. We describe the PIPER instrument and progress towards its first flight.

Simulations of Titan's atmospheric transmission and surface reflectivity have been developed in order to estimate how Titan's atmosphere and surface properties could affect performances of the Cassini radar experiment. In this paper we present a selection of models for Titan's haze, vertical rain distribution, and surface composition implemented in our simulations. We collected dielectric constant values for the Cassini radar wavelength ($\sim 2.2$ cm) for materials of interest for Titan: liquid methane, liquid mixture of methane-ethane, water ice and light hydrocarbon ices. Due to the lack of permittivity values for Titan's haze particles in the microwave range, we performed dielectric constant ($\varepsilon_r$) measurements around 2.2 cm on tholins synthesized in laboratory. We obtained a real part of $\varepsilon_r$ in the range of 2-2.5 and a loss tangent between $10^{-3}$ and $5.10^{-2}$. By combining aerosol distribution models (with hypothetical condensation at low altitudes) to surface models, we find the following results: (1) Aerosol-only atmospheres should cause no loss and are essentially transparent for Cassini radar, as expected by former analysis. (2) However, if clouds are present, some atmospheric models generate significant attenuation that can reach $-50 dB$, well below the sensitivity threshold of the receiver. In such cases, a $13.78 GHz$ radar would not be able to measure echoes coming from the surface. We thus warn about possible risks of misinterpretation if a \textquotedblleft wet atmosphere\textquotedblright $ $is not taken into account. (3) Rough surface scattering leads to a typical response of $\sim -17 dB$. These results will have important implications on future Cassini radar data analysis.

Clouds on Titan result from the condensation of methane and ethane and, as on other planets, are primarily structured by circulation of the atmosphere. At present, cloud activity mainly occurs in the southern (summer) hemisphere, arising near the pole and at mid-latitudes from cumulus updrafts triggered by surface heating and/or local methane sources, and at the north (winter) pole, resulting from the subsidence and condensation of ethane-rich air into the colder troposphere. General circulation models predict that this distribution should change with the seasons on a 15-year timescale, and that clouds should develop under certain circumstances at temperate latitudes (~40\degree) in the winter hemisphere. The models, however, have hitherto been poorly constrained and their long-term predictions have not yet been observationally verified. Here we report that the global spatial cloud coverage on Titan is in general agreement with the models, confirming that cloud activity is mainly controlled by the global circulation. The non-detection of clouds at latitude ~40\degree N and the persistence of the southern clouds while the southern summer is ending are, however, both contrary to predictions. This suggests that Titan's equator-to-pole thermal contrast is overestimated in the models and that its atmosphere responds to the seasonal forcing with a greater inertia than expected.

Titan is one of the primary scientific objectives of the NASA ESA ASI Cassini Huygens mission. Scattering by haze particles in Titan's atmosphere and numerous methane absorptions dramatically veil Titan's surface in the visible range, though it can be studied more easily in some narrow infrared windows. The Visual and Infrared Mapping Spectrometer (VIMS) instrument onboard the Cassini spacecraft successfully imaged its surface in the atmospheric windows, taking hyperspectral images in the range 0.4 5.2 ?m. On 26 October (TA flyby) and 13 December 2004 (TB flyby), the Cassini Huygens mission flew over Titan at an altitude lower than 1200 km at closest approach. We report here on the analysis of VIMS images of the Huygens landing site acquired at TA and TB, with a spatial resolution ranging from 16 to14.4 km/pixel. The pure atmospheric backscattering component is corrected by using both an empirical method and a first-order theoretical model. Both approaches provide consistent results. After the removal of scattering, ratio images reveal subtle surface heterogeneities. A particularly contrasted structure appears in ratio images involving the 1.59 and 2.03 ?m images north of the Huygens landing site. Although pure water ice cannot be the only component exposed at Titan's surface, this area is consistent with a local enrichment in exposed water ice and seems to be consistent with DISR/Huygens images and spectra interpretations. The images show also a morphological structure that can be interpreted as a 150 km diameter impact crater with a central peak.

In this paper, we characterize the morphology of the disk-integrated phase functions of satellites and rings around the giant planets of our Solar System. We find that the shape of the phase function is accurately represented by a logarithmic model (Bobrov, 1970, in Surfaces and Interiors of Planets and Satellites, Academic, edited by A. Dollfus). For practical purposes, we also parametrize the phase curves by a linear-exponential model (Kaasalainen et al., 2001, Journal of Quantitative Spectroscopy and Radiative Transfer, 70, 529-543) and a simple linear-by-parts model (Lumme and Irvine, 1976, Astronomical Journal, 81, 865-893), which provides three morphological parameters : the amplitude A and the Half-Width at Half-Maximum (HWHM) of the opposition surge, and the slope S of the linear part of the phase function at larger phase angles. Our analysis demonstrates that all of these morphological parameters are correlated with the single scattering albedos of the surfaces. By taking more accurately into consideration the finite angular size of the Sun, we find that the Galilean, Saturnian, Uranian and Neptunian satellites have similar HWHMs (0.5 degrees), whereas they have a wide range of amplitudes A. The Moon has the largest HWHM (2 degrees). We interpret that as a consequence of the solar size bias, via the finite size of the Sun which varies dramatically from the Earth to Neptune. By applying a new method that attempts to morphologically deconvolve the phase function to the solar angular size, we find that icy and young surfaces, with active resurfacing, have the smallest values of A and HWHM, whereas dark objects (and perhaps older surfaces) such as the Moon, Nereid and Saturn C ring have the largest A and HWHM.