<figure><img src="image_1.png" /></figure>

HAL
open science

# The high-energy environment and atmospheric escape of the mini-Neptune K2-18 b

Leonardo A dos Santos, David Ehrenreich, Vincent Bourrier, Nicola

Astudillo-Defru, Xavier Bonfils, François Forget, Christophe Lovis, Francesco

Pepe, Stéphane Udry

► To cite this version:

Leonardo A dos Santos, David Ehrenreich, Vincent Bourrier, Nicola Astudillo-Defru, Xavier Bonfils, et al.. The high-energy environment and atmospheric escape of the mini-Neptune K2-18 b. Astronomy and Astrophysics - A&A, 2020, 634, pp.L4. 10.1051/0004-6361/201937327. hal-03104121

HAL Id: hal-03104121

https://hal.science/hal-03104121

Submitted on 8 Jan 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
LETTER TO THE EDITOR

# The high-energy environment and atmospheric escape of the mini-Neptune K2-18 b*

Leonardo A. dos Santos<fnref n="1" />, David Ehrenreich<fnref n="1" />, Vincent Bourrier<fnref n="1" />, Nicola Astudillo-Defru<fnref n="2" />, Xavier Bonfils<fnref n="3" />, François Forget<fnref n="4" />, Christophe Lovis<fnref n="1" />, Francesco Pepe<fnref n="1" />, and Stéphane Udry<fnref n="1" />

Received 16 December 2019 / Accepted 13 January 2020

## ABSTRACT

K2-18 b is a transiting mini-Neptune that orbits a nearby (38 pc), cool M3 dwarf and is located inside its region of temperate irradiation. We report on the search for hydrogen escape from the atmosphere K2-18 b using Lyman-α transit spectroscopy with the Space Telescope Imaging Spectrograph instrument installed on the Hubble Space Telescope. We analyzed the time-series of fluxes of the stellar Lyman-α emission of K2-18 in both its blue- and redshifted wings. We found that the average blueshifted emission of K2-18 decreases by 67% ± 18% during the transit of the planet compared to the pre-transit emission, tentatively indicating the presence of H atoms escaping vigorously and being blown away by radiation pressure. This interpretation is not definitive because it relies on one partial transit. Based on the reconstructed Lyman-α emission of K2-18, we estimate an EUV irradiation in the range 10¹ – 10² erg s⁻<fnref n="1" /> cm⁻<fnref n="2" /> and a total escape rate on the order of 10⁸ g s⁻<fnref n="1" />. The inferred escape rate suggests that the planet will lose only a small fraction (< 1%) of its mass and retain its volatile-rich atmosphere during its lifetime. More observations are needed to rule out stellar variability effects, confirm the in-transit absorption, and better assess the atmospheric escape and high-energy environment of K2-18 b.

**Key words.** stars: individual: K2-18 – stars: chromospheres – planets and satellites: atmospheres – ISM: kinematics and dynamics

## 1. Introduction

Short-period exoplanets orbiting nearby, cool M dwarfs are prime targets for the search and characterization of atmospheres of low-mass, sub-Neptune-sized worlds. One particular target that falls in this category is K2-18 b, which was first pointed out as a transiting planet candidate by Montet et al. (2015) and later confirmed with Spitzer photometry (Benneke et al. 2017) and Doppler velocity measurements (Cloutier et al. 2017). This planet has a radius of $R_P = 2.711 \pm 0.065 R_⊕$, a mass of $M_P = 8.64 \pm 1.35 M_⊕$ and an orbital period of $T_{orb} = 32.9 d$ (Cloutier et al. 2019). At an average distance of 0.14 au from its host star, K2-18 b receives a similar amount of bolometric irradiation to that received by the Earth from the Sun. Nevertheless, its high-energy environment is unconstrained, and its density is consistent with either a significant H₂/He envelope or a 100% H₂O composition (Sarkis et al. 2018; Cloutier et al. 2019). The host star is a nearby M2.8-type dwarf located at 38 pc (Gaia Collaboration 2018), rendering K2-18 b one of the best mini-Neptunes suitable for atmospheric follow-up using the

Hubble Space Telescope (HST), the James Webb Space Telescope (JWST), and high-resolution infrared spectrographs.

