# Monitoring the large-scale magnetic field of AD Leo with SPIRou, ESPaDOnS, and Narval

S. Bellotti, Julien Morin, L. T. Lehmann, C. P. Folsom, G. A. J. Hussain, P.

Petit, J-F Donati, A. Lavail, A. Carmona, E. Martioli, et al.

► To cite this version:

S. Bellotti, Julien Morin, L. T. Lehmann, C. P. Folsom, G. A. J. Hussain, et al. Monitoring the large-scale magnetic field of AD Leo with SPIRou, ESPaDOnS, and Narval. Astronomy & Astrophysics - A&A, 2023, 676, pp.A56. 10.1051/0004-6361/202346845. hal-04187255

HAL Id: hal-04187255

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

Submitted on 24 Aug 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
# Monitoring the large-scale magnetic field of AD Leo with SPIRou, ESPaDOnS, and Narval

## Towards a magnetic polarity reversal?

S. Bellotti<sup>1,2</sup>, J. Morin<sup>3</sup>, L. T. Lehmann<sup>1</sup>, C. P. Folsom<sup>4</sup>, G. A. J. Hussain<sup>2</sup>, P. Petit<sup>1</sup>, J-F. Donati<sup>1</sup>, A. Lavail<sup>1,5</sup>, A. Carmona<sup>6</sup>, E. Martioli<sup>7,8</sup>, B. Romano Zaire<sup>9</sup>, E. Alecian<sup>6</sup>, C. Moutou<sup>1</sup>, P. Fouqué<sup>1</sup>, S. Alencar<sup>9</sup>, E. Artigau<sup>10</sup>, I. Boisse<sup>11</sup>, F. Bouchy<sup>12</sup>, C. Cadieux<sup>10</sup>, R. Cloutier<sup>15</sup>, N. J. Cook<sup>10</sup>, X. Delfosse<sup>6</sup>, R. Doyon<sup>10</sup>, G. Hébrard<sup>8,13</sup>, O. Kochukhov<sup>5</sup>, and G. A. Wade<sup>14</sup>

Received; accepted

## ABSTRACT

**Context.** One clear manifestation of dynamo action on the Sun is the 22-yr magnetic cycle, exhibiting a polarity reversal and a periodic conversion between poloidal and toroidal fields. For M dwarfs, several authors claim evidence of activity cycles from photometry and analyses of spectroscopic indices, but no clear polarity reversal has been identified from spectropolarimetric observations. These stars are excellent laboratories to investigate dynamo-powered magnetic fields under different stellar interior conditions, that is partly or fully convective.

**Aims.** Our aim is to monitor the evolution of the large-scale field of AD Leo, which has shown hints of a secular evolution from past dedicated spectropolarimetric campaigns. This is of central interest to inform distinct dynamo theories, contextualise the evolution of the solar magnetic field, and explain the variety of magnetic field geometries observed in the past.

**Methods.** We analysed near-infrared spectropolarimetric observations of the active M dwarf AD Leo taken with SPIRou between 2019 and 2020 and archival optical data collected with ESPaDOnS and Narval between 2006 and 2019. We searched for long-term variability in the longitudinal field, the width of unpolarised Stokes profiles, the unsigned magnetic flux derived from Zeeman broadening, and the geometry of the large-scale magnetic field using both Zeeman-Doppler imaging and principal component analysis.

**Results.** We found evidence of a long-term evolution of the magnetic field, featuring a decrease in axisymmetry (from 99% to 60%). This is accompanied by a weakening of the longitudinal field (−300 to −50 G) and a correlated increase in the unsigned magnetic flux (2.8 to 3.6 kG). Likewise, the width of the mean profile computed with selected near-infrared lines manifests a long-term evolution corresponding to field strength changes over the full time series, but does not exhibit modulation with the stellar rotation of AD Leo in individual epochs.

**Conclusions.** The large-scale magnetic field of AD Leo manifested first hints of a polarity reversal in late 2020 in the form of a substantially increased dipole obliquity, while the topology remained predominantly poloidal and dipolar for 14 yr. This suggests that low-mass M dwarfs with a dipole-dominated magnetic field can undergo magnetic cycles.

**Key words.** Stars: individual: AD Leo, Stars: magnetic field – Stars: activity – Techniques: polarimetric

<footnote n="1">¹ Institut de Recherche en Astrophysique et Planétologie, Université de Toulouse, CNRS, IRAP/UMR 5277, 14 avenue Edouard Belin, F-31400, Toulouse, France
e-mail: stefano.bellotti@irap.omp.eu</footnote>

<footnote n="2">² Science Division, Directorate of Science, European Space Research and Technology Centre (ESA/ESTEC), Keplerlaan 1, 2201 AZ, Noordwijk, The Netherlands</footnote>

<footnote n="3">³ Laboratoire Univers et Particules de Montpellier, Université de Montpellier, CNRS, F-34095, Montpellier, France</footnote>

<footnote n="4">⁴ Tartu Observatory, University of Tartu, Observatooriumi 1, Tõravere, 61602 Tartumaa, Estonia</footnote>

<footnote n="5">⁵ Department of Physics and Astronomy, Uppsala University, Box 516, SE-75120 Uppsala, Sweden</footnote>

<footnote n="6">⁶ Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote n="7">⁷ Laboratório Nacional de Astrofísica, Rua Estados Unidos 154, 37504-364, Itajubá - MG, Brazil</footnote>

<footnote n="8">⁸ Institut d'Astrophysique de Paris, CNRS, UMR 7095, Sorbonne Université, 98 bis bd Arago, 75014 Paris, France</footnote>

<footnote n="9">⁹ Universidade Federal de Minas Gerais, Belo Horizonte, MG, 31270-901, Brazil</footnote>

<footnote n="10">¹⁰ Université de Montréal, Département de Physique, IREX, Montréal, QC H3C 3J7, Canada</footnote>

<footnote n="11">¹¹ Aix Marseille Univ, CNRS, CNES, LAM, Marseille, France</footnote>

<footnote n="12">¹² Observatoire de Genève, Université de Genève, Chemin Pegasi, 51, 1290 Sauverny, Switzerland</footnote>

<footnote n="13">¹³ Observatoire de Haute Provence, St Michel l'Observatoire, France</footnote>

<footnote n="14">¹⁴ Department of Physics & Space Science, Royal Military College of Canada, PO Box 17000 Station Forces, Kingston, ON, Canada K7K 0C6</footnote>

<footnote n="15">¹⁵ Department of Physics & Astronomy, McMaster University, 1280 Main St West, Hamilton, ON, L8S 4L8, Canada</footnote>
# 1. Introduction

Studying stellar surface magnetic fields yields relevant insights into the internal structure of stars, as well as their essential role in stellar formation, evolution, and activity (Donati & Landstreet 2009). For cool stars, monitoring secular changes of the field's configuration provides useful feedback on the dynamo processes operating in the stellar interior and constraints on stellar wind models. The latter is fundamental to understanding atmospheric hydrodynamic escape of embedded planets since magnetic cycles modulate the star's activity level and thus its radiation output (Vidotto & Cleary 2020; Hazra et al. 2020).

The Sun is an important benchmark in this context: its long-term monitoring revealed a periodic variation in sunspot number, size, and latitude (Schwabe 1844; Maunder 1904; Hathaway 2010), and a polarity reversal of the large-scale magnetic field over a timescale of 11 yr (Hale et al. 1919). The proposed mechanism to reproduce these phenomena theoretically is the $\alpha\Omega$ dynamo (Parker 1955; Charbonneau 2010), namely the combination of differential rotation and cyclonic turbulence at the interface between the radiative and convective zones, known as tachocline. A different model is the Babcock-Leighton mechanism, which describes the conversion from a toroidal to poloidal field via a poleward migration of bipolar magnetic regions (Babcock 1961; Leighton 1969). However, there is still no model that can account for all the solar magnetic processes (Petrovay 2020).

For other stars, magnetic field measurements can be performed with two complementary approaches (Morin 2012; Reiners 2012). One is to model the Zeeman splitting in individual unpolarised spectral lines and estimate the total unsigned magnetic field, which is insensitive to polarity cancellation. The other is to apply tomographic techniques that use the polarisation properties of the Zeeman-split components to recover the orientation of the local field. In addition to these well-established methods, Lehmann & Donati (2022) show that fundamental properties of the large-scale field topology can be derived directly from the circularly polarised Stokes V time series using principal component analysis (PCA), without prior assumptions. This method allows us to qualitatively infer the predominant component of the field topology, as well as its complexity, axisymmetry, and evolution. Altogether, these observational constraints guide dynamo theories to a comprehensive description of the magnetic field generation and dynamic nature in the form of magnetic cycles (Reiners et al. 2010; Gregory et al. 2012; See et al. 2016).

Over the last three decades, Zeeman-Doppler imaging (ZDI, Semel 1989; Donati & Brown 1997) has been applied to reconstruct the poloidal and toroidal components of stellar magnetic fields, providing evidence of a wide variety of the large-scale magnetic topologies (e.g., Morin et al. 2016). Among rapidly rotating cool stars, the partly convective ones with masses above $0.5\,\text{M}_\odot$ tend to have moderate, predominantly toroidal large-scale fields generally featuring a non-axisymmetric poloidal component (Petit et al. 2008; Donati et al. 2008a; See et al. 2015). Those with masses between $0.2\,\text{M}_\odot$ and $0.5\,\text{M}_\odot$ – close to the fully convective boundary at $0.35\,\text{M}_\odot$ (Chabrier & Baraffe 1997) – generate stronger large-scale magnetic fields, dominated by a poloidal and axisymmetric component. For fully convective stars with $M<0.2\,\text{M}_\odot$, spectropolarimetric analyses have revealed a dichotomy of field geometries: either strong, mostly axisymmetric dipole-dominated or weak, non-axisymmetric multipole-dominated large-scale fields are observed (Morin et al. 2010). The latter findings could be understood either as a manifestation of dynamo bistability (Morin et al. 2011; Gastine et al. 2013), that is two dynamo branches that

coexist over a range of stellar rotation periods and masses, or of long magnetic cycles, implying that different topologies correspond to different phases of the cycle (Kitchatinov et al. 2014). Yet, no firm conclusion has been reached. In parallel, studies relying on the analysis of unpolarised spectra have shown that the average (unsigned) surface magnetic field of cool stars follows a classical rotation-activity relation including a non-saturated and a saturated (or quasi-saturated) regime, without a simple relation with the large-scale magnetic geometry (Reiners & Basri 2009; Shulyak et al. 2019; Kochukhov 2021; Reiners et al. 2022). Similarly, recent dynamo simulations conducted by, for instance, Zaire et al. (2022) confirm that the influence of rotation on convective motions alone could not explain the observed variety of magnetic geometry. Only in the case of fully convective very fast rotators, Shulyak et al. (2017) found that the strongest average fields were measured for stars with large-scale dipole-dominated fields. Kochukhov (2021) show that the fraction of magnetic energy contained in the large-scale field component is also the highest for these stars.

Cyclic trends for Sun-like stars were found via photometric and chromospheric activity (i.e. Ca II H&K lines) monitoring, and timescales shorter (e.g., 120 d for $\tau$ Boo, Mittag et al. 2017) or longer (≈ 20 yr for HD 1835; Boro Saikia et al. 2018) than the solar magnetic cycle were reported (Wilson 1968; Baliunas et al. 1995; Boro Saikia et al. 2018). Moreover, polarity flips of the large-scale field were detected for a handful of stars based on optical spectropolarimetric observations (Donati et al. 2008b; Petit et al. 2009; Fares et al. 2009; Morgenthaler et al. 2011; Boro Saikia et al. 2016; Rosén et al. 2016; Jeffers et al. 2018; Boro Saikia et al. 2022; Jeffers et al. 2022). For M dwarfs, numerous studies relying on photometry and spectroscopic indices claimed evidence of activity cycles (e.g., Gomes da Silva et al. 2012; Robertson et al. 2013; Mignon et al. 2023), and radio observations suggest the occurrence of polarity reversal at the end of the main sequence (Route 2016), but no polarity reversal has been directly observed with spectropolarimetry so far. This motivates long-term spectropolarimetric surveys, to reveal secular changes in the field topology and shed more light on the dynamo processes in action.

A well-known active M dwarf is AD Leo (GJ 388), whose mass ($0.42\,\text{M}_\odot$) falls at the boundary between the domains where toroidal- and dipole-dominated magnetic topologies have previously been identified, and thus represents an interesting laboratory to study stellar dynamos. Morin et al. (2008b) analysed the large-scale magnetic field from spectropolarimetric data sets collected with Narval at Télescope Bernard-Lyot in 2007 and 2008 and reported a stable, axisymmetric, dipole-dominated geometry. Later, Lavail et al. (2018) examined data collected with ESPaDOnS at Canada-France-Hawaii Telescope (CFHT) from 2012 and 2016, and showed an evolution of the field in the form of a global weakening (about 20%) and small-scale enhancement. The latter was quantitatively expressed by a decrease in the magnetic filling factor (from 13% to 7%), meaning that the field was more intense on local scales. No polarity reversal was reported on AD Leo (Lavail et al. 2018). The large-scale magnetic topology has remained stable since spectropolarimetric observations of AD Leo have been initiated (2007–2016): dominated by a strong axial dipole, the visible pole corresponding to negative radial field (magnetic field vector directed towards the star).

Here, we extend the magnetic analysis of AD Leo using both new optical ESPaDOnS observations collected in 2019 and near-infrared spectropolarimetric time series collected with SPIRou at CFHT in 2019 and 2020 under the SPIRou Legacy Survey (SLS), which adds to the previous optical data sets collected with
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Example Stokes profiles for AD Leo for the 9th December 2019 observation collected with SPIRou. From the top: Stokes *I* profile (unpolarised); Stokes *V* profile (circularly polarised) with a clear Zeeman signature of S/N = 6500, and Null *N* profile, used for quality check of the Stokes profiles (Donati et al. 1997; Bagnulo et al. 2009). The LSD profiles were obtained using a mask containing 838 lines. In all panels, the units are relative to the unpolarised continuum.</figcaption></figure>

ESPaDOnS and Narval between 2006 and 2016. The aim is to
apply distinct techniques to search for long-term variations that
may or may not resemble the solar behaviour.

The paper is structured as follows: in Sec. 2 we describe the
observations performed in the near-infrared and optical domains,
in Sec. 3 we outline the temporal analysis of the longitudinal
magnetic field, the Full-Width at Half Maximum (FWHM) of
the Stokes *I* profile, and the total magnetic flux inferred from
Zeeman broadening modelling. Then, we describe the magnetic
geometry reconstructions by means of ZDI and PCA. In Sec. 4
we discuss the wavelength dependence of magnetic field mea-
surements and in Sec. 5 we present our conclusions.

**2. Observations**

AD Leo is an M3.5 dwarf with a *V* and *H* band magnitude of
9.52 and 4.84, respectively (Zacharias et al. 2013), at a distance
of 4.9651±0.0007 pc (Gaia Collaboration et al. 2021). Its age
was estimated to be within 25 and 300 Myr by Shkolnik et al.
(2009). AD Leo has a rotation period of 2.23 days (Morin et al.
2008b; Carmona et al. 2023) and an inclination *i* = 20°, implying
an almost pole-on view (Morin et al. 2008b). Its high activity
level is seen in frequent flares (Muheki et al. 2020; Namekata
et al. 2020) and quantified by an X-ray-to-bolometric luminosity
ratio (log(L\*X/L\*bol)) of -3.62 (Wright et al. 2011) and a mean
CaII H&K index (logR\'*HK) of -4.00 (Boro Saikia et al. 2018).

AD Leo's mass is 0.42 M☉ (Mann et al. 2015; Cristofari
et al. 2023), which places it above the theoretical fully convec-
tive boundary at 0.35 M☉ (Chabrier & Baraffe 1997). The latter
value is in agreement with observations, as it has been invoked
to explain the dearth of stars with M_G ~10.2, known as Gaia
magnitude gap (Feiden et al. 2021). However, it is not an ab-
solute limit: age (Maeder & Meynet 2000) and metallicity affect
the depth of the convective envelope (van Saders & Pinsonneault
2012; Tanner et al. 2013), and the presence of strong magnetic
fields quenches convection and could push the theoretical bound-
ary towards later spectral type (Mullan & MacDonald 2001).

2.1. Near-infrared

A total of 77 spectropolarimetric observations in the near-
infrared were collected with the SpectroPolarimètre InfraRouge
(SPIRou) within the SLS. SPIRou is a stabilised high-resolution
near-infrared spectropolarimeter (Donati et al. 2020) mounted
on the 3.6 m CFHT atop Maunakea, Hawaii. It provides a full
coverage of the near-infrared spectrum from 0.96 to 2µm at a
spectral resolving power of R ~ 70,000. Optimal extraction of
SPIRou spectra was carried out with *A PipelinE to Reduce Ob-
servations* (APER0 v0.6.132), a fully automatic reduction pack-
age installed at CFHT (Cook et al. 2022). The same data set was
used in Carmona et al. (2023) to perform a velocimetric study
and reject the hypothesis of a planetary companion by Tuomi
et al. (2018) in favour of activity-induced variations, in agree-
ment with Carleo et al. (2020).

