We present the discovery of three sub-Neptune-sized planets transiting the nearby and bright Sun-like star HD 191939 (TIC 269701147, TOI 1339), a $K_{s}=7.18$ magnitude G8 V dwarf at a distance of only 54 parsecs. We validate the planetary nature of the transit signals by combining five months of data from the Transiting Exoplanet Survey Satellite with follow-up ground-based photometry, archival optical images, radial velocities, and high angular resolution observations. The three sub-Neptunes have similar radii ($R_{b} = 3.42^{+0.11}_{-0.11}\,R_{\oplus}$, $R_{c}=3.23_{-0.11}^{+0.11}\,R_{\oplus}$, and $R_{d}=3.16_{-0.11}^{+0.11}\,R_{\oplus}$) and their orbits are consistent with a stable, circular, and co-planar architecture near mean motion resonances of 1:3 and 3:4 ($P_{b}=8.88$ days, $P_{c}=28.58$ days, and $P_{d}=38.35$ days). The HD~191939 system is an excellent candidate for precise mass determinations of the planets with high-resolution spectroscopy due to the host star's brightness and low chromospheric activity. Moreover, the system's compact and near-resonant nature can provide an independent way to measure planetary masses via transit timing variations while also enabling dynamical and evolutionary studies. Finally, as a promising target for multi-wavelength transmission spectroscopy of all three planets' atmospheres, HD 191939 can offer valuable insight into multiple sub-Neptunes born from a proto-planetary disk that may have resembled that of the early Sun.

We report on the discovery of a transiting Earth-sized (0.95$R_\oplus$) planet around an M3.5 dwarf star at 57\,pc, EPIC~249631677. The planet has a period of $\sim$3.14 days, i.e. ${\sim}\pi$, with an instellation of 7.5\,S$_{\oplus}$. The detection was made using publicly available data from \textit{K2}'s Campaign 15. We observed three additional transits with SPECULOOS Southern and Northern Observatories, and a stellar spectrum from Keck/HIRES, which allowed us to validate the planetary nature of the signal. The confirmed planet is well suited for comparative terrestrial exoplanetology. While exoplanets transiting ultracool dwarfs present the best opportunity for atmospheric studies of terrestrial exoplanets with the \textit{James Webb Space Telescope}, those orbiting mid-M dwarfs within 100\,pc such as EPIC~249631677b will become increasingly accessible with the next generation of observatories (e.g., \textit{HabEx, LUVOIR, OST}).

{We investigate the atmospheric composition of the long period ($P_{\rm orb}=$~10 days), eccentric exo-Saturn WASP-117b. WASP-117b could be in atmospheric temperature and chemistry similar to WASP-107b. In mass and radius WASP-117b is similar to WASP-39b, which allows a comparative study of these planets.} % methods heading (mandatory) { We analyze a near-infrared transmission spectrum of WASP-117b taken with Hubble Space Telescope/WFC3 G141, which was reduced with two independent pipelines to ensure a robust detection of the $3 \sigma$-water spectrum. High resolution measurements were taken with VLT/ESPRESSO in the optical.} % results heading (mandatory) {Using a 1D atmosphere model with isothermal temperature, uniform cloud deck and equilibrium chemistry, the Bayesian evidence of a retrieval analysis of the transmission spectrum indicates a preference for a high atmospheric metallicity ${\rm [Fe/H]}=2.58^{+0.26}_{-0.37}$ and clear skies. The data are also consistent with a lower-metallicity composition ${\rm [Fe/H]}<1.75$ and a cloud deck between $10^{-2.2} - 10^{-5.1}$~bar, but with weaker Bayesian preference. We retrieve a low CH$_4$ abundance of $<10^{-4}$ volume fraction within $1 \sigma$. We cannot constrain the temperature between theoretically imposed limits of 700 and 1000~K. Further observations are needed to confirm quenching of CH$_4$ with $K_{zz}\geq 10^8$~cm$^2$/s. We report indications of Na and K in the VLT/ESPRESSO high resolution spectrum with substantial Bayesian evidence in combination with HST data.}