The 0.4–5.0 µm transmission spectrum of K2-18 b measured with data from the *Kepler* satellite (K2 mission), the Wide-Field Camera 3 (WFC3/HST), and the Infrared Array Camera (IRAC/Spitzer) revealed the presence of water vapor in its lower atmosphere (Benneke et al. 2019; Tsiaras et al. 2019). By comparing atmospheric models with the data, Benneke et al. (2019) concludes that the best match is a H₂-dominated atmosphere with water vapor absorbing above the cloud deck below the 10–1000 mbar pressure level. While these observations provided us with some initial information regarding the composition of its lower atmosphere, they do not constrain the abundances of molecular species.

Models predict that the deposition of high-energy photons (X-rays and far-ultraviolet-FUV) produced by the host star leads to an expansion of the planetary upper atmosphere, as well as the production of H atoms due to photodissociation of H₂O (e.g., Ip 1983; Wu & Chen 1993). This expansion populates the outer layers of the planetary atmosphere where the gas is collisionless, also known as the exosphere. It is therefore likely that the atmosphere of K2-18 b, which is rich in H₂ and H₂O, possesses a H-rich exosphere. Previous HST observations have shown evidence for the presence of large-scale,

<footnote n="1">¹ Observatoire Astronomique de l'Université de Genève, 51 Chemin des Maillettes, 1290 Versoix, Switzerland
e-mail: Leonardo.dosSantos@unige.ch</footnote>

<footnote n="2">² Departamento de Matemática y Física Aplicadas, Universidad Católica de la Santísima Concepción, Alonso de Rivera, 2850 Concepción, Chile</footnote>

<footnote n="3">³ Université Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote n="4">⁴ Laboratoire de Météorologie Dynamique, Institut Pierre Simon Laplace, Université Paris 6 Boite Postale 99, 75252 Paris Cedex 05, France</footnote>

<footnote>* The HST Lyman-α spectra are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/634/L4</footnote>
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** HST/STIS spectra of K2-18 in Visits A (**left panel**) and B (**right panel**) obtained after data reduction (lines). The geocoronal and background contamination is shown as vertical bars; in Visit A, the geocoronal emission overwhelms the stellar fluxes, preventing us from reliably measuring the latter. The Doppler velocities are in the stellar rest frame.</figcaption></figure>

H-rich exospheres around the warm Neptunes GJ 436 b (Ehrenreich et al. 2015; Lavie et al. 2017; dos Santos et al. 2019) and GJ 3470 b (Bourrier et al. 2018). However, to date, no evidence for extended atmospheres has been found for planets smaller than Neptune, as non-detections were reported for the super-Earths 55 Cnc e (Ehrenreich et al. 2012), HD 97658 b (Bourrier et al. 2017a), GJ 1132 b (Waalkes et al. 2019), and π Men c (García Muñoz et al. 2020), and marginal detections were reported for the small rocky-planet systems in TRAPPIST-1 (Bourrier et al. 2017b) and Kepler-444 (Bourrier et al. 2017c).

In this Letter, we report the results of a series of FUV observations of two transits of K2-18 b using the Space Telescope Imaging Spectrograph (STIS) installed on HST. The aim of these observations was to perform Lyman-α transmission spectroscopy of K2-18 b in order to detect its H-rich exosphere, which should produce an excess absorption in the stellar Lyman-α emission during the transit of the planet.

## 2. Observations and data reduction

K2-18 was observed with HST/STIS and the grating G140M (resolving power $R \approx 10\,000$) during two transits on June 18 2017 and March 9 2018 (Program GO-14221, PI: D. Ehrenreich). The first visit (A) allowed five exposures, of which only the first two were successful but were severely contaminated by strong geocoronal emission; this contamination completely swamped the stellar Lyman-α emission and rendered this visit unsuitable for analysis (see Fig. 1). The second visit (B) contained four successful exposures, of which two were performed before the transit ingress and the other two in-transit.