Observations were performed in circular polarisation mode
between February 2019 and June 2020, spanning 482 days in to-
tal; the journal of observations is available in Table D.1. The
mean airmass is 1.32 and the signal-to-noise ratio (S/N) at
1,650 nm per spectral element ranges from 68 to 218, with an
average of 168. We applied least-squares deconvolution (LSD)
to atomic spectral lines to derive averaged-line Stokes *I* (unpo-
larised) and *V* (circularly polarised) profiles (Donati et al. 1997;
Kochukhov et al. 2010). This numerical technique assumes the
spectrum to be the convolution between a mean line profile and
a line mask, that is to say a series of Dirac delta functions cen-
tred at each absorption line in the stellar spectrum, with corre-
sponding depths and Landé factors (i.e. sensitivities to the Zee-
man effect at a given wavelength). The output mean line pro-
file gathers the information of thousands of spectral lines and,
because of the consequent high S/N, enables the extraction of
polarimetric information from the spectrum. The adopted line
mask was generated using the Vienna Atomic Line Database<fnref n="1" />
(VALD, Ryabchikova et al. 2015) and a MARCS atmosphere
model (Gustafsson et al. 2008) with T<sub>eff</sub> = 3,500 K, log g =
5.0 [cm s<sup>-2</sup>] and v<sub>micro</sub> = 1 km s<sup>-1</sup>. It contains 1,400 atomic lines
between 950–2,600 nm and with known Landé factor (ranging
from 0 to 3) and with depth larger than 3 % of the continuum
level.

We discarded six observations in February 2019 since one
optical component of the instrument was not working nominally,
one observation in November 2019 because likely affected by a
flare (the corresponding radial velocity is >8 sigma lower than
the bulk of the measurements) and two observations in 2020 as
they led to noisier (by a factor of 10) LSD profiles. Therefore,
the data set analysed in this work comprises 68 polarimetric se-
quences, whose characteristics are reported in Table D.1.

The near-infrared observations were performed monthly be-
tween 2019 and 2020, except for two gaps of approximately two
and three months. There is also a gap of 1.5 month between
the end of 2019 and beginning of 2020. We thus split the time
series in four epochs to maintain coherency of magnetic activ-
ity over short time scales and for clearer visualisation: 2019a
(15th April 2019 to 21st June 2019, i.e. 2019.29 to 2019.47),
2019b (16th October 2019 to 12th December 2019, i.e. 2019.79
to 2019.95) 2020a (26th January 2020 to 12th March 2020, i.e.
2020.07 to 2020.19), and 2020b (8th May 2020 to 10th June
2020, i.e. 2020.35 to 2020.44).

The near-infrared domain covered by SPIRou is polluted by
strong and wide telluric bands due to Earth's atmospheric ab-
sorption. Their contribution to the stellar spectra is corrected

<footnote n="1">¹ http://vald.astro.uu.se/</footnote>
using a telluric transmission model (which is built from observations of standard stars since the start of SPIRou operations, and using Transmissions of the AtmosPhere for AStromomical data (TAPAS) atmospheric model Bertaux et al. 2014) and a PCA method implemented in the APERO pipeline (Artigau et al. 2014). To account for potential residuals in the telluric correction, we ignored the following intervals of the spectrum when computing the LSD profiles: [950, 979], [1116, 1163], [1331, 1490], [1790, 1985], [1995, 2029], [2250, 2500] nm. These intervals correspond to H₂O absorption regions, with transmission typically smaller than 40%. We assessed whether removing these telluric intervals optimises the quality of the Stokes V profiles. In a first test, we searched for stellar absorption lines deeper than 75% of the continuum level and within ±100 km s⁻¹ from telluric lines included in the transmission model of APERO. This approach allowed us to identify stellar lines that are contaminated by the telluric lines throughout the year. When the telluric-affected spectral lines were removed, no significant improvement was reported in the final LSD profiles, indicating a robust telluric correction as already reported in (Carmona et al. 2023). In a second test, we extended the intervals by 25 and 50 nm or reduced them by 10 nm and noticed an increment of the noise level in LSD profiles up to 20%, so we proceeded with the previous intervals.

Accounting for the ignored telluric intervals, the number of spectral lines used in LSD is 838. We show the LSD Stokes profiles for one example observation in Fig. 1. The average noise level in Stokes V for the entire time series is 1.6 · 10⁻⁴ relative to the unpolarised continuum, similar to the optical domain (Morin et al. 2008b). We also note that the profiles are broader than in the optical by more than 10 km s⁻¹, owing to a stronger Zeeman effect in the near-infrared domain (Zeeman 1897).

## 2.2. Optical

For most of the analyses presented here, we considered all archival observations collected with ESPaDOnS and Narval, and studied previously in Morin et al. (2008b) and Lavail et al. (2018). We also included six new observations taken in November 2019 (from 2019.87 to 2019.89) with ESPaDOnS for CFHT programme 19BC06, PI A. Lavail (reported in Table D.2). They are contemporaneous to the SPIRou ones for the same period, hence enabling us to study the dependence of the measured magnetic field strength on the wavelength domain employed (see Sec. 4).

ESPaDOnS is the optical spectropolarimeter on the 3.6 m CFHT located atop Mauna Kea in Hawaii, and Narval is the twin instrument on the 2 m TBL at the Pic du Midi Observatory in France (Donati 2003). The data reduction was performed with the LIBRE-ESPRIT pipeline (Donati et al. 1997), and the reduced spectra were retrieved from the PolarBase archive (Petit et al. 2014).

The LSD profiles were computed similar to the near-infrared, but using an optical VALD mask containing 3330 lines in range 350-1080 nm and with depths larger than 40% the continuum level, similar to Morin et al. (2008b); Bellotti et al. (2022). The number of lines used is 3240 and accounts for the removal of the following wavelength intervals, which are affected by telluric lines or in the vicinity of Hα: [627,632], [655.5,657], [686,697], [716,734], [759,770], [813,835], and [895,986] nm. For the 2019 observations, the average noise in Stokes V is 3 · 10⁻⁴ relative to the unpolarised continuum.

<figure><img src="image_3.png" /><figcaption>Fig. 2. Temporal evolution of the longitudinal magnetic field (Bᵢ). Top: full time series of measurements between 2006 and 2020 with ESPaDOnS, Narval and SPIRou. The shape of the data points corresponds to the wavelength domain, optical (circles) or near-infrared (pluses), and the colour to the epoch in which the observations were performed. Bottom: SPIRou time series split in four epochs phase-folded according to Eq. 1. The least-square sine fit corresponding to Eq. 3 (Stibbs 1950), to assess the change in magnetic obliquity, is shown as dashed lines. The six ESPaDOnS observations taken in 2019 are plotted as pink circles.</figcaption></figure>

In the next sections, the near-infrared and optical observations will be phased with the following ephemeris:

$$ \text{HJD} = 2458588.7573 + P_{\text{rot}} \cdot n_{\text{cyc}}, \quad (1) $$

where we used the first SPIRou observation taken in April 2019 as reference, P<sub>rot</sub>=2.23 days is the stellar rotation period (Morin et al. 2008b), and n<sub>cyc</sub> corresponds to the rotation cycle (see Table D.1).

# 3. Magnetic analysis

## 3.1. Longitudinal magnetic field

We measured the line-of-sight component of the magnetic field integrated over the stellar disk (Bᵢ) for all the available observations, in optical (2006–2019) and near-infrared (2019–2020). Since Bᵢ traces magnetic features present on the visible hemisphere, its temporal variations are modulated at the stellar rotation period and can be therefore used as a robust magnetic activity proxy (Folsom et al. 2016; Hébrard et al. 2016; Fouqué et al. 2023). Formally, it is computed as the first-order moment of Stokes V (Donati et al. 1997):

$$ B_{i} [G] = \frac{-2.14 \cdot 10^{11}}{\lambda_{0} g_{\text{eff}} c} \frac{\int vV(v)dv}{\int (I_{c} - I)dv}, \quad (2) $$
<table><caption>**Table 1.** Best fit values of the polar field and obliquity of Eq. 3 obtained for optical and near-infrared epochs.</caption><thead><tr><th></th><th>2006</th><th>2007</th><th>2008</th><th>2012</th><th>2016</th><th>2019a</th><th>2019b</th><th>2020a</th><th>2020b</th></tr></thead><tbody><tr><th>B<sub>p</sub> [G]</th><td>-868 ± 20</td><td>-842 ± 21</td><td>-873 ± 15</td><td>-736 ± 19</td><td>-555 ± 10</td><td>-814 ± 17</td><td>-882 ± 22</td><td>-828 ± 55</td><td>-923 ± 70</td></tr><tr><th>β [°]</th><td>13 ± 6</td><td>12 ± 4</td><td>26 ± 1</td><td>7 ± 4</td><td>12 ± 1</td><td>3 ± 4</td><td>23 ± 2</td><td>37 ± 4</td><td>59 ± 2</td></tr></tbody></table>

<table><caption>**Table 2.** Comparison of a constant line vs a sine fit for the FWHM phase variations.</caption><thead><tr><th>Epoch</th><th>Mask</th><th>Mean [km s<sup>-1</sup>]</th><th>Mean Error [km s<sup>-1</sup>]</th><th>STD [km s<sup>-1</sup>]</th><th>RMS<sub>const</sub> [km s<sup>-1</sup>]</th><th>χ<sup>2</sup><sub>r,const</sub></th><th>RMS<sub>sine</sub> [km s<sup>-1</sup>]</th><th>χ<sup>2</sup><sub>r,sine</sub></th></tr></thead><tbody><tr><td rowspan="2">2019a</td><td>default</td><td>19</td><td>0.29</td><td>0.40</td><td>0.40</td><td>1.8</td><td>0.39</td><td>2.0</td></tr><tr><td>g<sub>eff</sub> > 1.2</td><td>23</td><td>0.46</td><td>1.16</td><td>1.16</td><td>6.9</td><td>1.13</td><td>8.0</td></tr><tr><td rowspan="2">2019b</td><td>default</td><td>21</td><td>0.41</td><td>0.63</td><td>0.63</td><td>2.3</td><td>0.62</td><td>2.5</td></tr><tr><td>g<sub>eff</sub> > 1.2</td><td>27</td><td>0.84</td><td>1.77</td><td>1.77</td><td>6.7</td><td>1.71</td><td>6.9</td></tr><tr><td rowspan="2">2020a</td><td>default</td><td>21</td><td>0.41</td><td>0.68</td><td>0.68</td><td>3.2</td><td>0.61</td><td>3.4</td></tr><tr><td>g<sub>eff</sub> > 1.2</td><td>27</td><td>0.82</td><td>2.47</td><td>2.47</td><td>12.2</td><td>2.29</td><td>13.3</td></tr><tr><td rowspan="2">2020b</td><td>default</td><td>19</td><td>0.30</td><td>0.47</td><td>0.47</td><td>3.0</td><td>0.47</td><td>3.4</td></tr><tr><td>g<sub>eff</sub> > 1.2</td><td>22</td><td>0.46</td><td>0.83</td><td>0.83</td><td>5.3</td><td>0.79</td><td>7.3</td></tr></tbody></table>

**Notes.** The columns are: 1) subset of the time series, 2) line list used in LSD computation, 3) mean value of FWHM, 4) mean error bar on FWHM, 5) dispersion of the data set, 6) RMS (Root Mean Square) residual of a constant line fit equal to the average of the data set, 7) reduced $\chi^2$ of the constant fit, 8) RMS residual of a sine fit at the stellar rotation period, 9) reduced $\chi^2$ of the sine fit.

where $\lambda_0$ and $g_{eff}$ are the normalisation wavelength and Landé factor of the LSD profiles, $I_c$ is the continuum level, $v$ is the radial velocity associated to a point in the spectral line profile in the star's rest frame and $c$ the speed of light in vacuum. For the near-infrared and optical Stokes profiles, the normalisation wavelength and Landé factor are 1700 nm and 1.2144, and 700 nm and 1.1420, respectively. In accordance with the fact that near-infrared lines are broader than optical ones, the integration was carried out within $\pm 50$ km s$^{-1}$ from line centre in the former case and $\pm 30$ km s$^{-1}$ in the latter case, to include the absorption ranges of both Stokes I and V profiles.

The list of measurements is reported in Table D.3. The values are of constant sign (negative), which is expected when observing one polarity of a dipole almost aligned with the stellar rotation axis over the entire stellar rotation, especially for a star observed nearly pole-on as AD Leo (Morin et al. 2008b; Lavail et al. 2018). The near-infrared measurements range between -263 and -46 G, with an average of -179 G and a median error bar of 15 G. The optical measurements range between -297 and -155 G, with an average of -233 G and a median error bar of 10 G. The lower error bar is likely due to the narrower velocity range over which the optical measurements are performed, since less noise is introduced in Eq. 2. A discussion about chromatic differences in the longitudinal field measurements is presented in Sec. 4.

We plot the temporal evolution of $B_t$ in Fig. 2. In general, we note a secular weakening of the field strength over 14 yr, with an oscillation between 2016 and 2019 followed by a rapid decrease in strength (in absolute value). We also note that the intra-epoch dispersion increases for the last two epochs.

By phase-folding the near-infrared data at $P_{rot}$, we observe a systematic increase in the rotational modulation towards 2020b, meaning that the axisymmetry level of the field has likely decreased (see Fig. 2). For a first quantitative evaluation, we followed Stibbs (1950) and Preston (1967) to model the phase variations of the longitudinal field for a predominantly-dipolar magnetic configuration. Formally,

$$B_t[G] = \frac{15 + \varepsilon}{20(3 - \varepsilon)} B_p(\cos\beta \cos i + \sin\beta \sin i \cos(2\pi p)), \quad (3)$$

with $\varepsilon$ the limb darkening coefficient (set to 0.3; Claret & Bloemen 2011), $p$ the rotational phase, $B_p$ the longitudinal field of the dipole, $i$ the stellar inclination and $\beta$ the obliquity between magnetic and rotation axes. The results are listed in Table 1, for both near-infrared and optical time series for completeness. The six optical observations in November 2019 have poor coverage (three of them are clustered around phase 0.9) and lead to a less reliable sine fit. Nevertheless, they are compatible with the 2019b fit curve. These clues clearly indicate that the magnetic field of AD Leo is evolving, in agreement with Lavail et al. (2018), and demonstrate the interest of long-term spectropolarimetric monitoring of active M dwarf stars.

## 3.2. The mean line width

The width of near-infrared spectral lines of stars with intense fields and low equatorial velocity such as AD Leo ($v_e \sin(i) = 3$ km s$^{-1}$; Morin et al. 2008b) are sensitive to the Zeeman effect, given its proportionality to wavelength, field strength, and Landé factor. The rotationally-modulated line broadening correlates with the azimuthal distribution of the unsigned small-scale magnetic flux, a useful diagnostic for stellar activity radial velocity contamination, as shown for the Sun by Haywood et al. (2022). In this context, Klein et al. (2021) adopted a selection of magnetically sensitive lines for the young star AU Mic and saw a correspondence in the variations of RV and FWHM of the Stokes I profiles at the stellar rotation period. This confirmed the sensitivity of the FWHM to the distortions induced by magnetic regions on the stellar surface.

Here, we proceeded analogously in an attempt to connect modulations of the FWHM with variations of the large-scale field. We applied LSD on the near-infrared data using a mask of 417 lines characterised by $g_{eff} > 1.2$, following Bellotti et al. (2022). The near-infrared time series was divided in four epochs as in Sec. 3.1 for consistency.

In Table 2, we compare the phase variations of the FWHM when adopting the default and high-$g_{eff}$ masks, and we inspect whether they are more compatible with a sine fit or a constant fit equal to the mean of the data set. In all cases, there is no clear rotational modulation of the data points, as the sine fit does
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Phase variations of the FWHM, folded at a stellar rotation period of 2.23 days. The data sets are computed with the default (black), high-g<sub>eff</sub> (yellow), and low-g<sub>eff</sub> (purple) mask for three epochs: 2019a (top left), 2019b (top right), 2020a (bottom left), 2020b (bottom right). Although no evident rotational modulation is observed, we notice a long-term oscillation of the mean value in a similar way to the longitudinal field (absolute) values.</figcaption></figure>

not provide a better description (i.e. lower $\chi_r^2$) of the variations than the constant fit. This is confirmed by a quick inspection of the periodogram applied to the FWHM data for each individual epoch. The $\chi_r^2$ increase when using a sine model rather than a constant is not statistically significant. The observed variations are attributable to dispersion, as illustrated in Fig. 3. We observed that the FWHM is systematically larger in all epochs for the high-g<sub>eff</sub> mask, as expected given the linear dependence of Zeeman effect to g<sub>eff</sub>, and the dispersion is between 1.8 and 3.0 times larger. The lack of rotational modulation prevents us from searching for correlations with other quantities such as RV and B<sub>I</sub> as done in Klein et al. (2021).

