# Stable accretion and episodic outflows in the young transition disk system GM Aurigae

J. Bouvier, A. Sousa, K. Pouilly, J. Almenara, J.-F. Donati, S. Alencar, A.

Frasca, K. Grankin, A. Carmona, G. Pantolmos, et al.

► To cite this version:

J. Bouvier, A. Sousa, K. Pouilly, J. Almenara, J.-F. Donati, et al. Stable accretion and episodic outflows in the young transition disk system GM Aurigae. Astronomy & Astrophysics - A&A, 2023, 672, pp.A5. 10.1051/0004-6361/202245342. hal-04252213

HAL Id: hal-04252213

https://hal.science/hal-04252213v1

Submitted on 26 Oct 2023

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.

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

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Stable accretion and episodic outflows in the young transition disk system GM Aurigae

## A semester-long optical and near-infrared spectrophotometric monitoring campaign*,***

J. Bouvier<fnref n="1" />, A. Sousa<fnref n="1" />, K. Pouilly<fnref n="2" />, J. M. Almenara<fnref n="1" />, J.-F. Donati<fnref n="3" />, S. H. P. Alencar<fnref n="4" />, A. Frasca<fnref n="5" />, K. Grankin<fnref n="6" />, A. Carmona<fnref n="1" />, G. Pantolmos<fnref n="1" />, B. Zaire<fnref n="4" />, X. Bonfils<fnref n="1" />, A. Bayo<fnref n="7" />,<fnref n="8" />, L. M. Rebull<fnref n="9" />, J. Alonso-Santiago<fnref n="5" />, J. F. Gameiro<fnref n="10" />,<fnref n="11" />, N. J. Cook<fnref n="12" />, E. Artigau<fnref n="12" />, and the Spirou Legacy Survey (SLS) consortium

Received 31 October 2022 / Accepted 3 January 2023

## ABSTRACT

**Context.** Young stellar systems actively accrete from their circumstellar disk and simultaneously launch outflows. The physical link between accretion and ejection processes remains to be fully understood.

**Aims.** We investigate the structure and dynamics of magnetospheric accretion and associated outflows on a scale smaller than 0.1 au around the young transitional disk system GM Aur.

**Methods.** We devised a coordinated observing campaign to monitor the variability of the system on timescales ranging from days to months, including partly simultaneous high-resolution optical and near-infrared spectroscopy, multiwavelength photometry, and low-resolution near-infrared spectroscopy, over a total duration of six months, covering 30 rotational cycles. We analyzed the photometric and line profile variability to characterize the accretion and ejection processes.

**Results.** The optical and near-infrared light curves indicate that the luminosity of the system is modulated by surface spots at the stellar rotation period of 6.04 ± 0.15 days. Part of the Balmer, Paschen, and Brackett hydrogen line profiles as well as the HeI 5876 Å and HeI 10830 Å line profiles are modulated on the same period. The Paβ line flux correlates with the photometric excess in the *u'* band, which suggests that most of the line emission originates from the accretion process. High-velocity redshifted absorptions reaching below the continuum periodically appear in the near-infrared line profiles at the rotational phase in which the veiling and line fluxes are the largest. These are signatures of a stable accretion funnel flow and associated accretion shock at the stellar surface. This large-scale magnetospheric accretion structure appears fairly stable over at least 15 and possibly up to 30 rotational periods. In contrast, outflow signatures randomly appear as blueshifted absorption components in the Balmer and HeI 10830 Å line profiles. They are not rotationally modulated and disappear on a timescale of a few days. The coexistence of a stable, large-scale accretion pattern and episodic outflows supports magnetospheric ejections as the main process occurring at the star-disk interface.

**Conclusions.** Long-term monitoring of the variability of the GM Aur transitional disk system provides clues to the accretion and ejection structure and dynamics close to the star. Stable magnetospheric accretion and episodic outflows appear to be physically linked on a scale of a few stellar radii in this system.

**Key words.** accretion, accretion disks – stars: magnetic field – stars: individual: GM Aur – stars: pre-main sequence – stars: variables: T Tauri, Herbig Ae/Be – stars: winds, outflows

<footnote n="1">¹ Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France
e-mail: jerome.bouvier@univ-grenoble-alpes.fr</footnote>

<footnote n="2">² Department of Physics and Astronomy, Uppsala University, Box 516, 75120, Sweden</footnote>

<footnote n="3">³ Univ. de Toulouse, CNRS, IRAP, 14 avenue Belin, 31400 Toulouse, France</footnote>

<footnote n="4">⁴ Departamento de Física – ICEx – UFMG, Av. Antonio Carlos 6627, 30270-901 Belo Horizonte, MG, Brazil</footnote>

<footnote n="5">⁵ INAF – Osservatorio Astrofisico di Catania, via S. Sofia 78, 95123 Catania, Italy</footnote>

<footnote n="6">⁶ Crimean Astrophysical Observatory, Nauchny 298409, Crimea</footnote>

<footnote n="7">⁷ Instituto de Física y Astronomía, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile</footnote>

<footnote n="8">⁸ European Southern Observatory, Karl-Schwarzschild-Strasse 2, 85748 Garching bei München, Germany</footnote>

<footnote n="9">⁹ Infrared Science Archive (IRSA), IPAC, 1200 E. California Blvd., California Institute of Technology, Pasadena, CA 91125, USA</footnote>

<footnote n="10">¹⁰ Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal</footnote>

<footnote n="11">¹¹ Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto, rua do Campo Alegre 687, 4169-007 Porto, Portugal</footnote>

<footnote n="12">¹² Research on Exoplanets, Université de Montréal, Département de Physique, Montréal, QC H3C 3J7, Canada</footnote>

<footnote>* Tables containing the *u'*g'*r'*i'*' LCOGT photometric measurements are only available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/672/A5</footnote>

<footnote>** Based on observations obtained at the Canada-France-Hawaii Telescope (CFHT), at the Observatoire de Haute-Provence (OHP), at the European Organisation for Astronomical Research in the Southern Hemisphere (ESO), and at the Las Cumbres Observatory global telescope network (LCOGT).</footnote>
# 1. Introduction

Accretion and ejection processes are at the origin of most of the peculiar properties of young stellar systems. The structure and dynamics of the accretion flows within the disk and from the inner disk to the star, as well as the properties of the multiple outflows arising from the disk, from the star-disk interface, and from the stellar surface, remain to be fully deciphered, however. Low-mass pre-main-sequence stars, the so-called T Tauri stars (TTS), accrete from their circumstellar disks for a few million years, while contemporaneous planet formation impacts the disk structure and evolution. In the inner regions of the system, the disk is disrupted by the strong stellar magnetosphere that channels the accretion flow toward the star along magnetic field lines (see, e.g., the review by Hartmann et al. 2016). Thus, accretion funnel flows develop that connect the inner disk to the stellar surface, where the material is accreted at nearly free-fall velocity and is eventually halted in a strong accretion shock. Simultaneously, outflows are produced at the star-disk interface close to the magnetospheric truncation radius through the inflation and reconnection of magnetic field lines that are twisted by differential rotation (e.g., Zanni & Ferreira 2013). Ultimately, the release of gravitational energy delivered by the accretion process may trigger accretion-powered stellar winds (Matt & Pudritz 2005). The torque balance between accretion and ejection processes is a central issue for understanding the spin evolution of young stars (e.g., Pantolmos et al. 2020; Ireland et al. 2021).

The star-disk interaction takes place on a distance of a few stellar radii (e.g., Bessolaz et al. 2008), that is, on a scale of about 0.1 au or smaller. MHD models developed by several groups predict the structure and dynamics of the magnetospheric accretion region and associated outflows (see, e.g., the review by Romanova & Owocki 2015). Observationally, two main directions have been explored so far to investigate the properties of this region. On one hand, monitoring the spectroscopic and photometric variability of the system over a few rotational periods, that is, typically over a few weeks, allows identifying the signature of funnel flows, hot spots, and outflows, and relating them to the strength and topology of the surface magnetic field that is measured from spectropolarimetry (e.g., Pouilly et al. 2020, 2021; Bouvier et al. 2020a; Donati et al. 2019, 2020a; Alencar et al. 2018). On the other hand, a direct approach attempts to spatially resolve the star-disk interaction region on a scale of a few milliarcsecond on the sky, using long-baseline near-infrared interferometry (e.g., Eisner et al. 2014; Gravity Collaboration 2020; Bouvier et al. 2020b). Both approaches have been successful in mapping the inner region of accreting systems and have provided strong support to the magnetospheric accretion scenario and its MHD modeling. Following previous studies, we report here the results of a new observing campaign devoted to the young stellar system GM Aur.

GM Aur (RA = 04h55, Dec = +30°21, V = 12.1 mag) is a solar-type pre-main-sequence star located in the Taurus Auriga molecular cloud at a distance of 157.9 ± 1.2 pc (Gaia Collaboration 2021). This classical T Tauri star (cTTS) has a spectral type K6 (Herczeg & Hillenbrand 2014) and is surrounded by a circumstellar disk from which it actively accretes material at a rate of 0.6–2.0 × 10⁻⁸ M☉ yr⁻¹ (Robinson & Espaillat 2019). Based on its spectral energy distribution, which exhibits a small near-infrared excess compared to a significant mid-infrared one, the system has long been suspected to be in a transitional stage, that is, that it is surrounded by a disk whose inner regions are relatively devoid of matter

(Strom et al. 1989). High-resolution ALMA images of the circumstellar disk indeed reveal that it is highly structured. The large-scale disk, inclined at ~53° on the line of sight, features a large inner dust cavity extending over ~35–40 au and a succession of annular gaps and dusty rings on a wider scale up to 200 au (Macías et al. 2018; Huang et al. 2020). Much closer to the central star, long-baseline VLTI/GRAVITY interferometric observations unveil a compact dusty disk, whose inner edge was recently reported to be located at $r_{in} = 0.013^{+0.015}_{-0.008}$ au from the central star (Bohn et al. 2022) and that extends over at least a few 0.1 au (Akeson et al. 2005) and possibly up to 6.6 au (Varga et al. 2018; Woitke et al. 2019). The gaseous component of the inner disk has been detected from CO 4.7 micron emission down to 0.5 ± 0.2 au (Salyk et al. 2009). The inclination and position angle of the major axis of the inner dusty disk ($i = 68°^{+16}_{-28}$, PA = 37°$^{+31}_{-22}$) are found to be consistent with those of the outer disk, which suggests that the inner and outer disks are aligned (Bohn et al. 2022).

In an attempt to decipher the physical processes at work at the heart of the system, GM Aur has been the subject of several multiwavelength monitoring campaigns. The long-term light curve presented by Grankin et al. (2007) over the period 1986–1995 exhibits relatively low-level variability, with a V-band magnitude ranging from 11.74 to 12.35 mag. Photometric variations are modulated by surface spots at the stellar rotation period of 6.0–6.1 days (Percy et al. 2010; Artemenko et al. 2012). Ingleby et al. (2015) reported variability over the full wavelength range from the far-UV to the near-infrared, which they attributed in part to an accretion rate that varies by about a factor of 2–3 on a timescale of months, and for another part to dust inhomogeneities that are located in the inner disk close to the truncation radius. Variations in the mass accretion rate of similar amplitude have also been reported on a shorter timescale of about a week by Robinson & Espaillat (2019), and a connection between mass loss and mass accretion has been further suggested by Espaillat et al. (2019). McGinnis et al. (2020) presented the results of a high-resolution optical spectroscopic monitoring campaign performed on a timescale of a week that illustrated the variability of the $H\alpha$, $H\beta$, and HeI emission line profiles of the system. From the measured radial velocity variations of the HeI 5876 Å line profile, whose narrow component traces the accretion shock, they deduced that GM Aur accretes material from its circumstellar disk through an inclined magnetosphere, whose axis is tilted by about 13° relative to the stellar rotational axis. GM Aur indeed harbors a strong surface magnetic field, with a mean value of 2.2 kG (Johns-Krull 2007; Symington et al. 2005). Finally, from a recent multiwavelength X-UV-optical campaign, Espaillat et al. (2021) reported evidence for a transverse density stratification within the accretion shock at the base of the magnetic funnel flow.

We report here the results of a new coordinated monitoring campaign on GM Aur that combines high-resolution optical spectroscopy and near-infrared spectropolarimetry, multiwavelength optical and near-infrared photometry, and long-term low-resolution near-infrared spectroscopy. Part of the observations have been obtained simultaneously over a timescale of a few weeks, while the total duration of the campaign amounted to six months. The goal of the campaign was to investigate the physical processes that cause variability in GM Aur on a scale of a few stellar radii, and in particular, to constrain the structure and dynamics of the magnetospheric accretion flow from the inner disk to the star. We devised a long-term campaign in order to be able to probe various timescales, from days to months, and obtain a sufficiently long temporal baseline to investigate the
relation between accretion and ejection processes on small spa-
tial scales from the stellar surface to the inner disk regions.

The campaign whose results are reported here took place in
the framework of a larger project led by the ODYSSEUS team<fnref n="1" />
(see Espaillat et al. 2022), which uses the Hubble UV Legacy
Library of Young Stars as Essential Standards program (ULL-
YSES<fnref n="2" />, Roman-Duval et al. 2020), on HST Director’s Discre-
tionary time, to monitor a sample of T Tauri stars in the UV
domain, which includes GM Aur. Additional follow-up observa-
tions were acquired for this project at ESO in the framework of
the PENELLOPE Large Program<fnref n="3" /> (Manara et al. 2021).

In Sect. 2 we describe the observational techniques we
implemented to perform the campaign. In Sect. 3 we derive
the properties of the system and analyze its photometric and
spectroscopic variability over timescales from days to months,
including veiling measurements and emission line profiles. We
infer the global structure of the magnetospheric accretion flow
from the observed variability and characterize associated out-
flows. In Sect. 4 we discuss the dynamics of the accretion and
ejection structure and show that short-lived episodic outflows
coexist with a stable magnetospheric accretion pattern. In Sect. 5
we conclude on the ability of multiwavelength, multi-technique
coordinated observational campaigns to unveil the physical pro-
cesses at work in young stellar systems at the sub-au scale.

## 2. Observations

In this section, we describe the acquisition and data-reduction
processes of photometric, spectroscopic, and spectropolarimet-
ric datasets obtained during the large-scale campaign we per-
formed on the cTTS GM Aur from September 6, 2021, to
March 8, 2022, using CFHT/SPIRou, OHP/SOPHIE, ESO/
ExTrA, LCOGT, and ESO/REM. A summary plot of the GM
Aur observing campaign reported here is provided in Fig. 1.

2.1. LCOGT: Multiwavelength optical photometry

GM Aur was observed at Las Cumbres Observatory Global
Network (LCOGT, Brown et al. 2013) from September 6 to
December 30, 2021. We acquired 850 images in the Sloan
u'g'r'i' filters over two runs with a sub-day cadence
(LCO2021B-001, PI L. Rebull; CLN2021B-003, PI A. Bayo).
The u' images were obtained with the Sinistro 1m telescopes of
the network using an exposure time of 180 seconds and reading
the 2K × 2K central window of the detector with a 2 × 2 binning,
resulting in a 13 × 13 arcmin field of view on the sky. The g'r'i'
images were obtained with the 0.4m SBIG telescopes, offering a
field of view of 29.2 × 19.5 arcmin, with exposure times of 60,
20, and 20 seconds, respectively. We retrieved the BANZAI-
reduced images from the LCOGT archive service and the non-
calibrated photometric catalogs provided in the image headers
for all detected stars in the field.

In order to compute differential photometry, we considered
two stars, HD 282625 and HD 282626, both located within
3 arcmin of GM Aur. The first star was used as a reference
star to calibrate the differential light curve, and the second was
used as a check star to assess that these are nonvariable sources.
These field stars have spectral types F2 and F5, respectively,
and are only slightly brighter than GM Aur. We confirmed that
the two stars are nonvariable from their differential light curves,

and we deduced a mean rms photometric error of 0.025 mag in
the *u'*g'*r'* filters and 0.033 mag in the *i'* filter. We proceeded to
compute the differential light curve between GM Aur and HD
282625 in all four filters. We adopted the mean magnitude of HD
282625 listed in the APASS and Pan-STARRS surveys, namely
g' = 11.331 mag, r' = 10.916 mag, and i' = 10.756 mag to cal-
ibrate the GM Aur light curve in the g'*r'*i' filters to within an
accuracy of 0.02 mag.

We were not able to find an estimate of the *u'*-band mag-
nitude for the comparison stars in the literature. Instead, we
assumed the intrinsic (*u'* − *g'* ) colors of an F2 star (Covey et al.
2007; Kraus & Hillenbrand 2007) for HD 282625, to which we
applied interstellar reddening. From the observed versus intrin-
sic color indices of the HD 282625 (*g'* − *r'* ) and (*r'* − *i'* ) bands,
we derived *A*<sub>V</sub> = 0.8 ± 0.1fmag, using the *R* = 3.1 reddening
law from Fiorucci & Munari (2003). This procedure yielded an
estimate of the reddened (*u'* − *g'* ) color of the comparison star
from which we derived its *u'*-band magnitude. The photometric
calibration of the GM Aur *u'*-band light curve is thus relatively
indirect and probably not accurate to better than 0.1 mag<sup>4</sup>.

## 2.2. REM: Optical and near infrared photometry

Observations were performed with the 60 cm robotic REM tele-
scope located at the ESO La Silla Observatory (Chile), on
15 nights from JD 2 459 497 to JD 2 459 520 (October 9 to
November 1, 2021). By means of a dichroic, REM simultane-
ously feeds two cameras at the two Nasmyth focal stations, one
camera for the near-infrared (REMIR), and the other for the opti-
cal (ROSS2). The cameras have nearly the same field of view
of about 10' × 10' and use wide-band filters (J, H, and K' for
REMIR and Sloan/SDSS g', r', i', and z' for ROSS2).

Exposure times were 60 s for ROSS2, which simultaneously acquires images in the four Sloan bands, and five ditherings of 3 s each were adopted for each filter of REMIR. For the ROSS2 camera, we generated master flats using the twilight flat-fields taken during the observing run, which are available in the REM archive. The latter were used to correct for pixel-to-pixel sensitivity variations, as well as for the vignetting and illumination of the field of view. After subtracting the dark-frame, each scientific image was divided by the proper master-flat, depending on the filter. The prereduction of the REMIR images is automatically done by the AQuA pipeline (Testa et al. 2004), and the coadded and sky-subtracted frames, resulting from five individual ditherings, are made available to the observer.

The adopted comparison stars are reported in Table A.1
along with their *griz* (Tonry et al. 2018) and *JHK* (Cutri et al.
2003) magnitudes. Aperture photometry for all the stars listed
in Table A.1 was performed with DAOPHOT by using the IDL<fnref n="5" />
routine APER. For each frame and filter, we used the instrumen-
tal magnitudes of the stars listed in Table A.1 to generate an
artificial comparison, weighting them with the flux correspond-
ing to their standard magnitude in a way similar to the ensemble
photometry. This procedure also allowed us to evaluate a stan-
dard error based on the differences between the magnitudes cal-
culated with different comparison stars.

The optical photometry gathered at REM is listed in
Table A.2. In the common *g'*r'*i'* bands, we found it to agree well
with that obtained at LCOGT with a tighter sampling rate, and
therefore, we did not use it further in the analysis below. The

<footnote n="1">¹ https://sites.bu.edu/odysseus/</footnote>

<footnote n="2">² https://ullyses.stsci.edu/</footnote>

<footnote n="3">³ https://sites.google.com/view/cfmanara/penellope</footnote>

<footnote n="4">⁴ The table of photometric measurements is available at CDS.</footnote>

<footnote n="5">⁵ Interactive Data Language (IDL) is a registered trademark of Harris Corporation.</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Temporal sampling of the GM Aur campaign from September 6, 2021, to March 8, 2022. Bottom: *g'*-band light curve from LCOGT (black dots) and from REM (magenta crosses). The REM *g'*-band magnitudes are offset in this figure by +0.1 mag to match the LCOGT measurements. The mean photometric error on both the LCOGT and REM *g'*-band measurements is 0.025 mag. Vertical arrows: the vertical arrows show the dates of OHP/SOPHIE (red), CFHT/SPIRou (blue), and ESO/ExTrA (black) observations. The core of the campaign took place during October 2021 with contemporaneous measurements from the five instruments.</figcaption></figure>