The data were reduced using the standard STIS pipeline, except for the spectral extraction. Since the star is faint, the automated extraction is unable to accurately find the stellar spectrum in the flat-fielded frames. Furthermore, the dark current background of the FUV-MAMA detector of STIS reaches levels high enough to be comparable with the stellar spectrum. In order to correctly extract the spectrum and remove the dark current background, we use the x1d method of `stistools`<fnref n="1" /> with user-defined values for the location of: (i) the spectrum in the cross-dispersion direction, and (ii) the regions where the dark current background near the spectrum can be accurately estimated.

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** HST/STIS spectra of K2-18 in Visit B. The Doppler velocities are in the stellar rest frame. The shaded interval is the region with geocoronal contamination.</figcaption></figure>

Using visual inspection, we determined the location of the spectrum in the cross-dispersion direction to be $y = 389$ px (parameter `a2center` in the pipeline). Determining the best location of the background is not as straightforward; normally, the pipeline extracts the background from regions far from the spectrum, but these regions have discrepant levels of dark current compared to the region of the spectrum. Therefore, we chose to use regions immediately near the location of the spectrum to determine the background (parameters `bk1offst` and `bk2offst` in the pipeline), namely at a distance of $\Delta = 20$ px from `a2center`. We combine the two out-of-transit spectra and the two in-transit spectra separately in order to isolate potential signals of an atmospheric signal of K2-18 b, and the resulting spectra are shown in Fig. 2.

## 3. Tentative detection of a H-rich exosphere in K2-18 b

Since the interstellar medium (ISM) absorbs the core of the Lyman-α emission line, we can only observe the attenuated fluxes in the blue and red wings of the stellar line. We integrated the flux densities in wavelength space between Doppler velocities [-160, -50] km s⁻<fnref n="1" /> and [+50, +160] km s⁻<fnref n="1" />, respectively, to produce the Lyman-α light curves in the blue and red wings.

<footnote n="1">¹ Software freely available at https://stistools.readthedocs.io/en/latest/</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Light curves of the blue (left panel) and red (right panel) wings of the Lyman-α emission of K2-18 during the transit of planet b. Time = 0 corresponds to the transit center according to the ephemeris of Benneke et al. (2019). The vertical red lines represent the ingress and egress of K2-18 b. Visit A was affected by strong geocoronal contamination, and therefore the measured stellar fluxes are likely inaccurate.</figcaption></figure>

During Visit B, we observe a steep decrease in the blueshifted Lyman-α fluxes during the transit of K2-18 b (left panel in Fig. 3), reaching almost zero emission near the planetary egress. This decrease in flux is also seen when we compare the combined out-of-transit and in-transit spectra (see Fig. 2). The redshifted Lyman-α emission of the combined in-transit orbits varies by 14%±23%, which is consistent with a stable redshifted flux during the visit; this indicates that the variation in the blue wing is likely to be astrophysical in nature. Similar results are obtained even when we extract the background emission at different positions in the detector.

The blue wing flux measured in the combined in-transit spectra decreases by 67% ± 18% in relation to the combined out-of-transit spectra; in particular, the last orbit displays an absorption of 93% ± 18% in relation to the combined out-of-transit spectra. Although statistically significant, we conservatively deem this result tentative until it is repeated in future observations; for a reference, the intrinsic stellar variability of the Lyman-α emission of HD 97658 b is on the order of a few tens of percent at ~2σ confidence (Bourrier et al. 2017a). If confirmed to be linked to the transit of K2-18 b, the variation in Lyman-α flux can be interpreted as the absorption caused by an extended, H-rich exosphere of the planet. Such a large absorption signal can be explained by a combination of large atmospheric escape rate and long photoionization lifetime of the H atoms in the exosphere. This result gives further support to the hypothesis that K2-18 b possesses a H₂-dominated envelope (as in the conclusions of Benneke et al. 2019), unlike the super-Earths 55 Cnc e and HD 219134 b. Previous results for the super-Earths HD 97658 b and GJ 1132 b were inconclusive due to stellar variability for the first (Bourrier et al. 2017a) and lack of stellar blue wing emission in the second (Waalkes et al. 2019).

## 4. The high-energy environment of K2-18 b