From Fig. 3, we also noticed an evident long-term evolution of the mean FWHM. Such evolution has a moderate correlation (Pearson R coefficient of 0.5) with the variations of the mean B<sub>I</sub> for the same epochs, meaning that the FWHM is a reasonable proxy to trace long-term evolution of the field. This is consistent with the recent Donati et al. (2023) analysis analysis of AU Mic. When using the default mask, the mean FWHM oscillated from 19 km s<sup>-1</sup> in 2019a to 21 km s<sup>-1</sup> in 2019b and 2020a, and back to 19 km s<sup>-1</sup> in 2020b. As expected, such oscillation is enhanced when considering the magnetically sensitive lines and goes from 23 km s<sup>-1</sup> in 2019a to 27 km s<sup>-1</sup> in 2019b and 2020a, and back to 22 km s<sup>-1</sup> in 2020b. We performed the same analysis with low-Landé factor lines (i.e. g<sub>eff</sub> < 1.2 and 406 lines) and noticed no appreciable variation of the mean FWHM, since it remained stable at ~15 km s<sup>-1</sup>. A view of the Stokes I profiles computed with the three different line lists can be found in Appendix A.

The FWHM analysis was also carried out on the ESPaDOnS and Narval data between 2006 and 2019. When using low-g<sub>eff</sub> lines, the mean width of Stokes I is reasonably stable around 9.7 km s<sup>-1</sup>, stressing their potential for precise radial velocity measurements. The full (high-g<sub>eff</sub>) mask yields a mean value at 10 km s<sup>-1</sup> (12 km s<sup>-1</sup>) between 2006 and 2012, which then increases to 11 km s<sup>-1</sup> (13 km s<sup>-1</sup>) in 2016 and 2019. Such long-term evolution is only moderate compared to the one seen in the near-infrared time series. The entire evolution is illustrated in Fig. 4.

The difference between the mean FWHM of low-g<sub>eff</sub> lines in optical (~9.5 km s<sup>-1</sup>) and near-infrared (~16 km s<sup>-1</sup>) can be attributed to lines that have non-zero Landé factor. Indeed, the quadratic differential broadening between the two domains is 11.4 km s<sup>-1</sup>, corresponding to a total magnetic field of 2.5 kG for a line at 1700 nm with g<sub>eff</sub>=0.96 (the normalisation values of the low-g<sub>eff</sub> mask). Although we assumed that the Zeeman effect for low-g<sub>eff</sub> lines is negligible in the optical with this exercise, the inferred value of total magnetic field is reasonably consistent with what is reported in the literature (Saar 1994; Shulyak et al. 2017, 2019), indicating that the magnetic field accounts mostly for the difference in width between optical and near-infrared low-g<sub>eff</sub> lines.

Our analysis confirms that the FWHM is capable of tracing secular changes in the total, unsigned magnetic field, which could be used to better understand stellar activity jitter. Activity-mitigating techniques would benefit from this information even for low-inclination stars such as AD Leo, for which the phase modulation of the radial velocity jitter is more difficult to con-
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Secular evolution of the epoch-averaged FWHM of Stokes *I*. The scale on the left refers to optical data, while on the right to the near-infrared, and the colours indicate the use of different masks similar to Fig. 3. Both optical (circles) and near-infrared (pluses) observations feature a long-term variation, which is enhanced (quenched) when using magnetically sensitive (insensitive) lines. We note that the optical data point of 2019 falls behind the near-infrared ones when using magnetically sensitive lines.</figcaption></figure>

strain. At the same time, the analysis highlights the presence of short-term variability producing scatter and that is not rotationally modulated.

### 3.3. Modelling Zeeman broadening

To further investigate the small-scale magnetic field of AD Leo we conducted a Zeeman broadening analysis. For this analysis we used the full set of new and archival data, both in the near-infrared from SPIRou and in the optical from ESPaDOnS and Narval. All the data sets require a telluric correction, since telluric lines are present in much of the SPIRou wavelength range, and the red end of the ESPaDOnS and Narval range. For the SPIRou data we relied on the telluric correction from the APERO pipeline (Sec 2.1, Cook et al. 2022 for more detail). For the ESPaDOnS and Narval data, we made a telluric correction using the molecfit<fnref n="2" /> pipeline, originally designed for handling spectra from ESO instruments (Smette et al. 2015; Kausch et al. 2015). molecfit retrieves weather conditions and other relevant information at the time of observation and models the atmosphere in the line of sight. It performs radiative transfer and iteratively models the telluric component in the input spectrum while also fitting the continuum and the wavelength scale of the spectrum. It finally corrects telluric lines and provides a telluric-corrected output spectrum. After telluric correction the spectra were renormalised in the regions of interest using a low order polynomial fit through carefully selected continuum regions. A few ESPaDOnS and Narval spectra were affected by fringing effects, hence we adopted a higher-order polynomial fit to normalise to a flatter continuum. Finally, we discarded any observations where the telluric correction left a noticeable residual feature that was blended with the stellar lines of interest.

To characterise the magnetic field, we fitted synthetic spectra to the observed Stokes *I* spectra, incorporating both the Zeeman broadening and intensification effects. Synthetic spectra were calculated with ZEEMAN (Landstreet 1988; Wade et al. 2001; Folsom et al. 2016), using model atmospheres from MARCS (Gustafsson et al. 2008). ZEEMAN performs polarised radiative transfer including the Zeeman effect. However, a major limitation for M-dwarfs is that the programme does not currently include molecular lines, which are not typically used in Zeeman broadening analyses. Weak molecular lines are blended with many atomic lines in the spectra of M-dwarfs. With careful attention we identified a set of atomic lines suitable for AD Leo, with no evident distortion in the line shape by molecular blends. Thus the systematic error from this limitation is expected to be negligible, but the inclusion of molecular lines in the future would substantially simplify the selection of lines for Zeeman broadening analyses, as shown in the recent work of Cristofari et al. (2023) and previously applied to AD Leo by Shulyak et al. (2014, 2017). To check the validity of the analysis presented here, a second analysis of the ESPaDOnS and Narval spectra was carried out with the SYNAST code (Kochukhov et al. 2010). The analyses used nearly the same set of Ti I lines, and the results we obtained were consistent within uncertainty.

For the ESPaDOnS and Narval observations, we used the Ti I lines at 9675.54 Å (g<sub>eff</sub> = 1.35), 9688.87 Å (g<sub>eff</sub> = 1.50), 9705.66 Å (g<sub>eff</sub> = 1.26), 9728.40 Å (g<sub>eff</sub> = 1.00), 9743.61 Å (g<sub>eff</sub> = 0.00), and 9770.30 Å (g<sub>eff</sub> = 1.55). These lines have been used extensively for Zeeman broadening analysis (e.g., Kochukhov & Lavail 2017; Hill et al. 2019; Kochukhov 2021, and references therein) and have reliable oscillator strengths and Landé factors in VALD. These lines have relatively weak telluric blending, very little molecular blending, and a wide range of effective Landé factors.

For the SPIRou observations, we selected a set of lines using similar criteria, but also avoided lines with large pressure broadened wings, since small errors in the pressure broadening could cause larger errors in the Zeeman broadening estimation. In order to maximise the range of available effective Landé factors we used the Fe I lines at 11422.32 Å (g<sub>eff</sub> = 1.98), 11593.59 Å (g<sub>eff</sub> = 2.50), 11607.57 Å (g<sub>eff</sub> = 1.66), 11638.26 Å (g<sub>eff</sub> = 1.58), and 11783.26 Å (g<sub>eff</sub> = 1.14) and the Ti I lines at 11892.88 Å (g<sub>eff</sub> = 0.75), 12821.67 Å (g<sub>eff</sub> = 1.26), 12831.44 Å (g<sub>eff</sub> = 0.67), 12847.03 Å (g<sub>eff</sub> = 1.08), 22232.84 Å (g<sub>eff</sub> = 1.66), and 22310.61 Å (g<sub>eff</sub> = 2.50). This provides multiple lines with both high and low effective Landé factors, but uses lines from two different ions, which we compensated for by using the Ti and Fe abundances as independent free parameters in our analysis. There are a few other Ti I lines near 22000 Å with large effective Landé factors, but there is a relatively severe blending by many weak molecular lines in this region, hence we did not include these lines. Line data were extracted from VALD. In these line lists, experimental oscillator strengths for Ti I lines were from Lawler et al. (2013) (consistent with Blackwell-Whitehead et al. 2006), except for 22232.84 Å from Blackwell-Whitehead et al. (2006), and a theoretical value for 22310.61 Å from the compilation of R. L. Kurucz<fnref n="3" />. Oscillator strengths for the Fe I lines were taken from O'Brian et al. (1991).

The total magnetic field was modelled with a grid of field strengths and filling factors for the fraction of the surface area with the corresponding field strength (e.g. Johns-Krull et al. 1999). A uniform radial orientation was assumed for the magnetic field, since Stokes *I* spectra have little sensitivity to magnetic field orientation. This is also a reasonable assumption given the magnetic field maps reconstructed in Sec. 3.4. For the optical spectra, we adopted magnetic fields of 0, 2, 4, 6, 8, and 10 kG, and derived their filling factors. For the SPIRou spectra, we used

<footnote n="2">² https://www.eso.org/sci/software/pipelines/</footnote>

<footnote n="3">³ http://kurucz.harvard.edu</footnote>
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Example fits including Zeeman broadening resulting from the MCMC-based approach of radiative transfer modelling. Top: ESPaDOnS observation from 24 February 2016. Bottom: SPIRou observation from 3 November 2020. The panels have the same wavelength scale with minor ticks at 0.1 nm.</figcaption></figure>

a finer grid of 1 kG from 0 to 10 kG, since the sensitivity to Zeeman effect is larger at longer wavelengths, and a finer grid is needed to produce smooth line profiles.

To derive the magnetic filling factors we applied an MCMC-based approach, using the `emcee` package (Foreman-Mackey et al. 2013) integrated with `Zeeman`. The filling factors for $B > 0$ were treated as free parameters, with the filling factor for $B = 0$ ($f_{B=0}$) calculated from $1 - \sum_{B>0} f_B$. Proposed steps in the chain where $\sum_{B>0} f_B > 1$ were rejected to ensure that the filling factors sum to unity. The projected rotational velocity $v_e \sin(i)$ and the abundance of Ti (and Fe for SPIRou) were included as free parameters in the MCMC process. The modelling used $T_{\text{eff}} = 3500$ K, $\log g = 5.0$, a microturbulence of $1 \text{ km s}^{-1}$. The chemical abundances may be unreliable since they do not account for elements bound in molecules, making them effectively nuisance parameters in this study. However, this provides the code with flexibility for fitting line strength and width in the absence of a magnetic field, reducing the sensitivity of the results to small errors in non-magnetic parameters.

Example fits resulting from the MCMC-based approach are shown in Figs. 5 for ESPaDOnS and SPIRou. The shapes of the posterior distributions are generally similar for all observations using the same sets of lines, and are illustrated in Appendix C. There are important anti-correlations between filling factors with adjacent magnetic field strengths, and weak correlations between filling factors spaced by two bins in field strength. Therefore,

some caution should be taken in interpreting the uncertainties from this and similar analyses. The filling factor for $B = 0$ and the quantity $\sum_i B_i f_i$ summed over magnetic field bins (abbreviated to $\sum Bf$), were calculated from samples in the MCMC chain. The resulting distribution was used to provide the median value from the 50th percentile, with uncertainties from the 16th and 84th percentile.

The results for all observations, and averages for each epoch, are presented in Fig. 6, and values for each epoch are provided in Tables C.1 and C.2. The quantity $\sum Bf$ (sometimes called the magnetic flux, and analogous to a magnetic flux density) ranges between 2.6 kG and 3.7 kG, which is consistent with previous measurements (e.g., Saar 1994; Kochukhov et al. 2009; Shulyak et al. 2014, 2017, 2019; Cristofari et al. 2023). We observe a long-term increase of the average $\Sigma Bf$ from 2.8 kG in 2007 to 3.6 kG in 2016, followed by a weakening towards 3.4 kG with the latest SPIRou observations. Such behaviour correlates with the long-term decrease (in absolute value) of the longitudinal field (Pearson coefficient $R=0.6$, excluding the 2006 data point). Likewise, the average $\Sigma Bf$ time series correlates with the average FWHM of Stokes *I*, demonstrating its capability at tracing the evolution of the total, unsigned magnetic field (Donati et al. 2023).

The $\Sigma Bf$ values for the ESPaDOnS optical data acquired in 2006 fall out from this trend. This could stem from residuals of the telluric correction blending with the lines used in the
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Long-term evolution of the magnetic flux density measured from Zeeman broadening modelling. The data point format is the same as Fig. 2, and for each epoch we overplotted the mean value and standard deviation as error bar. There is an overall trend that reflects the variations of $B_I$ across both optical and near-infrared data. In addition, for the near-infrared SPIRou time series, we also observe the same oscillating behaviour found for the FWHM of Stokes *I*.</figcaption></figure>

modelling and/or instrumental effects such as fringing, for which the results are sensitive to the choice of continuum normalisation. Attempts were made to correct for these potential systematic errors: rejecting observations where the telluric correction left residual features in the used portion of the spectrum, and careful continuum normalisation to remove any weak fringing. However, it is possible these attempts were not fully successful, and thus the departure from the general trend of the 2006 result should be treated with caution.

### 3.4. Magnetic imaging

We applied ZDI to the SPIRou and 2019 ESPaDOnS time series of Stokes *V* profiles to recover the large-scale magnetic field at the surface of AD Leo. The magnetic geometry is modelled as the sum of a poloidal and a toroidal component, which are both expressed through spherical harmonics decomposition (Donati et al. 2006; Lehmann & Donati 2022). The algorithm compares observed and synthetic Stokes *V* profiles iteratively, fitting the spherical harmonics coefficients $\alpha_{\ell,m}$, $\beta_{\ell,m}$, and $\gamma_{\ell,m}$ (with $\ell$ and *m* the degree and order of the mode, respectively), until they match within a target reduced $\chi^2$. Because the inversion problem is ill-posed, a maximum-entropy regularisation scheme is applied to obtain the field map compatible with the data and with the lowest information content (for more details see Skilling & Bryan 1984; Donati & Brown 1997; Folsom et al. 2018).