We present the confirmation of two new planets transiting the nearby mid-M dwarf LTT 3780 (TIC 36724087, TOI-732, $V=13.07$, $K_s=8.204$, $R_s$=0.374 R$_{\odot}$, $M_s$=0.401 M$_{\odot}$, d=22 pc). The two planet candidates are identified in a single TESS sector and are validated with reconnaissance spectroscopy, ground-based photometric follow-up, and high-resolution imaging. With measured orbital periods of $P_b=0.77$ days, $P_c=12.25$ days and sizes $r_{p,b}=1.33\pm 0.07$ R$_{\oplus}$, $r_{p,c}=2.30\pm 0.16$ R$_{\oplus}$, the two planets span the radius valley in period-radius space around low mass stars thus making the system a laboratory to test competing theories of the emergence of the radius valley in that stellar mass regime. By combining 63 precise radial-velocity measurements from HARPS and HARPS-N, we measure planet masses of $m_{p,b}=2.62^{+0.48}_{-0.46}$ M$_{\oplus}$ and $m_{p,c}=8.6^{+1.6}_{-1.3}$ M$_{\oplus}$, which indicates that LTT 3780b has a bulk composition consistent with being Earth-like, while LTT 3780c likely hosts an extended H/He envelope. We show that the recovered planetary masses are consistent with predictions from both photoevaporation and from core-powered mass loss models. The brightness and small size of LTT 3780, along with the measured planetary parameters, render LTT 3780b and c as accessible targets for atmospheric characterization of planets within the same planetary system and spanning the radius valley.

We report the discovery and validation of four extrasolar planets hosted by the nearby, bright, Sun-like (G3V) star HD~108236 using data from the Transiting Exoplanet Survey Satellite (TESS). We present transit photometry, reconnaissance and precise Doppler spectroscopy as well as high-resolution imaging, to validate the planetary nature of the objects transiting HD~108236, also known as the TESS Object of Interest (TOI) 1233. The innermost planet is a possibly-rocky super-Earth with a period of $3.79523_{-0.00044}^{+0.00047}$ days and has a radius of $1.586\pm0.098$ $R_\oplus$. The outer planets are sub-Neptunes, with potential gaseous envelopes, having radii of $2.068_{-0.091}^{+0.10}$ $R_\oplus$, $2.72\pm0.11$ $R_\oplus$, and $3.12_{-0.12}^{+0.13}$ $R_\oplus$ and periods of $6.20370_{-0.00052}^{+0.00064}$ days, $14.17555_{-0.0011}^{+0.00099}$ days, and $19.5917_{-0.0020}^{+0.0022}$ days, respectively. With V and K$_{\rm s}$ magnitudes of 9.2 and 7.6, respectively, the bright host star makes the transiting planets favorable targets for mass measurements and, potentially, for atmospheric characterization via transmission spectroscopy. HD~108236 is the brightest Sun-like star in the visual (V) band known to host four or more transiting exoplanets. The discovered planets span a broad range of planetary radii and equilibrium temperatures, and share a common history of insolation from a Sun-like star ($R_\star = 0.888 \pm 0.017$ R$_\odot$, $T_{\rm eff} = 5730 \pm 50$ K), making HD 108236 an exciting, opportune cosmic laboratory for testing models of planet formation and evolution.

Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches.

The upcoming launch of the James Webb Space Telescope (JWST) combined with the unique features of the TRAPPIST-1 planetary system should enable the young field of exoplanetology to enter into the realm of temperate Earth-sized worlds. Indeed, the proximity of the system (12pc) and the small size (0.12 Rsun) and luminosity (0.05 Lsun) of its host star should make the comparative atmospheric characterization of its seven transiting planets within reach of an ambitious JWST program. Given the limited lifetime of JWST, the ecliptic location of the star that limits its visibility to 100d per year, the large number of observational time required by this study, and the numerous observational and theoretical challenges awaiting it, its full success will critically depend on a large level of coordination between the involved teams and on the support of a large community. In this context, we present here a community initiative aiming to develop a well-defined sequential structure for the study of the system with JWST and to coordinate on every aspect of its preparation and implementation, both on the observational (e.g. study of the instrumental limitations, data analysis techniques, complementary space-based and ground-based observations) and theoretical levels (e.g. model developments and comparison, retrieval techniques, inferences). Depending on the outcome of the first phase of JWST observations of the planets, this initiative could become the seed of a major JWST Legacy Program devoted to the study of TRAPPIST-1.