Determining the high-energy environment of K2-18 b provides a critical piece of information to interpret the evolution and current state of its atmosphere. To that end, we used the STIS observations of K2-18 to reconstruct its intrinsic Lyman-α spectrum and estimate the high-energy irradiation received by planet b. The reconstruction process follows the standard method used in, for example, Bourrier et al. (2017b, 2018). In short, we fit the observed spectrum to a model of the intrinsic emission line attenuated by ISM absorption, scaled for distance and convolved with

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Reconstructed intrinsic Lyman-α spectrum of K2-18 (green) and the average Lyman-α profile as observed with HST/STIS (black bars). The red (blue) curve shows the inferred observable spectrum with (without) instrumental convolution and the shaded region represents the 1σ uncertainty for the red curve. The dotted curve shows the inferred ISM absorption profile.</figcaption></figure>

the instrumental response; the fit yields an estimate of the intrinsic emission and certain properties of the ISM in the line of sight. In this process, we assume that the intrinsic Lyman-α emission of K2-18 possesses a Gaussian profile (applicable for M dwarfs and for the quality of the available spectra; see Bourrier et al. 2017c, 2018), and fix the temperature and turbulent velocity of the ISM to 8000 K and 1.23 km s⁻¹ (for the NGP cloud, as estimated by the LISM calculator Redfield & Linsky 2008), respectively. We also set the deuterium-to-hydrogen ratio (D/H) to 1.5 × 10⁻⁵ and the systemic velocity to 0.6537 km s⁻¹ (measured with high-resolution spectra by Cloutier et al. 2017).

The result of the Lyman-α line reconstruction is shown in Fig. 4. We simultaneously fit each exposure to a global Lyman-α line model within ±300 km s⁻¹ in the stellar rest frame, excluding the band that corresponds to significant airglow contamination (approximately the range ±50 km s⁻¹, shown in orange in Fig. 4). When fitting the in-transit orbits we excluded the pixels corresponding to the range potentially absorbed by the planet. We explore the parameter space using the Markov chain Monte Carlo ensemble sampler implementation of emcee (Foreman-Mackey et al. 2013); in addition, we report the uncertainties based on the highest density interval (HDI), which contains 68.3% of the distribution mass such that no point outside the interval has a higher density than any point within it. We
<table><caption>**Table 1.** Properties of the high-energy environment of K2-18 b.</caption><tr><td>Lyman-α flux (erg s<sup>-1</sup> cm<sup>-2</sup>)</td><td>100.7<sup>+96.1</sup><sub>-82.4</sub></td></tr><tr><td>EUV 10–91.2 nm flux (erg s<sup>-1</sup> cm<sup>-2</sup>)</td><td>107.9<sup>+124.7</sup><sub>-90.8</sub></td></tr><tr><td>Photoionization rate (×10<sup>-6</sup> s<sup>-1</sup>)</td><td>3.7<sup>+7.5</sup><sub>-3.5</sub></td></tr><tr><td>Photoionization lifetime (h)</td><td>3600<sup>+49900</sup><sub>-2400</sub></td></tr><tr><td>Escape rate at 100% efficiency (×10<sup>8</sup> g s<sup>-1</sup>)</td><td>3.5<sup>+4.0</sup><sub>-2.9</sub></td></tr></table>

fit for four free parameters in total: the temperature (assuming a Gaussian thermal broadening) and amplitude of the intrinsic stellar Lyman-α line, and the radial velocity and H I density of the ISM. We estimated the EUV flux in K2-18 b using the relation from Linsky et al. (2014).

We determined that the heliocentric radial velocity of the ISM is $V_R = 11.39_{-4.73}^{+4.04}$ km s⁻¹, and the H I column density in the line of sight is $log_{10} \eta = 18.16_{-0.34}^{+0.44}$ cm⁻². While $\eta$ is consistent with the results from Wood et al. (2005), the value we determined for $V_R$ differs by $2\sigma$ with the value predicted by LISM calculator (Redfield & Linsky 2008).