In practice, we used the zdipy code described in Folsom et al. (2018). In its initial version, the code performed tomographic inversion under weak-field approximation, for which Stokes *V* is proportional to the first derivative of Stokes *I* over velocity (Landi Degl'Innocenti & Landolfi 2004). For the present study, we have implemented the Unno-Rachkovsky's solutions to polarised radiative transfer equations in a Milne-Eddington atmosphere (Unno 1956; Rachkovsky 1967; Landi Degl'Innocenti & Landolfi 2004) and incorporated the filling factor formalism outlined in Morin et al. (2008b) and Donati et al. (2023). The implementation of Unno-Rachkovsky's solutions was motivated by the need of a more general model for the observed Stokes *V* profiles. Near-infrared observations of stars

with intense magnetic fields are indeed more susceptible to distortions and broadening due to an enhanced Zeeman effect.

As input parameters for ZDI, we assumed $i = 20^\circ$, $v_e \sin(i) = 3$ km s$^{-1}$, $P_{\text{rot}} = 2.23$ days, and solid body rotation. We adopted a linear limb darkening coefficient in H band of 0.3 and V band of 0.7 (Claret & Bloemen 2011). We set the maximum degree of the harmonic expansion $\ell_{\text{max}} = 8$ (considering the low $v_e \sin(i)$) and allowed an entropy weighting scheme proportional to $\ell$ during ZDI inversion, to favour simple geometries as in Lavail et al. (2018). The SPIRou near-infrared time series was split similarly to Sec. 3.1: 2019a (21 observations over 30 cycles), 2019b (21 observations over 26 cycles) 2020a (30 observations over 20 cycles), and 2020b (18 observations over 15 cycles). The Stokes *V* time series of SPIRou and 2019 ESPaDOnS data are shown in Fig. 7.

For the 2019a, 2019b, 2020a, and 2020b epochs, we fitted the Stokes *V* profiles to a $\chi_r^2$ level of 1.2, 1.0, 1.1, and 1.1 from an initial value of 10.3, 15.8, 14.7, and 8.5, respectively. For the ESPaDOnS 2019 time series, we fitted down to $\chi_r^2 = 2.5$ from an initial value of 156.3. We attempted to merge the 2019b and 2020a epochs and reconstruct a single map, since they are separated by the shortest time gap. The quality of the final model is deteriorated ($\chi_r^2=1.3$) with respect to the two epochs separately, but the corresponding map and magnetic energies are consistently recovered. We therefore kept these two epochs separate.

From Fig. 7, it is evident that the near-infrared Stokes *V* profiles manifest structures and stochastic variability in both lobes. This is not extended in the continuum, since the residuals with respect to the mean profile are compatible with the noise level, it is not rotationally-modulated, and it is not exhibited by Stokes *N*. The presence of such variability was already suggested by the phase-folded variations in $B_I$ (Fig. 2), as some data points featured a departure from a pure rotational modulation. Likewise, the residuals of the Stokes *I* profiles show clear variability, but the application of a 2D periodogram does not reveal any significant periodicity. While our ZDI model is capable of describing the general shape of Stokes *V* profiles, it is limited at reproducing the structures and at capturing all the information present. These considerations are also valid for optical observations in 2019, as the amplitude of Stokes *V* is not matched exactly by our ZDI model, and overall translate into an underestimate of the field strength. This motivates further the use of the PCA method described by Lehmann & Donati (2022), which is a data-driven approach offering a complementary view on the magnetic field evolution, as outlined in the next section.

We are able to constrain the filling factor $f_V$ following a $\chi^2$ minimisation prescription similar to (Petit et al. 2002). We found $f_V$ values oscillating between 9% in 2019a, 16% in 2019b and back to ~11% in the remaining epochs, compatible with Morin et al. (2008b), and larger by a factor of 1.7 than Lavail et al. (2018). This would indicate a weakening of the local small-scale field since 2016, on top of a decrease in large-scale field intensity as seen in the reconstructions (Fig. 8).

The filling factor $f_I$ was inspected by considering a grid of values between 0% and 100%; for each $f_I$ value, we synthesised a time series of model Stokes *I* profiles, computed the corresponding time series of $\chi_r^2$ with the observations, and phase-folded the $\chi_r^2$ curve at $P_{\text{rot}}$. We then assessed at what value of $f_I$ the $\chi_r^2$ curve would start manifesting rotational modulation, because it would indicate that certain model profiles deviate from the observations. We noticed that values above 30% deteriorate the fit of the profile core progressively, yielding variability and rotational modulation of the Stokes *I* profiles, which is not observed otherwise (see Sec. 3.2). Values of $f_I = 30\%$ are three
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Full SPIRou and ESPaDOnS 2019 time series of Stokes V profiles of AD Leo normalised by the unpolarised continuum intensity. Two panels are present for 2019a (top left), 2019b (top right), one panel for 2019 in optical (bottom left) and 2020a (bottom middle) and two panels for 2020b (bottom right). In all panels, both the observed profiles (black) and Unno-Rachkovsky models (red) are shown. They are shifted vertically for visualisation purposes, and their associated rotational cycle is reported on the right (see Eq. 1). A remarkable feature is the increased intermittency of the profiles amplitude for the most recent observations, for example at cycle 174.44 and 174.95 of 2020b against cycle 1.40 and 1.88 in 2019a, suggesting an evolution of the magnetic field obliquity.</figcaption></figure>

times larger than $f_V$, in agreement with (Morin et al. 2008b).

Since the plausible $f_I$ values are consistent with 0%, we adopted $f_I = 0\%$ in the ZDI modelling.
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Reconstructed ZDI maps in flattened polar view of AD Leo for the SPIRou and 2019 ESPaDOnS time series. The columns correspond to the epochs examined, from the left: 2019a, 2019b, 2019 optical, 2020a, and 2020b. In each column, the radial (top), azimuthal (middle) and meridional (bottom) components of the magnetic field vector are displayed. The radial ticks are located at the rotational phases when the observations were collected, while the concentric circles represent different stellar latitudes: -30°, +30° and +60° (dashed lines) and equator (solid line). The colour bar range is set by the maximum (in absolute value) of the magnetic field and illustrates the positive (red) and negative (blue) polarity for each epoch.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Evolution of the magnetic field obliquity with respect to the stellar rotation axis. Values estimated from Eq. 3 (Stibbs 1950) using optical and near-infrared data are illustrated as circles and pluses, colour-coded by epoch, whereas values from ZDI are shown as diamonds and crosses. The y axis is expressed in degrees of colatitude.</figcaption></figure>

The five maps of surface magnetic flux (one for each SPIRou epoch, and one for the ESPaDOnS 2019 epoch) are shown in Figure 8 and their properties are reported in Table 3. In all cases, the configuration is predominantly poloidal, storing >95% of the magnetic energy. The main modes are dipolar and quadrupolar, as they account for 70-90% and 15-20% of the magnetic energy. We report a weakening of the mean field strength ($\langle\langle B\rangle\rangle$) of factor of 1.5 and 2.4 relative to the optical maps reconstructed by Morin et al. (2008b) and Lavail et al. (2018), respectively. The most remarkable feature is the reduction of magnetic energy contained in the axisymmetric mode, going from > 99% in 2019a to 60% in

2020b, translating into an increase of the dipole obliquity relative to the rotation axis, from 3° to 38°.

We note that the maximum field strength reconstructed with ZDI is between 1.2 and 2.4 times smaller than obtained via Eq. 3 (Stibbs 1950). Likewise, the magnetic field obliquity is underestimated, as illustrated in Fig. 9. On one side, this difference stems from the limitation of the Stokes V ZDI model, since it does not encompass the full amplitude of the two lobes for some observations, and on the other side Eq. 3 assumes a purely dipolar field, contrarily to our reconstructions (the dipole accounts for 70-90% of the energy). Nevertheless, both approaches allow us to observe an evident evolution of the obliquity, featuring a rapid increase in the most recent epochs.

Finally, we merged the 2019b, 2020a and 2020b data sets and attempted a joint rotation period and differential rotation search following Petit et al. (2002). The results were inconclusive, likely due to the significant evolution of the surface magnetic field between each epoch.

The summary of the magnetic field's evolution is illustrated in Fig. 10. We performed ZDI reconstructions also for the archival ESPaDOnS and Narval data for consistency, finding reasonably compatible results with previous studies (Morin et al. 2008b; Lavail et al. 2018). We observe a globally simple geometry (i.e. predominantly poloidal and dipolar) over 14 yr, with a decreasing strength. Our latest SPIRou observations revealed a clear evolution of the dipole obliquity in the form of a reduced axisymmetry, suggesting a potential dynamo magnetic cycle. These features are indeed compatible with the variations observed by Sanderson et al. (2003) and Lehmann et al. (2021) for the solar cycle.
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Evolution of the magnetic topology of AD Leo over 14 yr. The time series includes optical ESPaDOnS and Narval data between 2007 and 2019, and near-infrared SPIRou data collected between 2019 and 2020. Field strength, level of axisymmetry and dominant field component (poloidal or toroidal) are encoded with the symbol size, shape, and colour, respectively, and were computed using the Unno-Rachkovsky implementation described in this work for consistency. The topology remained predominantly poloidal and dipolar, while the field strength has decreased over time. The most striking change is a decrease in axisymmetry in the most recent epoch.</figcaption></figure>

<table><caption>**Table 3.** Properties of the magnetic maps for the 2019a, 2019b, 2020a, and 2020b SPIRou epochs and the 2019 ESPaDOnS data set.</caption><thead><tr><th></th><th>2019a<br>NIR</th><th>2019b<br>NIR</th><th>2019<br>VIS</th><th>2020a<br>NIR</th><th>2020b<br>NIR</th></tr></thead><tbody><tr><td>f<sub>V</sub></td><td>9%</td><td>16%</td><td>12%</td><td>12%</td><td>11%</td></tr><tr><td>⟨B⟩ [G]</td><td>111.2</td><td>132.3</td><td>158.0</td><td>115.3</td><td>93.4</td></tr><tr><td>B<sub>max</sub> [G]</td><td>481.2</td><td>764.0</td><td>577.2</td><td>555.1</td><td>434.3</td></tr><tr><td>B<sub>pol</sub> [%]</td><td>100.0</td><td>99.9</td><td>95.0</td><td>99.3</td><td>98.7</td></tr><tr><td>B<sub>tor</sub> [%]</td><td>0.0</td><td>0.1</td><td>5.0</td><td>0.7</td><td>1.3</td></tr><tr><td>B<sub>dip</sub> [%]</td><td>81.3</td><td>71.1</td><td>81.7</td><td>75.7</td><td>70.1</td></tr><tr><td>B<sub>quad</sub> [%]</td><td>14.9</td><td>19.0</td><td>14.6</td><td>17.7</td><td>21.2</td></tr><tr><td>B<sub>oct</sub> [%]</td><td>2.8</td><td>6.2</td><td>2.7</td><td>4.4</td><td>5.9</td></tr><tr><td>B<sub>axisym</sub> [%]</td><td>99.8</td><td>94.5</td><td>77.0</td><td>85.8</td><td>58.3</td></tr><tr><td>Obliquity [°]</td><td>2.5</td><td>12.5</td><td>21.5</td><td>19.5</td><td>38.0</td></tr></tbody></table>

**Notes.** The following quantities are listed: filling factor, mean magnetic strength, maximum magnetic strength, poloidal and toroidal magnetic energy as a fraction of the total energy, dipolar, quadrupolar and octupolar magnetic energy as a fraction of the poloidal energy, axisymmetric magnetic energy as a fraction of the total energy, and tilt of the magnetic axis relative to the rotation axis. The time span of the 2019 optical epoch is encompassed in the 2019b near-infrared epoch.

### 3.5. Diagnosing the large-scale field using PCA

AD Leo is an ideal target for analysing large-scale field evolution with the data-driven PCA method recently presented by Lehmann & Donati (2022), given its magnetic field strength and $v_e \sin(i)$. Principal component analysis allows us to uncover details about the stellar large-scale field directly from the LSD Stokes V profiles and to trace its magnetic field evolution across the observation run, without prior assumptions. Here, we analyse only the near-infrared time series, because the number of optical 2019 observations is not sufficient.

First, we can get insights about the star's axisymmetric large-scale field by analysing the mean Stokes V profile determined over all Stokes V LSD profiles (see Lehmann & Donati (2022) for further details). Fig. 11 displays the mean profile and the decomposition into its antisymmetric and symmetric parts, denoting the poloidal and toroidal axisymmetric components, respectively. We clearly see that the mean profile is antisymmetric, which indicates a poloidal-dominated axisymmetric large-scale

field. The amplitude of the symmetric part is comparable to the noise, and likely due to an artefact of uneven phase coverage rather than a true toroidal field signal (Lehmann & Donati 2022). Compared to the mean-subtracted Stokes V profiles, the amplitude of the mean profile is generally strong, marking a dominant axisymmetric field. However, we observe an increase in the amplitude of the mean-subtracted Stokes V in the last two epochs 2020a and 2020b, which provides a first hint towards a less axisymmetric configuration.

Second, the application of PCA to the mean-subtracted Stokes V profiles yields insights on the non-axisymmetric field, (Lehmann & Donati 2022). For the mean-subtracted Stokes V profiles, we applied the mean profile computed across all epochs, which allows a direct reflection of the epoch-to-epoch variations in PCA coefficients (e.g. in amplitude and mean value). If the mean Stokes V profile were computed per epoch, we would miss such information, that is to say the mean value of the coefficients would be centred for each epoch, and the amplitudes could no longer be compared to each other. Fig. 12 presents the first three eigenvectors and their corresponding coefficients for the mean-subtracted Stokes V profiles separated by epoch and colour-coded by rotation cycle. The first eigenvector displays an antisymmetric shape proportional to the first derivative of the Stokes I profile and, together with the associated coefficient, scales mainly with the longitudinal magnetic field (Landi Degl'Innocenti & Landolfi 2004). The second eigenvector shows a more symmetric shape, more closely related to the second derivative of the Stokes I profile, and describes the temporal evolution of the Stokes V profiles between the maxima of the longitudinal field. According to Lehmann & Donati (2022), a strongly antisymmetric eigenvector traces the radial component and a symmetric eigenvector the azimuthal component for a dipole dominated field that is strongly poloidal, which is the case for AD Leo. The third eigenvector features a signal as well (antisymmetric, and related to the third derivative of the Stokes I profile), which is detectable due to the high S/N of the data set, while the further eigenvectors are dominated by noise. Seeing three eigenvectors indicates that even if the axisymmetric field is likely to be dominant, we are able to detect and to analyse the non-axisymmetric field in great detail.

The coefficients of the eigenvectors suggest an evolving large-scale field as their trend changes for every epoch, see
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Decomposition of the Stokes V mean profile of AD Leo SPIRou observations. The mean profile (solid red line) and its antisymmetric (dash dotted blue line) and symmetric (dashed yellow line) parts, related to the poloidal and toroidal components of the axisymmetric large-scale field, are shown.</figcaption></figure>

Fig. 12 2nd-5th row. In 2019a, the coefficients related to the first eigenvector show only a flat distribution around zero, implying a predominantly axisymmetric field. For the following epochs, 2019b, 2020a and 2020b, we see sine-like trends of the first two coefficients with rotational phase. The amplitude increases from epoch to epoch, indicating a growing obliquity of the dipole-dominated large-scale field. For the 2020b epoch, the obliquity becomes so large that the coefficients of the third eigenvector start to show a sine-like trend as well, which translates into a significant non-axisymmetric field.

Furthermore, the extremes of the coefficients associated to the antisymmetric and symmetric profile (first and second eigenvector) for the same epoch feature an apparent phase shift of ≈ 0.25, which demonstrates that the dipolar component is poloidal dominated with little to no toroidal contribution (Lehmann & Donati 2022). The extremes of the coefficients related to the antisymmetric eigenvector locate the pointing phase of the dipole (Lehmann & Donati 2022). For the last three epochs, the maximum of this coefficient occurs at a pointing phase of ≈ 0.3 for the northern pole of the dipole, and the sign of the eigenvector implies a negative polarity. The extremes of the coefficients occur at the same rotational phase throughout the whole observation run, designating a stable pointing phase of the dipole, in agreement with the $B_\phi$ measurements (see middle panel of Fig. 2).

By applying the PCA method on the time series of Stokes V (Lehmann & Donati 2022), we confirm that AD Leo features a dipolar large-scale field, whose obliquity increased during the latest epochs (2020a and 2020b). As the large-scale field became more non-axisymmetric, the pointing phase of the dipole remained stable.

## 4. Achromaticity of the magnetic field

The impact of stellar magnetic activity on radial velocity measurements features a chromatic dependence stemming from a combination of magnetic field and spot temperature contrast (Reiners et al. 2013; Baroch et al. 2020). Indeed, at near-infrared wavelengths the Zeeman broadening is expected to be stronger, while starspots contribute less owing to a lower contrast with the photosphere. The situation is reversed in the optical domain. For

<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Principal component analysis applied to the mean-subtracted Stokes V profiles. The first three eigenvectors (top row) and their corresponding coefficients (2nd-4th row) are shown. The vertical line in the eigenvector panels indicates the radial velocity offset for AD Leo, of approximately 12.4 km s⁻¹, while the horizontal dotted line in all panels represents the null line. The coefficients are displayed by epoch: 2019a (2nd row), 2019b (3rd row), 2020a (4th row) and 2020b (5th row). They are phase folded using Eq. 1 and colour-coded by rotation cycle.</figcaption></figure>

AD Leo, recent work by Carmona et al. (2023) demonstrated the strong chromatic behaviour of radial velocity jitter, the latter being significantly weaker in the near-infrared domain than in optical. The combination of these effects becomes increasingly important with the activity level of the star, since the number of spots would be correspondingly larger (Reinhold et al. 2019), and it could possibly result in distinct contributions to the magnetic field strength, which can then be used to facilitate the modelling of stellar activity.

Fast-rotating stars are expected to feature high active latitudes and large polar spots (e.g., Cang et al. 2020), because the Coriolis force would overcome the buoyancy force, making the flux tubes ascend parallel to the stellar rotation axis (Schuessler & Solanki 1992; Granzer et al. 2000). There are some cases, however, in which fast rotation does not correlate with the presence of a polar spot (Barnes et al. 2004; Morin et al. 2008a). The fact that AD Leo is a moderate rotator observed nearly pole-on makes it an interesting case to investigate whether longitudinal
field measurements are chromatic, reflecting the behaviour of an underlying spot.

Previous studies dedicated to the Sun have shown that the magnetic field strength measured in individual lines varies significantly (Demidov et al. 2008; Demidov & Balthasar 2012), and differences between optical and near-infrared domains have unveiled a dependence of the field strength on atmospheric height: the field increases while going towards deeper internal layers (Zayer et al. 1989; Solanki 1993). For other stars, Valenti et al. (1995) reported a chromatic difference in magnetic field strength for the moderately active K dwarf ε Eri, but attributed its origin to incomplete modelling of the spectral lines used for the Zeeman broadening analysis. No wavelength dependence of the field strength was reported more recently, neither for ε Eri (Petit et al. 2021) nor for T Tauri stars (Finociety et al. 2021). The same conclusion was reached by Bellotti et al. (2022) when computing longitudinal field values for the active M dwarf EV Lac using blue (< 550 nm) and red (> 550 nm) lines of an optical line list.

To investigate the longitudinal field chromaticity, we analyse the contemporaneous observations taken with SPIRou and ESPaDOnS in November 2019. We limit the LSD computation within successive wavelength bins of the line mask, and evaluate the longitudinal field for each case. Including both optical and near-infrared domains, we considered 11 subsets of lines in the following ranges: [350,390], [390,430], [430,480], [480,550], [550,650], [650,1100], [950,1100], [1100,1400], [1400,1600], [1600,1800], [1800,2500] nm. The [650,1100] and [950, 1100] nm ranges represent the red end of ESPaDOnS spectra and the blue end of SPIRou spectra, respectively. We adopt more wavelength regions than those presented in Bellotti et al. (2022), allowing a finer search of chromatic trends. The number of lines used varies between 100 and 1000 in the optical, and between 120 and 300 in the near-infrared (see Fig. B.1). In addition, we compute LSD using a 50-lines mask in the overlapping wavelength region of ESPaDOnS and SPIRou spectra ([950,1050] nm).

Stokes *I* and *V* profiles were computed for the simultaneous SPIRou and ESPaDOnS epochs, namely 2019b and 2019, respectively. To increase the S/N and allow a more precise estimate of $B_l$, the profiles obtained with a specific line list subset and belonging to the same epoch were co-added. This is reasonable considering the marginal amplitude variation over the epochs examined and the unchanged polarity of Stokes *V*. The longitudinal field was then computed with Eq. 2 using the specific normalisation wavelength and Landé factor of each line subset, and adapting the velocity integration range according to the width of the co-added Stokes *V* profile.

From Fig. 13, we observe no clear chromaticity of $B_l$. The distribution of field strength is flat around -200 G with a total scatter of 20 G. Such dispersion is mainly due to LSD computations with a low number of lines, implying Stokes shapes more sensitive to variations in individual lines, blends and residuals of telluric correction. For the same reason, some profiles appear deformed and lead to evident outliers (see Fig.B.1). For instance, the $B_l$ value obtained from ESPaDOnS data in the spectral region overlapping with SPIRou is 100 G weaker (in absolute value) than the $B_l$ value obtained from SPIRou data in the same wavelength region. This could be due to the low S/N at the very red edge of ESPaDOnS.

The case of [390,430] nm leads to a field value of -750 G, despite the Stokes profiles do not show a particular deformation. We attribute this behaviour to an imprecise continuum normalisation of the spectra, likely due to a challenging identification

<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Test on $B_l$ measurements chromaticity. Shown are $B_l$ measurements obtained with line mask subsets based on wavelength of the optical (circles) and near-infrared (pluses) domain. The horizontal black lines indicate the values of $B_l$ computed with the full masks and the vertical dashed line separates arbitrarily optical from near-infrared measurements. The error bars are smaller than the symbol size. For visualisation purposes, the [350, 390] nm wavelength bin is not shown since it leads to an outlier data point at around -1 kG, and the [390, 430] nm bin yields a value at -750 G so it is indicated with a downward arrow. Overall, no chromatic trend emerges from the data.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Dependence of the magnetic field strength on limb darkening. A monotonic decreasing (in absolute value) trend is observed from large (optical) to small (near-infrared) coefficient values. This effect is small, and negligible compared to noise in real observations.</figcaption></figure>

of the continuum level in the blue part of the spectrum, where M dwarfs feature forests of spectral lines. The effect is a smaller depth (and equivalent width) of the Stokes *I* profile relative to the other cases, which artificially increases the value of the field (in absolute value). Overall, although the [350,390] and [390,430] nm bins contain more than half of the lines in the optical mask, their weight in the LSD computation is small (Kochukhov et al. 2010) making their effect in the computation of $B_l$ with the full mask negligible.

We repeated the same exercise for the other SPIRou epochs and found a similar behaviour, the only difference being 2020b data points shifting upwards because of the field global weakening. A possible implication of the lack of a chromatic trend may be the absence of a polar spot for AD Leo. This would be justified considering that other faster-rotating M dwarfs like V374 Peg (Morin et al. 2008a) and HK Aqr (Barnes et al. 2004) do not show polar spots.
A potential source of chromaticity for $B_l$ values may come from limb darkening. This radial gradient in stellar brightness over the visible disk can be expressed as a linear function of the angle between the line of sight and the normal to a surface element ($\theta$)

$$ \frac{I}{I_0} = 1 - \varepsilon(1 - \cos \theta), \qquad (4) $$

where $I_0$ is the brightness at disk centre ($\theta = 0^\circ$) and $\varepsilon$ is the limb darkening coefficient. Claret & Bloemen (2011) show that $\varepsilon$ decreases with wavelength, being 0.7 in V band and 0.3 in H band. The linear limb darkening law in Eq. 4 is the one implemented in the ZDI reconstruction (Folsom et al. 2018).

Owing to the stronger limb darkening in optical than in near-infrared, there is the possibility of additional polarity cancellation in the latter domain, which would lead to weaker field measurements. For the specific case of AD Leo, the low stellar inclination makes the equator appear at the limb and near-infrared observations would be more sensitive to this region. In particular, the sign of large-scale dipolar magnetic field lines exiting the pole would cancel out more with those at equator, compared to optical observations.

To verify this, we 1) linearly interpolated the limb darkening coefficients in Claret & Bloemen (2011) at the wavelengths examined for the thermal contrast test (see Fig. 13), 2) synthesised Stokes profiles for the same coefficients assuming an axisymmetric dipole of 1 kG seen pole-on (akin to AD Leo in 2019a) and infinite S/N, and 3) computed the associated field values with Eq. 2. The results are illustrated in Fig. 14. We observe a small (7%) weakening of the field from optical to near-infrared, which is overwhelmed by noise in real observations.

