2018ApJ...853..122R

Transmission spectra are differential measurements that utilize stellar illumination to probe transiting exoplanet atmospheres. Any spectral difference between the illuminating light source and the disk-integrated stellar spectrum due to starspots and faculae will be imprinted in the observed transmission spectrum. However, few constraints exist for the extent of photospheric heterogeneities in M dwarfs. Here we model spot and faculae covering fractions consistent with observed photometric variabilities for M dwarfs and the associated 0.3-5.5 µm stellar contamination spectra. We find that large ranges of spot and faculae covering fractions are consistent with observations and corrections assuming a linear relation between variability amplitude, and covering fractions generally underestimate the stellar contamination. Using realistic estimates for spot and faculae covering fractions, we find that stellar contamination can be more than 10x larger than the transit depth changes expected for atmospheric features in rocky exoplanets. We also find that stellar spectral contamination can lead to systematic errors in radius and therefore the derived density of small planets. In the case of the TRAPPIST-1 system, we show that TRAPPIST-1's rotational variability is consistent with spot covering fractions f_spot=8_–7⁺¹⁸ % and faculae covering fractions f_fac=54_–46⁺¹⁶ % . The associated stellar contamination signals alter the transit depths of the TRAPPIST-1 planets at wavelengths of interest for planetary atmospheric species by roughly 1-15x the strength of planetary features, significantly complicating JWST follow-up observations of this system. Similarly, we find that stellar contamination can lead to underestimates of the bulk densities of the TRAPPIST-1 planets of Δ(ρ)=-8_–20⁺⁷ % , thus leading to overestimates of their volatile contents.