We estimated the properties of the high-energy environment of K2-18 b based on the reconstructed Lyman-α emission and the results are shown in Table 1; the large uncertainties in these estimates are due to the low signal-to-noise ratio (S/N) of the spectra and the uncertainties in the semi-empirical relations used to estimate the EUV flux. The best fit model of the intrinsic Lyman-α emission of K2-18 results in a $\beta$ (ratio between radiation pressure and stellar gravity; see right axis of Fig. 4) of about 2.2. Although more data would be necessary to confirm this line shape, and the corresponding high $\beta$, it nonetheless suggests that radiation pressure could blow away the escaping hydrogen atoms more strongly in K2-18 b than in GJ 436 b (Bourrier et al. 2016) or the TRAPPIST-1 planets (Bourrier et al. 2017b). This value of $\beta$ is similar to that inferred for GJ 3470 b (Bourrier et al. 2018), which could indicate that the exosphere of K2-18 b possesses a similar shape. Using the energy-limited escape (e.g., Salz et al. 2015) as an initial estimate, we predict that the total escape rate in K2-18 b is two orders of magnitude lower when compared to the values inferred for GJ 436 b (~2.2 × 10¹⁰ g s⁻¹; Bourrier et al. 2016) and GJ 3470 b (~8.5 × 10¹⁰ g s⁻¹; Bourrier et al. 2018). We estimated the photoionization rate and lifetime as in Bourrier et al. (2017a), which yields a value likely above 1200 h for the latter; this long lifetime means that the H atoms in the exosphere of K2-18 b could stay neutral for much longer than in GJ 436 (12 h) and GJ 3470 b (3.5 h). These results are not surprising because K2-18 b is subject to lower irradiation levels than the aforementioned warm Neptunes.

The inferred escape rate of K2-18 b is likely underestimated because we do not take into account the stellar X-ray flux, which is currently unknown. For a similar star like GJ 436, the ratio between X-ray and EUV emission is ~0.23 (Bourrier et al. 2016), which provides an upper limit on the level of underestimation contained in our calculation; however, the contribution of X-ray flux may be less important for the ionization conditions because of the smaller ionization cross-section when compared to EUV wavelengths.

## 5. Conclusions

K2-18 b is currently one of the best targets for transit spectroscopy among sub-Neptune planets due to its large scale

height, its short distance from the Sun, and the infrared brightness of the host star. Previous results have shown evidence that the atmosphere of the planet is dominated by H₂/He and contains water vapor. In these atmospheric conditions and under the expected high levels of EUV irradiation, K2-18 b is prone to efficiently losing its atmosphere and producing a detectable excess absorption of H in Lyman-α caused by a H-rich exosphere during transit. In this study we analyzed four HST orbits before and during the transit of K2-18 b with the STIS instrument to search for this feature.

We analyzed the flux time series of both the blue- and redshifted wings of the stellar Lyman-α emission. The blue wing displays a significant excess absorption during the transit; in particular, near the egress of K2-18 b, the flux in the blue wing is consistent with 100% absorption. The in-transit red wing fluxes vary by 14% ± 23% and are significantly more stable than the blue wing fluxes. A blueshifted absorption could indicate the presence of a H-rich exosphere around K2-18 b being swept away by radiation pressure from its host star towards the direction of the observer, similar to the exospheres of GJ 436 b and GJ 3470 b.

Despite the low S/N of the observed spectra, we were able to reconstruct the intrinsic stellar emission (without the ISM absorption) to assess the high-energy environment of K2-18 b. Our first estimate for the expected total escape rate of K2-18 b leads to a value on the order of 10⁸ g s⁻¹. The ratio between radiation pressure and gravity ($\beta$) suggests that the exosphere of K2-18 b is in a similar state to that observed for GJ 3470 b. We estimate that the EUV (10–91.2 nm) flux in the planet is on the order of 10¹⁻¹⁰ erg s⁻¹ cm⁻². At the estimated escape rate, it is likely that the planet will lose only a small fraction (1% or less) of its mass during its remaining lifetime, and therefore it is probably not an archetypal planet crossing the radius valley to become a bare rock (Fulton et al. 2017; Van Eylen et al. 2018; Fulton & Petigura 2018); as such, the planet will likely retain its volatile-rich atmosphere due to the more amenable EUV irradiation flux than for example GJ 3470 b, which is at least ten times more EUV irradiated than K2-18 b.

Since we observed only one partial transit of K2-18 b, we conclude that the H-rich exosphere detection is only tentative for now, and more observations are needed to rule out stellar activity effects and confirm the reported feature. Furthermore, additional observations of the Lyman-α spectrum of K2-18 b will help in better constraining the high-energy environment of the planet and its atmospheric escape history.