We present a new ground-based visual transmission spectrum of the hot Jupiter WASP-43b, obtained as part of the ACCESS Survey. The spectrum was derived from four transits observed between 2015 and 2018, with combined wavelength coverage between 5,300 \r{A}-9,000 \r{A} and an average photometric precision of 708 ppm in 230 \r{A} bins. We perform an atmospheric retrieval of our transmission spectrum combined with literature HST/WFC3 observations to search for the presence of clouds/hazes as well as Na, K, H$\alpha$, and $\rm H_2O$ planetary absorption and stellar spot contamination over a combined spectral range of 5,318 \r{A}-16,420 \r{A}. We do not detect a statistically significant presence of Na I or K I alkali lines, or H$\alpha$ in the atmosphere of WASP-43b. We find that the observed transmission spectrum can be best explained by a combination of heterogeneities on the photosphere of the host star and a clear planetary atmosphere with $\rm H_2O$. This model yields a log-evidence of $8.26\pm0.42$ higher than a flat (featureless) spectrum. In particular, the observations marginally favor the presence of large, low-contrast spots over the four ACCESS transit epochs with an average covering fraction $f_\text{het} = 0.27^{+0.42}_{-0.16}$ and temperature contrast $\Delta T = 132\text{ K} \pm 132\text{ K}$. Within the planet's atmosphere, we recover a log $\rm H_2O$ volume mixing ratio of $-2.78^{+1.38}_{-1.47}$, which is consistent with previous $\rm H_2O$ abundance determinations for this planet.

As we begin to discover rocky planets in the habitable zone of nearby stars with missions like TESS and CHEOPS, we will need quick advancements on instrumentation and observational techniques that will enable us to answer key science questions, such as What are the atmospheric characteristics of habitable zone rocky planets? How common are Earth-like biosignatures in rocky planets?} How similar or dissimilar are those planets to Earth? Over the next decade we expect to have discovered several Earth-analog candidates, but we will not have the tools to study the atmospheres of all of them in detail. Ground-based ELTs can identify biosignatures in the spectra of these candidate exo-Earths and understand how the planets' atmospheres compare to the Earth at different epochs. Transit spectroscopy, high-resolution spectroscopy, and reflected-light direct imaging on ELTs can identify multiple biosignatures for habitable zone, rocky planets around M stars at optical to near-infrared wavelengths. Thermal infrared direct imaging can detect habitable zone, rocky planets around AFGK stars, identifying ozone and motivating reflected-light follow-up observations with NASA missions like HabEx/LUVOIR. Therefore, we recommend that the Astro2020 Decadal Survey Committee support: (1) the search for Earth-like biosignatures on rocky planets around nearby stars as a key science case; (2) the construction over the next decade of ground-based Extremely Large Telecopes (ELTs), which will provide the large aperture and spatial resolution necessary to start revealing the atmospheres of Earth-analogues around nearby stars; (3) the development of instrumentation that optimizes the detection of biosignatures; and (4) the generation of accurate line lists for potential biosignature gases, which are needed as model templates to detect those molecules.

Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify and discuss four primary opportunities: (1) refining stellar magnetic active region properties through exoplanet crossing events; (2) spectral decomposition of active exoplanet host stars; (3) joint retrievals of stellar photospheric and planetary atmospheric properties with studies of transmission spectra; and (4) continued visual transmission spectroscopy studies to complement longer-wavelength studies from $\textit{JWST}$. We make five recommendations to the Astro2020 Decadal Survey Committee: (1) identify the transit light source (TLS) effect as a challenge to precise exoplanet transmission spectroscopy and an opportunity ripe for scientific advancement in the coming decade; (2) include characterization of host star photospheric heterogeneity as part of a comprehensive research strategy for studying transiting exoplanets; (3) support the construction of ground-based extremely large telescopes (ELTs); (4) support multi-disciplinary research teams that bring together the heliophysics, stellar physics, and exoplanet communities to further exploit transiting exoplanets as spatial probes of stellar photospheres; and (5) support visual transmission spectroscopy efforts as complements to longer-wavelength observational campaigns with $\textit{JWST}$.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