## 5. Discussion and conclusions

In this paper, we presented the results of an extended spectropolarimetric monitoring of the active M dwarf AD Leo, using near-infrared observations collected with SPIRou between 2019 and 2020 as part of the SLS survey. They add to the previous optical data obtained with ESPaDOnS and Narval between 2006 and 2019, making the entire time series encompass approximately 14 yr. To carry out our magnetic analysis, we computed the longitudinal magnetic field, tracked the variations of the Stokes I FWHM, modelled Zeeman broadening on individual selected lines, reconstructed the large-scale field topology via ZDI, and assessed axisymmetry variations by means of a novel PCA method.

Initially, Morin et al. (2008b) reported an axisymmetric, dipole-dominated structure that was stable over one year; later, Lavail et al. (2018) pointed out a large-scale weakening and small-scale enhancement of the field but no variation in the geometry. We found strong evidence of a large-scale field evolution, that is summarised as follows:

1. The longitudinal magnetic field has weakened between 2006 and 2020, from -300 to -50 G, with a rapid decrease of 100 G in the 2020b. The dipolar longitudinal magnetic field evolved in the same time frame, starting from -850 G in 2006, reaching -560 G in 2016 and restoring back to -900 G in 2020.

2. The FWHM of Stokes I profiles does not show rotational modulation, but a dispersion that may partly be due to short-term variability. The epoch-averaged FWHM manifests a long-term variation both in optical and near-infrared, being wider in 2019b and 2020a, and narrower in 2019a and 2020b.

The variations are enhanced when the Stokes profiles are computed with magnetically-sensitive lines, as opposed to the insensitive ones. The near-infrared data in particular feature a trend moderately correlated with $B_l$ (in absolute value).

3. The magnetic flux estimated from the modelling of Zeeman broadening exhibits a global increase over time, which is also correlated to the long-term trend of the longitudinal magnetic field (in absolute value). Moreover the epoch-averaged magnetic flux obtained for the near-infrared SPIRou time series oscillates in a similar manner to the FWHM of Stokes I, demonstrating that the latter is capable of tracing secular evolution of the total, unsigned magnetic field.

4. Zeeman-Doppler imaging reconstructions confirmed the same kind of topological evolution, with the axisymmetric level decreasing to 60% and the obliquity between magnetic and rotation axis increasing to 38°. This already found support by the enhanced intermittency of the amplitude of Stokes V profiles in late 2020.

5. The PCA method confirmed the predominantly poloidal and dipolar geometry of the large-scale field, as well as a lower axisymmetry in 2020a and 2020b. In addition, the pointing phase of the dipole remained stable during the evolution.

6. Measurements of the magnetic field strength are overall achromatic, since they manifest only a marginal wavelength dependence due to limb darkening.

Our results altogether suggest that AD Leo may be entering a polarity reversal phase of a long-term magnetic cycle, analogous to the solar one. The combination of chromospheric activity studies and spectropolarimetric campaigns show that some Sun-like stars may manifest magnetic cycles and polarity reversals in phase with chromospheric cycles (Boro Saikia et al. 2016; Jeffers et al. 2017), while others have a more complex behaviour where very regular chromospheric oscillations have no straightforward polarimetric counterpart (Boro Saikia et al. 2022).

Predicting when the polarity reversal may occur for AD Leo is not a trivial task, as the $B_l$ data set does not feature a clear minimum or maximum. Recently, Fuhrmeister et al. (2023) did not report any evident trends from a long-term campaign of chromospheric indexes, whereas previous studies based on photometric observations reported either two co-existing timescales for cycles, namely 7 yr and 2 yr (Buccino et al. 2014), or an individual one of about 11 yr (Tuomi et al. 2018). However, these time scales are not compatible with the variations in $B_l$ observed over 14 yr. The axisymmetric level of the large-scale topology is a more suitable proxy to track the cycle (Lehmann et al. 2021), but we recorded its change only in the most recent observations.

A comparison between the magnetic field evolution described here and that of the radial velocity jitter obtained in Carmona et al. (2023) leads to a puzzling situation. Carmona et al. (2023) show that radial velocity variations in optical are essentially due to the presence of a spot and that this signal has changed only slightly (in phase and amplitude, the latter varies from $25.6 \pm 0.3$ m s$^{-1}$ to $23.6 \pm 0.5$ m s$^{-1}$) between 2005 and 2021. Such radial velocity signal is not detected in infrared with SPIRou, corroborating its strong chromaticity and therefore its origin due to stellar activity Carmona et al. (2023). The fact that the dipolar field evolution is disjointed from a surface brightness evolution is not a surprise: Morin et al. (2008a) show that the mainly-dipolar topology of V374 Peg did not correlate with the complex brightness map reconstructed via Doppler imaging.

These considerations motivate long-term spectropolarimetric and velocimetric campaigns of active M dwarfs. For AD Leo in
particular, additional monitoring is required to observe the polar-
ity reversal and the cycle's extremes, to constrain a precise time
scale. An extended temporal baseline could also give more in-
sight on the link between topological variations and high-energy
flaring events (e.g., Stelzer et al. 2022). At the same time, we
could shed more light on the relation between the evolution of
the large-scale magnetic field topology and the stability of the
radial velocity jitter.

An additional detail we could infer about AD Leo's mag-
netic field is the helicity, which quantifies the linkage between
poloidal and toroidal field lines and thus describes the complex-
ity of the magnetic topology (Lund et al. 2020, 2021). For the
Sun, Pipin et al. (2019) reported a temporal variation of the value
correlated to the magnetic cycle. Indeed, helicity maxima and
minima occur when the axis of symmetry of the poloidal and
toroidal field components are aligned and orthogonal, respec-
tively.

For AD Leo, the fraction of toroidal energy is only a neg-
ligible fraction of the total one, hence we should exert caution
when deriving quantities from it. Over time, we observe that the
poloidal axisymmetric (m = 0) mode maintains >80% of the
magnetic energy and features a drop to 45% in 2020b, while the
energy in the toroidal axisymmetric mode decreases from 30% to
6%. As a result, the two components maintained an overall mis-
aligned configuration, but in the most recent epoch, the poloidal
component became more aligned with the toroidal one due to
the axisymmetry decrease. Following the practical visualisation
of Lund et al. (2021), this evolution would correspond to an in-
crease in field helicity.

The existence of a magnetic cycle for AD Leo is in agree-
ment with the observational evidence of such phenomena for
M dwarfs from radial velocity exoplanet searches (Gomes da
Silva et al. 2012; Lopez-Santiago et al. 2020). In general, stud-
ies have shown that magnetic cycles introduce long-term sig-
nals in radial velocity data sets that can dominate over planetary
signatures (Meunier et al. 2010; Meunier & Lagrange 2019), as
they modulate the appearance and number of heterogeneities on
the stellar surface. It is therefore necessary to have an accurate
constraint on the temporal variations of the cycle, in order to
remove its contamination and allow a more reliable planetary
detection and characterisation (Lovis et al. 2011; Costes et al.
2021; Sairam & Triaud 2022).

Furthermore, activity cycles modulate the stellar radiation
output and winds in which close-in planets are immersed (Yeo
et al. 2014; Hazra et al. 2020). This leads to a temporal varia-
tion in the planetary atmospheric stripping with consequent al-
teration of the chemical properties and habitability (Lanza 2013;
McCann et al. 2019; Louca et al. 2023; Konings et al. 2022).
Details on the occurrence of the cycle extremes can thus inform
the most suitable interpretation framework and observing plans
for missions dedicated to transmission spectroscopy like *Ariel*
(Tinetti et al. 2021). At the same time, periodic variations in the
large-scale field geometry need to be considered for an accurate
and updated modelling of the low-frequency radio emission dis-
covered for M dwarfs (Callingham et al. 2021), which has been
recently proposed to potentially reveal the presence of close-in
magnetic planets (Vedantham et al. 2020; Kavanagh et al. 2021,
2022).

Finally, AD Leo may not be an isolated case. To verify this, it
is essential to explore the possibility for such cycles over a wider
area of the stellar parameter space, namely mass and rotation
period.

***Acknowledgements.*** We acknowledge funding from the French National Research Agency (ANR) under contract number ANR-18-CE31-0019 (SPiLaSH).

