We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST, and K2 transit-time measurements, and revisit an N-body dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer light curves to derive the density of the host star and the planet densities. We find that all seven planets’ densities may be described with a single rocky mass–radius relation which is depleted in iron relative to Earth, with Fe 21 wt% versus 32 wt% for Earth, and otherwise Earth-like in composition. Alternatively, the planets may have an Earth-like composition but enhanced in light elements, such as a surface water layer or a core-free structure with oxidized iron in the mantle. We measure planet masses to a precision of 3%–5%, equivalent to a radial-velocity (RV) precision of 2.5 cm s−1, or two orders of magnitude more precise than current RV capabilities. We find the eccentricities of the planets are very small, the orbits are extremely coplanar, and the system is stable on 10 Myr timescales. We find evidence of infrequent timing outliers, which we cannot explain with an eighth planet; we instead account for the outliers using a robust likelihood function. We forecast JWST timing observations and speculate on possible implications of the planet densities for the formation, migration, and evolution of the planet system.

The TRAPPIST-1 planetary system represents an exceptional opportunity for the atmospheric characterization of temperate terrestrial exoplanets with the upcoming James Webb Space Telescope (JWST). Assessing the potential impact of stellar contamination on the planets' transit transmission spectra is an essential precursor step to this characterization. Planetary transits themselves can be used to scan the stellar photosphere and to constrain its heterogeneity through transit depth variations in time and wavelength. In this context, we present our analysis of 169 transits observed in the optical from space with K2 and from the ground with the SPECULOOS and Liverpool telescopes. Combining our measured transit depths with literature results gathered in the mid/near-IR with Spitzer/IRAC and HST/WFC3, we construct the broadband transmission spectra of the TRAPPIST-1 planets over the 0.8-4.5 µm spectral range. While planets b, d, and f spectra show some structures at the 200-300ppm level, the four others are globally flat. Even if we cannot discard their instrumental origins, two scenarios seem to be favored by the data: a stellar photosphere dominated by a few high-latitude giant (cold) spots, or, alternatively, by a few small and hot (3500-4000K) faculae. In both cases, the stellar contamination of the transit transmission spectra is expected to be less dramatic than predicted in recent papers. Nevertheless, based on our results, stellar contamination can still be of comparable or greater order than planetary atmospheric signals at certain wavelengths. Understanding and correcting the effects of stellar heterogeneity therefore appears essential to prepare the exploration of TRAPPIST-1's with JWST.

Context. The Rosetta spacecraft provided us with a unique opportunity to study comet 67P/Churyumov–Gerasimenko (67P) from a close perspective and over a 2-yr time period. Comet 67P is a weakly active comet. It was therefore unexpected to find an active and dynamic ionosphere where the cometary ions were largely dominant over the solar wind ions, even at large heliocentric distances. Aims. Our goal is to understand the different drivers of the cometary ionosphere and assess their variability over time and over the different conditions encountered by the comet during the Rosetta mission. Methods. We used a multi-instrument data-based ionospheric model to compute the total ion number density at the position of Rosetta. In-situ measurements from the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) and the Rosetta Plasma Consortium (RPC)–Ion and Electron Sensor (IES), together with the RPC–LAngmuir Probe instrument (LAP) were used to compute the local ion total number density. The results are compared to the electron densities measured by RPC–Mutual Impedance Probe (MIP) and RPC–LAP. Results. We were able to disentangle the physical processes responsible for the formation of the cometary ions throughout the 2-yr escort phase and we evaluated their respective magnitudes. The main processes are photo-ionization and electron-impact ionization. The latter is a significant source of ionization at large heliocentric distance (>2 au) and was predominant during the last 4 months of the mission. The ionosphere was occasionally subject to singular solar events, temporarily increasing the ambient energetic electron population. Solar photons were the main ionizer near perihelion at 1.3 au from the Sun, during summer 2015.

Context. With more than 1000 h of observation from Feb. 2016 to Oct. 2019, the Spitzer Exploration Program Red Worlds (ID: 13067, 13175 and 14223) exclusively targeted TRAPPIST-1, a nearby (12 pc) ultracool dwarf star, finding that it is orbited by seven transiting Earth-sized planets. At least three of these planets orbit within the classical habitable zone of the star, and all of them are well-suited for a detailed atmospheric characterization with the upcoming JWST. Aims. The main goals of the Spitzer Red Worlds program were (1) to explore the system for new transiting planets, (2) to intensively monitor the planets’ transits to yield the strongest possible constraints on their masses, sizes, compositions, and dynamics, and (3) to assess the infrared variability of the host star. In this paper, we present the global results of the project. Methods. We analyzed 88 new transits and combined them with 100 previously analyzed transits, for a total of 188 transits observed at 3.6 or 4.5 μm. For a comprehensive study, we analyzed all light curves both individually and globally. We also analyzed 29 occultations (secondary eclipses) of planet b and eight occultations of planet c observed at 4.5 μm to constrain the brightness temperatures of their daysides. Results. We identify several orphan transit-like structures in our Spitzer photometry, but all of them are of low significance. We do not confirm any new transiting planets. We do not detect any significant variation of the transit depths of the planets throughout the different campaigns. Comparing our individual and global analyses of the transits, we estimate for TRAPPIST-1 transit depth measurements mean noise floors of ~35 and 25 ppm in channels 1 and 2 of Spitzer/IRAC, respectively. We estimate that most of this noise floor is of instrumental origins and due to the large inter-pixel inhomogeneity of IRAC InSb arrays, and that the much better interpixel homogeneity of JWST instruments should result in noise floors as low as 10 ppm, which is low enough to enable the atmospheric characterization of the planets by transit transmission spectroscopy. Our analysis reveals a few outlier transits, but we cannot conclude whether or not they correspond to spot or faculae crossing events. We construct updated broadband transmission spectra for all seven planets which show consistent transit depths between the two Spitzer channels. Although we are limited by instrumental precision, the combined transmission spectrum of planet b to g tells us that their atmospheres seem unlikely to be CH4-dominated. We identify and model five distinct high energy flares in the whole dataset, and discuss our results in the context of habitability. Finally, we fail to detect occultation signals of planets b and c at 4.5 μm, and can only set 3-σ upper limits on their dayside brightness temperatures (611 K for b 586 K for c).