We present an optical transmission spectrum of the atmosphere of WASP-4b obtained through observations of four transits with Magellan/IMACS. Using a Bayesian approach to atmospheric retrieval, we find no evidence for scattering or absorption features in our transit spectrum. Our models include a component to model the transit light source effect (spectral contamination from unocculted spots on the stellar photosphere), which we show can have a marked impact on the observed transmission spectrum for reasonable spot covering fractions (< 5%); this is the first such analysis for WASP-4b. We are also able to fit for the size and temperature contrast of spots observed during the second and third transits, finding evidence for both small, cool and large, warm spot-like features on the photosphere. Finally, we compare our results to those published by Huitson et al. (2017) using Gemini/GMOS and May et al. (2018) using IMACS, and find that our data are in agreement.

Transmission spectra probe exoplanetary atmospheres, but they can also be strongly affected by heterogeneities in host star photospheres through the transit light source effect. Here we build upon our recent study of the effects of unocculted spots and faculae on M-dwarf transmission spectra, extending the analysis to FGK dwarfs. Using a suite of rotating model photospheres, we explore spot and facula covering fractions for varying activity levels and the associated stellar contamination spectra. Relative to M dwarfs, we find that the typical variabilities of FGK dwarfs imply lower spot covering fractions, though they generally increase with later spectral types, from $\sim 0.1\%$ for F dwarfs to 2-4$\%$ for late-K dwarfs. While the stellar contamination spectra are considerably weaker than those for typical M dwarfs, we find that typically active G and K dwarfs produce visual slopes that are detectable in high-precision transmission spectra. We examine line offsets at H$\alpha$ and the Na and K doublets and find that unocculted faculae in K dwarfs can appreciably alter transit depths around the Na D doublet. We find that band-averaged transit depth offsets at molecular bands for CH$_{4}$, CO, CO$_{2}$, H$_{2}$O, N$_{2}$O, O$_{2}$, and O$_{3}$ are not detectable for typically active FGK dwarfs, though stellar TiO/VO features are potentially detectable for typically active late-K dwarfs. Generally, this analysis shows that inactive FGK dwarfs do not produce detectable stellar contamination features in transmission spectra, though active FGK host stars can produce such features and care is warranted in interpreting transmission spectra from these systems.

KIC 12557548 b is first of a growing class of intriguing disintegrating planet candidates, which lose mass in the form of a metal rich vapor that condenses into dust particles. Here, we follow up two perplexing observations of the system: 1) the transits appeared shallower than average in 2013 and 2014 and 2) the parameters derived from a high resolution spectrum of the star differed from other results using photometry and low resolution spectroscopy. We observe 5 transits of the system with the 61-inch Kuiper telescope in 2016 and show that they are consistent with photometry from the Kepler spacecraft in 2009-2013, suggesting that the dusty tail has returned to normal length and mass. We also evaluate high resolution archival spectra from the Subaru HDS spectrograph and find them to be consistent with a main-sequence Teff=4440 +/- 70 K star in agreement with the photometry and low resolution spectroscopy. This disfavors the hypothesis that planet disintegration affected the analysis of prior high resolution spectra of this star. We apply Principal Component Analysis to the Kepler long cadence data to understand the modes of disintegration. There is a tentative 491 day periodicity of the second principal component, which corresponds to possible long-term evolution of the dust grain sizes, though the mechanism on such long timescales remains unclear.