SB acknowledges funding from the European Space Agency (ESA), under the visiting researcher programme. LTL acknowledges funding from the European Research Council under the H2020 research & innovation programme (grant #740651 NewWorlds). XD and AC acknowledge for funding in the framework of the Investissements d'Avenir programme (ANR-15-IDEX-02), through the funding of the 'Origin of Life' project of the Univ. Grenoble-Alpes. OK acknowledges support by the Swedish Research Council (grant agreement no. 2019-03548), the Swedish National Space Agency, and the Royal Swedish Academy of Sciences. 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. We gratefully acknowledge the CFHT QSO observers who made this project possible. This work has made use of the VALD database, operated at Uppsala University, the Institute of Astronomy RAS in Moscow, and the University of Vienna; Astropy, 12 a community-developed core Python package for Astronomy (Astropy Collaboration et al. 2013, 2018); NumPy (van der Walt et al. 2011); Matplotlib: Visualization with Python (Hunter 2007); SciPy (Virtanen et al. 2020).

References

Artigau, É., Astudillo-Defru, N., Delfosse, X., et al. 2014, in Society of Photo-
Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9149, Ob-
servatory Operations: Strategies, Processes, and Systems V, ed. A. B. Peck,
C. R. Benn, & R. L. Seaman, 914905

Astropy Collaboration, Price-Whelan, A. M., Sipőcz, B. M., et al. 2018, AJ, 156, 123

Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33

Babcock, H. W. 1961, ApJ, 133, 572

Bagnulo, S., Landolfi, M., Landstreet, J. D., et al. 2009, PASP, 121, 993

Baliunas, S. L., Donahue, R. A., Soon, W. H., et al. 1995, ApJ, 438, 269

Barnes, J. R., James, D. J., & Collier Cameron, A. 2004, MNRAS, 352, 589

Baroch, D., Morales, J. C., Ribas, I., et al. 2020, A&A, 641, A69

Bellotti, S., Petit, P., Morin, J., et al. 2022, A&A, 657, A107

Bertaux, J. L., Lallement, R., Ferron, S., Boonne, C., & Bodichon, R. 2014, A&A, 564, A46

Blackwell-Whitehead, R. J., Lundberg, H., Nave, G., et al. 2006, MNRAS, 373, 1603

Boro Saikia, S., Jeffers, S. V., Morin, J., et al. 2016, A&A, 594, A29

Boro Saikia, S., Lüftinger, T., Folsom, C. P., et al. 2022, A&A, 658, A16

Boro Saikia, S., Marvin, C. J., Jeffers, S. V., et al. 2018, A&A, 616, A108

Buccino, A. P., Petrucci, R., Jofré, E., & Mauas, P. J. D. 2014, ApJ, 781, L9

Callingham, J. R., Vedantham, H. K., Shimwell, T. W., et al. 2021, Nature As-
tronomy, 5, 1233

Cang, T. Q., Petit, P., Donati, J. F., et al. 2020, A&A, 643, A39

Carleo, I., Malavolta, L., Lanza, A. F., et al. 2020, A&A, 638, A5

Carmona, A., Delfosse, X., Bellotti, S., et al. 2023, A&A, 674, A110

Chabrier, G. & Baraffe, I. 1997, A&A, 327, 1039

Charbonneau, P. 2010, Living Reviews in Solar Physics, 7, 3

Claret, A. & Bloemen, S. 2011, A&A, 529, A75

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

Costes, J. C., Watson, C. A., de Mooij, E., et al. 2021, MNRAS, 505, 830

Cristofari, P. I., Donati, J. F., Folsom, C. P., et al. 2023, MNRAS, 522, 1342

Demidov, M. L. & Balthasar, H. 2012, Sol. Phys., 276, 43

Demidov, M. L., Golubeva, E. M., Balthasar, H., Staude, J., & Grigoryev, V. M. 2008, Sol. Phys., 250, 279

Donati, J. F. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 307, Solar Polarization, ed. J. Trujillo-Bueno & J. Sanchez Almeida, 41

Donati, J. F. & Brown, S. F. 1997, A&A, 326, 1135

Donati, J. F., Cristofari, P. I., Finociety, B., et al. 2023, MNRAS[arXiv:2304.09642]

Donati, J. F., Howarth, I. D., Jardine, M. M., et al. 2006, MNRAS, 370, 629

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

Donati, J. F. & Landstreet, J. D. 2009, ARA&A, 47, 333

Donati, J. F., Morin, J., Petit, P., et al. 2008a, MNRAS, 390, 545

Donati, J. F., Moutou, C., Farès, R., et al. 2008b, MNRAS, 385, 1179

Donati, J. F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, MNRAS, 291, 658

Fares, R., Donati, J. F., Moutou, C., et al. 2009, MNRAS, 398, 1383

Feiden, G. A., Skidmore, K., & Jao, W.-C. 2021, ApJ, 907, 53

Finociety, B., Donati, J. F., Klein, B., et al. 2021, MNRAS, 508, 3427

Folsom, C. P., Bouvier, J., Petit, P., et al. 2018, MNRAS, 474, 4956

Folsom, C. P., Petit, P., Bouvier, J., et al. 2016, MNRAS, 457, 580

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306
Fouqué, P., Martioli, E., Donati, J. F., et al. 2023, A&A, 672, A52

Fuhrmeister, B., Czesla, S., Perdelwitz, V., et al. 2023, A&A, 670, A71

Gaia Collaboration, Smart, R. L., Sarro, L. M., et al. 2021, A&A, 649, A6

Gastine, T., Morin, J., Duarte, L., et al. 2013, A&A, 549, L5

Gomes da Silva, J., Santos, N. C., Bonfils, X., et al. 2012, A&A, 541, A9

Granzer, T., Schüssler, M., Caligari, P., & Strassmeier, K. G. 2000, A&A, 355, 1087

Gregory, S. G., Donati, J. F., Morin, J., et al. 2012, ApJ, 755, 97

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

Hale, G. E., Ellerman, F., Nicholson, S. B., & Joy, A. H. 1919, ApJ, 49, 153

Hathaway, D. H. 2010, Living Reviews in Solar Physics, 7, 1

Haywood, R. D., Milbourne, T. W., Saar, S. H., et al. 2022, ApJ, 935, 6

Hazra, G., Vidotto, A. A., & D'Angelo, C. V. 2020, MNRAS, 496, 4017

Hébrard, É. M., Donati, J. F., Delfosse, X., et al. 2016, MNRAS, 461, 1465

Hill, C. A., Folsom, C. P., Donati, J. F., et al. 2019, MNRAS, 484, 5810

Hunter, J. D. 2007, Computing in Science and Engineering, 9, 90

Jeffers, S. V., Boro Saikia, S., Barnes, J. R., et al. 2017, MNRAS, 471, L96

Jeffers, S. V., Cameron, R. H., Marsden, S. C., et al. 2022, A&A, 661, A152

Jeffers, S. V., Mengel, M., Moutou, C., et al. 2018, MNRAS, 479, 5266

Johns-Krull, C. M., Valenti, J. A., & Koresko, C. 1999, ApJ, 516, 900

Kausch, W., Noll, S., Smette, A., et al. 2015, A&A, 576, A78

Kavanagh, R. D., Vidotto, A. A., Klein, B., et al. 2021, MNRAS, 504, 1511

Kavanagh, R. D., Vidotto, A. A., Vedantham, H. K., et al. 2022, MNRAS, 514, 675

Kitchatinov, L. L., Moss, D., & Sokoloff, D. 2014, MNRAS, 442, L1

Klein, B., Donati, J.-F., Moutou, C., et al. 2021, MNRAS, 502, 188

Kochukhov, O. 2021, A&A Rev., 29, 1

Kochukhov, O., Heiter, U., Piskunov, N., et al. 2009, in American Institute of Physics Conference Series, Vol. 1094, 15th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, ed. E. Stempels, 124–129

Kochukhov, O. & Lavail, A. 2017, ApJ, 835, L4

Kochukhov, O., Makaganiuk, V., & Piskunov, N. 2010, A&A, 524, A5

Konings, T., Baeyens, R., & Decin, L. 2022, A&A, 667, A15

Landi Degl'Innocenti, E. & Landolfi, M. 2004, Polarization in Spectral Lines, Vol. 307

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

Lanza, A. F. 2013, A&A, 557, A31

Lavail, A., Kochukhov, O., & Wade, G. A. 2018, MNRAS, 479, 4836

Lawler, J. E., Guzman, A., Wood, M. P., Sneden, C., & Cowan, J. J. 2013, ApJS, 205, 11

Lehmann, L. T. & Donati, J. F. 2022, MNRAS, 514, 2333

Lehmann, L. T., Hussain, G. A. J., Vidotto, A. A., Jardine, M. M., & Mackey, D. H. 2021, MNRAS, 500, 1243

Leighton, R. B. 1969, ApJ, 156, 1

Lopez-Santiago, J., Martino, L., Míguez, J., & Vázquez, M. A. 2020, AJ, 160, 273

Louca, A. J., Miguel, Y., Tsai, S.-M., et al. 2023, MNRAS, 521, 3333

Lovis, C., Dumusque, X., Santos, N. C., et al. 2011, arXiv e-prints, arXiv:1107.5325

Lund, K., Jardine, M., Lehmann, L. T., et al. 2020, MNRAS, 493, 1003

Lund, K., Jardine, M., Russell, A. J. B., et al. 2021, MNRAS, 502, 4903

Maeder, A. & Meynet, G. 2000, ARA&A, 38, 143

Mann, A. W., Feiden, G. A., Gaidos, E., Boyajian, T., & von Braun, K. 2015, ApJ, 804, 64

Maunder, E. W. 1904, MNRAS, 64, 747

McCann, J., Murray-Clay, R. A., Kratter, K., & Krumholz, M. R. 2019, ApJ, 873, 89

Meunier, N., Desort, M., & Lagrange, A. M. 2010, A&A, 512, A39

Meunier, N. & Lagrange, A. M. 2019, A&A, 628, A125

Mignon, L., Meunier, N., Delfosse, X., et al. 2023, arXiv e-prints, arXiv:2303.03998

Mittag, M., Robrade, J., Schmitt, J. H. M. M., et al. 2017, A&A, 600, A119

Morgenthaler, A., Petit, P., Morin, J., et al. 2011, Astronomische Nachrichten, 332, 866

Morin, J. 2012, in EAS Publications Series, Vol. 57, EAS Publications Series, ed. C. Reylé, C. Charbonnel, & M. Schultheis, 165–191

Morin, J., Delfosse, X., Donati, J. F., et al. 2011, in SF2A-2011: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, ed. G. Alecian, K. Belkacem, R. Samadi, & D. Valls-Gabaud, 503–508

Morin, J., Donati, J. F., Forveille, T., et al. 2008a, MNRAS, 384, 77

Morin, J., Donati, J. F., Petit, P., et al. 2008b, MNRAS, 390, 567

Morin, J., Donati, J. F., Petit, P., et al. 2010, MNRAS, 407, 2269

Morin, J., Hill, C. A., & Watson, C. A. 2016, in Astrophysics and Space Science Library, Vol. 439, Astronomy at High Angular Resolution, ed. H. M. J. Boffin, G. Hussain, J.-P. Berger, & L. Schmidtobreick, 223

Muheki, P., Guenther, E. W., Mutabazi, T., & Jurua, E. 2020, A&A, 637, A13

Mullan, D. J. & MacDonald, J. 2001, ApJ, 559, 353

Namekata, K., Maehara, H., Sasaki, R., et al. 2020, PASJ, 72, 68

O'Brian, T. R., Wickliffe, M. E., Lawler, J. E., Whaling, W., & Brault, J. W. 1991, Journal of the Optical Society of America B Optical Physics, 8, 1185

Parker, E. N. 1955, ApJ, 122, 293

Petit, P., Dintrans, B., Morgenthaler, A., et al. 2009, A&A, 508, L9

Petit, P., Dintrans, B., Solanki, S. K., et al. 2008, MNRAS, 388, 80

Petit, P., Donati, J. F., & Collier Cameron, A. 2002, MNRAS, 334, 374

Petit, P., Folsom, C. P., Donati, J. F., et al. 2021, A&A, 648, A55

Petit, P., Louge, T., Théado, S., et al. 2014, PASP, 126, 469

Petrovay, K. 2020, Living Reviews in Solar Physics, 17, 2

Pipin, V. V., Pevtsov, A. A., Liu, Y., & Kosovichev, A. G. 2019, ApJ, 877, L36

Preston, G. W. 1967, ApJ, 150, 547

Rachkovsky, D. N. 1967, Izvestiya Ordena Trudovogo Krasnogo Znameni Krymskoj Astrofizicheskoj Observatorii, 37, 56

Reiners, A. 2012, Living Reviews in Solar Physics, 9, 1

Reiners, A. & Basri, G. 2009, A&A, 496, 787

Reiners, A., Bean, J. L., Huber, K. F., et al. 2010, ApJ, 710, 432

Reiners, A., Shulyak, D., Anglada-Escudé, G., et al. 2013, A&A, 552, A103

Reiners, A., Shulyak, D., Käpylä, P. J., et al. 2022, A&A, 662, A41

Reinhold, T., Bell, K. J., Kuszlewicz, J., Hekker, S., & Shapiro, A. I. 2019, A&A, 621, A21

Robertson, P., Endl, M., Cochran, W. D., & Dodson-Robinson, S. E. 2013, in American Astronomical Society Meeting Abstracts, Vol. 221, American Astronomical Society Meeting Abstracts #221, 423.05

Rosén, L., Kochukhov, O., Hackman, T., & Lehtinen, J. 2016, A&A, 593, A35

Route, M. 2016, ApJ, 830, L27

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

Saar, S. H. 1994, in Infrared Solar Physics, ed. D. M. Rabin, J. T. Jefferies, & C. Lindsey, Vol. 154, 493

Sairam, L. & Triaud, A. H. M. J. 2022, MNRAS, 514, 2259

Sanderson, T. R., Appourchaux, T., Hoeksema, J. T., & Harvey, K. L. 2003, Journal of Geophysical Research (Space Physics), 108, 1035

Schuessler, M. & Solanki, S. K. 1992, A&A, 264, L13

Schwabe, H. 1844, Astronomische Nachrichten, 21, 233

See, V., Jardine, M., Vidotto, A. A., et al. 2016, MNRAS, 462, 4442

See, V., Jardine, M., Vidotto, A. A., et al. 2015, MNRAS, 453, 4301

Semel, M. 1989, A&A, 225, 456

Shkolnik, E., Liu, M. C., & Reid, I. N. 2009, ApJ, 699, 649

Shulyak, D., Reiners, A., Engeln, A., et al. 2017, Nature Astronomy, 1, 0184

Shulyak, D., Reiners, A., Nagel, E., et al. 2019, A&A, 626, A86

Shulyak, D., Reiners, A., Seemann, U., Kochukhov, O., & Piskunov, N. 2014, A&A, 563, A35

Skilling, J. & Bryan, R. K. 1984, MNRAS, 211, 111

Smette, A., Sana, H., Noll, S., et al. 2015, A&A, 576, A77

Solanki, S. K. 1993, Space Sci. Rev., 63, 1

Stelzer, B., Caramazza, M., Raetz, S., Argiroffi, C., & Coffaro, M. 2022, A&A, 667, L9

Stibbs, D. W. N. 1950, MNRAS, 110, 395

Tanner, J. D., Basu, S., & Demarque, P. 2013, ApJ, 767, 78

Tinetti, G., Eccleston, P., Haswell, C., et al. 2021, arXiv e-prints, arXiv:2104.04824

Tuomi, M., Jones, H. R. A., Barnes, J. R., et al. 2018, AJ, 155, 192

Unno, W. 1956, PASJ, 8, 108

Valenti, J. A., Marcy, G. W., & Basri, G. 1995, ApJ, 439, 939

van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Computing in Science and Engineering, 13, 22

van Saders, J. L. & Pinsonneault, M. H. 2012, ApJ, 746, 16

Vedantham, H. K., Callingham, J. R., Shimwell, T. W., et al. 2020, Nature Astronomy, 4, 577

Vidotto, A. A. & Cleary, A. 2020, MNRAS, 494, 2417

Virtanen, P., Gommers, R., Burovski, E., et al. 2020, scipy/scipy: SciPy 1.5.3

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

Wilson, O. C. 1968, ApJ, 153, 221

Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48

Yeo, K. L., Krivova, N. A., & Solanki, S. K. 2014, Space Sci. Rev., 186, 137

Zacharias, N., Finch, C. T., Girard, T. M., et al. 2013, AJ, 145, 44

Zaire, B., Jouve, L., Gastine, T., et al. 2022, MNRAS, 517, 3392

Zayer, I., Solanki, S. K., & Stenflo, J. O. 1989, A&A, 211, 463

Zeeman, P. 1897, Nature, 55, 347
## Appendix A: FWHM of Stokes I

In this appendix, the Stokes *I* profiles computed using different line selections are shown. We compare the width of the profiles obtained with the full atomic mask, high-Landé factor ($g_{eff} > 1.2$) lines and low-Landé factor ($g_{eff} < 1.2$) lines.

## Appendix B: Chromatic Stokes profiles

The various Stokes profiles computed with different wavelength-based line lists for LSD are reported. The wavelength intervals of the line subsets are: [350,390], [390,430], [430,480], [480,550], [550,650], [650,1100], [950,1100], [1100,1400], [1400,1600], [1600,1800], [1800,2500] nm.

## Appendix C: Zeeman broadening examples

Example plots of the posterior distributions from the Zeeman broadening MCMC analysis are shown in Fig. C.1 for a ES-PaDOnS observation and in Fig. C.2 for a SPIRou observation. Summaries of the results of the MCMC analysis for each epoch are provided in Tables C.1 and C.2.

## Appendix D: Observing log

This appendix contains the journal of observations of AD Leo, for both optical and near-infrared observations. It also includes all measurements of longitudinal magnetic field and $\Sigma Bf$.

<figure><img src="image_16.png" /><figcaption>Fig. A.1. Median Stokes *I* profiles computed for the four near-infrared epochs: 2019a, 2019b, 2020a, and 2020b. The profiles are obtained using the full mask (top) magnetically sensitive (middle) and insensitive (bottom) line lists in LSD. Following the linear dependence of the Zeeman effect on the Landé factor, high-*g*<sub>eff</sub> lines are broader than low-*g*<sub>eff</sub> lines on average. The width of the profiles changes over time in the high-*g*<sub>eff</sub> case, while it remains reasonably constant in the low-*g*<sub>eff</sub> case.</figcaption></figure>
Optical

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

Near-Infrared

<figure><img src="image_18.png" /><figcaption>**Fig. B.1.** Series of Stokes profiles computed with subsets of the line list used for LSD, as described in Sec.4. Every panel displays Stokes *I* (top), *V* (middle) and *N* (bottom) for each line subset and vertical dashed lines indicate the adjusted velocity range over which *B*<sub>I</sub> was estimated. The optical profiles are obtained by stacking the ESPaDOnS observations from 2019, whereas the near-infrared ones from SPIRou 2019b time series, for each panel. The number of lines used in the LSD computation as well as the mean wavelength are displayed.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Fig. C.1.** Posterior distribution for parameters from the Zeeman broadening analysis of the ESPaDOoS observation on 24 February 2016. The filling factor for a given magnetic field strength is *f*, and abundances are in log *n*<sub>X</sub>/n<sub>H</sub> units.</figcaption></figure>

<table><caption>**Table C.1.** Parameters from the Zeeman broadening analysis, averaged for each epoch, for ESPaDOoS and Narval optical observations. Uncertainties are the standard deviations for each epoch.</caption><thead><tr><th>epoch</th><th>2006</th><th>2007</th><th>2008</th><th>2011</th><th>2012</th><th>2016</th><th>2019</th></tr></thead><tbody><tr><td>ν<sub>e</sub> sin(i) [km s<sup>-1</sup>]</td><td>4.27±0.38</td><td>4.25±0.17</td><td>4.16±0.29</td><td>3.46±0.50</td><td>2.65±0.48</td><td>3.57±0.21</td><td>3.39±0.45</td></tr><tr><td>[Ti/H]</td><td>-6.97±0.01</td><td>-6.90±0.01</td><td>-6.93±0.02</td><td>-6.99±0.01</td><td>-7.00±0.01</td><td>-6.96±0.01</td><td>-6.95±0.01</td></tr><tr><td>f<sub>0kG</sub></td><td>0.06±0.02</td><td>0.24±0.06</td><td>0.18±0.06</td><td>0.05±0.06</td><td>0.02±0.01</td><td>0.05±0.02</td><td>0.06±0.01</td></tr><tr><td>f<sub>2kG</sub></td><td>0.46±0.01</td><td>0.32±0.07</td><td>0.39±0.07</td><td>0.60±0.06</td><td>0.52±0.02</td><td>0.44±0.01</td><td>0.41±0.01</td></tr><tr><td>f<sub>4kG</sub></td><td>0.33±0.02</td><td>0.31±0.02</td><td>0.28±0.04</td><td>0.19±0.07</td><td>0.29±0.01</td><td>0.33±0.01</td><td>0.36±0.01</td></tr><tr><td>f<sub>6kG</sub></td><td>0.09±0.01</td><td>0.06±0.01</td><td>0.08±0.01</td><td>0.11±0.05</td><td>0.09±0.01</td><td>0.11±0.01</td><td>0.09±0.01</td></tr><tr><td>f<sub>8kG</sub></td><td>0.05±0.01</td><td>0.03±0.01</td><td>0.04±0.01</td><td>0.04±0.03</td><td>0.05±0.01</td><td>0.05±0.01</td><td>0.05±0.01</td></tr><tr><td>f<sub>10kG</sub></td><td>0.04±0.01</td><td>0.02±0.01</td><td>0.02±0.01</td><td>0.03±0.03</td><td>0.02±0.01</td><td>0.04±0.01</td><td>0.04±0.01</td></tr><tr><td>Σ Bf [kG]</td><td>3.59±0.04</td><td>2.82±0.13</td><td>2.97±0.13</td><td>3.22±0.15</td><td>3.37±0.06</td><td>3.65±0.04</td><td>3.62±0.03</td></tr></tbody></table>
<figure><img src="image_20.png" /><figcaption>Fig. C.2. Posterior distribution for parameters from the Zeeman broadening analysis of the SPIRou observation on 3 November 2020.</figcaption></figure>

<table><caption>**Table C.2.** Parameters from the Zeeman broadening analysis, averaged for each epoch, for SPIRou observations. Uncertainties are the standard deviations for each epoch.</caption><thead><tr><th>epoch</th><th>2019a</th><th>2019b</th><th>2020a</th><th>2020b</th></tr></thead><tbody><tr><td>v<sub>eq</sub> sin(i) [km s<sup>-1</sup>]</td><td>3.42±0.23</td><td>3.34±0.18</td><td>3.48±0.17</td><td>3.48±0.21</td></tr><tr><td>[Ti/H]</td><td>-6.78±0.02</td><td>-6.78±0.02</td><td>-6.79±0.02</td><td>-6.78±0.01</td></tr><tr><td>[Fe/H]</td><td>-4.20±0.01</td><td>-4.22±0.01</td><td>-4.20±0.03</td><td>-4.21±0.02</td></tr><tr><td>f<sub>0kG</sub></td><td>0.06±0.02</td><td>0.05±0.02</td><td>0.05±0.02</td><td>0.06±0.01</td></tr><tr><td>f<sub>1kG</sub></td><td>0.15±0.02</td><td>0.12±0.03</td><td>0.15±0.02</td><td>0.17±0.02</td></tr><tr><td>f<sub>2kG</sub></td><td>0.20±0.04</td><td>0.23±0.05</td><td>0.19±0.04</td><td>0.19±0.04</td></tr><tr><td>f<sub>3kG</sub></td><td>0.24±0.04</td><td>0.19±0.05</td><td>0.23±0.05</td><td>0.24±0.04</td></tr><tr><td>f<sub>4kG</sub></td><td>0.12±0.04</td><td>0.16±0.04</td><td>0.14±0.03</td><td>0.13±0.03</td></tr><tr><td>f<sub>5kG</sub></td><td>0.06±0.02</td><td>0.08±0.02</td><td>0.06±0.02</td><td>0.06±0.01</td></tr><tr><td>f<sub>6kG</sub></td><td>0.05±0.01</td><td>0.05±0.02</td><td>0.05±0.01</td><td>0.04±0.01</td></tr><tr><td>f<sub>7kG</sub></td><td>0.03±0.01</td><td>0.02±0.01</td><td>0.03±0.01</td><td>0.03±0.01</td></tr><tr><td>f<sub>8kG</sub></td><td>0.02±0.01</td><td>0.02±0.01</td><td>0.02±0.01</td><td>0.02±0.01</td></tr><tr><td>f<sub>9kG</sub></td><td>0.03±0.01</td><td>0.03±0.01</td><td>0.03±0.01</td><td>0.03±0.01</td></tr><tr><td>f<sub>10kG</sub></td><td>0.03±0.01</td><td>0.03±0.01</td><td>0.04±0.01</td><td>0.03±0.01</td></tr><tr><td>Σ Bf [kG]</td><td>3.40±0.09</td><td>3.49±0.09</td><td>3.44±0.04</td><td>3.29±0.07</td></tr></tbody></table>
<table><caption>**Table D.1.** List of AD Leo observations collected with SPIRou. The columns are: (1 and 2) date and universal time of the observations, (3) rotational cycle of the observations found using Eq. 1, (4) exposure time of a polarimetric sequence, (5) signal-to-noise ratio at 1650 nm per spectral element, (6) RMS noise level of Stokes V signal in units of unpolarised continuum.</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-4</sup>I<sub>c</sub>]</th></tr></thead><tbody><tr><td colspan="6" style="text-align: center; font-weight: bold;">2019</td></tr><tr><td>April 15</td><td>06:11:02.35</td><td>0.00</td><td>4x61</td><td>151</td><td>1.9</td></tr><tr><td>April 16</td><td>12:03:07.61</td><td>0.56</td><td>4x61</td><td>130</td><td>1.8</td></tr><tr><td>April 18</td><td>09:00:59.77</td><td>1.40</td><td>4x61</td><td>138</td><td>1.9</td></tr><tr><td>April 19</td><td>10:36:14.55</td><td>1.88</td><td>4x61</td><td>143</td><td>1.8</td></tr><tr><td>April 20</td><td>08:56:19.44</td><td>2.29</td><td>4x61</td><td>143</td><td>1.8</td></tr><tr><td>April 21</td><td>05:42:43.43</td><td>2.68</td><td>4x61</td><td>165</td><td>1.5</td></tr><tr><td>April 22</td><td>08:58:10.67</td><td>3.19</td><td>4x61</td><td>147</td><td>1.9</td></tr><tr><td>April 23</td><td>07:01:28.73</td><td>3.60</td><td>4x61</td><td>154</td><td>1.9</td></tr><tr><td>April 24</td><td>06:04:06.74</td><td>4.03</td><td>4x61</td><td>147</td><td>1.9</td></tr><tr><td>April 25</td><td>09:41:45.90</td><td>4.55</td><td>4x61</td><td>152</td><td>1.9</td></tr><tr><td>April 26</td><td>08:13:21.75</td><td>4.97</td><td>4x61</td><td>162</td><td>1.9</td></tr><tr><td>April 27</td><td>08:17:48.22</td><td>5.42</td><td>4x61</td><td>139</td><td>2.1</td></tr><tr><td>May 01</td><td>08:59:11.25</td><td>7.23</td><td>4x61</td><td>153</td><td>1.8</td></tr><tr><td>May 15</td><td>06:11:38.01</td><td>13.45</td><td>4x61</td><td>165</td><td>1.6</td></tr><tr><td>June 13</td><td>06:30:59.49</td><td>26.46</td><td>4x61</td><td>186</td><td>2.1</td></tr><tr><td>June 14</td><td>05:44:29.79</td><td>26.89</td><td>4x61</td><td>193</td><td>2.0</td></tr><tr><td>June 15</td><td>06:04:31.16</td><td>27.35</td><td>4x61</td><td>192</td><td>2.1</td></tr><tr><td>June 16</td><td>05:44:52.21</td><td>27.79</td><td>4x61</td><td>173</td><td>1.8</td></tr><tr><td>June 17</td><td>06:08:46.48</td><td>28.25</td><td>4x61</td><td>150</td><td>1.8</td></tr><tr><td>June 19</td><td>05:47:00.42</td><td>29.14</td><td>4x61</td><td>175</td><td>1.5</td></tr><tr><td>June 21</td><td>06:18:29.14</td><td>30.05</td><td>4x61</td><td>169</td><td>1.7</td></tr><tr><td>October 16</td><td>15:31:29.14</td><td>82.69</td><td>4x61</td><td>180</td><td>1.4</td></tr><tr><td>October 31</td><td>15:32:56.51</td><td>89.41</td><td>4x61</td><td>169</td><td>1.6</td></tr><tr><td>November 01</td><td>15:22:10.08</td><td>89.86</td><td>4x61</td><td>159</td><td>1.5</td></tr><tr><td>November 02</td><td>15:36:52.90</td><td>90.31</td><td>4x61</td><td>170</td><td>1.4</td></tr><tr><td>November 03</td><td>14:53:22.07</td><td>90.75</td><td>4x61</td><td>151</td><td>1.5</td></tr><tr><td>November 04</td><td>15:43:18.06</td><td>91.21</td><td>4x61</td><td>137</td><td>1.9</td></tr><tr><td>November 05</td><td>15:33:34.93</td><td>91.66</td><td>4x61</td><td>164</td><td>1.3</td></tr><tr><td>November 06</td><td>15:37:27.47</td><td>92.10</td><td>4x61</td><td>190</td><td>1.5</td></tr><tr><td>November 07</td><td>15:01:11.73</td><td>92.54</td><td>4x61</td><td>165</td><td>1.5</td></tr><tr><td>November 09</td><td>14:03:02.38</td><td>93.42</td><td>4x61</td><td>116</td><td>1.6</td></tr><tr><td>November 10</td><td>15:51:32.68</td><td>93.90</td><td>4x61</td><td>151</td><td>1.5</td></tr><tr><td>November 13</td><td>14:47:09.89</td><td>95.22</td><td>4x61</td><td>201</td><td>1.8</td></tr><tr><td>November 14</td><td>14:13:58.20</td><td>95.66</td><td>4x61</td><td>197</td><td>1.4</td></tr><tr><td>December 05</td><td>15:12:55.68</td><td>105.10</td><td>4x61</td><td>116</td><td>1.4</td></tr><tr><td>December 05</td><td>15:35:03.10</td><td>105.10</td><td>4x61</td><td>68</td><td>1.5</td></tr><tr><td>December 07</td><td>15:08:39.83</td><td>105.99</td><td>4x61</td><td>133</td><td>2.2</td></tr><tr><td>December 08</td><td>14:29:19.93</td><td>106.43</td><td>4x61</td><td>181</td><td>1.8</td></tr><tr><td>December 09</td><td>14:21:48.74</td><td>106.88</td><td>4x61</td><td>193</td><td>1.5</td></tr><tr><td>December 10</td><td>13:15:39.48</td><td>107.31</td><td>4x61</td><td>194</td><td>1.5</td></tr><tr><td>December 11</td><td>14:58:07.44</td><td>107.79</td><td>4x61</td><td>189</td><td>3.0</td></tr><tr><td>December 12</td><td>14:33:27.45</td><td>108.23</td><td>4x61</td><td>190</td><td>3.6</td></tr><tr><td colspan="6" style="text-align: center; font-weight: bold;">2020</td></tr><tr><td>January 26</td><td>12:11:48.90</td><td>128.36</td><td>4x61</td><td>218</td><td>2.0</td></tr><tr><td>January 27</td><td>12:07:16.45</td><td>128.81</td><td>4x61</td><td>166</td><td>1.5</td></tr><tr><td>January 28</td><td>12:06:25.46</td><td>129.26</td><td>4x61</td><td>193</td><td>1.5</td></tr><tr><td>February 05</td><td>08:18:51.61</td><td>132.78</td><td>4x61</td><td>174</td><td>1.4</td></tr><tr><td>February 16</td><td>07:39:02.87</td><td>137.70</td><td>4x61</td><td>193</td><td>1.8</td></tr><tr><td>February 17</td><td>09:20:45.46</td><td>138.18</td><td>4x61</td><td>171</td><td>1.9</td></tr><tr><td>February 18</td><td>07:39:13.56</td><td>138.59</td><td>4x61</td><td>191</td><td>1.2</td></tr><tr><td>February 19</td><td>08:37:09.84</td><td>139.06</td><td>4x61</td><td>210</td><td>1.5</td></tr><tr><td>March 12</td><td>09:16:56.32</td><td>148.94</td><td>4x61</td><td>217</td><td>2.0</td></tr><tr><td>May 08</td><td>05:59:01.58</td><td>174.43</td><td>4x61</td><td>194</td><td>1.4</td></tr><tr><td>May 09</td><td>09:42:30.22</td><td>174.95</td><td>4x61</td><td>206</td><td>1.5</td></tr><tr><td>May 12</td><td>09:39:45.19</td><td>176.30</td><td>4x61</td><td>177</td><td>1.4</td></tr></tbody></table>
<table><caption>**Table D.1.** continued.</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-4</sup>I<sub>c</sub>]</th></tr></thead><tbody><tr><td>May 13</td><td>09:43:17.10</td><td>176.75</td><td>4x61</td><td>204</td><td>1.4</td></tr><tr><td>May 14</td><td>07:46:35.66</td><td>177.16</td><td>4x61</td><td>208</td><td>1.5</td></tr><tr><td>May 15</td><td>09:50:31.27</td><td>177.65</td><td>4x61</td><td>164</td><td>1.6</td></tr><tr><td>May 31</td><td>06:22:44.29</td><td>184.76</td><td>4x61</td><td>215</td><td>1.2</td></tr><tr><td>June 01</td><td>07:29:57.32</td><td>185.23</td><td>4x61</td><td>189</td><td>1.3</td></tr><tr><td>June 02</td><td>06:23:31.70</td><td>185.65</td><td>4x61</td><td>202</td><td>1.4</td></tr><tr><td>June 03</td><td>07:38:43.10</td><td>186.13</td><td>4x61</td><td>196</td><td>1.3</td></tr><tr><td>June 04</td><td>06:32:43.13</td><td>186.55</td><td>4x61</td><td>197</td><td>1.3</td></tr><tr><td>June 05</td><td>06:42:42.88</td><td>187.01</td><td>4x61</td><td>139</td><td>1.8</td></tr><tr><td>June 06</td><td>07:53:54.53</td><td>187.48</td><td>4x61</td><td>167</td><td>1.2</td></tr><tr><td>June 07</td><td>06:58:07.88</td><td>187.91</td><td>4x61</td><td>154</td><td>1.6</td></tr><tr><td>June 08</td><td>06:03:47.85</td><td>188.34</td><td>4x61</td><td>135</td><td>1.3</td></tr><tr><td>June 08</td><td>06:10:09.48</td><td>188.34</td><td>4x61</td><td>143</td><td>1.4</td></tr><tr><td>June 09</td><td>06:31:50.17</td><td>188.79</td><td>4x61</td><td>120</td><td>1.3</td></tr><tr><td>June 10</td><td>06:57:13.44</td><td>189.25</td><td>4x61</td><td>200</td><td>1.8</td></tr><tr><td>October 31</td><td>15:06:44.49</td><td>253.75</td><td>4x61</td><td>205</td><td>1.5</td></tr><tr><td>November 03</td><td>15:29:52.78</td><td>255.11</td><td>4x61</td><td>220</td><td>1.5</td></tr></tbody></table>