**Acknowledgements.** LADs thanks M. Turbet and M. López-Morales for the insightful discussions about K2-18 b. This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (project FOUR ACES; grant agreement No 724427), and it has been carried out in the frame of the National Centre for Competence in Research PlanetS supported by the Swiss National Science Foundation (SNSF). This research is based on observations made with the NASA/ESA Hubble Space Telescope. NA-D acknowledges the support of FONDECYT project 3180063. The authors are grateful to the anonymous referee for the quick and helpful review. The data are openly available in the Mikulski Archive for Space Telescopes (MAST), which is maintained by the Space Telescope Science Institute (STScI). STScI is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 5-26555. This research made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. We used the open source software SciPy (Jones et al. 2001), Jupyter (Kluyver et al. 2016), Astropy (Astropy Collaboration 2013), Matplotlib (Hunter 2007) and emcee (Foreman-Mackey et al. 2013).
References

Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33

Benneke, B., Werner, M., Petigura, E., et al. 2017, ApJ, 834, 187

Benneke, B., Wong, I., Piaulet, C., et al. 2019, ApJ, 887, L14

Bourrier, V., Lecavelier des Etangs, A., Ehrenreich, D., Tanaka, Y. A., & Vidotto, A. A. 2016, A&A, 591, A121

Bourrier, V., Ehrenreich, D., King, G., et al. 2017a, A&A, 597, A26

Bourrier, V., de Wit, J., Bolmont, E., et al. 2017b, AJ, 154, 121

Bourrier, V., Ehrenreich, D., Allart, R., et al. 2017c, A&A, 602, A106

Bourrier, V., Lecavelier des Etangs, A., Ehrenreich, D., et al. 2018, A&A, 620, A147

Cloutier, R., Astudillo-Defru, N., Doyon, R., et al. 2017, A&A, 608, A35

Cloutier, R., Astudillo-Defru, N., Doyon, R., et al. 2019, A&A, 621, A49

dos Santos, L. A., Ehrenreich, D., Bourrier, V., et al. 2019, A&A, 629, A47

Ehrenreich, D., Bourrier, V., Bonfils, X., et al. 2012, A&A, 547, A18

Ehrenreich, D., Bourrier, V., Wheatley, P. J., et al. 2015, Nature, 522, 459

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306

Fulton, B. J., & Petigura, E. A. 2018, AJ, 156, 264

Fulton, B. J., Petigura, E. A., Howard, A. W., et al. 2017, AJ, 154, 109

Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1

García Muñoz, A., Youngblood, A., Fossati, L., et al. 2020, ApJ, 888, L21

Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90

Ip, W. H. 1983, ApJ, 264, 726

Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy: Open Source Scientific Tools for Python

Kluyver, T., Ragan-Kelley, B., Pérez, F., et al. 2016, in Positioning and Power in Academic Publishing: Players, Agents and Agendas, eds. F. Loizides, & B. Schmidt (IOS Press), 87

Lavie, B., Ehrenreich, D., Bourrier, V., et al. 2017, A&A, 605, L7

Linsky, J. L., Fontenla, J., & France, K. 2014, ApJ, 780, 61

Montet, B. T., Morton, T. D., Foreman-Mackey, D., et al. 2015, ApJ, 809, 25

Redfield, S., & Linsky, J. L. 2008, ApJ, 673, 283

Salz, M., Schneider, P. C., Czesla, S., & Schmitt, J. H. M. M. 2015, A&A, 576, A42

Sarkis, P., Henning, T., Kürster, M., et al. 2018, AJ, 155, 257

Tsiaras, A., Waldmann, I. P., Tinetti, G., Tennyson, J., & Yurchenko, S. N. 2019, Nat. Astron., 3, 1086

Van Eylen, V., Agentoft, C., Lundkvist, M. S., et al. 2018, MNRAS, 479, 4786

Waalkes, W. C., Berta-Thompson, Z., Bourrier, V., et al. 2019, AJ, 158, 50

Wood, B. E., Redfield, S., Linsky, J. L., Müller, H.-R., & Zank, G. P. 2005, ApJS, 159, 118

Wu, C. Y. R., & Chen, F. Z. 1993, J. Geophys. Res., 98, 7415