The short period ($0.94$-day) transiting exoplanet WASP-19b is an exceptional target for transmission spectroscopy studies, due to its relatively large atmospheric scale-height ($\sim 500$ km) and equilibrium temperature ($\sim 2100$ K). Here we report on six precise spectroscopic Magellan/IMACS observations, five of which target the full optical window from $0.45-0.9\mu$m and one targeting the $0.4-0.55\mu$m blue-optical range. Five of these datasets are consistent with a transmission spectrum without any significant spectral features, while one shows a significant slope as a function of wavelength, which we interpret as arising from photospheric heterogeneities in the star. Coupled with HST/WFC3 infrared observations, our optical/near-infrared measurements point to the presence of high altitude clouds in WASP-19b's atmosphere in agreement with previous studies. Using a semi-analytical retrieval approach, considering both planetary and stellar spectral features, we find a water abundance consistent with solar for WASP-19b and strong evidence for sub-solar abundances for optical absorbers such as TiO and Na; no strong optical slope is detected, which suggests that if hazes are present, they are much weaker than previously suggested. In addition, two spot-crossing events are observed in our datasets and analyzed, including one of the first unambiguously detected bright spot-crossing events on an exoplanet host star.

The seven approximately Earth-sized transiting planets in the \object{TRAPPIST-1} system provide a unique opportunity to explore habitable zone and non-habitable zone small planets within the same system. Its habitable zone exoplanets -- due to their favorable transit depths -- are also worlds for which atmospheric transmission spectroscopy is within reach with the Hubble Space Telescope (HST) and with the James Webb Space Telescope (JWST). We present here an independent reduction and analysis of two \textit{HST} Wide Field Camera 3 (WFC3) near-infrared transit spectroscopy datasets for six planets (b through g). Utilizing our physically-motivated detector charge trap correction and a custom cosmic ray correction routine, we confirm the general shape of the transmission spectra presented by \textbf{\citet{deWit2016, deWit2018}}. Our data reduction approach leads to a 25\% increase in the usable data and reduces the risk of confusing astrophysical brightness variations (e.g., flares) with instrumental systematics. No prominent absorption features are detected in any individual planet's transmission spectra; by contrast, the combined spectrum of the planets shows a suggestive decrease around 1.4\,$\micron$ similar to an inverted water absorption feature. Including transit depths from \textit{K2}, the SPECULOOS-South Observatory, and \textit{Spitzer}, we find that the complete transmission spectrum is fully consistent with stellar contamination owing to the transit light source effect. These spectra demonstrate how stellar contamination can overwhelm planetary absorption features in low-resolution exoplanet transit spectra obtained by \textit{HST} and \textit{JWST} and also highlight the challenges in combining multi epoch observations for planets around rapidly rotating spotted stars.

Transmission spectroscopy provides a powerful probe of the atmospheric properties of transiting exoplanets. To date, studies of exoplanets in transit have focused on inferring their atmospheric properties such as chemical compositions, cloud/haze properties, and temperature structures. However, surface inhomogeneities in the host stars of exoplanets in the form of cool spots and hot faculae can in principle imprint signatures on the observed planetary transit spectrum. Here we present Aura, a new retrieval paradigm for inferring both planetary and stellar properties from a transmission spectrum. We apply our retrieval framework to a sample of hot giant exoplanets to determine the significance of stellar heterogeneity and clouds/hazes in their spectra. The retrieval analyses distinguish four groups of planets. First, the spectra of WASP-6b and WASP-39b are best characterised by imprints of stellar heterogeneity and hazes and/or clouds. HD 209458b and HAT-P-12b comprise the second group for which there is weak evidence for stellar heterogeneity and a high significance of hazes and/or clouds. The third group constitutes HAT-P-1b and WASP-31b and shows weak evidence against stellar heterogeneity but weak to substantial indications of clouds/hazes. The fourth group -- WASP-19b, WASP-17b, and WASP-12b -- is fit best by molecular and alkali absorbers with H$_2$ scattering without evidence for stellar heterogeneity and weak to no evidence for clouds/hazes. Our retrieval methodology paves the way to simultaneous information on the star and planet from higher resolution spectra using future facilities such as the James Webb Space Telescope and large ground-based facilities.