<table><caption>**Table D.2.** List of AD Leo observations collected with ESPaDOnS in 2019. The columns are: (1 and 2) date and universal time of the observations, (3) rotational cycle of the observations found using Eq. 1, (4) exposure time of a polarimetric sequence, (5) signal-to-noise ratio at 650 nm per spectral element, (6) RMS noise level of Stokes V signal in units of unpolarised continuum.</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-4</sup>I<sub>c</sub>]</th></tr></thead><tbody><tr><td>November 15</td><td>13:24:01.40</td><td>96.32</td><td>4x300</td><td>234</td><td>1.7</td></tr><tr><td>November 16</td><td>14:20:01.30</td><td>96.79</td><td>4x300</td><td>230</td><td>1.9</td></tr><tr><td>November 19</td><td>14:08:04.00</td><td>98.13</td><td>4x300</td><td>151</td><td>3.6</td></tr><tr><td>November 19</td><td>14:33:03.00</td><td>98.14</td><td>4x300</td><td>159</td><td>3.2</td></tr><tr><td>November 19</td><td>14:56:04.00</td><td>98.15</td><td>4x300</td><td>180</td><td>2.8</td></tr><tr><td>November 21</td><td>15:43:01.00</td><td>99.06</td><td>4x300</td><td>265</td><td>1.8</td></tr></tbody></table>