individual *JHK'* measurements are listed in Table 1. The median errors on *JHK'* measurements are 0.05, 0.05, and 0.06 mag, respectively. However, some measurements were affected by the nearby bright moon around JD 2459512, and had an error larger than 0.1 mag. We chose to discard these measurements. We eventually derived the following values for the median near-infrared magnitudes of the GM Aur system and their rms variations at the time of the observations: *J* = 9.41 ± 0.10 mag, *H* = 8.71 ± 0.03 mag, and *K* = 8.40 ± 0.07 mag.

## 2.3. OHP SOPHIE: High-resolution optical spectroscopy

Observations were carried out from October 12 to 29, 2021, at Observatoire de Haute-Provence using the fiber-fed SOPHIE spectrograph (Perruchot et al. 2008) in high-efficiency mode, which delivers a spectral resolution of *R* ~ 40 000 over the wavelength range 387–694 nm. We obtained 15 spectra over 18 nights, with an exposure time of 3600 s, yielding a signal-to-noise ratio ranging from 42 to 67 at 600 nm. The raw spectra were fully reduced at the telescope by the SOPHIE real-time pipeline (Bouchy et al. 2009). The data products include a resampled 1D spectrum with a constant wavelength step of 0.01 Å, corrected for barycentric radial velocity, an order-by-order estimate of the signal-to-noise ratio, and a measurement of the source radial velocity, $V_r$, and projected rotational velocity, $v \sin i$. The latter two quantities are derived from a

<table><caption>**Table 1.** REM *JHK* photometry.</caption><thead><tr><th>Julian date<br>(2450 000+)</th><th>J<br>(mag)</th><th>Julian date<br>(2450 000+)</th><th>H<br>(mag)</th><th>Julian date<br>(2450 000+)</th><th>K<br>(mag)</th></tr></thead><tbody><tr><td>9497.73017</td><td>9.4</td><td>9497.73222</td><td>8.73</td><td>9497.73431</td><td>8.43</td></tr><tr><td>9497.73084</td><td>9.41</td><td>9497.73289</td><td>8.69</td><td>9497.73497</td><td>8.39</td></tr><tr><td>9500.80130</td><td>9.39</td><td>9500.80335</td><td>8.62</td><td>9500.80543</td><td>8.37</td></tr><tr><td>9500.80197</td><td>9.36</td><td>9500.80402</td><td>8.67</td><td>9500.80609</td><td>8.36</td></tr><tr><td>9501.80968</td><td>9.41</td><td>9501.81171</td><td>8.67</td><td>9501.81379</td><td>8.39</td></tr><tr><td>9501.81034</td><td>9.38</td><td>9501.81238</td><td>8.69</td><td>9501.81445</td><td>8.4</td></tr><tr><td>9502.83416</td><td>9.57</td><td>9502.83554</td><td>8.69</td><td>9502.83764</td><td>8.56</td></tr><tr><td>9503.83843</td><td>9.66</td><td>9502.83621</td><td>8.71</td><td>9502.83830</td><td>8.49</td></tr><tr><td>9504.84267</td><td>9.72</td><td>9503.83980</td><td>8.72</td><td>9507.79811</td><td>8.34</td></tr><tr><td>9507.79393</td><td>9.43</td><td>9503.84047</td><td>8.72</td><td>9507.79877</td><td>8.35</td></tr><tr><td>9507.79459</td><td>9.4</td><td>9504.84406</td><td>8.75</td><td>9509.81133</td><td>8.23</td></tr><tr><td>9509.80714</td><td>9.28</td><td>9504.84473</td><td>8.74</td><td>9509.81198</td><td>8.27</td></tr><tr><td>9513.71118</td><td>9.36</td><td>9506.79183</td><td>8.75</td><td>9513.71545</td><td>8.26</td></tr><tr><td>9515.79210</td><td>9.35</td><td>9506.79250</td><td>8.73</td><td>9513.71611</td><td>8.3</td></tr><tr><td>9515.79276</td><td>9.41</td><td>9517.72696</td><td>8.72</td><td>9515.79695</td><td>8.33</td></tr><tr><td>9517.72494</td><td>9.42</td><td>9517.72762</td><td>8.7</td><td>9517.72906</td><td>8.42</td></tr><tr><td>9517.72559</td><td>9.42</td><td>9518.73394</td><td>8.7</td><td>9517.72972</td><td>8.4</td></tr><tr><td>9518.73191</td><td>9.43</td><td>9518.73461</td><td>8.69</td><td>9518.73604</td><td>8.4</td></tr><tr><td>9518.73256</td><td>9.44</td><td>9519.74207</td><td>8.71</td><td>9518.73670</td><td>8.43</td></tr><tr><td>9519.74005</td><td>9.41</td><td>9519.74274</td><td>8.71</td><td>9519.74416</td><td>8.43</td></tr><tr><td>9519.74071</td><td>9.41</td><td>9520.74628</td><td>8.68</td><td>9519.74482</td><td>8.41</td></tr><tr><td>9520.74424</td><td>9.38</td><td>9520.74695</td><td>8.68</td><td>9520.74837</td><td>8.42</td></tr><tr><td>9520.74490</td><td>9.37</td><td>—</td><td>—</td><td>9520.74903</td><td>8.39</td></tr></tbody></table>
<table><caption>**Table 2.** Journal of OHP/SOPHIE observations.</caption><thead><tr><th>Julian date<br>(2 450 000+)</th><th>S/N</th><th>V<sub>r</sub><br>km s<sup>-1</sup></th><th>v sin i<br>km s<sup>-1</sup></th><th>CCF span<br>km s<sup>-1</sup></th></tr></thead><tbody><tr><td>9499.5168</td><td>52</td><td>14.03</td><td>12.12</td><td>0.39</td></tr><tr><td>9500.5044</td><td>44</td><td>15.17</td><td>12.04</td><td>-0.25</td></tr><tr><td>9501.6560</td><td>44</td><td>14.94</td><td>11.99</td><td>0.36</td></tr><tr><td>9502.5564</td><td>65</td><td>15.14</td><td>12.3</td><td>-0.77</td></tr><tr><td>9503.6296</td><td>59</td><td>14.91</td><td>12.24</td><td>-0.44</td></tr><tr><td>9504.6045</td><td>67</td><td>14.5</td><td>12.78</td><td>0.21</td></tr><tr><td>9505.5450</td><td>63</td><td>14.89</td><td>12.49</td><td>-0.22</td></tr><tr><td>9506.4838</td><td>59</td><td>15.59</td><td>12.46</td><td>-0.67</td></tr><tr><td>9510.5691</td><td>61</td><td>15.07</td><td>12.96</td><td>-0.33</td></tr><tr><td>9511.5284</td><td>67</td><td>15.23</td><td>12.08</td><td>-0.75</td></tr><tr><td>9512.5162</td><td>50</td><td>16.85<sup>(†)</sup></td><td>9.74<sup>(†)</sup></td><td>-2.19<sup>(†)</sup></td></tr><tr><td>9513.6282</td><td>52</td><td>15.13</td><td>11.76</td><td>0.09</td></tr><tr><td>9514.5918</td><td>59</td><td>14.98</td><td>12.06</td><td>-0.85</td></tr><tr><td>9515.6107</td><td>42</td><td>15.19</td><td>11.53</td><td>-1.02</td></tr><tr><td>9516.5321</td><td>49</td><td>14.46</td><td>12.7</td><td>0.24</td></tr></tbody></table>

**Notes.** (†)Uncertain value due to lunar contamination (see text).

cross-correlation analysis of nearly 7000 spectral lines between the observed spectrum and a K5 spectral mask template (e.g., Boisse et al. 2010). We list the values of these parameters in Table 2. The mean formal error provided by the SOPHIE pipeline on the $V_r$ measurement is 0.013 km s$^{-1}$. This accuracy is well suited to investigating the significantly larger amplitude of photospheric line profile variability induced by surface spots and/or accretion flows in young stars (e.g., Petrov et al. 2001). No error is provided by the pipeline for the $v \sin i$ measurements, and we assumed that an upper limit is given by the rms deviation of the individual measurements, excluding JD 2 459 512 (see below), namely 0.38 km s$^{-1}$.

The SOPHIE spectrograph includes a second fiber that simultaneously records the spectrum of the nearby sky. Inspection of the cross-correlation function (CCF) of the sky fiber with a synthetic mask of spectral type G2 revealed that the signature of the moon becomes apparent at the expected barycentric Earth radial velocity from JD 2 459 505 onward because the growing moon approaches the target. The lunar contamination culminates on JD 2 459 512, as the bright moon is located about 10 degrees away from GM Aur, which explains the discrepant values measured for $V_r$ and $v \sin i$ on this date. Table 2 suggests that except for JD 2 459 512, the contamination of the CCF by the moon only marginally impacts the $V_r$ and $v \sin i$ measurements. However, to be conservative, we only considered the $V_r$ and $v \sin i$ measurements obtained from the first six spectra of the observing run for the subsequent analysis, from JD 2 459 499 to JD 2 459 504, where no lunar contamination is present.

The journal of observations is given in Table 2. It lists the Julian date, the signal-to-noise ratio of individual spectra at 600 nm, the radial and rotational velocities derived from each spectrum, and the bisector span computed from the cross-correlation function (Queloz et al. 2001).

## 2.4. CFHT SPIRou: Near-infrared spectropolarimetry

Near-infrared spectropolarimetric observations of GM Aur were performed at CFHT using the SPIRou near-infrared spectropolarimeter. It has a spectral range covering from 0.95 to 2.50 µm in a single exposure at a spectral resolution of 70 000 (Donati et al.

2020b). The observations were completed in the framework of the CFHT SPIRou Legacy Survey over four observing runs extending from September 15 to December 18, 2021. Each monthly run was scheduled around the full moon, and we aimed at obtaining one spectrum per night during the run. An additional single spectrum was obtained on January 6, 2022. We thus gathered 34 spectra, whose temporal sampling is shown in Fig. 1.

Each spectrum consists of four polarimetric sub-exposures<fnref n="6" /> for a total integration time of 2,200 s. Individual exposures were combined to yield a single spectrum with a signal-to-noise ratio ranging from ~60 to 100. In one instance, on JD 2 459 503.08, the polarimetric sequence was aborted, and the spectrum consists of a single sub-exposure. The raw data were reduced within the SPIRou consortium, using version V6.132 of the APERO pipeline (Cook et al. 2022). Spectra were cross-correlated with a K2 spectral mask template over about 6,700 spectral lines, and the radial velocity of the object was derived with sub-km s$^{-1}$ accuracy (0.08 km s$^{-1}$ mean rms uncertainty) by fitting a Gaussian to the resulting CCF. The median $V_r$ amounts to 14.65 ± 0.27 km s$^{-1}$. Using TWA 9A, a WTTS of spectral type K5, as a template, the veiling was derived in the JHK bands following the procedure described in Sousa et al. (2023). The median rms errors on veiling measurements in the JHK bands are 0.011, 0.025, and 0.027, respectively.

The journal of observations is presented in Table 3. It lists the Julian date, the signal-to-noise ratio at 2.16 µm, the photospheric radial velocity and its uncertainty, and the veiling in the JHK bands.

## 2.5. ESO ExTrA: Low-resolution near-infrared spectroscopy

The ExTrA facility (Bonfils et al. 2015), located at La Silla Observatory in Chile, consists of three 60 cm telescopes and a single near-infrared (0.88–1.55 µm) fibre-fed spectrograph. We observed GM Aur on 89 nights between October 13, 2021, and March 8, 2022, using either one telescope (21 nights) or two telescopes simultaneously (68 nights). Five fiber units are located at the focal plane of each telescope, each consisting of two 8" aperture fibers. One fiber is used to observe a star and the other is used to observe the nearby sky background. We observed GM Aur with one fiber unit and used another fiber unit to simultaneously observe 2MASS J04535474+3021441 ($J = 8.450 \pm 0.027$ mag) as a comparison star to compute differential photometry. We used the higher-resolution mode of the spectrograph ($R \sim 200$) and 300-second exposures. We obtained between 1 and 30 exposures per night for a total of 1898 spectra with a median signal-to-noise ratio of 105 at 1.05 µm for GM Aur. The ExTrA data were corrected for dark current, extracted using the flat-field, corrected for sky background emission, and were wavelength calibrated using custom data reduction software. Median spectra of GM Aur were computed for each night and telescope, yielding a total of 157 spectra with a median signal-to-noise ratio of 179 and a standard deviation of 62.

We computed differential photometry of GM Aur relative to the comparison star by integrating the individual ExTrA spectra over the J-filter passband. The UKIRT-WFCAM filter transmission curves were retrieved from the SVO Filter Profile Service<fnref n="7" /> (Rodrigo et al. 2012; Rodrigo & Solano 2020). We multiplied each corrected individual spectrum of GM Aur and the comparison star by the filter transmission curve, integrated

<footnote n="6">⁶ The polarimetric analysis of the dataset will be published in a companion paper (Zaire et al., in prep.).</footnote>

<footnote n="7">⁷ http://svo2.cab.inta-csic.es/theory/fps/</footnote>
<table><caption>**Table 3.** Journal of CFHT/SPIRou observations.</caption><thead><tr><th>Julian date<br>(2 450 000+)</th><th>S/N</th><th>V<sub>r</sub><br>km s<sup>-1</sup></th><th>σV<sub>r</sub></th><th>r<sub>J</sub></th><th>r<sub>H</sub></th><th>r<sub>K</sub></th></tr></thead><tbody><tr><td colspan="7" style="text-align: center; font-weight: bold;">September 2021</td></tr><tr><td>9473.066</td><td>101</td><td>14.77</td><td>0.06</td><td>0.05</td><td>0.05</td><td>0.27</td></tr><tr><td>9475.043</td><td>109</td><td>14.09</td><td>0.17</td><td>0.05</td><td>0.06</td><td>0.22</td></tr><tr><td>9476.969</td><td>122</td><td>14.88</td><td>0.16</td><td>0.12</td><td>0.09</td><td>0.27</td></tr><tr><td>9478.027</td><td>119</td><td>14.85</td><td>0.04</td><td>0.17</td><td>0.14</td><td>0.38</td></tr><tr><td>9480.082</td><td>121</td><td>14.62</td><td>0.09</td><td>0.05</td><td>0.09</td><td>0.29</td></tr><tr><td>9481.086</td><td>113</td><td>14.42</td><td>0.08</td><td>0.06</td><td>0.10</td><td>0.23</td></tr><tr><td>9482.086</td><td>120</td><td>14.70</td><td>0.06</td><td>0.05</td><td>0.07</td><td>0.23</td></tr><tr><td colspan="7" style="text-align: center; font-weight: bold;">October 2021</td></tr><tr><td>9502.074<sup>(†)</sup></td><td>112</td><td>14.97</td><td>0.06</td><td>0.15</td><td>0.14</td><td>0.32</td></tr><tr><td>9503.074</td><td>64</td><td>14.48</td><td>–</td><td>0.26</td><td>0.15</td><td>0.07</td></tr><tr><td>9504.098</td><td>94</td><td>14.61</td><td>0.05</td><td>0.10</td><td>0.11</td><td>0.30</td></tr><tr><td>9506.082</td><td>112</td><td>14.65</td><td>0.01</td><td>0.04</td><td>0.12</td><td>0.29</td></tr><tr><td>9508.082</td><td>113</td><td>14.68</td><td>0.08</td><td>0.18</td><td>0.24</td><td>0.43</td></tr><tr><td>9509.086</td><td>110</td><td>14.85</td><td>0.06</td><td>0.32</td><td>0.26</td><td>0.51</td></tr><tr><td>9510.086</td><td>92</td><td>14.40</td><td>0.12</td><td>0.41</td><td>0.29</td><td>0.58</td></tr><tr><td>9511.074</td><td>116</td><td>14.33</td><td>0.05</td><td>0.02</td><td>-0.04</td><td>0.28</td></tr><tr><td>9513.090</td><td>98</td><td>14.70</td><td>0.25</td><td>0.09</td><td>0.13</td><td>0.29</td></tr><tr><td>9514.051</td><td>105</td><td>14.66</td><td>0.16</td><td>0.18</td><td>0.12</td><td>0.36</td></tr><tr><td>9515.078</td><td>108</td><td>14.66</td><td>0.03</td><td>0.10</td><td>0.11</td><td>0.28</td></tr><tr><td>9516.102</td><td>104</td><td>14.73</td><td>0.16</td><td>0.10</td><td>0.09</td><td>0.26</td></tr><tr><td colspan="7" style="text-align: center; font-weight: bold;">November 2021</td></tr><tr><td>9535.102</td><td>102</td><td>14.27</td><td>0.09</td><td>0.13</td><td>0.07</td><td>0.25</td></tr><tr><td>9537.098</td><td>101</td><td>14.95</td><td>0.15</td><td>0.11</td><td>0.08</td><td>0.26</td></tr><tr><td>9538.039</td><td>104</td><td>14.22</td><td>0.10</td><td>0.05</td><td>0.11</td><td>0.22</td></tr><tr><td>9538.984</td><td>95</td><td>14.40</td><td>0.10</td><td>0.12</td><td>0.08</td><td>0.28</td></tr><tr><td>9539.988</td><td>107</td><td>14.51</td><td>0.03</td><td>0.14</td><td>0.13</td><td>0.27</td></tr><tr><td>9541.012</td><td>110</td><td>13.96</td><td>0.10</td><td>0.14</td><td>0.14</td><td>0.38</td></tr><tr><td colspan="7" style="text-align: center; font-weight: bold;">December 2021</td></tr><tr><td>9557.980</td><td>110</td><td>14.42</td><td>0.11</td><td>0.22</td><td>0.19</td><td>0.50</td></tr><tr><td>9558.949</td><td>110</td><td>14.12</td><td>0.03</td><td>0.21</td><td>0.19</td><td>0.45</td></tr><tr><td>9559.961</td><td>95</td><td>14.71</td><td>0.11</td><td>0.09</td><td>0.11</td><td>0.35</td></tr><tr><td>9560.965</td><td>94</td><td>15.04</td><td>0.06</td><td>0.07</td><td>0.12</td><td>0.35</td></tr><tr><td>9563.078</td><td>106</td><td>14.73</td><td>0.06</td><td>0.21</td><td>0.06</td><td>0.43</td></tr><tr><td>9564.059</td><td>102</td><td>14.34</td><td>0.10</td><td>0.32</td><td>0.11</td><td>0.43</td></tr><tr><td>9566.051</td><td>96</td><td>14.37</td><td>0.13</td><td>0.07</td><td>0.03</td><td>0.34</td></tr><tr><td>9566.953</td><td>95</td><td>14.93</td><td>0.03</td><td>0.06</td><td>0.02</td><td>0.25</td></tr><tr><td colspan="7" style="text-align: center; font-weight: bold;">January 2022</td></tr><tr><td>9585.941</td><td>112</td><td>14.30</td><td>0.06</td><td>0.11</td><td>0.12</td><td>0.25</td></tr></tbody></table>

Notes. (†) Single sub-exposure.

the flux, and obtained a magnitude difference from the flux ratio. We computed a differential magnitude measurement for each night as the mean and standard deviation of the individual measurements taken on that night. We then derived the J-band magnitude of GM Aur from the 2MASS magnitude of the comparison star. We obtained a median value and 68.3% confidence interval of $J = 9.417 \pm 0.061$ mag for the ExTrA observations. The results are listed in Table C.1.

# 3. Results

## 3.1. Multicolor photometry

The LCOGT $u'g'r'i'$ light curves of GM Aur span nearly four months and are shown in Fig. 2. A TESS light curve extending over 50 days, from JD 2459474 to JD 2459524, is also

shown. We derive a mean magnitude and rms variability of $u' = 14.32 \pm 0.35$ mag, $g' = 12.77 \pm 0.17$ mag, $r' = 11.64 \pm 0.11$ mag, and $i' = 11.21 \pm 0.09$ mag, respectively. The variability amplitude is much larger than the photometric errors in all bands; it amounts to 0.025 mag in the $g'r'i'$ bands and 0.033 mag in the $u'$ band, and it decreases toward longer wavelengths, from $u'$ to $i'$. We ran two period-search algorithms on the light curves: a CLEAN periodogram analysis (Roberts et al. 1987), which is conceptually similar to a Fourier transform, and the String-Length method (Dworetsky 1983), which finds the period that minimizes the average distance between consecutive points in the phased light curve. Both yielded the same period at all wavelengths, namely $P = 6.04 \pm 0.15$ days, with the uncertainty being estimated from the standard deviation of a Gaussian fitted to the periodogram peak. The results of the period search are shown in Fig. 3. The light curves folded in phase at the 6.04 d period are shown in Fig. 2. They clearly show the modulation of the brightness level at this period, particularly in the $u'$ band, where the amplitude is largest. We ascribe this low-level modulation to surface spots and the $P = 6.04$ d period to stellar rotation. This estimate agrees within the error bars with the previously reported periods for GM Aur, namely $P = 6.1$ d by Percy et al. (2010) and $P = 6.02$ d by Artemenko et al. (2012). We therefore adopt the following ephemeris:

$$ \text{HJD}(d) = 2459460.80 + 6.04 \times E, \quad (1) $$

which defines the rotational phase $\Phi_{\text{rot}} = (\text{HJD} - 2459460.80)/P_{\text{rot}}$ (modulo $P_{\text{rot}}$), where phase zero ($\Phi_{\text{rot}} = 0$) is chosen as the epoch of maximum optical brightness in the rotational cycle.

Superimposed onto spot modulation, additional signs of intrinsic variability are visible. This is the case, for instance, of an apparent brightness event that is centered on JD 2459509 and lasted for a few days, as well as more pronounced wide dips toward the end of the observing period, from JD 2459544 on. Interestingly, the first half of the light curve exhibits several day-long brightening events, while the last third is dominated by wide dimming events. We note that most of the brightening events, with $g' \le 12.5$ mag, tend to occur at rotational phases shorter than 0.15 or longer than 0.85, that is, close to the maximum brightness of the spot modulation. If the photometric variations of the system are modulated by the visibility of a bright accretion spot at the stellar surface, these brightening events seen around the time of maximum accretion shock visibility most likely reflect varying accretion on the stellar surface. The wide dips in the last part of the light curve exhibit the same periodicity and phase as the spot modulation. They reach a minimum brightness close to phase 0.5, with an amplitude that steadily decreases from 0.6 mag to 0.3 mag in the $g'$ band over a timescale of a few weeks, from JD 2459550 to JD 2459570.