Helium is the second-most abundant element in the Universe after hydrogen and is one of the main constituents of gas-giant planets in our Solar System. Early theoretical models predicted helium to be among the most readily detectable species in the atmospheres of exoplanets, especially in extended and escaping atmospheres. Searches for helium, however, have hitherto been unsuccessful. Here we report observations of helium on an exoplanet, at a confidence level of 4.5 standard deviations. We measured the near- infrared transmission spectrum of the warm gas giant WASP-107b and identified the narrow absorption feature of excited metastable helium at 10,833 angstroms. The amplitude of the feature, in transit depth, is 0.049 +/- 0.011 per cent in a bandpass of 98 angstroms, which is more than five times greater than what could be caused by nominal stellar chromospheric activity. This large absorption signal suggests that WASP-107b has an extended atmosphere that is eroding at a total rate of 10^10 to 3 x 10^11 grams per second (0.1-4 per cent of its total mass per billion years), and may have a comet-like tail of gas shaped by radiation pressure.

Transmission spectroscopy during planetary transits is expected to be a major source of information on the atmospheres of small (approximately Earth-sized) exoplanets in the next two decades. This technique, however, is intrinsically affected by stellar spectral contamination caused by the fact that stellar photo- and chromospheres are not perfectly homogeneous. Such stellar contamination will often reach or exceed the signal introduced by the planetary spectral features. Finding effective methods to correct stellar contamination -- or at least to quantify its possible range -- for the most important exoplanets is a necessary step for our understanding of exoplanet atmospheres. This will require significantly deepening our understanding of stellar heterogeneity, which is currently limited by the available data.

The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2.

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 $\mu$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 stellar contamination can be more than $10 \times$ larger than 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^{+18}_{-7}\%$ and faculae covering fractions $f_{fac} = 54^{+16}_{-46}\%$. The associated stellar contamination signals alter transit depths of the TRAPPIST-1 planets at wavelengths of interest for planetary atmospheric species by roughly 1-15 $\times$ the strength of planetary features, significantly complicating $JWST$ follow-up observations of this system. Similarly, we find stellar contamination can lead to underestimates of bulk densities of the TRAPPIST-1 planets of $\Delta(\rho) = -3^{+3}_{-8} \%$, thus leading to overestimates of their volatile contents.

GJ 1214b is the most studied sub-Neptune exoplanet to date. Recent measurements have shown its near-infrared transmission spectrum to be flat, pointing to a high-altitude opacity source in the exoplanet's atmosphere, either equilibrium condensate clouds or photochemical hazes. Many photometric observations have been reported in the optical by different groups, though simultaneous measurements spanning the entire optical regime are lacking. We present an optical transmission spectrum ($4,500-9,260$\AA) of GJ 1214b in 14 bins measured with Magellan/IMACS repeatedly over three transits. We measure a mean planet-to-star radius ratio of ${R_{p}/R_{s} = 0.1146\pm{2\times10^{-4}}}$ and mean uncertainty of $\sigma(R_{p}/R_{s})=8.7\times10^{-4}$ in the spectral bins. The optical transit depths are shallower on average than observed in the near-infrared. We present a model for jointly incorporating the effects of a composite photosphere and atmospheric transmission (CPAT) through the exoplanet's limb, and use it to examine the cases of absorber and temperature heterogeneities in the stellar photosphere. We find the optical and near-infrared measurements are best explained by the combination of (1) photochemical haze in the exoplanetary atmosphere with a mode particle size $r=0.1~\mu$m and haze-forming efficiency $f_{haze}=10 \%$ and (2) faculae in the unocculted stellar disk with a temperature contrast $\Delta T=354^{+46}_{-46}$ K, assuming 3.2% surface coverage. The CPAT model can be used to assess potential contributions of heterogeneous stellar photospheres to observations of exoplanet transmission spectra, which will be important for searches for spectral features in the optical.