<table><caption>**Table D.3.** List of optical and near-infrared measurements of longitudinal magnetic field and magnetic flux. The columns are: (1) Heliocentric Julian date of the observation, (2) $B_l$ with formal error bar (see Eq. 2), and (3) magnetic flux from Zeeman broadening modelling, when a reliable measurement was possible, and (4) the instrument used.</caption><thead><tr><th>HJD<br>[-2450000]</th><th>B<sub>l</sub><br>[G]</th><th>Bf<br>[kG]</th><th>Instrument</th></tr></thead><tbody><tr><td>3747.0876</td><td>-269.4 ± 26.4</td><td>3.65<sup>+0.07</sup><sub>-0.08</sub></td><td>ESPaDOnS</td></tr><tr><td>3748.8868</td><td>-272.4 ± 15.1</td><td>3.52<sup>+0.09</sup><sub>-0.10</sub></td><td>ESPaDOnS</td></tr><tr><td>3780.0705</td><td>-280.7 ± 13.1</td><td>3.55<sup>+0.10</sup><sub>-0.10</sub></td><td>ESPaDOnS</td></tr><tr><td>3895.8047</td><td>-291.2 ± 10.7</td><td>3.62<sup>+0.08</sup><sub>-0.09</sub></td><td>ESPaDOnS</td></tr><tr><td>3896.8124</td><td>-260.1 ± 7.6</td><td>3.61<sup>+0.08</sup><sub>-0.09</sub></td><td>ESPaDOnS</td></tr><tr><td>3897.8005</td><td>-266.7 ± 7.4</td><td>3.63<sup>+0.08</sup><sub>-0.09</sub></td><td>ESPaDOnS</td></tr><tr><td>3898.7785</td><td>-294.2 ± 8.1</td><td>3.58<sup>+0.08</sup><sub>-0.10</sub></td><td>ESPaDOnS</td></tr><tr><td>4127.5975</td><td>-294.6 ± 16.5</td><td>...</td><td>Narval</td></tr><tr><td>4128.6088</td><td>-248.6 ± 10.8</td><td>2.95<sup>+0.16</sup><sub>-0.16</sub></td><td>Narval</td></tr><tr><td>4129.5717</td><td>-296.2 ± 11.8</td><td>2.85<sup>+0.21</sup><sub>-0.22</sub></td><td>Narval</td></tr><tr><td>4130.6084</td><td>-253.8 ± 9.4</td><td>2.65<sup>+0.21</sup><sub>-0.20</sub></td><td>Narval</td></tr><tr><td>4133.6312</td><td>-274.8 ± 10.4</td><td>...</td><td>Narval</td></tr></tbody></table>
<table><caption>**Table D.3. continued.**</caption><thead><tr><th>HJD<br/>[-2450000]</th><th>B<sub>i</sub><br/>[G]</th><th>B<sub>f</sub><br/>[kG]</th><th>Instrument</th></tr></thead><tbody><tr><td>4134.6112</td><td>-271.7 ± 11.4</td><td>...</td><td>Narval</td></tr><tr><td>4135.6217</td><td>-231.4 ± 10.4</td><td>...</td><td>Narval</td></tr><tr><td>4136.5925</td><td>-294.9 ± 13.0</td><td>...</td><td>Narval</td></tr><tr><td>4276.7715</td><td>-261.0 ± 8.2</td><td>...</td><td>ESPaDOns</td></tr><tr><td>4485.5177</td><td>-290.1 ± 13.2</td><td>3.18<sup>+0.17</sup><sub>-0.20</sub></td><td>Narval</td></tr><tr><td>4489.5683</td><td>-248.6 ± 10.3</td><td>3.10<sup>+0.17</sup><sub>-0.18</sub></td><td>Narval</td></tr><tr><td>4492.5379</td><td>-285.3 ± 10.7</td><td>...</td><td>Narval</td></tr><tr><td>4493.5486</td><td>-204.6 ± 10.6</td><td>3.03<sup>+0.18</sup><sub>-0.19</sub></td><td>Narval</td></tr><tr><td>4495.5611</td><td>-227.2 ± 12.6</td><td>2.81<sup>+0.19</sup><sub>-0.19</sub></td><td>Narval</td></tr><tr><td>4499.5675</td><td>-256.8 ± 11.1</td><td>...</td><td>Narval</td></tr><tr><td>4501.5473</td><td>-288.6 ± 12.1</td><td>...</td><td>Narval</td></tr><tr><td>4502.5475</td><td>-202.8 ± 9.7</td><td>...</td><td>Narval</td></tr><tr><td>4506.5576</td><td>-218.9 ± 10.0</td><td>2.84<sup>+0.19</sup><sub>-0.20</sub></td><td>Narval</td></tr><tr><td>4508.5516</td><td>-265.5 ± 11.0</td><td>...</td><td>Narval</td></tr><tr><td>4509.5564</td><td>-219.6 ± 12.2</td><td>2.90<sup>+0.18</sup><sub>-0.20</sub></td><td>Narval</td></tr><tr><td>4510.5523</td><td>-286.6 ± 15.3</td><td>...</td><td>Narval</td></tr><tr><td>4511.5694</td><td>-200.2 ± 10.5</td><td>2.86<sup>+0.17</sup><sub>-0.17</sub></td><td>Narval</td></tr><tr><td>4512.5537</td><td>-296.7 ± 10.8</td><td>2.99<sup>+0.16</sup><sub>-0.16</sub></td><td>Narval</td></tr><tr><td>5896.7560</td><td>-249.3 ± 11.5</td><td>3.22<sup>+0.15</sup><sub>-0.16</sub></td><td>Narval</td></tr><tr><td>5934.6407</td><td>-247.0 ± 9.5</td><td>3.39<sup>+0.11</sup><sub>-0.12</sub></td><td>Narval</td></tr><tr><td>5935.6765</td><td>-225.9 ± 9.3</td><td>3.50<sup>+0.09</sup><sub>-0.11</sub></td><td>Narval</td></tr><tr><td>5936.6050</td><td>-241.7 ± 10.9</td><td>3.38<sup>+0.09</sup><sub>-0.10</sub></td><td>Narval</td></tr><tr><td>5937.7575</td><td>-254.6 ± 11.2</td><td>3.44<sup>+0.08</sup><sub>-0.10</sub></td><td>Narval</td></tr><tr><td>5938.6659</td><td>-243.9 ± 8.8</td><td>3.31<sup>+0.08</sup><sub>-0.10</sub></td><td>Narval</td></tr><tr><td>5939.6031</td><td>-234.2 ± 10.6</td><td>3.34<sup>+0.08</sup><sub>-0.10</sub></td><td>Narval</td></tr><tr><td>5940.6416</td><td>-203.8 ± 9.1</td><td>3.31<sup>+0.09</sup><sub>-0.09</sub></td><td>Narval</td></tr><tr><td>5941.6411</td><td>-232.7 ± 9.7</td><td>3.36<sup>+0.11</sup><sub>-0.13</sub></td><td>Narval</td></tr><tr><td>5942.6256</td><td>-220.5 ± 10.2</td><td>3.30<sup>+0.09</sup><sub>-0.10</sub></td><td>Narval</td></tr><tr><td>7435.7957</td><td>-177.0 ± 5.5</td><td>3.68<sup>+0.07</sup><sub>-0.08</sub></td><td>ESPaDOnS</td></tr><tr><td>7436.8831</td><td>-165.8 ± 5.3</td><td>3.67<sup>+0.08</sup><sub>-0.09</sub></td><td>ESPaDOnS</td></tr><tr><td>7441.8954</td><td>-181.4 ± 5.3</td><td>3.57<sup>+0.06</sup><sub>-0.07</sub></td><td>ESPaDOnS</td></tr><tr><td>7443.0074</td><td>-166.1 ± 5.3</td><td>3.60<sup>+0.09</sup><sub>-0.10</sub></td><td>ESPaDOnS</td></tr><tr><td>7447.9103</td><td>-160.7 ± 6.4</td><td>3.65<sup>+0.14</sup><sub>-0.16</sub></td><td>ESPaDOnS</td></tr><tr><td>7449.0011</td><td>-182.9 ± 7.7</td><td>3.66<sup>+0.11</sup><sub>-0.13</sub></td><td>ESPaDOnS</td></tr><tr><td>7449.9154</td><td>-154.6 ± 6.5</td><td>3.66<sup>+0.12</sup><sub>-0.15</sub></td><td>ESPaDOnS</td></tr><tr><td>7450.8328</td><td>-190.5 ± 6.9</td><td>3.70<sup>+0.09</sup><sub>-0.11</sub></td><td>ESPaDOnS</td></tr><tr><td>7495.8253</td><td>-192.9 ± 8.8</td><td>3.66<sup>+0.12</sup><sub>-0.15</sub></td><td>ESPaDOnS</td></tr><tr><td>7498.7442</td><td>-170.6 ± 5.8</td><td>3.68<sup>+0.09</sup><sub>-0.09</sub></td><td>ESPaDOnS</td></tr><tr><td>8588.7573</td><td>-194.7 ± 17.2</td><td>3.40<sup>+0.08</sup><sub>-0.09</sub></td><td>SPIRou</td></tr></tbody></table>
<table><caption>**Table D.3.** continued.</caption><thead><tr><th>HJD<br/>[-2450000]</th><th>B<sub>l</sub><br/>[G]</th><th>B<sub>f</sub><br/>[kG]</th><th>Instrument</th></tr></thead><tbody><tr><td>8590.0016</td><td>-212.3±19.9</td><td>...</td><td>SPIRou</td></tr><tr><td>8592.9410</td><td>-206.5±17.4</td><td>3.33<sup>+0.09</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8593.8715</td><td>-202.5±16.5</td><td>3.37<sup>+0.11</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8594.7370</td><td>-234.3±16.8</td><td>3.42<sup>+0.10</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8595.8726</td><td>-211.7±16.4</td><td>3.47<sup>+0.10</sup><sub>-0.13</sub></td><td>SPIRou</td></tr><tr><td>8596.7915</td><td>-217.9±16.0</td><td>3.37<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8597.7516</td><td>-219.1±18.6</td><td>3.33<sup>+0.10</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8598.9026</td><td>-240.0±16.9</td><td>3.49<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8599.8411</td><td>-263.6±17.6</td><td>3.36<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8600.8441</td><td>-182.9±16.7</td><td>3.36<sup>+0.10</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8604.8725</td><td>-190.5±16.3</td><td>3.27<sup>+0.11</sup><sub>-0.12</sub></td><td>SPIRou</td></tr><tr><td>8618.7549</td><td>-216.5±15.1</td><td>3.27<sup>+0.08</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8647.7665</td><td>-217.0±12.9</td><td>3.26<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8648.7341</td><td>-221.3±13.1</td><td>3.38<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8649.7480</td><td>-230.0±15.3</td><td>3.23<sup>+0.10</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8650.7343</td><td>-194.7±13.6</td><td>3.49<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8651.7509</td><td>-209.4±16.0</td><td>3.49<sup>+0.08</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8653.7357</td><td>-227.2±13.6</td><td>3.38<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8655.7575</td><td>-250.3±20.3</td><td>3.42<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8773.1472</td><td>-223.6±14.4</td><td>3.39<sup>+0.10</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8788.1497</td><td>-223.7±13.4</td><td>3.32<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8789.1423</td><td>-184.8±13.7</td><td>3.43<sup>+0.10</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8790.1526</td><td>-245.2±15.4</td><td>3.50<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8791.1225</td><td>-192.7±13.8</td><td>3.50<sup>+0.12</sup><sub>-0.12</sub></td><td>SPIRou</td></tr><tr><td>8792.1572</td><td>-256.3±16.3</td><td>3.46<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8793.1506</td><td>-204.8±13.6</td><td>3.50<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8794.1534</td><td>-246.5±13.5</td><td>3.53<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8795.1283</td><td>-209.1±15.4</td><td>3.48<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8797.0881</td><td>-189.3±25.7</td><td>3.51<sup>+0.10</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8798.1635</td><td>-202.8±17.8</td><td>3.63<sup>+0.10</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8801.1190</td><td>-254.3±14.6</td><td>3.56<sup>+0.11</sup><sub>-0.12</sub></td><td>SPIRou</td></tr><tr><td>8802.0961</td><td>-186.1±14.1</td><td>3.58<sup>+0.08</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8803.0577</td><td>-208.3±7.0</td><td>3.67<sup>+0.13</sup><sub>-0.14</sub></td><td>ESPaDOnS</td></tr><tr><td>8804.0967</td><td>-189.1±7.2</td><td>3.58<sup>+0.10</sup><sub>-0.11</sub></td><td>ESPaDOnS</td></tr><tr><td>8807.0890</td><td>-178.9±12.9</td><td>3.65<sup>+0.07</sup><sub>-0.09</sub></td><td>ESPaDOnS</td></tr><tr><td>8807.1062</td><td>-198.8±11.5</td><td>3.61<sup>+0.08</sup><sub>-0.09</sub></td><td>ESPaDOnS</td></tr><tr><td>8807.1223</td><td>-169.3±10.1</td><td>3.63<sup>+0.07</sup><sub>-0.07</sub></td><td>ESPaDOnS</td></tr><tr><td>8809.1547</td><td>-161.1±6.4</td><td>3.57<sup>+0.08</sup><sub>-0.09</sub></td><td>ESPaDOnS</td></tr><tr><td>8823.1385</td><td>-254.6±30.7</td><td>3.48<sup>+0.08</sup><sub>-0.09</sub></td><td>SPIRou</td></tr></tbody></table>
<table><caption>**Table D.3.** continued.</caption><thead><tr><th>HJD<br/>[-2450000]</th><th>B<sub>i</sub><br/>[G]</th><th>B<sub>f</sub><br/>[kG]</th><th>Instrument</th></tr></thead><tbody><tr><td>8823.1539</td><td>-242.9±42.2</td><td>3.69<sup>+0.10</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8825.1357</td><td>-193.2±22.1</td><td>3.60<sup>+0.08</sup><sub>-0.08</sub></td><td>SPIRou</td></tr><tr><td>8826.1084</td><td>-236.2±14.4</td><td>3.41<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8827.1033</td><td>-178.9±13.7</td><td>3.42<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8828.0574</td><td>-249.2±13.9</td><td>3.40<sup>+0.11</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8829.1286</td><td>-152.6±14.4</td><td>3.44<sup>+0.10</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8830.1115</td><td>-248.1±14.8</td><td>3.35<sup>+0.10</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8875.0136</td><td>-176.0±12.2</td><td>3.43<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8876.0104</td><td>-118.6±13.3</td><td>3.54<sup>+0.08</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8877.0098</td><td>-238.7±14.5</td><td>3.40<sup>+0.10</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8884.8515</td><td>-122.9±12.6</td><td>3.42<sup>+0.10</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8895.8233</td><td>-132.6±11.9</td><td>3.40<sup>+0.10</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8896.8939</td><td>-219.3±13.7</td><td>3.44<sup>+0.11</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8897.8233</td><td>-157.3±12.7</td><td>3.43<sup>+0.09</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8898.8635</td><td>-216.3±13.3</td><td>3.44<sup>+0.10</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8920.8895</td><td>-176.4±12.6</td><td>3.42<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8977.7467</td><td>-177.7±11.8</td><td>3.27<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8978.9018</td><td>-114.9±11.5</td><td>3.33<sup>+0.09</sup><sub>-0.08</sub></td><td>SPIRou</td></tr><tr><td>8981.8996</td><td>-214.6±14.0</td><td>3.28<sup>+0.10</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>8982.9020</td><td>-81.2±11.5</td><td>3.43<sup>+0.11</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>8983.8209</td><td>-190.0±11.2</td><td>3.27<sup>+0.08</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>8984.9069</td><td>-100.5±16.8</td><td>3.44<sup>+0.11</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>9000.7614</td><td>-50.4±10.5</td><td>3.26<sup>+0.12</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>9001.8080</td><td>-217.7±17.7</td><td>3.32<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>9002.7618</td><td>-61.7±11.9</td><td>3.20<sup>+0.09</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>9003.8140</td><td>-199.2±12.7</td><td>3.18<sup>+0.10</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>9004.7681</td><td>-110.4±11.8</td><td>3.26<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>9005.7750</td><td>-96.4±15.5</td><td>3.30<sup>+0.09</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>9006.8244</td><td>-119.8±14.9</td><td>3.30<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>9007.7856</td><td>-70.0±15.1</td><td>3.30<sup>+0.10</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>9008.7478</td><td>-178.4±16.2</td><td>3.28<sup>+0.10</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>9008.7522</td><td>-185.2±16.6</td><td>3.30<sup>+0.09</sup><sub>-0.10</sub></td><td>SPIRou</td></tr><tr><td>9009.7672</td><td>-64.5±18.5</td><td>3.35<sup>+0.10</sup><sub>-0.11</sub></td><td>SPIRou</td></tr><tr><td>9010.7848</td><td>-194.7±13.2</td><td>3.17<sup>+0.10</sup><sub>-0.09</sub></td><td>SPIRou</td></tr><tr><td>9154.1315</td><td>-62.2±10.1</td><td>3.31<sup>+0.07</sup><sub>-0.07</sub></td><td>SPIRou</td></tr><tr><td>9157.1479</td><td>-46.1±9.9</td><td>3.28<sup>+0.08</sup><sub>-0.09</sub></td><td>SPIRou</td></tr></tbody></table>