Figure 4 shows the color behavior of the system. As the system fades, it becomes redder, with a color slope larger than expected from ISM-like extinction at least in the ($u' - g'$) and ($g' - r'$) color indices. According to Venuti et al. (2015), a large color slope at short wavelengths is characteristics of accretion-driven photometric variability. However, both the large scatter seen in the ($u' - g'$) color index at low brightness levels and the changing shape of the light curve during the semester suggest that several sources of variability might be present, such as a combination of accretion and obscuration events.
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** GM Aur light curves. **Top:** GM Aur's *u'*-*g'*-*r'*-*i'* light curves from LCOGT observations that extended over nearly four months. The mean photometric error is 0.025 mag in the *g'*-*r'*-*i'* bands and 0.033 mag in the *u'* band. A scaled TESS light curve (black) obtained contemporaneously is overplotted onto the LCOGT *r'*-band light curve. **Middle:** LCOGT light curves folded in phase with a period of 6.04 days with the ephemeris of Eq. (1). **Bottom:** same as above, with a color scale for data points that reflects the Julian date. The amplitude of the light variations appears to increase slightly toward the end of the observing run. In all panels, the error bars on the measurements are smaller than the symbol size.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Period search results. Top: CLEAN periodogram analysis of the *u'*g'*r'*i'*' light curves. A peak occurs at the frequency of 0.166 day⁻¹, which corresponds to a period of 6.04 ± 0.15 days. Bottom: string-length analysis of the *u'*g'*r'*i'*' light curves. A clear minimum appears for a period of 6.03 ± 0.15 days.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Color-magnitude relation. The (*u'* − *g'*), (*g'* − *r'*), and (*r'* − *i'* colors are plotted as a function of the system brightness in the *g'*-band. The dashed lines indicate the expected ISM reddening slope computed from Fiorucci & Munari (2003) with *R* = 3.1.</figcaption></figure>

## 3.2. High-resolution optical spectroscopy

We took advantage of the wide wavelength range covered by the SOPHIE spectrograph at high spectral resolution to derive the stellar parameters of GM Aur, to measure optical veiling, and to investigate the emission line profiles and their variability across the optical range. These results are described in the next subsections.

### 3.2.1. Stellar parameters

To derive the stellar parameters (*T*<sub>eff</sub>, v sin i, *V*<sub>r</sub>, and *V*<sub>mic</sub>), we fit synthetic ZEEMAN spectra (Landstreet 1988; Wade et al. 2001; Folsom et al. 2012) to the average of the first six SOPHIE spectra, which are not contaminated by the moon. Synthetic spectra were computed from MARCS atmosphere models (Gustafsson et al. 2008) and the VALD line list database (Ryabchikova et al. 2015). We applied a χ² minimization procedure based on a Levenberg-Marquart algorithm over seven wavelength windows ranging from 455 to 649 nm, excluding the regions affected by tellurics, emission, or molecular lines. We set the macroturbulent velocity to 2.0 km s⁻¹ and the surface gravity log *g* to 4.0, and we assumed solar metallicity. These are typical parameters for low-mass TTSs. During the fitting procedure, we applied the mean veiling value derived below for the GM Aur mean optical spectrum (*r*<sub>0.55</sub> = 0.3, see Sect. 3.2.2) to the synthetic spectra. Finally, we averaged the results obtained from the various wavelength windows, except for one window that yielded values higher than 2σ from the mean. We thus derived *T*<sub>eff</sub> = 4287 ± 35 K, *V*<sub>r</sub> = 14.94 ± 0.14 km s⁻¹, *V*<sub>mic</sub> = 3.3 ± 0.4 km s⁻¹, and v sin i = 14.9 ± 0.3 km s⁻¹. We also used the ROTFIT package (Frasca et al. 2003, 2006) applied to the mean GM Aur's SOPHIE spectrum, from which we derive *T*<sub>eff</sub> = 4505 ± 53 K, *V*<sub>r</sub> = 15.37 ± 0.26 km s⁻¹, and v sin i = 13.0 ± 0.7 km s⁻¹.

While all these values are consistent within 3σ, the large *T*<sub>eff</sub> difference derived from ZEEMAN and ROTFIT illustrates model-dependent uncertainties that are likely related to the use of different model templates (MARCs for ZEEMAN vs. BTSetll for ROTFIT) and possibly to wavelength-dependent systematics induced by starspots (Gangi et al. 2022; Flores et al. 2022). For consistency with similar observing campaigns that we previously performed on young stars (e.g., Pouilly et al. 2020, 2021; Bouvier et al. 2020a; Alencar et al. 2018), we adopted the results of the ZEEMAN fitting. This estimate also agrees better with the K6 spectral type derived by Herczeg & Hillenbrand (2014), from which they deduced *T*<sub>eff</sub> = 4115 K. When the *T*<sub>eff</sub>-SpT conversion tables from Herczeg & Hillenbrand (2014) and Pecaut & Mamajek (2013) are used, the *T*<sub>eff</sub> value derived above corresponds to a K4-K5 spectral type, which agrees fairly well with previous estimates obtained from optical and near-infrared spectroscopy (K5-K6; e.g., Espaillat et al. 2010; Herczeg & Hillenbrand 2014).

The *V*<sub>r</sub> and v sin i values can be compared to those derived from the uncontaminated CCF of the first six SOPHIE spectra, namely ⟨*V*<sub>r</sub⟩ = 14.79 ± 0.40 km s⁻¹ and ⟨v sin i⟩ = 12.25 ± 0.26 km s⁻¹. The quoted uncertainties are the rms of the six measurements and therefore include intrinsic variability of the CCF profiles due to spot modulation, for example. The two estimates of *V*<sub>r</sub> agree within the errors as well as with the median *V*<sub>r</sub> value derived from the SPIRou spectra (see Sect. 2.4), while the v sin i value derived from the CCF is significantly lower than the value deduced from spectral fitting. We suspect that the discrepancy may arise from the color-dependent FWHM-v sin i relation used for SOPHIE CCFs, which is calibrated on main-sequence stars (Boisse et al. 2010) and may not be fully adequate for pre-main-sequence objects.

From the average *g'*r'*i'* colors reported for GM Aur in Sect. 2.1, the intrinsic colors of a K4-K5 dwarf listed by Kraus & Hillenbrand (2007) and Covey et al. (2007), and using ISM extinction coefficients from Fiorucci & Munari (2003), we derive a visual extinction on the line of sight *A*<sub>V</sub> = 0.3 ± 0.3 mag. This agrees with previous determinations
<table><caption>**Table 4.** Properties of the GM Aur system.</caption><thead><tr><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td>SpT</td><td>K4–K5</td></tr><tr><td>A<sub>V</sub></td><td>0.3 ± 0.3 mag</td></tr><tr><td>T<sub>eff</sub></td><td>4287 ± 35 K</td></tr><tr><td>L<sub>★</sub>/L<sub>☉</sub></td><td>0.9 ± 0.2</td></tr><tr><td>R<sub>★</sub>/R<sub>☉</sub></td><td>1.7 ± 0.2</td></tr><tr><td>M<sub>★</sub>/M<sub>☉</sub></td><td>0.95 ± 0.13</td></tr><tr><td>M&#775;<sub>acc</sub></td><td>0.7 ± 0.3 10<sup>-8</sup> M<sub>☉</sub> yr<sup>-1</sup></td></tr><tr><td>EW(LiI)</td><td>420 ± 23 mÅ</td></tr><tr><td>v sin i</td><td>14.9 ± 0.3 km s<sup>-1</sup></td></tr><tr><td>P<sub>rot</sub></td><td>6.04 ± 0.15 d</td></tr><tr><td>r<sub>cor</sub></td><td>8.3 ± 0.5R<sub>★</sub> (0.064 au)</td></tr></tbody></table>

(e.g., Herczeg & Hillenbrand 2014). The REM photometry reported above yields a median J-band magnitude of 9.42 ± 0.11 mag, which is close to that of 2MASS (J = 9.34 mag). We dereddened the median J-band magnitude with $A_j = 0.28 \times A_V = 0.084 \pm 0.084$ mag and used the J-band bolo- metric correction listed in Pecaut & Mamajek (2013) for a K4-K5 dwarf, $BC_J = 1.55 \pm 0.03$ mag, to derive the stellar luminosity, $L_★ = 0.9 \pm 0.2 L_☉$, assuming the Gaia distance of 157.9 ± 1.2 pc (Gaia Collaboration 2021). We thus derive a stellar radius $R_★ = 1.7 \pm 0.2 R_☉$ from Stefan's law, and a stellar mass $M_★ = 1.05 \pm 0.05 M_☉$ from the Siess et al. (2000) pre-main-sequence evolution models, while CESAM models (Marques et al. 2013) yield $M_★ = 0.88 \pm 0.12 M_☉$ (E. Alécian, priv. comm.). We therefore adopt $M_★ = 0.95 \pm 0.13 M_☉$, in agreement with Baraffe et al. (2015) models. This estimate is also consistent within the errors with the dynamical mass estimate, $M_{dyn} = 1.00 \pm 0.02 M_☉$, reported by Guilloteau et al. (2014), which was later revised to $M_{dyn} = 1.14 \pm 0.02 M_☉$ by Simon et al. (2019). Table 4 summarizes the derived stellar parameters.

Finally, we combined the $v \sin i$ and rotational period measurements with the stellar radius estimate to derive the stellar inclination $\sin i = P_{rot} \times v \sin i/(2 \pi R_★) = 1.05$ with the values listed in Table 4. Accounting for $1\sigma$ uncertainties on the stellar parameters, we derive a lower limit of $i_★ \ge 63^\circ$ for the stellar rotational axis onto the line of sight. This value is significantly higher than the inclination inferred from high-resolution ALMA images of the outer disk of GM Aur observed at millimeter wavelength, which yield $i_{disk} = 53.2$ deg (Huang et al. 2020). It agrees better, however, with the inclination value derived for the disk seen in scattered light with adaptive optics on a scale of a few dozen au, for which Oh et al. (2016) obtained $i_{disk} = 64 \pm 2$ deg. On the much smaller scale of 0.013 au, Bohn et al. (2022) measured an inner-disk inclination $i_{disk} = 68^{+16}_{-28}$ deg from long-baseline K-band interferometry using VLTI/GRAVITY. They did not find evidence for an inner and outer disk misalignment in this system. As outlined by Appenzeller & Bertout (2013), the uncertainty on the determination of stellar inclinations from rotation measurements rapidly increases at large angles and is prone to systematic errors. For GM Aur, inferring the stellar inclination from the disk inclination might therefore be more reliable than estimating it from rotation measurements, owing in particular to the significant uncertainty on the stellar radius. Nevertheless, all independent measurements indicate a moderate to high inclination for the system.

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Optical veiling measurements from the OHP/SOPHIE spectra. Top: part of a spectral window showing the template spectrum (gray), the observed spectrum (black), and the velocity broadened and veiled template spectrum (red) that fits the observed spectrum. Bottom: veiling measured in several spectral windows (see text) as a function of Julian date. The central wavelength of the spectral windows is indicated in the top left corner of the panel.</figcaption></figure>

### 3.2.2. Veiling

At optical wavelengths, accreting T Tauri stars exhibit an additional source of continuum flux, which presumably arises from the accretion shock on the stellar surface. This optical excess continuum partly veils the stellar photospheric spectrum and is therefore referred to as “veiling” (Hartigan et al. 1995). We measured the optical veiling from the high-resolution OHP/SOPHIE spectra by comparing the photospheric spectrum of GM Aur to that of the unveiled nonaccreting template V819 Tau, a WTTS with $T_{eff} = 4250 \pm 50$ K, $V_r = 16.6$ km s$^{-1}$, and $v \sin i = 9.5$ km s$^{-1}$ (Donati et al. 2015). We retrieved an archival CFHT/ESPADOnS spectrum of V819 Tau, which we resampled at the spectral resolution of OHP/SOPHIE spectra, translated into the radial velocity of GM Aur, and rotationally broadened using the rotational function from Gray (1973) to match the rotational velocity of GM Aur. The template spectrum was then fit to the GM Aur spectra over the same wavelength windows as discussed in Sect. 3.2.1 by adjusting the veiling using the following formula:

$$I = \frac{I_0 + r}{1 + r}, \qquad (2)$$

where I is the veiled spectrum, $I_0$ is the spectrum without veiling, and r is the continuum veiling. Rei et al. (2018) showed that an additional line-dependent veiling component may be present in the strongest photospheric lines. We therefore retained only weak to moderate lines with an EW between 0.01 and 0.1 Å to perform the fit. The veiling measured for each spectrum in each spectral window is shown in Fig. 5. Systematic offsets are clearly seen between the veiling values measured in the different spectral windows. These offsets may partly reflect a wavelength-dependent veiling, as veiling seems to decrease toward longer wavelengths for all but the bluest spectral window, but it may also result from systematic errors depending on the specific sample of lines included in each window. We therefore computed an average veiling value, $r_{0.55}$, over all spectral windows for individual GM Aur spectra. This value is listed in Table 5 with its
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Mean optical veiling plotted as a function of rotational phase. The color code indicates the Julian date.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** LiI 6707 Å line profile. Left: the 34 line profile measurements from SOPHIE spectra are shown superimposed. The color code corresponds to successive rotational cycles. Right: 2D periodogram across the line profile. The dotted horizontal line drawn at a frequency of 0.166 day⁻¹ indicates the stellar rotational period. The white curve displays the mean line profile. The color code reflects the periodogram power.</figcaption></figure>

associated uncertainty. The optical veiling is moderate, ranging from 0.17 to 0.55 at 5500 Å. Similar mean values were derived from ROTFIT, namely *r* = 0.5, 0.4, and 0.3 at 4500, 6000, and 6500 Å, respectively, using a library of 400 spectral templates with spectral types FGKM from the OHP/ELODIE database, which yielded a best χ² fit for spectral types ranging from K3.5 to K5. Figure 6 shows the mean optical veiling plotted along the rotational phase. A hint of higher veiling values towards phases 0 and 1, and minimum values around phase 0.5–0.6 appears. A periodogram analysis of the mean veiling variation did not yield significant results, however.

The lack of significant veiling variability is confirmed by the examination of the LiI 6707 Å photospheric line profile shown in Fig. 7. The shape and depth of the line appear to remain quite stable over the duration of the observing run, and a periodogram analysis across the line profile (Giampapa et al. 1993) revealed no signs of periodic modulation. This suggests that the source of optical veiling remains at least partly in view throughout the rotational cycle. This might arise from a high-latitude accretion shock.

### 3.2.3. Emission lines

The main emission lines seen in the optical spectrum of GM Aur, namely Hα, Hβ, and HeI 5876 Å, are displayed in Fig. 8.

The Balmer lines show a broad emission peak and a slightly blueshifted absorption component, whose depth varies from being hardly discernable to reaching below the continuum in the Hβ profile. The wide, slightly blueshifted absorption component peaks at −30 to −20 km s⁻¹ and covers a velocity range from about −90 to +40 km s⁻¹ in the Hβ line profile. Additional absorption components appear at higher blueshifted velocities, peaking from −110 to −90 km s⁻¹, and extending down to about −160 km s⁻¹. These blueshifted absorption components cause most of the variability in the Balmer line profile. The Hβ profile also exhibits significant variability over the red wing, up to velocities of +200 km s⁻¹. However, none of the profiles exhibits redshifted absorption components reaching below the continuum level.

Owing to the complex line shapes, equivalent widths (EW) of the Hα, Hβ, and HeI 5876 Å lines were computed by directly integrating below the line profile. The results are listed in Table 5. The measurement accuracy is 10% or better for EW(Hα) and EW(Hβ), and about 20% for EW(HeI) due to the more uncertain continuum location. We note that on JD 2459512, the EWs measurements are systematically lower than during the rest of the observations, which might be due to lunar contamination, as discussed in Sect. 2.3. The average and rms values we obtain are EW(Hα) = 83 ± 8 Å, EW(Hβ) = 9.3 ± 2.5 Å, and EW(HeI) = 0.40 ± 0.13 Å.

The HeI 5876 Å line profile is roughly symmetric and consists of a narrow component (FWHM ~ 30–40 km s⁻¹) superimposed on a broad pedestal, as previously reported by McGinnis et al. (2020). The peak intensity of the narrow component varies significantly, while the broad component appears relatively stable (see Fig. 8). We fit the HeI line profile with a two-component Gaussian model to extract the properties of the narrow (NC) and broad (BC) components. The EWs were derived from the Gaussian fit of the components. The radial velocity, FWHM, and EW of the NC and BC, as well as their uncertainty derived from the covariance matrix of the Levenberg-Marquart fitting algorithm, are listed in Table 6.

We derived the radial velocity of the narrow component and found it to be variable and redshifted by ~5–10 km s⁻¹ relative to the stellar velocity. Figure 9 shows the HeI NC radial velocity curve plotted as a function of Julian date and rotational phase. As previously reported by McGinnis et al. (2020), HeI NC $V_r$ appears to be rotationally modulated with a full amplitude of ~6 km s⁻¹, as expected if the NC component of the HeI 5876 Å line profile were produced in a high-latitude accretion shock at the stellar surface. We fit the observed NC $V_r$ curve with the geometrical accretion shock model described in Pouilly et al. (2021). The model computes the variation of the HeI NC radial velocity as the combination of the rotational modulation of the accretion shock and the intrinsic inflow velocity. The free parameters of the model are the inflow velocity, the latitude of the accretion shock, the phase at which it faces the observer, and the stellar inclination. The HeI NC $V_r$ curve is best reproduced with an accretion shock located at a latitude of 83 ± 1.5° that faces the observer at phase 0.2 ± 0.08 and has a radial post-shock velocity of 18.3 ± 1.0 km s⁻¹ in the stellar rest frame. Because the model now includes the inflow velocity, the HeI NC radial velocity curve does not reach the mean velocity at the time when the spot faces the observer, as it would for the case of static stellar spots. The stellar inclination we derive from the model is $i = 64 \pm 2.2°$, which is consistent with the inner disk inclination derived from K-band VLTI/GRAVITY data (see Sect. 3.2.1). According to the model, the accretion shock faces the observer

<footnote n="8">⁸ The [OI] 6300 Å line is also seen in emission in the spectrum. However, it is significantly contaminated by sky, which cannot be easily corrected for due to the different response of the object and sky fibers of the SOPHIE spectrograph. The other emission lines seen in the high signal-to-noise ratio mean SOPHIE spectrum are Ca II H&K, Hγ, Hδ, He, and He I 6678 Å, the latter being weak and affected by a deep photospheric FeI line.</footnote>
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Optical line profile variability. **Top:** series of optical line profiles Hα, Hβ, and HeI plotted as a function of Julian date (left subpanels) and rotational phase (right subpanels). The color code corresponds to successive rotational cycles. **Middle:** same profiles, superimposed to illustrate their variability. **Bottom:** 2D periodograms across the line profiles. The color code reflects the periodogram power, from zero (blue) to 0.6 (red). The dotted horizontal red line drawn at a frequency of 0.166 day⁻¹ indicates the stellar rotational period. The white curve is the mean line profile.</figcaption></figure>

close to the origin of the phase, which is consistent with the photometric behavior described above (see Sect. 3.1). The HeI line profile is also strongest over rotational phases ranging from 0.75 to 1.09 (excluding the probable accretion burst occurring at JD 2459510, see Sect. 3.1), and this is also when the highest veiling

values are measured in the JHK bands (see Fig. D.1). These two results further support a maximum visibility of the accretion shock close to $\Phi_{\text{rot}} = 0$.

A bidimensional periodogram analysis (Giampapa et al. 1993) of the Balmer and HeI 5876 Å line profiles reveals a
<table><caption>**Table 5.** Optical line EW and veiling measurements from the OHP/SOPHIE spectra.</caption><thead><tr><th rowspan="2">Julian date<br>(2 450 000+)</th><th colspan="3">EW</th><th colspan="2">Veiling</th></tr><tr><th>Hα<br>Å</th><th>Hβ<br>Å</th><th>HeI<br>Å</th><th>r<sub>0.55</sub></th><th>rms</th></tr></thead><tbody><tr><td>9499.51683266</td><td>92.3</td><td>9.2</td><td>0.31</td><td>0.35</td><td>0.19</td></tr><tr><td>9500.50445956</td><td>81.1</td><td>6.8</td><td>0.37</td><td>0.17</td><td>0.15</td></tr><tr><td>9501.65604543</td><td>86</td><td>10.3</td><td>0.5</td><td>0.26</td><td>0.16</td></tr><tr><td>9502.55640645</td><td>73.3</td><td>7.4</td><td>0.45</td><td>0.46</td><td>0.17</td></tr><tr><td>9503.62969684</td><td>79.6</td><td>8.5</td><td>0.4</td><td>0.4</td><td>0.15</td></tr><tr><td>9504.60450624</td><td>88.1</td><td>10.5</td><td>0.44</td><td>0.3</td><td>0.14</td></tr><tr><td>9505.54502964</td><td>85.8</td><td>10.2</td><td>0.23</td><td>0.27</td><td>0.13</td></tr><tr><td>9506.48382743</td><td>86.4</td><td>9</td><td>0.33</td><td>0.37</td><td>0.15</td></tr><tr><td>9510.56919423</td><td>102</td><td>15.1</td><td>0.74</td><td>0.51</td><td>0.18</td></tr><tr><td>9511.52849738</td><td>91.1</td><td>12.6</td><td>0.34</td><td>0.38</td><td>0.16</td></tr><tr><td>9512.51625045</td><td>68.3<sup>(†)</sup></td><td>3.7<sup>(†)</sup></td><td>0.18<sup>(†)</sup></td><td>0.29<sup>(†)</sup></td><td>0.13<sup>(†)</sup></td></tr><tr><td>9513.62825832</td><td>82.6</td><td>10.1</td><td>0.36</td><td>0.37</td><td>0.15</td></tr><tr><td>9514.59185217</td><td>81.5</td><td>9</td><td>0.5</td><td>0.55</td><td>0.23</td></tr><tr><td>9515.61072618</td><td>76.8</td><td>7.7</td><td>0.42</td><td>0.3</td><td>0.19</td></tr><tr><td>9516.53216217</td><td>74.7</td><td>9.6</td><td>0.36</td><td>0.27</td><td>0.15</td></tr></tbody></table>

Notes. (†)Uncertain value due to lunar contamination (see text).

periodic modulation of part of the profiles (see Fig. 8). The intensity of the narrow component of the HeI line profile is modulated at a frequency of ~0.15 d⁻¹, corresponding to a period of 6.7 d, with a large uncertainty due to the limited temporal sampling of the spectral series, however, that translates into a poor frequency resolution. As the HeI line profile NC component arises in the accretion shock (Beristain et al. 2001), it is expected to be modulated at the stellar rotational period or close to it in case of latitudinal differential rotation at the stellar surface. In the Balmer line profiles, a modulation close to the stellar period also appears over three distinct locations: in the highly redshifted wing over the velocity range ~200–400 km s⁻¹, at slightly redshifted velocities between 0 and ~50 km s⁻¹, and at highly blueshifted velocities from ~400 to ~200 km s⁻¹. While the maximum power of the periodogram in the blue wing appears at the stellar rotation period, it seems to drift to longer periods, similar to that seen in the HeI NC component, toward the red wing. If the red wing modulation of Balmer line profiles is caused by the absorption of shock emission by a funnel flow crossing the line of sight, this might be an indication of differential rotation along the funnel flow. Noticeably, no sign of periodic variability is seen in the Balmer line profiles in the velocity channels in which the variable blueshifted absorption components arise. This suggests that these components either result from sporadic ejection processes or that they vary on periods longer than ten days.

Correlation matrices (Johns & Basri 1995) between line profiles were computed and are presented in Appendix B. They display the degree to which temporal flux variations in a pair of spectral lines are correlated. Matrices can be computed for a single profile (autocorrelation), for instance, $Hα★Hβ$, to investigate how the different parts of the profile vary with respect to each other, or between two profiles (cross-correlation), for example, $Hα★Hβ$, to compare the intensity variations of different lines. The $Hα★Hα$ and $Hβ★Hβ$ matrices shown in Fig. B.1 are quite similar. The blue and red wings of the profiles vary in a correlated way. The high-velocity red wings are anticorrelated with the emission peak region. This may be a sign that high-velocity redshifted absorption components appear when the peak intensity is higher, although the absorptions do not reach below the continuum. Many absorption features are superimposed on the emission line component. These features do not present a periodicity and are not correlated with each other, nor with the emission part of the profile. The $Hβ★Hα$ matrix mimics the matrices of $Hα★Hα$ and the $Hβ★Hβ$, showing that both lines are formed in the same region, as expected from magnetospheric accretion models (e.g., Muzerolle et al. 2001).

### 3.3. High-resolution near-infrared spectroscopy

The main emission lines seen in the high-resolution near-infrared spectrum of GM Aur, namely HeI 10830 Å, Paβ, and Bry, are depicted in Figs. 10 and 11⁹. The near-infrared hydrogen lines exhibit broad emission profiles, with FWHM ~200 km s⁻¹, whose peaks are slightly blueshifted compared to the stellar velocity, and are located at ~30 and ~10 km s⁻¹ for Paβ and Bry, respectively. The profiles are roughly symmetric, but have pronounced high-velocity redshifted absorptions that extend up to +400 km s⁻¹ and reach below the continuum level. Paβ and Bry exhibit inverse P Cygni (IPC) profiles of varying depth in most SPIRou spectra, which is in stark contrast to the optical hydrogen line profiles, $Hα$ and $Hβ$, which do not exhibit such features (see Sect. 3.2.3).

Correlation matrices for the near-infrared hydrogen line profiles are shown in Fig. B.2. The Paβ ★ Paβ, Bry ★ Bry, and Paβ ★ Bry matrices are quite similar. The emission part of the profiles is overall correlated, and anti-correlated with the high-velocity redshifted absorption component. This indicates that the redshifted absorption is deeper when the emission line is more intense. Cross-correlation matrices between optical and near-infrared line profiles observed in the same night during the October runs are shown in Fig. B.3. The Bry and the Paβ emission components correlate well with the main $Hα$ and $Hβ$ emission components (~100 km s⁻¹ < v < 200 km s⁻¹), which suggests that at least part of the line emission forms in the same extended region. The Paβ redshifted absorption, and to a lesser extent, the Bry redshifted absorption, correlates well with the $Hα$ and $Hβ$ blue and red high-velocity wings. This indicates that they all form in the same region close to the star. The near-infrared profiles are strongly anticorrelated with the highly variable blueshifted absorption components that appear at velocities of about ~100 km s⁻¹ in the Balmer line profiles and presumably trace outflows. Except for this feature, the overall correlated variations of optical and near-infrared line profiles suggest that they are all good diagnostics of the accretion flow, as previously suggested by Alcalá et al. (2014), for example. Detailed modeling of the line profile shapes and variability is needed, however, to ascertain the exact location from which they arise (Tessore et al. 2023).

The HeI 10830 Å line profile is more complex. It is usually dominated by an emission peak at low velocities that becomes quite weak at specific times, however. The profile also often displays highly redshifted absorption features, similar to the IPC components seen in the Paβ and Bry lines. Unlike the near-infrared hydrogen profiles however, the HeI line additionally exhibits various absorption components in the blue wing of

<footnote n="9">⁹ Pay and Paδ also appear in SPIRou spectra. Their shape and variability behavior is similar to that of Paβ and Bry. Brδ is also included in the SPIRou wavelength range, but lies at the intersection of spectral orders. In all SPIRou spectra, we also detect a weak H2 2.12 µm line, with EW = 0.22 ± 0.04 Å, and FWHM = 1.40 ± 0.05 Å (≈20 km s⁻¹). This line is located at the stellar rest velocity ($V_r$ = 0.88 ± 0.75 km s⁻¹).</footnote>
<table><caption>**Table 6.** Properties of the narrow (NC) and broad (BC) components of the HeI 5876 Å line profile and their $1\sigma$ error.</caption><thead><tr><th rowspan="2">Julian date<br>(2 450 000+)</th><th colspan="4">V<sub>r</sub><br>km s<sup>-1</sup></th><th colspan="4">FWHM<br>km s<sup>-1</sup></th><th colspan="4">EW<br>Å</th></tr><tr><th>NC</th><th>err</th><th>BC</th><th>err</th><th>NC</th><th>err</th><th>BC</th><th>err</th><th>NC</th><th>err</th><th>BC</th><th>err</th></tr></thead><tbody><tr><td>9499.51683</td><td>11.8</td><td>1.0</td><td>14.8</td><td>1.8</td><td>27.7</td><td>3.6</td><td>100.7</td><td>6.7</td><td>0.07</td><td>0.02</td><td>0.30</td><td>0.05</td></tr><tr><td>9500.50445</td><td>3.3</td><td>1.3</td><td>14.5</td><td>3.1</td><td>38.7</td><td>4.5</td><td>106.6</td><td>9.6</td><td>0.13</td><td>0.03</td><td>0.30</td><td>0.07</td></tr><tr><td>9501.65604</td><td>4.6</td><td>0.8</td><td>20.0</td><td>3.7</td><td>40.7</td><td>2.7</td><td>91.0</td><td>7.4</td><td>0.28</td><td>0.04</td><td>0.29</td><td>0.08</td></tr><tr><td>9502.55640</td><td>7.7</td><td>0.7</td><td>20.0</td><td>2.5</td><td>36.2</td><td>2.4</td><td>81.2</td><td>5.5</td><td>0.22</td><td>0.03</td><td>0.27</td><td>0.06</td></tr><tr><td>9503.62969</td><td>8.6</td><td>0.7</td><td>12.9</td><td>2.3</td><td>37.6</td><td>2.8</td><td>113.2</td><td>9.4</td><td>0.18</td><td>0.03</td><td>0.30</td><td>0.06</td></tr><tr><td>9504.60450</td><td>9.2</td><td>0.8</td><td>10.1</td><td>1.8</td><td>34.0</td><td>2.9</td><td>111.4</td><td>7.3</td><td>0.13</td><td>0.02</td><td>0.36</td><td>0.06</td></tr><tr><td>9505.54503</td><td>11.9</td><td>1.4</td><td>8.7</td><td>2.5</td><td>28.2</td><td>4.9</td><td>101.5</td><td>9.3</td><td>0.06</td><td>0.02</td><td>0.22</td><td>0.05</td></tr><tr><td>9506.48382</td><td>9.1</td><td>1.3</td><td>11.9</td><td>1.9</td><td>31.3</td><td>4.6</td><td>104.3</td><td>7.7</td><td>0.08</td><td>0.02</td><td>0.31</td><td>0.06</td></tr><tr><td>9510.56919</td><td>9.0</td><td>0.4</td><td>11.0</td><td>1.3</td><td>36.3</td><td>1.4</td><td>120.7</td><td>5.4</td><td>0.29</td><td>0.02</td><td>0.51</td><td>0.06</td></tr><tr><td>9511.52849</td><td>10.0</td><td>0.8</td><td>14.5</td><td>1.8</td><td>30.9</td><td>3.0</td><td>97.3</td><td>7.0</td><td>0.11</td><td>0.02</td><td>0.27</td><td>0.05</td></tr><tr><td>9512.51624</td><td>6.7</td><td>1.8</td><td>7.6</td><td>4.5</td><td>38.6</td><td>7.7</td><td>99.3</td><td>21.6</td><td>0.08</td><td>0.04</td><td>0.14</td><td>0.09</td></tr><tr><td>9513.62825</td><td>4.6</td><td>0.8</td><td>20.0</td><td>4.6</td><td>38.3</td><td>2.6</td><td>89.3</td><td>9.2</td><td>0.21</td><td>0.03</td><td>0.18</td><td>0.06</td></tr><tr><td>9514.59185</td><td>9.5</td><td>0.6</td><td>20.0</td><td>1.8</td><td>33.6</td><td>2.3</td><td>90.1</td><td>5.0</td><td>0.19</td><td>0.03</td><td>0.36</td><td>0.06</td></tr><tr><td>9515.61072</td><td>7.7</td><td>0.7</td><td>20.0</td><td>2.0</td><td>31.4</td><td>2.5</td><td>84.4</td><td>4.9</td><td>0.15</td><td>0.03</td><td>0.32</td><td>0.05</td></tr><tr><td>9516.53216</td><td>8.9</td><td>1.0</td><td>15.5</td><td>2.0</td><td>29.4</td><td>3.3</td><td>102.5</td><td>7.1</td><td>0.10</td><td>0.02</td><td>0.32</td><td>0.05</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Radial velocity of the narrow component of the HeI 5876 Å line profile in the stellar rest frame plotted as a function of Julian date (top) and rotational phase computed from the ephemeris of Eq. (1), (bottom). The color code corresponds to successive rotational cycles. The fit by a geometrical accretion shock model (see text) is shown (dash-dotted curve) together with its $1\sigma$ uncertainty (gray area).</figcaption></figure>

the profile. At least two systems of blueshifted absorptions can be identified: a low-velocity system extending between -20 to -80 km s⁻¹, and a high-velocity system ranging from -100 to -250 km s⁻¹. Figures 10 and 11 show that the two systems are quite variable. These two absorption systems are reminiscent of those observed in the Hα and Hβ lines, although they occur at slightly bluer velocities in the HeI profile, with an offset of about -40 km s⁻¹ compared to the optical profiles.

Correlation matrices involving the HeI 10830 Å line are displayed in Fig. B.4. The HeI autocorrelation matrix shows several correlated components. Over the velocity channels extending from -100 km s⁻¹ up to 200 km s⁻¹, the line flux varies in a correlated fashion, but this region of the line does not correlate with the rest of the profile. Similarly, the redshifted absorption region, which extends from 200 km s⁻¹ to 400 km s⁻¹, varies as a whole

and does not correlate with the rest of the profile. The blue wing of the profile presents many short-duration absorption components superimposed on the emission component, and, probably due to this, each velocity bin varies independently of the other. Comparison of the HeI 10830 Å line profile to optical and near-infrared hydrogen profiles indicates that the emission core and peak intensity are correlated, as are the high-velocity redshifted absorptions seen in the HeI, Paβ, and Brγ line profiles. However, the HeI blue wing shows a strong anticorrelation with the core of the Hβ profile.

The periodogram analysis of the line profiles is shown in Fig. 11. In the three lines, the IPC components appear to be periodically modulated at the stellar rotation period. As the SPIRou dataset extends over four months, this suggests that a quite stable structure, presumably the accretion funnel flow, gives rise to this spectral feature<fnref n="10" />. The redshifted absorption component is deeper when the emission line is more intense, which agrees with the assumption that the two components trace the densest part of the accretion column close to the star. The low-velocity red wing of the HeI line profile also shows a periodic modulation at the same period, which might trace the visibility of the accretion shock, as seen in the optical HeI line profile. The periodogram power in the blue wing of the lines is weaker, except perhaps over the high-velocity channels of the Paβ line. In particular, the variable blueshifted absorption systems seen in the HeI line profile are not periodically modulated on this long timescale. We performed a similar analysis of each SPIRou run individually, and the results are shown in Appendix E. Figures E.1-E.4 reveal that the most conspicuous high-velocity redshifted absorptions in the HeI, Paβ, and Brγ line profiles occur preferentially around $\Phi_{rot} = 0$. During the SPIRou October run, which was simultaneous with the OHP/SOPHIE observations, we realized that the highly blueshifted velocity channels of the near-infrared lines appear to be modulated at the stellar rotation period (see Fig. E.2). Hence, while the modulation of high-velocity blueshifted absorptions seen in the HeI profile disappear on a timescale of several months, the

<footnote n="10">¹⁰ The long time coverage of the SPIRou observations enabled us to explore periods up to 100 days. However, we did not find significant periods longer than those reported here.</footnote>
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Near-infrared line profiles: HeI (**left**), Paβ (**center**), and Bry (**right**), plotted with arbitrary offsets as a function of Julian date. Each color represents a SPIRou run, namely September (blue), October (orange), November (green), December 2021 (red), and January 2022 (purple). The October SPIRou run is contemporaneous to the OHP/SOPHIE observations.</figcaption></figure>

modulation may survive over a few rotational periods. Finally, similar to what is seen in the Balmer line profiles over a much shorter time span (see Sect. 3.2.3), there is a marginal indication from the periodograms of the near-infrared lines that the modulation period drifts from the blue to the red wing of the profiles. This might be a sign of differential rotation in the source of the variability.

The EW of the HeI 10830 Å, Paβ, and Bry lines was computed on the residual profiles using two methods. The first method consists of adjusting a Gaussian fit to the line profile, and the second method is integrating below the line profile. The two methods were applied to the Paβ and Bry lines that most often exhibit a strong Gaussian-like emission and, at times, a pronounced redshifted absorption. The difference between the
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Near-infrared line profile variability. Top: series of residual near-infrared line profiles HeI (left), Paβ (center), and Bry (right), plotted superimposed to illustrate their variability. Each color represents a SPIRou run as in Fig. 10. Bottom: 2D periodograms across the line profiles. The color code reflects the periodogram power from zero (blue) to 0.5 (red). The dotted red horizontal line drawn at a frequency of 0.166 day⁻¹ indicates the stellar rotational period. The white curve displays the mean line profile.</figcaption></figure>

EW derived from the Gaussian fit and from the profile integration then provides an estimate of the strength of the redshifted absorption component. The HeI line exhibits a complex profile, with pronounced absorptions appearing in both the blue and red wings, and it cannot be fit by a Gaussian. We report here HeI EWs measured through profile integration only. EW measurements are usually accurate to within 5% because of the well-defined adjacent stellar continuum level.

The results are listed in Table 7, and the evolution of EW during the observing campaign is illustrated on Fig. 12. While the EW variations of the three lines are usually correlated, there are notable exceptions. For instance, during the first SPIRou run around JD 2 459 478, the HeI line goes into absorption, driven by a broad redshifted absorption component that reaches half the continuum value, while the Paβ and Bry lines reach a local intensity maximum, even though they also exhibit a pronounced redshifted absorption (see Fig. E.1). In contrast, during the second SPIRou run, the variations of the three lines are well correlated. Line EWs and near-infrared veiling measurements (listed in Table 3) are both shown as a function of Julian date and rotational phase in Fig. D.1. Both quantities appear to be rotationally modulated, with maximum values occurring close to Φ<sub>rot</sub> = 0.

Finally, the photospheric radial velocity curve derived from SPIRou spectra is shown in Fig. 13. $V_r$ is found to vary between 13.96 and 15.04 km s⁻¹, with a median value of 14.65 km s⁻¹. A CLEAN periodogram analysis reveals a period $P = 5.94 \pm 0.11$ d, and the String-Length method yields $P = 5.98 \pm 0.12$ d, both consistent with the $6.04 \pm 0.15$ d photometric period derived in Sect. 3.1. The photospheric radial velocity curve folded in phase at the stellar rotation period is shown in Fig. 13. $V_r$ exhibits a roughly sinusoidal variations in rotational phase, which suggests that it is modulated by a surface spot. It reaches the median velocity going blueward around $\Phi_{\text{rot}}=0$, as expected from a stellar spot facing the observer at this phase. Because this is also the phase of maximum brightness of the system (see Sect. 3.1), this suggests that a hot spot modulates the photospheric $V_r$ curve. The $V_r$ curve derived from the near-infrared photospheric lines (see Fig. 13) is inverted

compared to the $V_r$ curve derived for the HeI 5876 Å NC line profile (see Fig. 9). Similar antiphase radial velocity variations between absorption and emission lines have been reported for other accreting T Tauri stars and were interpreted as caused by the modulation of the line profiles by a hot spot at the stellar surface (Petrov et al. 2001, 2011; Gahm et al. 2013). Alternatively, the photospheric radial velocity variations may also be produced by a cool spot that coexists with the accretion shock at the same location on the stellar surface, as Doppler images of accreting T Tauri stars suggest (e.g., Donati et al. 2010, 2011, 2019, 2020a).

### 3.4. Low-resolution near-infrared spectroscopy

The median nightly spectra obtained from the ExTrA-T2 and ExTrA-T3 telescopes are shown in Fig. 14. We computed the EW of the HeI 10830 Å, Paβ, Pay, and Paδ lines from the ExTrA spectra using specutils<fnref n="11" />. We analyzed each telescope independently because the point spread function (and therefore the resolution) of the spectrograph depends on the position on the detector. First, we fit a Gaussian to each line on the median spectra to measure the line center and full width, the latter we defined as amounting to six times the standard deviation of the Gaussian fit in order to isolate the line profile from the nearby continuum. The local continuum around each line was modeled with a third-degree polynomial, adjusted on three line widths centered on the line, but excluding the line. Then, for each line in each of the 1898 individual spectra, we computed the EW by integrating the flux over the full line width using the parameters and the normalization region defined from the median spectra. Finally, we computed the median of the individual measurements for each night, regardless of the telescope. The median is less affected by outliers than the mean, and the differences between the mean and the median values are within the error bars. Table C.1 lists the results.

In order to estimate the reliability of the procedure, we compared EWs derived from ExTrA and from SPIRou spectra for the Paβ line using 20 measurements obtained on each

<footnote n="11">¹¹ https://github.com/astropy/specutils</footnote>
<table><caption>**Table 7.** Near-infrared line EW measurements from the CFHT/SPIRou spectra.</caption><thead><tr><th rowspan="2">Julian date<br>(2 450 000+)</th><th colspan="6">EW (Å)</th></tr><tr><th>HeI<sub>int</sub></th><th>Paβ<sub>g</sub></th><th>Paβ<sub>int</sub></th><th>Paβ<sub>diff</sub></th><th>Bryγ<sub>g</sub></th><th>Bryγ<sub>int</sub></th><th>Bryγ<sub>diff</sub></th></tr></thead><tbody><tr><td>9473.066</td><td>2.88</td><td>8.29</td><td>7.41</td><td>-0.87</td><td>3.26</td><td>2.94</td><td>-0.32</td></tr><tr><td>9475.043</td><td>0.98</td><td>7.58</td><td>7.40</td><td>-0.18</td><td>3.07</td><td>2.73</td><td>-0.34</td></tr><tr><td>9476.969</td><td>-0.76</td><td>6.05</td><td>5.82</td><td>-0.23</td><td>2.60</td><td>1.95</td><td>-0.65</td></tr><tr><td>9478.027</td><td>-2.45</td><td>11.58</td><td>9.64</td><td>-1.93</td><td>6.97</td><td>5.70</td><td>-1.26</td></tr><tr><td>9480.082</td><td>2.57</td><td>10.74</td><td>9.76</td><td>-0.97</td><td>4.56</td><td>4.05</td><td>-0.51</td></tr><tr><td>9481.086</td><td>1.86</td><td>9.90</td><td>9.65</td><td>-0.25</td><td>4.15</td><td>3.80</td><td>-0.35</td></tr><tr><td>9482.086</td><td>0.98</td><td>5.18</td><td>5.53</td><td>0.35</td><td>1.72</td><td>1.47</td><td>-0.25</td></tr><tr><td>9502.074</td><td>3.33</td><td>10.82</td><td>10.02</td><td>-0.80</td><td>5.37</td><td>4.73</td><td>-0.63</td></tr><tr><td>9503.074</td><td>2.43</td><td>9.02</td><td>8.12</td><td>-0.89</td><td>3.02</td><td>2.17</td><td>-0.85</td></tr><tr><td>9504.098</td><td>2.58</td><td>12.18</td><td>11.48</td><td>-0.70</td><td>5.09</td><td>4.63</td><td>-0.45</td></tr><tr><td>9506.082</td><td>2.19</td><td>11.36</td><td>11.80</td><td>0.44</td><td>6.54</td><td>6.29</td><td>-0.24</td></tr><tr><td>9508.082</td><td>7.33</td><td>16.93</td><td>16.12</td><td>-0.81</td><td>9.00</td><td>8.48</td><td>-0.51</td></tr><tr><td>9509.086</td><td>9.45</td><td>19.84</td><td>19.11</td><td>-0.73</td><td>9.68</td><td>9.16</td><td>-0.51</td></tr><tr><td>9510.086</td><td>10.53</td><td>17.07</td><td>16.80</td><td>-0.27</td><td>8.27</td><td>7.57</td><td>-0.70</td></tr><tr><td>9511.074</td><td>4.73</td><td>14.56</td><td>14.20</td><td>-0.36</td><td>7.48</td><td>7.48</td><td>0.00</td></tr><tr><td>9513.090</td><td>0.73</td><td>8.00</td><td>8.32</td><td>0.31</td><td>4.68</td><td>4.82</td><td>0.13</td></tr><tr><td>9514.051</td><td>4.82</td><td>13.19</td><td>12.92</td><td>-0.27</td><td>6.88</td><td>6.37</td><td>-0.51</td></tr><tr><td>9515.078</td><td>1.03</td><td>12.18</td><td>11.88</td><td>-0.30</td><td>5.81</td><td>5.71</td><td>-0.09</td></tr><tr><td>9516.102</td><td>0.51</td><td>11.26</td><td>10.78</td><td>-0.48</td><td>5.02</td><td>4.62</td><td>-0.39</td></tr><tr><td>9535.102</td><td>5.37</td><td>11.00</td><td>10.67</td><td>-0.33</td><td>6.03</td><td>5.83</td><td>-0.20</td></tr><tr><td>9537.098</td><td>0.96</td><td>9.65</td><td>9.96</td><td>0.30</td><td>4.93</td><td>5.04</td><td>0.11</td></tr><tr><td>9538.039</td><td>2.84</td><td>6.93</td><td>6.13</td><td>-0.80</td><td>2.74</td><td>2.31</td><td>-0.42</td></tr><tr><td>9538.984</td><td>4.10</td><td>7.99</td><td>6.49</td><td>-1.50</td><td>3.12</td><td>2.23</td><td>-0.88</td></tr><tr><td>9539.988</td><td>3.45</td><td>9.12</td><td>8.52</td><td>-0.59</td><td>3.08</td><td>2.63</td><td>-0.45</td></tr><tr><td>9541.012</td><td>5.32</td><td>12.48</td><td>12.53</td><td>0.05</td><td>5.46</td><td>5.46</td><td>0.00</td></tr><tr><td>9557.980</td><td>3.13</td><td>9.12</td><td>8.55</td><td>-0.57</td><td>3.96</td><td>3.42</td><td>-0.53</td></tr><tr><td>9558.949</td><td>4.88</td><td>13.61</td><td>13.41</td><td>-0.20</td><td>6.25</td><td>6.14</td><td>-0.10</td></tr><tr><td>9559.961</td><td>1.81</td><td>10.18</td><td>10.16</td><td>-0.02</td><td>5.77</td><td>5.63</td><td>-0.14</td></tr><tr><td>9560.965</td><td>1.29</td><td>6.34</td><td>6.67</td><td>0.33</td><td>3.85</td><td>3.82</td><td>-0.02</td></tr><tr><td>9563.078</td><td>7.11</td><td>12.10</td><td>11.03</td><td>-1.07</td><td>5.48</td><td>4.56</td><td>-0.92</td></tr><tr><td>9564.059</td><td>7.27</td><td>10.85</td><td>10.41</td><td>-0.44</td><td>5.09</td><td>4.72</td><td>-0.37</td></tr><tr><td>9566.051</td><td>2.98</td><td>6.62</td><td>6.78</td><td>0.15</td><td>2.62</td><td>2.44</td><td>-0.17</td></tr><tr><td>9566.953</td><td>1.05</td><td>7.17</td><td>7.57</td><td>0.39</td><td>4.39</td><td>4.47</td><td>0.08</td></tr><tr><td>9585.941</td><td>6.31</td><td>13.67</td><td>13.56</td><td>-0.11</td><td>6.47</td><td>6.50</td><td>0.02</td></tr></tbody></table>

**Notes.** EW$_g$ is obtained from a Gaussian fit to the line profile, while EW$_{int}$ is measured from line profile integration. EW$_{diff}$ is the difference between EW$_{int}$ and EW$_g$.

instrument less than one day apart. The results show an excellent correlation, with a slight tendency for the EW measured from ExTrA to exceed those measured from SPIRou, with a mean difference of $0.7 \pm 1.1$ Å. A similar result is obtained for the HeI line, with a mean difference of $0.31$ Å and an rms of $1.1$ Å between ExTrA and SPIRou estimates. The comparison with high-resolution measurements thus validates the EWs obtained from low-resolution spectra.

The line variability is illustrated on Fig. 14, where the median nightly spectra are superimposed. The median and extreme values of EWs measured for the HeI, Paβ, Paγ, and Paδ lines are listed in Table 8, and the night-by-night measurements are listed in Table C.1. We note that the HeI line at times appears to be in absorption at this low spectral resolution.

Figure 15 shows the EW measurements plotted as a function of time and its generalized Lomb-Scargle periodogram (GLS; Zechmeister & Kürster 2009). The EW of the four lines is found to be modulated with a period of $6.028 \pm 0.087$ days, consistent with the stellar rotation period, where the error estimate is the standard deviation of a Gaussian fit to the periodogram peak. Figure 16 shows the HeI and Paβ line EWs folded in phase at the stellar rotation period. The modulation of the line strength clearly appears, with maximum flux around phase zero. As shown in the same figure, these variations follow the modulated brightness level of the system in the J band.

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

**Fig. 12.** EW of near-infrared lines. Top: EW of the HeI, Paβ, and Bry lines measured on SPIRou spectra plotted as a function of Julian date. The measurement uncertainties are about the symbol size. Bottom: difference between the EWs measured by direct line profile integration (i.e., including the redshifted absorption component) and the EW derived from the Gaussian fitting of the emission component only plotted as a function of Julian date. This differential quantity measures the strength of the redshifted absorption component in the Paβ and Bry line profiles. The more negative the differential quantity, the deeper the redshifted absorption.

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

**Fig. 13.** Radial velocity variations measured in the near-infrared photospheric lines of the SPIRou spectra. Top: radial velocity as a function of Julian date. Bottom: radial velocity curve folded in phase at the stellar rotational period $P = 6.04$ days. The color code reflects the Julian date.

Longer-term EW variations of higher amplitudes are also clearly seen in Fig. 15. These variations seem to be correlated with the multicolor photometry presented in Sect. 3.1, in particular, during the brightness event centered on JD 2 459 509 and the wide dip around JD 2 459 549. In Fig. 17, a clear correlation appears between the photometric variations in the u' band and the Paβ line variations. If the brightening of the system in the
<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Median spectrum for each night from the ExTrA-T2 telescope (lower part, 88 spectra) and ExTrA-T3 telescope (upper part, 69 spectra). The main emission lines are indicated. The color is proportional to the EW of the Paβ line shown in the top panel as a function of Julian date.</figcaption></figure>

<table><caption>**Table 8.** Minimum, median, and maximum near-infrared line EW measurements from the ExTrA spectra.</caption><thead><tr><th rowspan="2"></th><th colspan="4">EW (Å)</th></tr><tr><th>HeI</th><th>Paβ</th><th>Paγ</th><th>Paδ</th></tr></thead><tbody><tr><th>Min.</th><td>-1.9</td><td>2.8</td><td>3.4</td><td>0.2</td></tr><tr><th>Med.</th><td>2.8</td><td>10.5</td><td>7.0</td><td>2.1</td></tr><tr><th>Max.</th><td>11.8</td><td>18.6</td><td>11.4</td><td>4.8</td></tr></tbody></table>

u' band is linked to accretion, this correlation suggests that most of the emission line flux is connected to the same process.

### 3.5. Mass-accretion rate

The combination of optical and near-infrared spectroscopy offers a number of emission lines from which we can estimate the mass accretion rate onto the star, using the empirical relations between line luminosity and accretion luminosity proposed by Alcalá et al. (2017). Combining the range of Hα and Hβ EWs reported above with the nearby continuum fluxes computed from the r' and g' magnitudes corrected for extinction, we obtain the line fluxes and luminosities as follows: $F_{line} = F_{\lambda}^{o} \times EW(line) \times 10^{-0.4(m_{\lambda}-A_{\lambda})}$ and $L_{line} = 4\pi d^2 F_{line}$, where $F_{\lambda}^{o}$ is the flux of a zero-magnitude star<fnref n="12" />, $m_{\lambda}$ and $A_{\lambda}$ are the magnitude and extinction in the photometric band of interest, and d is the distance to the star. The accretion luminosity was derived from the line luminosity using the relation reported by Alcalá et al. (2017), and $\dot{M}_{acc}$ was deduced from $L_{acc}$ assuming a magnetospheric radius of $5R_{\star}$

(see Alcalá et al. 2017, Eq.(1)). Taking the uncertainties on all involved quantities into account, we thus derive $\dot{M}_{acc} = 0.7 \pm 0.3$ and $0.5 \pm 0.4 \times 10^{-8} M_{\odot} \text{yr}^{-1}$ from the Hα and Hβ line fluxes, respectively.

We performed a similar analysis using the extensive measurements of Paβ EWs obtained from the ExTrA spectra. From the extreme values, namely $EW(Paβ) = 2.8-18.6$ Å, and the mean REM J-band magnitude of the system during the observing period, $J = 9.41 \pm 0.10$, we derive $\dot{M}_{acc} = 0.3-2.0 \times 10^{-8} M_{\odot} \text{yr}^{-1}$, with a median value of $\dot{M}_{acc} = 1.0 \times 10^{-8} M_{\odot} \text{yr}^{-1}$. Finally, we derived additional estimates of $\dot{M}_{acc}$ from the Bry line flux, using the median and extreme values of EW(Bry) measured on SPIRou spectra, namely 5.1, 1.7, and 9.7 Å, which we calibrated with the mean REM K'-band magnitude of the system, $K' = 8.40$. Using the relations reported by Alcalá et al. (2017) between line and accretion luminosity, we thus derived $\dot{M}_{acc} = 0.2-1.8 \times 10^{-8} M_{\odot} \text{yr}^{-1}$, with a median value of $\dot{M}_{acc} = 0.8 \times 10^{-8} M_{\odot} \text{yr}^{-1}$.

The dispersion in the $\dot{M}_{acc}$ estimates partly results from intrinsic $\dot{M}_{acc}$ variability. For instance, the SOPHIE spectra were obtained during a brightening of the system that occurred around JD 2459512, at a time at which all optical and infrared lines were relatively strong. We thus derive a relatively high value of $\dot{M}_{acc}$ from these spectra compared, for example, to the smaller $\dot{M}_{acc}$ estimate derived from the Bry line in the SPIRou spectra that were obtained during more quiescent phases of the system. We cannot exclude, however, that some of the dispersion may also arise from systematic uncertainties in the empirical line-to-accretion luminosity relations over the optical and near-infrared wavelength ranges. In any case, the various estimates agree globally, with $\dot{M}_{acc}$ typically varying between 0.2 and $2.0 \times 10^{-8} M_{\odot} \text{yr}^{-1}$.

### 4. Discussion

The monitoring campaign we performed on GM Aur reveals significant but relatively low-level temporal variability over a timescale of six months. The photometric variations are mild, as might be expected for this moderately accreting young system ($\dot{M}_{acc} \sim 0.8 \times 10^{-8} M_{\odot} \text{yr}^{-1}$), with amplitudes ranging from 1.5 mag in the u' band to 0.3 mag in the i' band, and 0.1 mag in the J band. The brightness of the system varies smoothly because it is modulated by the visibility of surface spots at the stellar rotational period of 6.04 d. Except for the HeI 10830 Å line profile, the spectral appearance of the system does not change drastically on a timescale of months, and the veiling is low and stable, amounting to about 0.3 in the optical range. The strength of the main emission lines (Hα, Hβ, HeI 5876 Å, Paβ, and Bry) varies by a factor of 2–3 over the course of the semester. In contrast, the HeI 10830 Å line profile exhibits extreme variability and is sometimes barely noticeable in emission. The line profile shapes are strongly variable on a timescale of days, and the development of both blueshifted and redshifted absorption components is superimposed onto a broad emission component. Blueshifted absorption components indicate outflows (e.g., Edwards et al. 2003; Kwan et al. 2007), and redshifted components probe funnel flows (e.g., Edwards et al. 2006; Fischer et al. 2008). Remarkably, the redshifted absorption features that reach below the continuum level, that is, inverse P Cygni profiles (IPC; see Calvet & Hartmann 1992), which are seen in the high-resolution near infrared line profiles (HeI 10830 Å, Paβ, Bry), are steadily modulated by stellar rotation over an extended observing time span of 3 months. Rotational modulation is also

<footnote n="12">¹² $F_{\lambda}^{o} = 2.43 \times 10^{-9}$ and $5.27 \times 10^{-9}$ erg s⁻¹ cm⁻² Å⁻¹ in the r' and g' bands, respectively.</footnote>
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** EW variability of the near-infrared line profiles. **Left:** EW measurements plotted as a function of Julian date for the HeI, Paβ, Pay, and Paδ lines from the ExTrA spectra. **Right:** GLS periodogram of the EW measurements. The period on the x-axis is displayed on a log scale.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Fig. 16.** EWs of near-infrared line profiles and J-band magnitude variability. **Top:** EW of the near-infrared HeI and Paβ lines folded in phase at the stellar rotational period. EW(HeI) is offset by -10 Å for clarity. **Bottom:** J-band light curve, deduced from ExTrA spectra, folded in phase at the stellar period. The brightness modulation is similar to that observed at optical wavelengths (see Fig. 2).</figcaption></figure>

clearly detected in the strength of the near-infrared lines and is continuously observed at low resolution for more than five months. Assuming that most of the line flux arises from the magnetospheric accretion region, as suggested by the periodic appearance of the IPCs and the correlation between line flux and *u'*-band excess, this indicates a stable large-scale accretion structure on this timescale.

In contrast, high-velocity blueshifted absorption components are neither periodic nor stable on this timescale. While they are ubiquitous in the optical line profiles, most notably Hα and Hβ, and are also quite conspicuous in the HeI 10830 Å line profile, their signature evolves on a timescale of a few days, sometimes drifting in velocity before disappearing altogether. As an example, Fig. E.3 shows the evolution of the blueshifted absorption

<figure><img src="image_18.png" /><figcaption>**Fig. 17.** Mid-term variation of EW(Paβ), (red) compared to the system *u'*-band light curve (blue) for measurements taken less than one day apart. To facilitate the comparison, the EW measurements are plotted on a magnitude scale and are offset, namely -2.5 log EW(Paβ) + 5.</figcaption></figure>

components in the HeI line profile over the course of the SPIRou November run. On JD 9537, a high-velocity blueshifted feature appears in the profile and drifts toward lower velocities over the next several days, from -200 km s⁻¹ to -120 km s⁻¹. Soon after, on JD 9540, a new high-velocity component appears and follows the same trend. A similar behavior is seen in the high-velocity blueshifted component appearing in the HeI 10830 Å line profile during the SPIRou September run, which drifted from -240 km s⁻¹ to -150 km s⁻¹ on a timescale of five days, from JD 9475 to JD 9480 (see Fig. E.1). This suggests episodic outflows lasting for a few days only. We see no sign of a steady, constant velocity wind in the HeI line profiles over the observing period, nor do we find evidence for a rotational modulation of the outflow signatures. The only stable outflow signature seen in the optical and near-infrared line profiles consists of a narrow, low-velocity blueshifted absorption in the Hα profile that peaks at -20 km s⁻¹ and remains visible over more than two weeks.

Observations thus suggest that we witness a globally stable accretion structure and a succession of short-lived episodic outflows. The contrasting behavior of accretion and outflow diagnostics observed on a timescale of months thus raises the question whether the two processes are physically connected on the scale of a few stellar radii that we probe here. To examine this issue, we investigated the aftermath of the brightening event GM Aur underwent around JD 2 459 509. On this date, the system exhibited a significant brightening at optical and near-infrared wavelengths, most notably in the *u'* band (~1 mag), as well as some of the strongest line fluxes and highest optical and
<figure><img src="image_19.png" /><figcaption>**Fig. 18.** Magnetospheric ejection model. Development of a magnetospheric ejection computed from star-disk interaction MHD simulations. The snapshots shown here are extracted from model 3 of Pantolmos et al. (2020), where the magnetospheric truncation radius amounts to 54% of the corotation radius. The three snapshots are shown at 0.2, 0.5, 0.8, and 1.1*P*<sub>rot</sub>. The white curves indicate expanding magnetic field lines that give rise to the ejection of plasmons. The color scale indicates density. The green arrows show the velocity field of the stellar and disk winds and of MEs at their interface.</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>**Fig. 19.** Magnetospheric ejection model. Velocity of the gas in the ejected plasmoid at a distance of 10 (magenta), 20 (blue), and 30 (cyan) stellar radii as a function of time in units of the stellar rotational period. Successive ejections of plasmons are featured. The vertical lines correspond to the snapshots shown in Fig. 18.</figcaption></figure>

near-infrared veiling values measured during the campaign. A simultaneous TESS light curve of the system recorded the brightening event (see Fig. 2), which started on JD 2 459 508 and ended on JD 2 459 511, and exhibited a flat peak lasting for two days with a 20% flux increase. We therefore interpret this episode as an accretion burst that occurred around the rotational phase $\Phi_{rot} = 0$, corresponding to the maximum visibility of the accretion shock, and lasted for several days. From the measured continuum level and Balmer line EW during the burst, we derive an increase of a factor of 2 in the accretion rate, reaching $\sim 2 \times 10^{-8} M_{\odot} \text{yr}^{-1}$. Inspection of the optical and near-infrared line profiles during and after this event reveals high-velocity blueshifted absorption components that appear in the Balmer and HeI 10830 Å line profiles. On JD 2 459 509 a new blueshifted absorption component appears in the HeI infrared line at a velocity of $-280 \text{ km s}^{-1}$ and extends down to $-350 \text{ km s}^{-1}$ before it reaches the continuum level again, and it drifts to lower velocities ($\sim -150 \text{ km s}^{-1}$) over the next few days. The closest opti-

cal spectrum was recorded only toward the end of the accretion burst, on JD 2 459 510.6. In this spectrum, the $H\alpha$ and $H\beta$ profiles feature a new blueshifted absorption component at a velocity of $-110 \text{ km s}^{-1}$ that lasts for a few days. The contemporaneous occurrence of an accretion burst rapidly followed by outflow signatures in the line profiles therefore suggests that the accretion and ejection processes are physically connected on small scales.

We propose that the most likely scenario that accounts for these episodic events is magnetospheric ejections, possibly triggered by magnetic reconnections in the accreting magnetosphere. Magnetospheric ejections (Zanni & Ferreira 2013; Sauty et al. 2022), or nonstationary conical winds (Romanova et al. 2009), are caused by the expansion and reconnection of the field lines that connect the star with the disk. The inflation process is the result of the star-disk differential rotation and the consequent build-up of toroidal magnetic field pressure (e.g., Goodson et al. 1997). Quasi-periodic ejections of plasmons are predicted to occur throughout the magnetospheric inflation cycle on a timescale of about the rotational period (Hayashi et al. 1996). The speed and variability of the outflows likely depend on various parameters, such as the magnetic field strength and topology, the thermal disk pressure, nonideal MHD effects, and the interaction of the magnetospheric-ejection region with the surrounding outflows, that is, stellar and disk winds (Miller & Stone 1997; Romanova et al. 2009; Zanni & Ferreira 2013). Previous monitoring campaigns on young stellar objects have reported evidence for magnetospheric inflation cycles (Bouvier et al. 2003; Alencar et al. 2018). We show a 2.5D MHD simulation of the interaction between an inner accretion disk and a dipolar magnetosphere from Pantolmos et al. (2020) in Fig. 18. The figure illustrates the ejection of plasmons along expanding field lines at the distance of a few stellar radii from the stellar surface. The timescale for successive ejections is about that of the stellar rotation period. The outflow speed, shown in Fig. 19, reaches more than $200 \text{ km s}^{-1}$ in the early phases of the ejection, then decelerates to about $150 \text{ km s}^{-1}$ on a timescale of days ($0.3 \times P_{rot}$), and finally vanishes altogether. A direct comparison of the model to observations is not straightforward. The
terminal speed we derive from observations is higher than predicted by the simulation, and its temporal evolution on a timescale of a few days might be dominated by projection effects as the system rotates and does not trace the evolution of the plasmoid velocity field. Moreover, the sporadic ejections we observe do not appear to have the quasi-periodic character of the model. Both the wind speed and ejection timescale may, however, depend on numerical effects. In any case, the behavior of the magnetospheric ejection model qualitatively matches the dynamics of the high-velocity blueshifted absorptions seen in the line profiles of GM Aur on a timescale of days to weeks. Moreover, in the scenario of a magnetospheric inflation cycle, magnetic reconnection leads to an accretion burst and simultaneously triggers an ejection episode (e.g., Goodson & Winglee 1999). This is quite reminiscent indeed of what is suggested by the variability of GM Aur.

We also noted a change in the GM Aur light curve whose first part is dominated by successive low-level brightening events, up to the major accretion burst described above, while it exhibits luminosity dips toward the end of the observations. It is unclear whether the contrasting behavior of the system luminosity observed before and after the main accretion burst is a consequence of the burst itself, perhaps inducing a structural change in the star-disk interaction process. It is conceivable that the rearrangement of the magnetic topology after the inflation or reconnection event that led to the burst has impacted the large-scale geometry of the star-disk interaction region. Either a modest increase in the inner disk scale-height (Nagel 2017) or a reduction of the extent of the truncation radius is prone to trigger a periodic obscuration of the central star by circumstellar dust, that is, a dipper phenomenon (Cody et al. 2014), especially in young systems seen at high inclination (McGinnis et al. 2015; Bodman et al. 2017). The magnetic topology and the mass accretion rate may thus have slightly evolved over the six-month time span of the campaign, as suggested by the long-term variations of the emission line EWs shown in Fig. 15. However, the results we report here clearly indicate that the large-scale geometry of the star-disk interaction was not drastically modified over this timescale. This is evidenced by the phase stability of the modulated light curve, the smooth variability of the emission line profiles around their mean shape, and the strictly periodic appearance of IPC profiles over at least three months. All these accretion diagnostics support a globally stable magnetospheric accretion structure during the campaign.

Finally, it is interesting to compare the line profile shape and variability reported here to those reported by McGinnis et al. (2020), which were obtained in 2011, that is, ten years prior to our observing campaign. The shapes of the $H\alpha$ and $H\beta$ profiles are quite different in the two studies. In McGinnis et al. (2020), these are pure emission profiles with a triangular shape, without any significant absorption components, and they exhibit little variability over the timescale of a week (see their Figs. S5 and S6). Here, the same profiles appear to be much more structured, with highly variable absorption components on the same timescale. The HeI 5876 Å line profile variability also differs between the two studies, but in the opposite direction. In both studies, the line profile consists of a broad and a narrow component. However, in McGinnis et al. (2020), the intensity of the broad component clearly varies, especially on the blue wing, while we found it to be quite stable here. It seems that the behavior of the system was different between the two epochs. The $H\alpha$ and $H\beta$ line EWs are indeed systematically higher in McGinnis et al. (2020) and were comparable to the highest values we measured here during the JD 2 459 509 accretion burst.

It is therefore likely that GM Aur was in a state of more active accretion during the 2011 observations. This is consistent with the more triangular shape and lack of structure of the Balmer emission line profiles, which are predicted to become more optical thick as the funnel flow density increases (Muzerolle et al. 2001). It is also consistent with the higher level of variability seen in the blue wing of the broad component of the HeI 5876 Å line profile that may betray the existence of a hot accretion-driven wind at times of enhanced accretion (Beristain et al. 2001). Long-term changes in the system behavior driven by a varying mass accretion rate and/or a change in the magnetic topology are therefore likely to occur. Over the six-month span of our observing campaign, the significant variation observed in the EW of the $Pa\beta$ line profile, which by the end of the campaign reaches similar levels to those measured during the JD 2 459 509 accretion burst (see Fig. 15), suggests that such changes may occur on a timescale of a few weeks to months.

## 5. Conclusion

By combining optical and near-infrared high-resolution spectroscopic time series, seconded by a long-term monitoring of the photometric variability of the system and low-resolution near-infrared spectrophotometry, we were able to characterize the accretion and ejection process occurring in the young system GM Aur on a timescale ranging from days to months. We report a stable accretion pattern according to which the large-scale magnetic field of the star controls the accretion of gas from the inner disk onto the central star along funnel flows. The appearance of inverse P Cygni profiles that signal the crossing of funnel flows on the line of sight is remarkably periodic at the stellar rotation period of 6.04 days. Similarly, the photometric and line flux variations, both driven by the visibility of the accretion shock located at the foot of the main accretion column, are modulated at the same period. While the amplitude varies, the phase of variability of all these accretion diagnostics remains stable over the 30 rotational periods covered by the campaign. This suggests that the underlying magnetic topology that controls the non-axisymmetric accretion flow, presumably an inclined dipole on the large scale, did not undergo major changes over a timescale of six months. In stark contrast, high-velocity blueshifted absorption components that indicate outflows appear at random times in the emission line profiles. They are not rotationally modulated, and their signatures last for a few days only. We argue that these transient outflows associated with a stable accretion pattern are best accounted for by magnetospheric ejection models, as predicted by MHD simulations. Thus, by probing the dynamics of the star-disk interaction region, these results show that the physical connection between accretion and ejection processes that has long been established on large scales also appears to be valid on the much smaller sub-au scales.

**Acknowledgements.** We thank the referee for a prompt and detailed report. This study is based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. The observations at the CFHT were performed with care and respect from the summit of Maunakea which is a significant cultural and historic site; based on observations made at Observatoire de Haute Provence (CNRS), France; based on data collected under the ExTrA project at the ESO La Silla Paranal Observatory. ExTrA is a project of Institut de Planétologie et d'Astrophysique de Grenoble (IPAG/CNRS/UGA), funded by the European Research Council under the ERC Grant Agreement n. 337591-ExTrA. We thank Ágnes Kóspál for providing a reduced TESS light curve of GM Aur. Funding for the TESS mission is provided by NASA's Science Mission directorate. This project has received
funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement no. 742095; SPIDI: Star-Planets-Inner Disk-Interactions, http://www.spidi-eu.org). We acknowledge funding from the French National Research Agency (ANR) under contract number ANR-18-CE31-0019 (SPlaSH). SHPA acknowledges financial support from CNPq, CAPES and Fapemig. JFD acknowledges funding from the European Research Council (ERC) under the H2020 research & innovation programme (grant agreement no. 740651 NewWorlds). AF acknowledges support by the PRIN-INAFF 2019 STRADE (Spectroscopically TRAcing the Disk dispersal Evolution) and by the Large Grant INAF YODA (YSOs Outflow, Disks and Accretion). JFG was supported by fundação para a Ciência e Tecnologia (FCT) through the research grants UIDB/04434/2020 and UIDP/04434/2020. This work benefited from discussions with the ODYSSEUS (HST AR-16129) and PENELLOPE teams. Some of the plots presented in this paper were built using TOPCAT (Taylor 2005).

## References

Akeson, R. L., Boden, A. F., Monnier, J. D., et al. 2005, ApJ, 635, 1173

Alcalá, J. M., Natta, A., Manara, C. F., et al. 2014, A&A, 561, A2

Alcalá, J. M., Manara, C. F., Natta, A., et al. 2017, A&A, 600, A20

Alencar, S. H. P., Bouvier, J., Donati, J. F., et al. 2018, A&A, 620, A195

Appenzeller, I., & Bertout, C. 2013, A&A, 558, A83

Artemenko, S. A., Grankin, K. N., & Petrov, P. P. 2012, Astron. Lett., 38, 783

Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42

Beristain, G., Edwards, S., & Kwan, J. 2001, ApJ, 551, 1037

Bessolaz, N., Zanni, C., Ferreira, J., Keppens, R., & Bouvier, J. 2008, A&A, 478, 155

Bodman, E. H. L., Quillen, A. C., Ansdell, M., et al. 2017, MNRAS, 470, 202

Bohn, A. J., Benisty, M., Perraut, K., et al. 2022, A&A, 658, A183

Boisse, I., Eggenberger, A., Santos, N. C., et al. 2010, A&A, 523, A88

Bonfils, X., Almenara, J. M., Jocou, L., et al. 2015, SPIE Conf. Ser., 9605, 96051L

Bouchy, F., Hébrard, G., Udry, S., et al. 2009, A&A, 505, 853

Bouvier, J., Grankin, K. N., Alencar, S. H. P., et al. 2003, A&A, 409, 169

Bouvier, J., Alecian, E., Alencar, S. H. P., et al. 2020a, A&A, 643, A99

Bouvier, J., Perraut, K., Le Bouquin, J. B., et al. 2020b, A&A, 636, A108

Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, PASP, 125, 1031

Calvet, N., & Hartmann, L. 1992, ApJ, 386, 239

Cody, A. M., Stauffer, J., Baglin, A., et al. 2014, AJ, 147, 82

Cook, N. J., Artigau, E., Doyon, R., et al. 2022, PASP, 134

Covey, K. R., Ivezić, Z., Schlegel, D., et al. 2007, AJ, 134, 2398

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR On-line Data Catalog: II/246

Donati, J. F., Skelly, M. B., Bouvier, J., et al. 2010, MNRAS, 409, 1347

Donati, J. F., Bouvier, J., Walter, F. M., et al. 2011, MNRAS, 412, 2454

Donati, J. F., Hébrard, E., Hussain, G. A. J., et al. 2015, MNRAS, 453, 3706

Donati, J. F., Bouvier, J., Alencar, S. H., et al. 2019, MNRAS, 483, L1

Donati, J. F., Bouvier, J., Alencar, S. H., et al. 2020a, MNRAS, 491, 5660

Donati, J. F., Kouach, D., Moutou, C., et al. 2020b, MNRAS, 498, 5684

Dworetsky, M. M. 1983, MNRAS, 203, 917

Edwards, S., Fischer, W., Kwan, J., Hillenbrand, L., & Dupree, A. K. 2003, ApJ, 599, L41

Edwards, S., Fischer, W., Hillenbrand, L., & Kwan, J. 2006, ApJ, 646, 319

Eisner, J. A., Hillenbrand, L. A., & Stone, J. M. 2014, MNRAS, 443, 1916

Espaillat, C., D'Alessio, P., Hernández, J., et al. 2010, ApJ, 717, 441

Espaillat, C. C., Macías, E., Hernández, J., & Robinson, C. 2019, ApJ, 877, L34

Espaillat, C. C., Robinson, C. E., Romanova, M. M., et al. 2021, Nature, 597, 41

Espaillat, C. C., Herczeg, G. J., Thanathibodee, T., et al. 2022, AJ, 163, 114

Fiorucci, M., & Munari, U. 2003, A&A, 401, 781

Fischer, W., Kwan, J., Edwards, S., & Hillenbrand, L. 2008, ApJ, 687, 1117

Flores, C., Connelley, M. S., Reipurth, B., & Duchêne, G. 2022, ApJ, 925, 21

Folsom, C. P., Bagnulo, S., Wade, G. A., et al. 2012, MNRAS, 422, 2072

Frasca, A., Alcalá, J. M., Covino, E., et al. 2003, A&A, 405, 149

Frasca, A., Guillout, P., Marilli, E., et al. 2006, A&A, 454, 301

Gahm, G. F., Stempels, H. C., Walter, F. M., Petrov, P. P., & Herczeg, G. J. 2013, A&A, 560, A57

Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1

Gangi, M., Antoniucci, S., Biazzo, K., et al. 2022, A&A, 667, A124

Giampapa, M. S., Basri, G. S., Johns, C. M., & Imhoff, C. 1993, ApJS, 89, 321

Goodson, A. P., & Winglee, R. M. 1999, ApJ, 524, 159

Goodson, A. P., Winglee, R. M., & Böhm, K.-H. 1997, ApJ, 489, 199

Grankin, K. N., Melnikov, S. Y., Bouvier, J., Herbst, W., & Shevchenko, V. S. 2007, A&A, 461, 183

Gravity Collaboration (Garcia Lopez, R., et al.) 2020, Nature, 584, 547

Gray, D. F. 1973, ApJ, 184, 461

Guilloteau, S., Simon, M., Piétu, V., et al. 2014, A&A, 567, A117

Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

Hartigan, P., Edwards, S., & Ghandour, L. 1995, ApJ, 452, 736

Hartmann, L., Herczeg, G., & Calvet, N. 2016, ARA&A, 54, 135

Hayashi, M. R., Shibata, K., & Matsumoto, R. 1996, ApJ, 468, L37

Herczeg, G. J., & Hillenbrand, L. A. 2014, ApJ, 786, 97

Huang, J., Andrews, S. M., Dullemond, C. P., et al. 2020, ApJ, 891, 48

Ingleby, L., Espaillat, C., Calvet, N., et al. 2015, ApJ, 805, 149

Ireland, L. G., Zanni, C., Matt, S. P., & Pantolmos, G. 2021, ApJ, 906, 4

Johns, C. M., & Basri, G. 1995, AJ, 109, 2800

Johns-Krull, C. M. 2007, ApJ, 664, 975

Kraus, A. L., & Hillenbrand, L. A. 2007, AJ, 134, 2340

Kwan, J., Edwards, S., & Fischer, W. 2007, ApJ, 657, 897

Landstreet, J. D. 1988, ApJ, 326, 967

Macías, E., Espaillat, C. C., Ribas, A., et al. 2018, ApJ, 865, 37

Manara, C. F., Frasca, A., Venuti, L., et al. 2021, A&A, 650, A196

Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74

Matt, S., & Pudritz, R. E. 2005, ApJ, 632, L135

McGinnis, P. T., Alencar, S. H. P., Guimarães, M. M., et al. 2015, A&A, 577, A11

McGinnis, P., Bouvier, J., & Gallet, F. 2020, MNRAS, 497, 2142

Miller, K. A., & Stone, J. M. 1997, ApJ, 489, 890

Muzerolle, J., Calvet, N., & Hartmann, L. 2001, ApJ, 550, 944

Nagel, E., & álvarez-Meraz, R., & Rendón, F., 2017, Rev. Mex. A&A, 53, 227

Oh, D., Hashimoto, J., Carson, J. C., et al. 2016, ApJ, 831, L7

Pantolmos, G., Zanni, C., & Bouvier, J. 2020, A&A, 643, A129

Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9

Percy, J. R., Grynkso, S., Seneviratne, R., & Herbst, W. 2010, PASP, 122, 753

Perruchot, S., Kohler, D., Bouchy, F., et al. 2008, SPIE Conf. Ser., 7014, 70140J

Petrov, P. P., Gahm, G. F., Gameiro, J. F., et al. 2001, A&A, 369, 993

Petrov, P. P., Gahm, G. F., Stempels, H. C., Walter, F. M., & Artemenko, S. A. 2011, A&A, 535, A6

Pouilly, K., Bouvier, J., Alecian, E., et al. 2020, A&A, 642, A99

Pouilly, K., Bouvier, J., Alecian, E., et al. 2021, A&A, 656, A50

Queloz, D., Henry, G. W., Sivan, J. P., et al. 2001, A&A, 379, 279

Rei, A. C. S., Petrov, P. P., & Gameiro, J. F. 2018, A&A, 610, A40

Roberts, D. H., Lehar, J., & Dreher, J. W. 1987, AJ, 93, 968

Robinson, C. E., & Espaillat, C. C. 2019, ApJ, 874, 129

Rodrigo, C., & Solano, E. 2020, in XIV.0 Scientific Meeting (virtual) of the Spanish Astronomical Society, 182

Rodrigo, C., Solano, E., & Bayo, A. 2012, IVOA Working Draft 15 October 2012

Roman-Duval, J., Proffitt, C. R., Taylor, J. M., et al. 2020, Res. Notes Am. Astron. Soc., 4, 205

Romanova, M. M., & Owocki, S. P. 2015, Space. Sec. Rev., 191, 339

Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & Lovelace, R. V. E. 2009, MNRAS, 399, 1802

Ryabchikova, T., Piskunov, N., Kurucz, R. L., et al. 2015, Phys. Scr., 90

Salyk, C., Blake, G. A., Boogert, A. C. A., & Brown, J. M. 2009, ApJ, 699, 330

Sauty, C., de Albuquerque, R. M. G., Cayatte, V., Lima, J. J. G., & Gameiro, J. F. 2022, A&A, 664, A176

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593

Simon, M., Guilloteau, S., Beck, T. L., et al. 2019, ApJ, 884, 42

Sousa, A. P., Bouvier, J., Alencar, S. H. P., et al. 2023, A&A, 670, A142

Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, M. F. 1989, AJ, 97, 1451

Symington, N. H., Harries, T. J., & Kurosawa, R. 2005, MNRAS, 356, 1489

Taylor, M. B. 2005, ASP Conf. Ser., 347, 29

Tessore, B., Soulain, A., Pantolmos, G., et al. 2023, A&A, 671, A129

Testa, V., Antonelli, L. A., Di Paola, A., et al. 2004, SPIE Conf. Ser., 5496, 729

Tonry, J. L., Denneau, L., Flewelling, H., et al. 2018, ApJ, 867, 105

Varga, J., Ábrahám, P., Chen, L., et al. 2018, A&A, 617, A83

Venuti, L., Bouvier, J., Irwin, J., et al. 2015, A&A, 581, A66

Wade, G. A., Bagnulo, S., Kochukhov, O., et al. 2001, A&A, 374, 265

Woitke, P., Kamp, I., Antonellini, S., et al. 2019, PASP, 131

Zanni, C., & Ferreira, J. 2013, A&A, 550, A99

Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577
## Appendix A: REM comparison stars and optical photometry

Table A.1 lists the comparison stars we used to calibrate the REM optical photometry (see Section 2.2). Table A.2 lists the derived g'r'i'z' magnitudes for GM Aur.

<table><caption>**Table A.1.** Literature data for comparison stars in the field of GM Aur on the REM cameras.</caption><thead><tr><th>Id</th><th>Name</th><th>2MASS</th><th>g</th><th>r</th><th>i</th><th>z</th><th>J</th><th>H</th><th>K'</th></tr><tr><td></td><td></td><td></td><td>(mag)</td><td>(mag)</td><td>(mag)</td><td>(mag)</td><td>(mag)</td><td>(mag)</td><td>(mag)</td></tr></thead><tbody><tr><td>*2</td><td>HD 282625</td><td>J04551650+3022369</td><td>11.310</td><td>10.914</td><td>10.772</td><td>10.707</td><td>9.781</td><td>9.554</td><td>9.487</td></tr><tr><td>*3</td><td>TIC 96533048</td><td>J04551015+3021333</td><td>12.361</td><td>11.399</td><td>10.936</td><td>10.668</td><td>9.295</td><td>8.716</td><td>8.563</td></tr><tr><td>*4</td><td>HD 282626</td><td>J04550536+3020382</td><td>11.747</td><td>11.375</td><td>11.238</td><td>11.174</td><td>10.239</td><td>10.052</td><td>9.987</td></tr><tr><td>*5</td><td>...</td><td>J04551400+3020168</td><td>13.279</td><td>12.637</td><td>12.312</td><td>12.109</td><td>10.986</td><td>10.601</td><td>10.501</td></tr><tr><td>*6</td><td>...</td><td>J04550078+3020090</td><td>13.351</td><td>12.288</td><td>11.604</td><td>11.135</td><td>9.564</td><td>8.828</td><td>8.608</td></tr><tr><td>*7</td><td>...</td><td>J04550253+3019228</td><td>14.529</td><td>13.494</td><td>12.892</td><td>12.496</td><td>10.994</td><td>10.264</td><td>10.086</td></tr></tbody></table>

Notes: griz magnitudes from Pan-STARRS (Tonry et al. 2018); JHK' magnitudes from 2MASS (Cutri et al. 2003).
<table><caption>**Table A.2.** REM g'r'i'z' photometry</caption><thead><tr><th>Julian date<br>(2,450,000+)</th><th>g'<br>(mag)</th><th>err<br>(mag)</th><th>r'<br>(mag)</th><th>err<br>(mag)</th><th>i'<br>(mag)</th><th>err<br>(mag)</th><th>z'<br>(mag)</th><th>err<br>(mag)</th></tr></thead><tbody><tr><td>9497.72986</td><td>12.525</td><td>0.045</td><td>11.566</td><td>0.008</td><td>11.161</td><td>0.006</td><td>10.861</td><td>0.012</td></tr><tr><td>9497.73190</td><td>12.559</td><td>0.042</td><td>11.587</td><td>0.011</td><td>11.171</td><td>0.007</td><td>10.872</td><td>0.017</td></tr><tr><td>9497.73396</td><td>12.566</td><td>0.051</td><td>11.601</td><td>0.019</td><td>11.193</td><td>0.004</td><td>10.892</td><td>0.007</td></tr><tr><td>9500.80099</td><td>12.67</td><td>0.016</td><td>11.589</td><td>0.005</td><td>11.174</td><td>0.016</td><td>10.827</td><td>0.004</td></tr><tr><td>9500.80302</td><td>12.667</td><td>0.02</td><td>11.592</td><td>0.001</td><td>11.175</td><td>0.013</td><td>10.838</td><td>0.006</td></tr><tr><td>9500.80505</td><td>12.682</td><td>0.01</td><td>11.608</td><td>0.007</td><td>11.182</td><td>0.019</td><td>10.853</td><td>0.005</td></tr><tr><td>9501.80935</td><td>12.617</td><td>0.02</td><td>11.58</td><td>0.01</td><td>11.168</td><td>0.01</td><td>10.871</td><td>0.007</td></tr><tr><td>9501.81138</td><td>12.626</td><td>0.024</td><td>11.582</td><td>0.01</td><td>11.169</td><td>0.007</td><td>10.871</td><td>0.006</td></tr><tr><td>9501.81344</td><td>12.635</td><td>0.019</td><td>11.585</td><td>0.013</td><td>11.173</td><td>0.008</td><td>10.868</td><td>0.004</td></tr><tr><td>9502.83317</td><td>12.5</td><td>0.022</td><td>11.544</td><td>0.003</td><td>11.144</td><td>0.012</td><td>10.825</td><td>0.01</td></tr><tr><td>9502.83523</td><td>12.52</td><td>0.024</td><td>11.548</td><td>0.003</td><td>11.149</td><td>0.01</td><td>10.834</td><td>0.011</td></tr><tr><td>9502.83729</td><td>12.519</td><td>0.021</td><td>11.55</td><td>0.004</td><td>11.149</td><td>0.012</td><td>10.843</td><td>0.011</td></tr><tr><td>9503.83744</td><td>12.494</td><td>0.035</td><td>11.522</td><td>0.007</td><td>11.129</td><td>0.014</td><td>10.812</td><td>0.009</td></tr><tr><td>9503.83947</td><td>12.494</td><td>0.032</td><td>11.515</td><td>0.004</td><td>11.127</td><td>0.013</td><td>10.802</td><td>0.009</td></tr><tr><td>9503.84153</td><td>12.505</td><td>0.033</td><td>11.521</td><td>0.004</td><td>11.121</td><td>0.013</td><td>10.816</td><td>0.01</td></tr><tr><td>9504.84169</td><td>12.75</td><td>0.028</td><td>11.656</td><td>0.008</td><td>11.223</td><td>0.004</td><td>10.902</td><td>0.02</td></tr><tr><td>9504.84375</td><td>12.749</td><td>0.029</td><td>11.653</td><td>0.008</td><td>11.222</td><td>0.007</td><td>10.897</td><td>0.019</td></tr><tr><td>9504.84578</td><td>12.75</td><td>0.027</td><td>11.654</td><td>0.009</td><td>11.223</td><td>0.006</td><td>10.89</td><td>0.02</td></tr><tr><td>9506.78944</td><td>12.81</td><td>0.041</td><td>11.67</td><td>0.01</td><td>11.239</td><td>0.009</td><td>10.907</td><td>0.006</td></tr><tr><td>9506.79150</td><td>12.812</td><td>0.031</td><td>11.671</td><td>0.006</td><td>11.243</td><td>0.012</td><td>10.906</td><td>0.01</td></tr><tr><td>9506.79353</td><td>12.797</td><td>0.034</td><td>11.679</td><td>0.008</td><td>11.233</td><td>0.011</td><td>10.886</td><td>0.009</td></tr><tr><td>9507.79360</td><td>12.503</td><td>0.045</td><td>11.53</td><td>0.014</td><td>11.134</td><td>0.008</td><td>10.814</td><td>0.01</td></tr><tr><td>9507.79566</td><td>12.495</td><td>0.043</td><td>11.534</td><td>0.011</td><td>11.131</td><td>0.009</td><td>10.822</td><td>0.015</td></tr><tr><td>9507.79772</td><td>12.477</td><td>0.044</td><td>11.519</td><td>0.013</td><td>11.127</td><td>0.009</td><td>10.82</td><td>0.012</td></tr><tr><td>9508.80257</td><td>12.135</td><td>0.048</td><td>11.284</td><td>0.006</td><td>10.957</td><td>0</td><td>10.661</td><td>0</td></tr><tr><td>9508.80461</td><td>12.141</td><td>0.042</td><td>11.277</td><td>0.001</td><td>10.953</td><td>0</td><td>10.65</td><td>0.005</td></tr><tr><td>9508.80664</td><td>12.13</td><td>0.037</td><td>11.278</td><td>0.004</td><td>10.952</td><td>0.002</td><td>10.653</td><td>0.01</td></tr><tr><td>9509.80682</td><td>12.039</td><td>0.042</td><td>11.218</td><td>0</td><td>10.929</td><td>0.013</td><td>10.627</td><td>0.008</td></tr><tr><td>9509.80888</td><td>12.021</td><td>0.035</td><td>11.209</td><td>0.001</td><td>10.924</td><td>0.017</td><td>10.63</td><td>0.013</td></tr><tr><td>9509.81092</td><td>12.027</td><td>0.033</td><td>11.217</td><td>0.001</td><td>10.929</td><td>0.014</td><td>10.618</td><td>0.012</td></tr><tr><td>9510.86553</td><td>12.498</td><td>0.025</td><td>11.493</td><td>0.001</td><td>11.105</td><td>0.016</td><td>10.788</td><td>0.008</td></tr><tr><td>9510.86759</td><td>12.511</td><td>0.024</td><td>11.494</td><td>0.002</td><td>11.123</td><td>0.013</td><td>10.787</td><td>0.008</td></tr><tr><td>9510.86963</td><td>12.517</td><td>0.022</td><td>11.495</td><td>0.001</td><td>11.109</td><td>0.008</td><td>10.801</td><td>0.019</td></tr><tr><td>9513.71086</td><td>12.66</td><td>0.017</td><td>11.618</td><td>0.002</td><td>11.202</td><td>0.011</td><td>10.866</td><td>0.013</td></tr><tr><td>9513.71300</td><td>12.659</td><td>0.017</td><td>11.615</td><td>0.001</td><td>11.209</td><td>0.013</td><td>10.867</td><td>0</td></tr><tr><td>9513.71503</td><td>12.653</td><td>0.016</td><td>11.601</td><td>0.001</td><td>11.216</td><td>0.009</td><td>10.868</td><td>0.011</td></tr><tr><td>9514.78750</td><td>12.531</td><td>0.034</td><td>11.536</td><td>0.001</td><td>11.169</td><td>0.013</td><td>10.841</td><td>0.004</td></tr><tr><td>9514.78954</td><td>12.543</td><td>0.028</td><td>11.534</td><td>0.001</td><td>11.158</td><td>0.013</td><td>10.832</td><td>0.013</td></tr><tr><td>9514.79157</td><td>12.544</td><td>0.026</td><td>11.535</td><td>0.003</td><td>11.17</td><td>0.01</td><td>10.819</td><td>0.008</td></tr><tr><td>9515.79177</td><td>12.677</td><td>0.016</td><td>11.617</td><td>0.003</td><td>11.184</td><td>0.002</td><td>10.861</td><td>0.025</td></tr><tr><td>9515.79383</td><td>12.693</td><td>0.017</td><td>11.622</td><td>0.006</td><td>11.191</td><td>0.007</td><td>10.868</td><td>0.017</td></tr><tr><td>9515.79586</td><td>12.689</td><td>0.019</td><td>11.621</td><td>0.008</td><td>11.184</td><td>0</td><td>10.854</td><td>0.02</td></tr><tr><td>9517.72459</td><td>12.762</td><td>0.019</td><td>11.641</td><td>0.007</td><td>11.252</td><td>0.008</td><td>10.898</td><td>0.009</td></tr><tr><td>9517.72665</td><td>12.774</td><td>0.014</td><td>11.642</td><td>0.01</td><td>11.249</td><td>0.002</td><td>10.904</td><td>0.004</td></tr><tr><td>9517.72869</td><td>12.769</td><td>0.016</td><td>11.638</td><td>0.01</td><td>11.247</td><td>0.006</td><td>10.9</td><td>0.009</td></tr><tr><td>9518.73156</td><td>12.774</td><td>0.021</td><td>11.663</td><td>0.002</td><td>11.228</td><td>0.012</td><td>10.898</td><td>0.008</td></tr><tr><td>9518.73362</td><td>12.778</td><td>0.019</td><td>11.659</td><td>0.002</td><td>11.222</td><td>0.012</td><td>10.887</td><td>0.008</td></tr><tr><td>9518.73568</td><td>12.771</td><td>0.017</td><td>11.659</td><td>0.001</td><td>11.227</td><td>0.013</td><td>10.892</td><td>0.011</td></tr><tr><td>9519.73968</td><td>12.76</td><td>0.051</td><td>11.693</td><td>0.012</td><td>11.223</td><td>0.01</td><td>10.879</td><td>0.006</td></tr><tr><td>9519.74174</td><td>12.758</td><td>0.051</td><td>11.683</td><td>0.011</td><td>11.218</td><td>0.008</td><td>10.877</td><td>0.006</td></tr><tr><td>9519.74378</td><td>12.766</td><td>0.051</td><td>11.682</td><td>0.011</td><td>11.214</td><td>0.01</td><td>10.882</td><td>0.011</td></tr><tr><td>9520.74391</td><td>12.641</td><td>0.016</td><td>11.597</td><td>0.001</td><td>11.176</td><td>0.012</td><td>10.846</td><td>0.008</td></tr><tr><td>9520.74595</td><td>12.651</td><td>0.016</td><td>11.596</td><td>0.001</td><td>11.177</td><td>0.015</td><td>10.866</td><td>0.007</td></tr><tr><td>9520.74801</td><td>12.662</td><td>0.014</td><td>11.603</td><td>0.001</td><td>11.177</td><td>0.011</td><td>10.867</td><td>0.008</td></tr></tbody></table>
**Appendix B: Line profile correlation matrices**

This section provides correlation matrices between line profiles (Johns & Basri 1995). Correlation coefficients are computed across the line profiles for every pair of velocity channels common to the two profiles. Fig. B.1 to B.4 present the correlation matrices of the optical and near-infrared line profiles studied here.

<figure><img src="image_21.png" /><figcaption>**Fig. B.1.** Correlation matrices for hydrogen optical line profiles computed from the 15 OHP/SOPHIE spectra obtained during the campaign: $H\alpha \star H\alpha$ (left), $H\beta \star H\beta$ (center), and $H\alpha \star H\beta$ (right).</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Fig. B.2.** Correlation matrices for the hydrogen near-infrared line profiles computed from the 34 CFHT/SPIRou spectra obtained during the campaign: $\text{Pa}\beta \star \text{Pa}\beta$ (left), $\text{Bry} \star \text{Bry}$ (center), and $\text{Pa}\beta \star \text{Bry}$ (right).</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>**Fig. B.3.** Correlation matrices between optical and near-infrared hydrogen line profiles computed from ten OHP/SOPHIE and ten CFHT/SPIRou spectra obtained over the same nights during the October runs: $H\alpha \star \text{Pa}\beta$ (left), $H\beta \star \text{Pa}\beta$ (center left), $H\alpha \star \text{Bry}$ (center right), and $H\beta \star \text{Bry}$ (right).</figcaption></figure>
<figure><img src="image_24.png" /><figcaption>Fig. B.4. Correlation matrices for the HeI 10830 Å and hydrogen lines: HeI ★ HeI (**left**), Paβ ★ HeI (**center**), and Hβ ★ HeI (**right**). The Bry ★ HeI and Hα ★ HeI matrices are not shown as they are similar to those of Paβ ★ HeI and Hβ ★ HeI, respectively.</figcaption></figure>

## Appendix C: EWs and J-band measurements from the ExTrA spectra

We provide the HeI 10830 Å, Paβ, Paγ, and Paδ line EWs and the J-band photometry measured from the ExTrA spectra in Table C.1. For each night, the table lists the mean observation time, the median EW measurement and its standard deviation for each spectral line, and the J-band magnitude and its error.

<table><caption>Table C.1. EW measurements and J-band photometry from the ExTrA spectra.</caption><thead><tr><th>Julian date<br>(2,450,000+)</th><th>EW(HeI)<br>(Å)</th><th>EW(Paβ)<br>(Å)</th><th>EW(Pay)<br>(Å)</th><th>EW(Paδ)<br>(Å)</th><th>J<br>(mag)</th><th>errJ<br>(mag)</th></tr></thead><tbody><tr><td>9501.89138</td><td>2.95 ± 0.80</td><td>10.45 ± 0.28</td><td>8.23 ± 0.88</td><td>2.37 ± 0.80</td><td>9.487</td><td>0.028</td></tr><tr><td>9504.85514</td><td>2.43 ± 0.50</td><td>9.33 ± 0.61</td><td>4.93 ± 0.70</td><td>1.48 ± 0.67</td><td>9.499</td><td>0.027</td></tr><tr><td>9506.86219</td><td>2.43 ± 0.58</td><td>10.91 ± 0.72</td><td>7.28 ± 0.60</td><td>2.80 ± 0.64</td><td>9.451</td><td>0.027</td></tr><tr><td>9507.81204</td><td>7.87 ± 0.70</td><td>14.65 ± 1.00</td><td>9.57 ± 0.62</td><td>3.92 ± 0.73</td><td>9.411</td><td>0.027</td></tr><tr><td>9509.84834</td><td>10.06 ± 0.61</td><td>16.56 ± 0.82</td><td>9.05 ± 0.60</td><td>4.23 ± 0.69</td><td>9.316</td><td>0.027</td></tr><tr><td>9510.88153†</td><td>5.91 ± 0.61</td><td>12.95 ± 0.93</td><td>7.62 ± 0.58</td><td>2.42 ± 0.59</td><td>—</td><td>—</td></tr><tr><td>9512.79907</td><td>1.87 ± 0.30</td><td>7.09 ± 1.20</td><td>7.14 ± 0.27</td><td>1.14 ± 0.02</td><td>9.461</td><td>0.027</td></tr><tr><td>9513.79425</td><td>5.50 ± 0.69</td><td>10.50 ± 1.07</td><td>8.73 ± 1.39</td><td>2.90 ± 1.29</td><td>9.407</td><td>0.027</td></tr><tr><td>9514.84366</td><td>2.16 ± 0.67</td><td>10.56 ± 0.83</td><td>7.42 ± 0.71</td><td>2.56 ± 0.68</td><td>9.444</td><td>0.027</td></tr><tr><td>9517.84528</td><td>4.26 ± 0.74</td><td>11.42 ± 0.78</td><td>7.85 ± 0.95</td><td>2.74 ± 0.83</td><td>9.495</td><td>0.037</td></tr><tr><td>9518.84178</td><td>1.08 ± 0.55</td><td>7.88 ± 0.73</td><td>7.77 ± 0.50</td><td>2.08 ± 0.82</td><td>9.482</td><td>0.027</td></tr><tr><td>9519.86259</td><td>4.16 ± 0.75</td><td>9.48 ± 0.71</td><td>7.76 ± 0.56</td><td>3.05 ± 1.01</td><td>9.464</td><td>0.027</td></tr><tr><td>9520.85926</td><td>4.42 ± 0.66</td><td>13.95 ± 0.57</td><td>8.53 ± 0.78</td><td>2.29 ± 0.35</td><td>9.444</td><td>0.027</td></tr><tr><td>9521.87770</td><td>5.36 ± 0.48</td><td>12.38 ± 0.25</td><td>9.63 ± 0.94</td><td>1.63 ± 0.23</td><td>9.424</td><td>0.027</td></tr><tr><td>9522.79321</td><td>3.04 ± 0.31</td><td>14.32 ± 0.41</td><td>7.54 ± 0.56</td><td>2.51 ± 0.48</td><td>9.448</td><td>0.027</td></tr><tr><td>9524.87753</td><td>1.01 ± 0.84</td><td>12.81 ± 1.07</td><td>8.01 ± 0.69</td><td>1.40 ± 0.27</td><td>9.467</td><td>0.027</td></tr><tr><td>9525.86383</td><td>2.47 ± 0.42</td><td>11.31 ± 0.86</td><td>6.61 ± 0.49</td><td>1.13 ± 0.76</td><td>9.473</td><td>0.027</td></tr><tr><td>9526.86452</td><td>2.24 ± 0.47</td><td>8.44 ± 0.63</td><td>5.00 ± 0.70</td><td>1.17 ± 1.18</td><td>9.476</td><td>0.027</td></tr><tr><td>9527.86505</td><td>2.71 ± 0.64</td><td>10.61 ± 0.76</td><td>5.91 ± 0.46</td><td>1.08 ± 0.74</td><td>9.442</td><td>0.027</td></tr><tr><td>9528.82389</td><td>1.59 ± 0.52</td><td>8.35 ± 1.12</td><td>4.64 ± 0.72</td><td>0.02 ± 0.51</td><td>9.475</td><td>0.027</td></tr><tr><td>9529.85920</td><td>0.03 ± 0.65</td><td>6.70 ± 0.70</td><td>3.97 ± 0.85</td><td>0.33 ± 0.59</td><td>9.463</td><td>0.027</td></tr><tr><td>9530.85025</td><td>0.30 ± 0.81</td><td>9.43 ± 0.96</td><td>5.83 ± 1.17</td><td>1.13 ± 0.75</td><td>9.477</td><td>0.029</td></tr><tr><td>9531.81805</td><td>2.70 ± 0.46</td><td>8.98 ± 0.62</td><td>6.50 ± 0.47</td><td>2.82 ± 0.46</td><td>9.460</td><td>0.027</td></tr><tr><td>9532.82519</td><td>3.27 ± 0.61</td><td>10.56 ± 0.82</td><td>5.81 ± 0.57</td><td>2.67 ± 0.65</td><td>9.429</td><td>0.035</td></tr><tr><td>9533.85615</td><td>4.04 ± 0.65</td><td>6.78 ± 0.60</td><td>4.30 ± 0.61</td><td>1.31 ± 0.63</td><td>9.428</td><td>0.027</td></tr><tr><td>9534.84045</td><td>4.57 ± 0.86</td><td>8.43 ± 0.95</td><td>5.70 ± 0.60</td><td>2.12 ± 0.74</td><td>9.425</td><td>0.090</td></tr><tr><td>9535.83697</td><td>1.85 ± 0.52</td><td>8.08 ± 0.71</td><td>5.85 ± 0.68</td><td>1.84 ± 0.62</td><td>9.456</td><td>0.053</td></tr><tr><td>9536.86516</td><td>0.15 ± 0.64</td><td>6.93 ± 2.51</td><td>6.94 ± 0.30</td><td>2.34 ± 0.03</td><td>9.358</td><td>0.027</td></tr><tr><td>9537.78315</td><td>3.45 ± 1.60</td><td>6.97 ± 0.98</td><td>5.43 ± 1.00</td><td>2.06 ± 0.69</td><td>9.545</td><td>0.027</td></tr><tr><td>9538.80149</td><td>2.83 ± 1.64</td><td>5.43 ± 1.42</td><td>4.09 ± 1.20</td><td>1.59 ± 0.86</td><td>9.531</td><td>0.027</td></tr><tr><td>9539.81823</td><td>3.43 ± 0.68</td><td>7.42 ± 1.27</td><td>5.50 ± 0.60</td><td>1.48 ± 0.65</td><td>9.396</td><td>0.028</td></tr><tr><td>9540.84522</td><td>3.61 ± 0.72</td><td>10.38 ± 2.17</td><td>6.80 ± 0.84</td><td>2.58 ± 0.70</td><td>9.313</td><td>0.037</td></tr><tr><td>9543.79258</td><td>2.94 ± 0.46</td><td>10.10 ± 0.79</td><td>6.48 ± 0.40</td><td>1.94 ± 0.61</td><td>9.419</td><td>0.027</td></tr><tr><td>9544.73410</td><td>2.34 ± 0.48</td><td>9.12 ± 0.59</td><td>5.92 ± 0.42</td><td>1.68 ± 0.52</td><td>9.376</td><td>0.027</td></tr><tr><td>9547.68435</td><td>0.08 ± 0.37</td><td>4.16 ± 0.69</td><td>4.52 ± 0.58</td><td>0.46 ± 0.36</td><td>9.443</td><td>0.027</td></tr><tr><td>9548.72798</td><td>-1.19 ± 0.62</td><td>2.79 ± 0.90</td><td>4.06 ± 0.62</td><td>0.30 ± 0.68</td><td>9.568</td><td>0.027</td></tr></tbody></table>

<table><caption>Table C.1. continued.</caption><thead><tr><th>Julian date<br>(2,450,000+)</th><th>EW(HeI)<br>(Å)</th><th>EW(Paβ)<br>(Å)</th><th>EW(Pay)<br>(Å)</th><th>EW(Paδ)<br>(Å)</th><th>J<br>(mag)</th><th>errJ<br>(mag)</th></tr></thead><tbody><tr><td>9549.73881</td><td>1.79 ± 0.63</td><td>6.06 ± 1.24</td><td>5.67 ± 0.82</td><td>1.21 ± 0.69</td><td>9.537</td><td>0.027</td></tr><tr><td>9550.69491</td><td>0.74 ± 0.41</td><td>3.22 ± 0.74</td><td>3.84 ± 0.59</td><td>1.16 ± 0.44</td><td>9.492</td><td>0.027</td></tr><tr><td>9551.70856</td><td>2.48 ± 0.64</td><td>5.54 ± 0.70</td><td>4.66 ± 0.48</td><td>1.22 ± 0.54</td><td>9.474</td><td>0.028</td></tr><tr><td>9553.69774</td><td>1.13 ± 0.64</td><td>6.05 ± 0.62</td><td>4.50 ± 0.63</td><td>0.87 ± 0.61</td><td>9.443</td><td>0.028</td></tr><tr><td>9555.72248</td><td>4.05 ± 0.75</td><td>12.83 ± 1.35</td><td>7.78 ± 0.94</td><td>2.04 ± 0.72</td><td>9.266</td><td>0.028</td></tr><tr><td>9556.67400</td><td>1.85 ± 0.57</td><td>11.64 ± 1.85</td><td>6.23 ± 1.62</td><td>1.48 ± 0.61</td><td>9.333</td><td>0.027</td></tr><tr><td>9560.71477</td><td>0.99 ± 0.56</td><td>6.39 ± 0.68</td><td>5.64 ± 0.63</td><td>1.31 ± 0.70</td><td>9.507</td><td>0.028</td></tr><tr><td>9561.73574</td><td>2.54 ± 0.44</td><td>7.87 ± 0.52</td><td>7.21 ± 0.44</td><td>2.16 ± 0.57</td><td>9.481</td><td>0.028</td></tr><tr><td>9564.76316</td><td>6.31 ± 0.66</td><td>11.98 ± 0.74</td><td>6.73 ± 0.67</td><td>1.83 ± 0.66</td><td>9.394</td><td>0.029</td></tr><tr><td>9566.66048</td><td>-1.91 ± 0.85</td><td>3.89 ± 0.75</td><td>3.44 ± 0.59</td><td>0.30 ± 0.51</td><td>9.505</td><td>0.028</td></tr><tr><td>9567.68335</td><td>2.96 ± 0.69</td><td>8.80 ± 0.92</td><td>5.48 ± 0.39</td><td>1.18 ± 0.64</td><td>9.470</td><td>0.027</td></tr><tr><td>9569.70386</td><td>2.77 ± 0.56</td><td>7.92 ± 0.78</td><td>5.88 ± 0.90</td><td>1.35 ± 0.75</td><td>9.397</td><td>0.029</td></tr><tr><td>9570.68277</td><td>4.73 ± 0.48</td><td>10.35 ± 1.11</td><td>7.26 ± 0.55</td><td>2.46 ± 0.79</td><td>9.399</td><td>0.027</td></tr><tr><td>9573.66353</td><td>4.33 ± 0.48</td><td>10.45 ± 0.41</td><td>9.01 ± 0.66</td><td>2.64 ± 0.41</td><td>9.432</td><td>0.027</td></tr><tr><td>9577.61918</td><td>2.11 ± 0.36</td><td>8.54 ± 0.65</td><td>5.62 ± 0.51</td><td>1.43 ± 0.57</td><td>9.401</td><td>0.027</td></tr><tr><td>9578.66562</td><td>7.11 ± 0.65</td><td>12.34 ± 0.76</td><td>7.62 ± 0.91</td><td>3.88 ± 0.55</td><td>9.368</td><td>0.027</td></tr><tr><td>9580.62157</td><td>4.88 ± 0.50</td><td>10.48 ± 0.72</td><td>7.06 ± 0.53</td><td>2.42 ± 0.56</td><td>9.410</td><td>0.027</td></tr><tr><td>9581.66843</td><td>6.99 ± 0.42</td><td>12.16 ± 0.54</td><td>7.38 ± 0.64</td><td>2.40 ± 0.46</td><td>9.331</td><td>0.027</td></tr><tr><td>9582.60702</td><td>3.56 ± 0.56</td><td>10.93 ± 0.72</td><td>6.81 ± 0.86</td><td>1.99 ± 0.52</td><td>9.390</td><td>0.027</td></tr><tr><td>9586.69987</td><td>6.29 ± 0.58</td><td>13.34 ± 1.26</td><td>9.30 ± 0.66</td><td>3.94 ± 0.51</td><td>9.406</td><td>0.029</td></tr><tr><td>9588.70981</td><td>2.08 ± 0.26</td><td>6.93 ± 0.69</td><td>6.47 ± 0.55</td><td>1.77 ± 0.67</td><td>9.426</td><td>0.027</td></tr><tr><td>9589.60234</td><td>2.20 ± 0.68</td><td>9.23 ± 1.17</td><td>7.62 ± 0.58</td><td>2.38 ± 0.52</td><td>9.429</td><td>0.027</td></tr><tr><td>9596.62015</td><td>2.49 ± 0.58</td><td>9.68 ± 1.33</td><td>6.58 ± 0.47</td><td>2.26 ± 0.62</td><td>9.458</td><td>0.028</td></tr><tr><td>9597.67444</td><td>4.58 ± 0.65</td><td>11.33 ± 0.95</td><td>9.61 ± 0.56</td><td>2.91 ± 0.64</td><td>9.450</td><td>0.028</td></tr><tr><td>9598.59529</td><td>4.30 ± 0.54</td><td>10.94 ± 0.86</td><td>7.80 ± 0.48</td><td>2.43 ± 0.63</td><td>9.397</td><td>0.028</td></tr><tr><td>9600.55490</td><td>7.14 ± 0.48</td><td>16.46 ± 0.86</td><td>10.65 ± 0.68</td><td>4.19 ± 0.58</td><td>9.424</td><td>0.028</td></tr><tr><td>9601.56718</td><td>2.41 ± 0.72</td><td>8.23 ± 1.12</td><td>6.02 ± 0.59</td><td>1.59 ± 0.60</td><td>9.456</td><td>0.028</td></tr><tr><td>9602.53523</td><td>1.35 ± 0.63</td><td>9.00 ± 0.75</td><td>6.44 ± 0.80</td><td>2.51 ± 0.61</td><td>9.486</td><td>0.027</td></tr><tr><td>9605.54455</td><td>10.34 ± 0.74</td><td>17.45 ± 1.23</td><td>11.39 ± 0.82</td><td>4.42 ± 0.90</td><td>9.372</td><td>0.027</td></tr><tr><td>9606.55663</td><td>4.71 ± 0.58</td><td>12.46 ± 1.33</td><td>7.92 ± 0.82</td><td>2.47 ± 0.82</td><td>9.409</td><td>0.028</td></tr><tr><td>9607.55053</td><td>2.61 ± 0.51</td><td>11.21 ± 1.02</td><td>7.28 ± 0.60</td><td>2.37 ± 0.57</td><td>9.443</td><td>0.027</td></tr><tr><td>9608.52952</td><td>-0.27 ± 0.75</td><td>7.99 ± 1.16</td><td>5.62 ± 1.00</td><td>1.60 ± 0.60</td><td>9.484</td><td>0.027</td></tr><tr><td>9615.52033</td><td>5.17 ± 0.56</td><td>13.40 ± 1.17</td><td>8.94 ± 0.97</td><td>3.34 ± 0.62</td><td>9.409</td><td>0.027</td></tr><tr><td>9616.55594</td><td>6.72 ± 0.64</td><td>13.72 ± 1.09</td><td>9.56 ± 0.88</td><td>3.19 ± 0.62</td><td>9.360</td><td>0.027</td></tr><tr><td>9617.55930</td><td>6.54 ± 0.92</td><td>15.62 ± 0.69</td><td>9.79 ± 0.72</td><td>3.73 ± 0.72</td><td>9.354</td><td>0.030</td></tr><tr><td>9618.55535</td><td>8.27 ± 1.11</td><td>18.41 ± 1.05</td><td>11.10 ± 0.65</td><td>4.79 ± 0.72</td><td>9.332</td><td>0.028</td></tr><tr><td>9619.56342</td><td>7.84 ± 0.76</td><td>14.54 ± 0.90</td><td>9.04 ± 0.72</td><td>3.59 ± 0.78</td><td>9.417</td><td>0.028</td></tr><tr><td>9623.52755</td><td>11.79 ± 0.81</td><td>17.61 ± 1.58</td><td>10.38 ± 0.73</td><td>4.42 ± 0.75</td><td>9.320</td><td>0.027</td></tr><tr><td>9624.55795</td><td>9.62 ± 0.78</td><td>14.33 ± 1.01</td><td>8.17 ± 0.75</td><td>3.10 ± 0.65</td><td>9.365</td><td>0.028</td></tr><tr><td>9625.53847</td><td>6.42 ± 0.92</td><td>10.27 ± 1.41</td><td>6.52 ± 0.85</td><td>1.64 ± 0.70</td><td>9.438</td><td>0.028</td></tr><tr><td>9626.54670</td><td>-0.82 ± 0.70</td><td>5.39 ± 1.01</td><td>4.91 ± 0.92</td><td>1.34 ± 0.86</td><td>9.527</td><td>0.030</td></tr><tr><td>9629.53317</td><td>7.30 ± 0.63</td><td>12.52 ± 0.87</td><td>7.25 ± 0.65</td><td>2.74 ± 0.60</td><td>9.336</td><td>0.028</td></tr><tr><td>9630.53061</td><td>11.17 ± 0.86</td><td>18.58 ± 0.89</td><td>10.69 ± 0.81</td><td>4.68 ± 0.87</td><td>9.342</td><td>0.027</td></tr><tr><td>9631.50698</td><td>8.20 ± 0.41</td><td>15.59 ± 1.21</td><td>9.15 ± 0.94</td><td>3.27 ± 0.30</td><td>9.409</td><td>0.027</td></tr><tr><td>9632.52380</td><td>1.15 ± 0.82</td><td>9.75 ± 1.99</td><td>5.78 ± 0.74</td><td>1.04 ± 0.73</td><td>9.627</td><td>0.028</td></tr><tr><td>9633.54300</td><td>1.99 ± 0.41</td><td>10.77 ± 1.32</td><td>6.97 ± 0.60</td><td>1.70 ± 0.31</td><td>9.445</td><td>0.027</td></tr><tr><td>9634.53461</td><td>2.68 ± 0.59</td><td>8.75 ± 1.34</td><td>5.18 ± 0.56</td><td>1.25 ± 0.63</td><td>9.413</td><td>0.027</td></tr><tr><td>9634.52095</td><td>1.44 ± 0.58</td><td>14.07 ± 1.00</td><td>10.88 ± 0.90</td><td>3.03 ± 0.62</td><td>9.393</td><td>0.028</td></tr><tr><td>9644.52110</td><td>1.53 ± 0.68</td><td>10.47 ± 1.09</td><td>8.49 ± 0.89</td><td>1.87 ± 0.83</td><td>9.467</td><td>0.029</td></tr><tr><td>9645.51054</td><td>2.58 ± 0.50</td><td>10.74 ± 2.45</td><td>9.04 ± 0.41</td><td>2.61 ± 0.71</td><td>9.416</td><td>0.027</td></tr><tr><td>9646.52604</td><td>6.67 ± 1.02</td><td>14.56 ± 1.53</td><td>9.61 ± 0.92</td><td>2.84 ± 0.83</td><td>9.375</td><td>0.028</td></tr><tr><td>9647.50937</td><td>2.58 ± 0.66</td><td>10.79 ± 1.33</td><td>7.04 ± 1.01</td><td>1.48 ± 0.64</td><td>9.392</td><td>0.028</td></tr></tbody></table>

<footnote>† No J-band photometry could be obtained for that night due to poor seeing.</footnote>
**Appendix D: nIR veiling and line EWs**

We present in Figure D.1 the evolution of near-infrared veiling
in the JHK bands and of the HeI, Paβ, and Bry line EWs over the
course of the campaign.

<figure><img src="image_25.png" /><figcaption>**Fig. D.1.** Veiling measured on SPIRou spectra in the JHK bands (*triangles*) and the EW of the HeI, Paβ, and Bry lines (*dots*) plotted as a function of Julian date (*top*) and rotational phase (*bottom*).</figcaption></figure>
**Appendix E: Line profile variability on successive SPIRou runs**

Figures E.1 to E.4 show the line profile variability of the HeI 10830 Å, Paβ, and Bry line profiles for the successive SPIRou runs in September, October, November, and December 2021. The figures also include 2D periodograms for each line. We

note, however, that only the October SPIRou run is long enough, extending over 14 days, to yield significant results when searching for a modulation of the line profiles around the rotational period of 6.04 days. The September, November, and December runs lasted for only 9, 6, and 9 days, respectively, which is too short to reliably investigate periods longer than 5 days.

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

**Fig. E.1.** Near-infrared HeI (**left**), Paβ (**center**), and Bry (**right**) line profiles obtained over nine days during the September 2021 SPIRou run. **Top:** Line profiles are plotted as a function of Julian date (**left subpanels**) and rotational phase (**right subpanels**). The colors represent successive rotational cycles. **Bottom:** 2D periodograms across the line profiles. The dotted horizontal red line drawn at a frequency of 0.166 day⁻¹ indicates the stellar rotational period. The white curve displays the mean line profile. The color code reflects the periodogram power from zero (**blue**) to 1 (**red**).
<figure><img src="image_27.png" /></figure>

**Fig. E.2.** Near-infrared HeI (**left**), Paβ (**center**), and Bry (**right**) line profiles obtained over 14 days during the October 2021 SPIRou run. **Top:** Line profiles are plotted as a function of Julian date (**left subpanels**) and rotational phase (**right subpanels**). The colors represent successive rotational cycles. **Bottom:** 2D periodograms across the line profiles. The dotted horizontal red line drawn at a frequency of 0.166 day⁻¹ indicates the stellar rotational period. The white curve displays the mean line profile. The color code reflects the periodogram power from zero (**blue**) to 0.9 (**red**).
<figure><img src="image_28.png" /></figure>

**Fig. E.3.** Near-infrared HeI (**left**), Paβ (**center**), and Bry (**right**) line profiles obtained over six days during the November 2021 SPIRou run. **Top:** Line profiles are plotted as a function of Julian date (**left subpanels**) and rotational phase (**right subpanels**). The colors represent successive rotational cycles. **Bottom:** 2D periodograms across the line profiles. The dotted horizontal red line drawn at a frequency of 0.166 day⁻¹ indicates the stellar rotational period. The white curve displays the mean line profile. The color code reflects the periodogram power from zero (**blue**) to 1 (**red**).
<figure><img src="image_29.png" /></figure>

**Fig. E.4.** Near-infrared HeI (**left**), Paβ (**center**), and Bry (**right**) line profiles obtained over nine days during the December 2021 SPIRou run. **Top:** Line profiles are plotted as a function of Julian date (**left subpanels**) and rotational phase (**right subpanels**). The colors represent successive rotational cycles. **Bottom:** 2D periodograms across the line profiles. The dotted horizontal red line drawn at a frequency of 0.166 day⁻¹ indicates the stellar rotational period. The white curve displays the mean line profile. The color code reflects the periodogram power from zero (**blue**) to 1 (**red**).