# Seismic differences between solar magnetic cycles 23 and 24 for low-degree modes

R. A. García, S. N. Breton, D. Salabert, S. C. Tripathy, K. Jain, S. Mathur, E. Panetier

► To cite this version:

R. A. García, S. N. Breton, D. Salabert, S. C. Tripathy, K. Jain, et al. Seismic differences between solar magnetic cycles 23 and 24 for low-degree modes. Astronomy & Astrophysics - A&A, 2024, 691, pp.L20. 10.1051/0004-6361/202451949 . cea-04799708

HAL Id: cea-04799708

https://cea.hal.science/cea-04799708v1

Submitted on 23 Nov 2024

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 Attribution 4.0 International License
LETTER TO THE EDITOR

Seismic differences between solar magnetic cycles 23 and 24 for
low-degree modes

R. A. García<fnref n="1" />,*, S. N. Breton<fnref n="2" />,<fnref n="3" />, D. Salabert<fnref n="4" />, S. C. Tripathy<fnref n="5" />, K. Jain<fnref n="5" />, S. Mathur<fnref n="6" />,<fnref n="7" />, and E. Panetier<fnref n="3" />

Received 21 August 2024 / Accepted 25 October 2024

ABSTRACT

Solar magnetic activity follows regular cycles of about 11 years with an inversion of polarity in the poles every ~22 years. This changing surface magnetism impacts the properties of the acoustic modes. The acoustic mode frequency shifts are a good proxy of the magnetic cycle. In this Letter we investigate solar magnetic activity cycles 23 and 24 through the evolution of the frequency shifts of low-degree modes (ℓ = 0, 1, and 2) in three frequency bands. These bands probe properties between 74 and 1575 km beneath the surface. The analysis was carried out using observations from the space instrument Global Oscillations at Low Frequency and the ground-based Birmingham Solar Oscillations Network and Global Oscillation Network Group. The frequency shifts of radial modes suggest that changes in the magnetic field amplitude and configuration likely occur near the Sun's surface rather than near its core. The maximum shifts of solar cycle 24 occurred earlier at mid and high latitudes (relative to the equator) and about 1550 km beneath the photosphere. At this depth but near the equator, this maximum aligns with the surface activity but has a stronger magnitude. At around 74 km deep, the behaviour near the equator mirrors the behaviour at the surface, while at higher latitudes, it matches the strength of cycle 23.

**Key words.** methods: data analysis – Sun: activity – Sun: helioseismology

**1. Introduction**

Main-sequence solar-like stars have a convective envelope where the interaction between convection, rotation, and magnetic fields results in a magnetic activity cycle (e.g. Brun & Browning 2017). However, the reasons why some stars with similar fundamental properties show very different magnetic activity levels are still unclear (e.g. Baliunas et al. 1995; Egeland 2017; See et al. 2021; Santos et al. 2024). Large spectroscopic surveys, such as the one led at the Mount Wilson Observatory (Wilson 1978) or at Lowell Observatory (Hall et al. 2007), and surveys that use spectropolarimetry (e.g. Marsden et al. 2014) or space-based photometric observations (García et al. 2014; McQuillan et al. 2014; Mathur et al. 2014; Reinhold et al. 2017; Santos et al. 2019, 2021) show a large diversity of magnetic variability (cyclic, flat, or variable stars; Baliunas & Vaughan 1985). Recent studies even suggest that the magnetic solar cycle might soon stop or change its nature (e.g. Metcalfe 2018; McIntosh et al. 2019; Noraz et al. 2022, 2024). Thus, it is necessary to keep monitoring the temporal evolution of the Sun with a broad set of observables at different heights, from the corona to the sub-surface layers, to impose tighter constraints on current physical models.

Over the last 50 years, helioseismology, and more recently
asteroseismology, has been able to probe the properties
of solar and stellar interiors with high precision (e.g.
García & Ballot 2019; Christensen-Dalsgaard 2021). Surface
magnetism perturbs the properties of the acoustic oscilla-
tion modes (p modes) in several ways. In particular, p-
mode properties are affected by magnetic activity; this was
first detected in the Sun (e.g. Pallé et al. 1989; Elsworth et al.
1990; Anguera Gubau et al. 1992; Jiménez-Reyes et al. 2003;
Tripathy et al. 2010; Howe et al. 2017, 2022a) and later in solar-
like stars (e.g. García et al. 2010; Salabert et al. 2011a, 2016;
Mathur et al. 2013, 2019; Kiefer et al. 2017; Karoff et al. 2018;
Santos et al. 2018). Mode frequencies are shifted and correlated
with the magnetic cycle, while the mode amplitudes are anti-
correlated. These frequency shifts can modify the age estimation
of a star like the Sun, with variations of up to 6.5% between solar
minima and maxima (Bétrisey et al. 2024). Studies focused on
the low-activity phase between cycles also suggest complicated
spatial distributions of magnetic fields beneath the solar surface
(Salabert et al. 2009; Jain et al. 2022).

Recently, a detailed seismic analysis that included
medium- and high-degree modes collected by the Solar
Oscillation Imager/Michelson Doppler Imager (SOI/MDI;
Scherrer et al. 1995), the Helioseismic and Magnetic Imager
(HMI; Scherrer et al. 2012), and the Global Oscillation Network

<footnote n="1">¹ Université Paris-Saclay, Université Paris Cité, CEA, CNRS, AIM, 91191 Gif-sur-Yvette, France</footnote>

<footnote n="2">² INAF – Osservatorio Astrofisico di Catania, Via S. Sofia, 78, 95123 Catania, Italy</footnote>

<footnote n="3">³ Université Paris Cité, Université Paris-Saclay, CEA, CNRS, AIM, 91191 Gif-sur-Yvette, France</footnote>

<footnote n="4">⁴ Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Laboratoire Lagrange, Nice, France</footnote>

<footnote n="5">⁵ National Solar Observatory, 3665 Discovery Drive, Boulder, CO 80303, USA</footnote>

<footnote n="6">⁶ Instituto de Astrofísica de Canarias (IAC), E-38205 La Laguna, Tenerife, Spain</footnote>

<footnote n="7">⁷ Universidad de La Laguna (ULL), Departamento de Astrofísica, E-38206 La Laguna, Tenerife, Spain</footnote>

<footnote>* Corresponding author; rgarcia@cea.fr</footnote>
Group (GONG; Harvey et al. 1996) by Basu (2021) has revealed significant changes in the solar convection zone (and even below), indicating possible small changes in the position of the tachocline. Moreover, Inceoglu et al. (2021) studied solar quasi-biennial oscillations (QBOs; Sakurai 1979) in the rotation rate residuals and conclude that QBO-like signals were present at different latitudinal bands, with their amplitudes increasing with depth. In addition, Kurtz et al. (2014) showed that the quasi-biennial variability in the solar interior depended on the cycle (see also Mehta et al. 2022). In addition, a recent study by Jain et al. (2023) found different QBO periods in cycles 23 and 24.

Unfortunately, in the case of the asteroseismology of solar-like stars, only low-degree modes can be detected, and the capabilities to diagnose the impact of the surface magnetic activity at different layers are limited. Nevertheless, using these low-degree *p* modes, it is still possible to evaluate structural perturbations in the ~3000 km beneath the surface (Basu et al. 2012; Salabert et al. 2015; Jain et al. 2018). In this work we investigated the solar seismic signatures of cycles 23 and 24 (from 1996 to 2020) through the evolution of averaged frequency shifts, $\langle\delta\nu\rangle$, of low-degree modes ($\ell = 0, 1$, and 2) in different frequency ranges. This allowed us to study, on one hand, the dependence on latitude because radial modes are more sensitive to mid and high latitudes, while dipolar and quadrupolar sectoral modes ($\ell = |m|$) are more sensitive towards the equator and with increasing $\ell$ (see Jiménez-Reyes et al. 2004a, and Fig. 3 in Chaplin & Basu 2014). On the other hand, by averaging the $\langle\delta\nu\rangle$ in different frequency bands, we were able to study the sensitivity of the magnetic perturbation to depth. In this context, we note that Simoniello et al. (2016) analysed intermediate- and high-degree modes of GONG data for a part of cycle 23 using different frequency ranges and harmonic degrees and conclude that the frequency shifts behave differently at low and high latitudes.

The layout of this Letter is as follows. In Sect. 2 we present the observations and the procedure for extracting the *p*-mode frequency shifts. In Sect. 3 we present the results and discuss them in Sect. 4.

## 2. Observations and data analysis

We used three different datasets of velocity measurements. The first one was collected by the Global Oscillations at Low Frequency (GOLF) instrument<fnref n="1" /> (Gabriel et al. 1995, 1997) on board the Solar and Heliospheric Observatory (SoHO; Domingo et al. 1995). We obtained 24.5 years of continuous data – from April 11, 1996, until September 4, 2020 – of Doppler velocity time series following García et al. (2005) and ensuring a proper timing of the measurements (for more details, see Appourchaux et al. 2018; Breton et al. 2022a). The second dataset was collected with the Birmingham Solar Oscillation Network (BiSON<fnref n="2" />; Hale et al. 2016) starting at the same time as GOLF but extending only to April 4, 2020 and corrected following the procedure described in Davies et al. (2014). The third dataset is composed of Doppler velocity measurements from the GONG network<fnref n="3" /> (Jain et al. 2021). The series in this third dataset start on May 7, 1995, and are available through October 9, 2020, with a median duty cycle value of 87% (varying between 78% and 92%).

Both the GOLF and BiSON datasets are divided into consecutive sub-series of 365 days shifted by 91.25 days, while the GONG dataset comprises 360 days shifted by 72 days. Series with a duty cycle below 80%, 70%, and 56% respectively for GOLF, GONG, and BiSON were removed. 15 consecutive radial orders between 1800 µHz and 3790 µHz were fitted (as explained in Appendix A) following a traditional, deterministic, maximum likelihood (ML) algorithm or through the apollinaire module<fnref n="4" /> (Breton et al. 2022b), which implements a Markov chain Monte Carlo (MCMC) Bayesian sampling of the model parameter posterior distribution (Foreman-Mackey et al. 2013). The lowest fitted radial order is *n* = 12 for $\ell$ = 0 and 1, and *n* = 11 for $\ell$ = 2. We did not use $\ell$ = 3 modes in our analysis because of their low signal-to-noise ratio in Sun-as-a-star data.

Temporal variations in mode frequencies, $\langle\delta\nu\rangle$, were defined as the differences between the frequencies observed at a given time and the average of the frequencies of years 1996–1997 (independent series 1 and 5), corresponding to the minimum of activity of cycles 22–23. The frequency-shift errors were computed as

$$ \sigma(\delta\nu) = \sqrt{\frac{1}{\sum \sigma(\nu)^2}} \qquad (1) $$

Once the frequency shifts were obtained, they were averaged in three different frequency bands following Basu et al. (2012): (i) the low-frequency band (LFB) from 1800 to 2450 µHz; (ii) the medium-frequency band (MFB) from 2450 to 3110 µHz; and (iii) the high-frequency band (HFB) from 3110 to 3790 µHz. The formal uncertainties resulting from the peak-fitting analysis were used as weights in the computation of the averages of the different degree modes and frequency bands. The resultant $\langle\delta\nu_{\ell=0,1,2}\rangle$ for the three instruments are given in Appendix B (Tables B.1, B.2, and B.3). Using model BS05 from Bahcall & Serenelli (2005), we computed the average sensitivity for each degree in each frequency band. The results are shown in Table 1 and Fig. 1. Since the kernels for individual ($n, \ell$) modes have different depth sensitivities, we calculated these values by averaging them over all the *n* and $\ell$ values in each selection taken from the observations. Our kernels are slightly different from those in Basu et al. (2012) and Salabert et al. (2015) because the modes averaged in each frequency band are not the same. For example, Basu et al. (2012) used the same number of overtones (i.e. four) in all three frequency ranges, while we used five consecutive radial orders in each category. This led to the use of a different set of modes in this study and hence a different depth sensitivity.

Finally, because we are interested in comparing solar cycles 23 and 24, we removed the contribution of the QBO (e.g. Fletcher et al. 2010; Jain et al. 2011; Simoniello et al. 2012) by smoothing the data with a boxcar of 2.5 years. To reduce any border effects associated with the filtering, we extended the $\langle\delta\nu\rangle$ by half of the filter width, assuming them to be symmetric with respect to the beginning and the end of the series.

In addition, mean values of daily measurements of the 10.7 cm radio flux<fnref n="5" />, $F_{10.7}$, were used as a proxy of the solar surface activity since it has been shown that this magnetic activity proxy is the best correlated with the mean frequency shifts (e.g.

<footnote n="1">¹ http://irfu.cea.fr/Phocea/Vie_des_labos/Ast/ast_technique.php?id_ast=3842</footnote>

<footnote n="2">² http://bison.ph.bham.ac.uk/portal/timeseries</footnote>

<footnote n="3">³ https://nispdata.nso.edu/ftp/TSERIES/vmt/</footnote>

<footnote n="4">⁴ Documentation is available at https://apollinaire.readthedocs.io/en/latest/</footnote>

<footnote n="5">⁵ The 10.7 cm radio flux data are available from the LASP Interactive Solar Irradiance Datacenter at https://lasp.colorado.edu/lisird/data/pentictonRadioFlux/</footnote>
<table><caption>**Table 1.** Depths from the surface of the sensitivity kernels.</caption><thead><tr><th>ν<br>band</th><th>Depth<br>(ℓ = 0)</th><th>Depth<br>(ℓ = 1)</th><th>Depth<br>(ℓ = 2)</th><th>Basu et al. (2012)</th></tr></thead><tbody><tr><td>LFB</td><td>0.99777 R<sub>☉</sub><br>(1550 km)</td><td>0.99798 R<sub>☉</sub><br>(1402 km)</td><td>0.99774 R<sub>☉</sub><br>(1575 km)</td><td>0.9963 R<sub>☉</sub><br>(2576 km)</td></tr><tr><td>MFB</td><td>0.99932 R<sub>☉</sub><br>(474 km)</td><td>0.99941 R<sub>☉</sub><br>(412 km)</td><td>0.99930 R<sub>☉</sub><br>(487 km)</td><td>0.9981 R<sub>☉</sub><br>(1323 km)</td></tr><tr><td>HFB</td><td>0.99989 R<sub>☉</sub><br>(74 Km)</td><td>0.99989 R<sub>☉</sub><br>(74 km)</td><td>0.99989 R<sub>☉</sub><br>(75 km)</td><td>0.9989 R<sub>☉</sub><br>(766 km)</td></tr></tbody></table>

**Notes.** The depths were computed with the BS05 model from Bahcall & Serenelli (2005) for modes $ℓ = 0, 1$, and $2$ and the three frequency bands. The last column indicates the depth computed by Basu et al. (2012).

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Integrated sensitivity, $\langle \mathcal{H} \rangle$, as a function of the relative radius for the LFB, MFB, and HFB (top, middle, and bottom panels, respectively) from the BS05 model in Bahcall & Serenelli (2005). Vertical dotted lines indicate the depth of the sensitivity kernels for $ℓ = 0$ modes.</figcaption></figure>

Jain et al. 2009; Howe et al. 2017). The same 2.5-a-boxcar filter was applied to remove the QBO contribution from this proxy. Using the filtered $F_{10.7}$, we determined the reference dates for the maximum of cycles 23 and 24 and the minimum in between, respectively March 2001, January 2014, and July 2008.

By comparing the mode parameters from MDI and HMI, which used two different spectral lines (Ni I 676.8 nm and Fe I 617.3 nm, respectively), Larson & Schou (2018) found that differences in frequencies and splitting coefficients were not significant. In addition, no differences were found when comparing the temporal evolution of p-mode parameters of GOLF and BiSON (Jiménez-Reyes et al. 2004b) except for the mode asymmetry (Jiménez-Reyes et al. 2007). Thus, we conclude that our results should not be affected by the fact that they are based on the study of the temporal evolution of $\langle \delta v \rangle$ involving different instruments observing the signals at slightly different heights.

### 3. Results: Temporal evolution of GOLF frequency shifts

The temporal evolution of the 2.5-a-boxcar smoothed $\langle \delta v \rangle$ from GOLF is shown in Fig. 2. Each column corresponds to a different mode ($ℓ = 0, 1$, and $2$), and each row corresponds to the aforementioned frequency bands. A figure containing the non-smoothed $\langle \delta v \rangle$ is shown in Appendix C. In each panel, the $F_{10.7}$ radio flux is scaled to the $\langle \delta v \rangle$ using the rising phase of cycle 23. For the sake of clarity and because the qualitative results are the same for the three instruments, only the GOLF results are presented. In Appendix D, results obtained with the other two datasets are compared and discussed.

Overall, the $\langle \delta v \rangle$ follow the evolution of the $F_{10.7}$, but there are some notable differences:

- The rising phase of cycle 24 in the LFB $\langle \delta v_{ℓ=0} \rangle$ is reached in 2012, two years earlier than $F_{10.7}$.

- Cycle 24 is known to have a lower amplitude than cycle 23, as shown for example by Howe et al. (2022a) and by the $F_{10.7}$. In this respect, $\langle \delta v \rangle$ modulations have surprisingly large amplitudes in cycle 24 for some frequency bands and degrees, $ℓ$. The amplitude of the LFB $\langle \delta v_{ℓ=1} \rangle$ of cycle 24 is nearly twice that of cycle 23, while for the MFB $\langle \delta v_{ℓ=2} \rangle$, the amplitude of cycle 24 is comparable to that of cycle 23. In all of these cases, the amplitude of cycle 24 is clearly higher than that of the $F_{10.7}$. The MFB $\langle \delta v_{ℓ=0} \rangle$ reaches its minimum between the two cycles around a year earlier than the $F_{10.7}$.

- The MFB $\langle \delta v_{ℓ=1} \rangle$ and HFB $\langle \delta v_{ℓ=0} \rangle$ are both greater than $F_{10.7}$ in cycle 24.

- With larger and larger error bars in the LFB $\langle \delta v_{ℓ=2} \rangle$ since 2005, due to the lower counting rates of the GOLF instrument, we cannot extract information on cycle 24.

- For the remaining three cases, the MFB $\langle \delta v_{ℓ=0} \rangle$ and the HFB $\langle \delta v_{ℓ=1} \rangle$ and $\langle \delta v_{ℓ=2} \rangle$, they closely follow $F_{10.7}$.

### 4. Discussion

Recent seismic analyses of averaged frequency shifts for both low-degree modes (Howe et al. 2022a) and medium- and high-degree modes (Jain et al. 2022) show that the influence of the magnetic perturbations on the acoustic modes are different during consecutive solar cycles. Although it is usually assumed that the main perturbations on the mode frequencies are located in the sub-surface layers, it has also been speculated that a fossil magnetic field in the core of the Sun could play a crucial role in the generation of the solar magnetic cycle and thus could also modify the properties of the acoustic modes (e.g. Sonett 1983; Mursula et al. 2001). Such magnetic fields have not been observed in main-sequence solar-like stars yet, but they could modify the structure of stellar oscillations (e.g. Goode & Thompson 1992; Loi & Papaloizou 2018; Bugnet et al. 2021). They have been observed in the core of red giants (Li et al. 2022, 2024; Deheuvels et al. 2023).

As discussed in the previous section, the maximum of the $\langle \delta v_{ℓ=0} \rangle$ in cycle 24 arrives two years earlier in the LFB. This means that something changed between the two cycles at a depth of ~1550 km below the photosphere (see Table 1 and Fig. 1). The dominant perturbation in the LFB $\langle \delta v_{ℓ=0} \rangle$ is weighted towards mid to high latitudes, while the maximum of $\langle \delta v_{ℓ=1} \rangle$ is synchronised with $F_{10.7}$ except in regions close to the equator, where $\langle \delta v_{ℓ=1} \rangle$ are the most sensitive. Therefore, in cycle 24, the magnetic perturbation modifying $\langle \delta v_{ℓ=0} \rangle$ reaches an extended 4-a maximum earlier (~2010) at high latitudes and at a depth of
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Temporal evolution of the averaged and smoothed GOLF frequency shifts ($\langle\delta\nu\rangle$, dark blue lines) of $\ell = 0, 1$, and 2 modes (left, central, and right panels respectively). The top, middle, and bottom rows correspond to the three frequency bands considered in the analysis: 1800–2450 (LFB), 2450–3110 (MFB), and 3110–3790 µHz (HFB), respectively. Light blue regions represent the errors in $\langle\delta\nu\rangle$. Orange lines depict the scaled smoothed $F_{10.7}$ magnetic proxy. Vertical dashed and dot-dashed lines depict the maxima of cycles 23 and 24 and the minimum in between. The horizontal blue arrow indicates the change in sensitivity from the pole to the equator, and the vertical green arrow depicts the maximum sensitivity to the depth of the subsurface layers.</figcaption></figure>

~1550 km beneath the surface. This perturbation reaches shallow layers (>1402 km) at the same time as the photosphere, with a maximum of around ~2014. At these mid and high latitudes, just below the photosphere (~74 km), the HFB of $\langle\delta\nu_{\ell=0}\rangle$ shows stronger variations than the surface activity represented by the $F_{10.7}$ in cycle 24. Based on our analysis of the radial modes, we do not find evidence supporting that the changes comes from an hypothetical magnetic field in the core of the Sun as the MFB $\langle\delta\nu_{\ell=0}\rangle$ seems to follow the evolution of $F_{10.7}$ very closely. The fact that all the radial modes reach the centre rules out a change coming from the core of the Sun. Only the upper turning points are different. Thus, the change in the perturbation should be located at the subsurface layers, which likely have a complicated vertical structure.

Near the solar equator, we see that the LFB of $\langle\delta\nu_{\ell=1}\rangle$ is stronger in cycle 24 than in cycle 23. The perturbation of the dipolar modes at ~1402 km is larger in the later cycle. The differences in the $\langle\delta\nu\rangle$ of the $\ell = 1$ and 2 modes compared to the $F_{10.7}$ in cycle 24 are also large in the MFB. In this cycle, the maximum of the dipolar $\langle\delta\nu\rangle$ is slightly smaller than in cycle 23, while it is larger for the quadrupolar modes. That means that the perturbation is larger very close to the equator since the difference in the depth is very small (~412 versus ~487 km). Finally, very close to the surface (~74–75 km), the behaviour of $\langle\delta\nu\rangle$ is almost identical to the surface activity. Despite the different frequency bands and mode degrees, we find similar qualitative results regarding the dependence of the progression of solar cycle 23 to those obtained by Simoniello et al. (2016).

# 5. Conclusions

Similar to previous studies using intermediate- and high-degree modes (e.g. Howe et al. 2022b, and references therein), the present study allowed us to trace how the solar-cycle-induced perturbations evolve with time, latitude, and depth using low-degree modes during solar cycles 23 and 24, and with three different instruments, GOLF, BiSON, and GONG. Our first conclusion is that the temporal evolution of the radial-mode frequency shifts suggests that these changes most likely did not occur in the solar core but rather close to the surface. The second conclusion is that the maximum shifts of cycle 24 seem to arrive earlier at mid and high latitudes compared to cycle 23 and at a depth of around 1550 km. At this depth and near the equator, the maximum of cycle 24 is synchronised with the surface, but the magnitude is stronger than at the surface. Finally, in a band around 74 km beneath the surface and near the equator, the behaviour is similar to that at the surface, while at higher latitudes it matches the strength of cycle 23. The results obtained with the three instruments are in close agreement in most of the cases. However, as shown in Appendix B, the main difference is the presence of a clear QBO pattern in the GONG LFB of $\langle\delta\nu_{\ell=0}\rangle$ that is not seen in the integrated-light instruments. The smoothed frequency shift also behaves differently, closely following the $F_{10.7}$. Our analysis can also be used to test new theories on the origin of the solar dynamo (Vasil et al. 2024).
## Data availability

Full Tables B.1, B.2, and B.3 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/691/L20

**Acknowledgements.** We thank J. Ballot and F. Pérez Hernández for useful comments and discussions. K.J. thanks Sarbani Basu for providing the standard solar model values. The GOLF instrument on board SoHO is a cooperative effort of many individuals, to whom we are indebted. SoHO is a project of international collaboration between ESA and NASA. BiSON is funded by the UK Science and Technology Facilities Council (STFC). This work also utilises GONG data obtained by the NSO Integrated Synoptic Program, managed by the National Solar Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA), Inc. under a cooperative agreement with the National Science Foundation and with contribution from the National Oceanic and Atmospheric Administration. The GONG network of instruments is hosted by the Big Bear Solar Observatory, High Altitude Observatory, Learmonth Solar Observatory, Udaipur Solar Observatory, Instituto de Astrofísica de Canarias, and Cerro Tololo Interamerican Observatory. R.A.G., S.N.B., and E.P. acknowledge the support from the GOLF and PLATO Centre National D'Études Spatiales grants. S.N.B acknowledges support from PLATO ASI-NAF agreement n. 2015-019-R.1-2018. S.M. acknowledges support by the Spanish Ministry of Science and Innovation with the Ramon y Cajal fellowship number RYC-2015-17697 and the grant number PID2019-107187GB-100, and through AEI under the Severo Ochoa Centres of Excellence Programme 2020-2023 (CEX2019-000920-S). S.C.T. and K.J. acknowledge partial funding support from the NASA DRIVE Science Center COFFIES Phase II CAN 80NSSC22M0162. Software: AstroPy (Astropy Collaboration 2013, 2018), Matplotlib (Hunter 2007), NumPy (van der Walt et al. 2011), SciPy (Jones et al. 2001), pandas (Wes McKinney 2010; The pandas development team 2020), emcee (Foreman-Mackey et al. 2013), apollinaire (Breton et al. 2022b), Interactive Data Language (IDL).

## References

Anderson, E. R., Duvall, T. L., Jr., & Jefferies, S. M. 1990, ApJ, 364, 699

Anguera Gubau, M., Palle, P. L., Perez Hernandez, F., Regulo, C., & Roca Cortes, T. 1992, A&A, 255, 363

Appourchaux, T., Boumier, P., Leibacher, J. W., & Corbard, T. 2018, A&A, 617, A108

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

Atropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123

Bahcall, J. N., & Serenelli, A. M. 2005, ApJ, 626, 530

Baliunas, S. L., & Vaughan, A. H. 1985, ARA&A, 23, 379

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

Basu, S. 2021, ApJ, 917, 45

Basu, S., Broomhall, A.-M., Chaplin, W. J., & Elsworth, Y. 2012, ApJ, 758, 43

Benomar, O., Appourchaux, T., & Baudin, F. 2009, A&A, 506, 15

Bétrisey, J., Farnir, M., Breton, S. N., et al. 2024, A&A, 688, L17

Breton, S. N., Pallé, P. L., García, R. A., et al. 2022a, A&A, 658, A27

Breton, S. N., García, R. A., Ballot, J., Delsanti, V., & Salabert, D. 2022b, A&A, 663, A118

Brun, A. S., & Browning, M. K. 2017, Liv. Rev. Sol. Phys., 14, 4

Bugnet, L., Prat, V., Mathis, S., et al. 2021, A&A, 650, A53

Chaplin, W. J., & Basu, S. 2014, Space Sci. Rev., 186, 437

Christensen-Dalsgaard, J. 2021, Liv. Rev. Sol. Phys., 18, 2

Davies, G. R., Chaplin, W. J., Elsworth, Y., & Hale, S. J. 2014, MNRAS, 441, 3009

Deheuvels, S., Li, G., Ballot, J., & Lignières, F. 2023, A&A, 670, L16

Domingo, V., Fleck, B., & Poland, A. I. 1995, Sol. Phys., 162, 1

Egeland, R. 2017, Ph.D. Thesis, Montana State University, Bozeman, Montana, USA

Elsworth, Y., Howe, R., Isaak, G. R., McLeod, C. P., & New, R. 1990, Nature, 345, 322

Fletcher, S. T., Broomhall, A., Salabert, D., et al. 2010, ApJ, 718, L19

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

Gabriel, A. H., Grec, G., Charra, J., et al. 1995, Sol. Phys., 162, 61

Gabriel, A. H., Charra, J., Grec, G., et al. 1997, Sol. Phys., 175, 207

García, R. A., & Ballot, J. 2019, Liv. Rev. Sol. Phys., 16, 4

García, R. A., Turck-Chièze, S., Boumier, P., et al. 2005, A&A, 442, 385

García, R. A., Mathur, S., Salabert, D., et al. 2010, Science, 329, 1032

García, R. A., Ceillier, T., Salabert, D., et al. 2014, A&A, 572, A34

Goode, P. R., & Thompson, M. J. 1992, ApJ, 395, 307

Hale, S. J., Howe, R., Chaplin, W. J., Davies, G. R., & Elsworth, Y. P. 2016, Sol. Phys., 291, 1

Hall, J. C., Lockwood, G. W., & Skiff, B. A. 2007, AJ, 133, 862

Harvey, J. W., Hill, F., Hubbard, R., et al. 1996, Science, 272, 1284

Howe, R., Davies, G. R., Chaplin, W. J., et al. 2017, MNRAS, 470, 1935

Howe, R., Chaplin, W. J., Elsworth, Y. P., Hale, S. J., & Nielsen, M. B. 2022a, MNRAS, 514, 3821

Howe, R., Chaplin, W. J., Christensen-Dalsgaard, J., Elsworth, Y. P., & Schou, J. 2022b, RNAAS, 6, 261

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

Inceoglu, F., Howe, R., & Loto'aniu, P. T. M. 2021, ApJ, 920, 49

Jain, K., Tripathy, S. C., & Hill, F. 2009, ApJ, 695, 1567

Jain, K., Tripathy, S. C., & Hill, F. 2011, ApJ, 739, 6

Jain, K., Tripathy, S. C., Hill, F., et al. 2018, IAU Symp., 340, 27

Jain, K., Tripathy, S. C., Hill, F., & Pevtsov, A. A. 2021, PASP, 133, 105001

Jain, K., Jain, N., Tripathy, S. C., & Dikpati, M. 2022, ApJ, 924, L20

Jain, K., Chowdhury, P., & Tripathy, S. C. 2023, ApJ, 959, 16

Jiménez-Reyes, S. J., García, R. A., Jiménez, A., & Chaplin, W. J. 2003, ApJ, 595, 446

Jiménez-Reyes, S. J., García, R. A., Chaplin, W. J., & Korzennik, S. G. 2004a, ApJ, 610, L65

Jiménez-Reyes, S. J., Chaplin, W. J., Elsworth, Y., & García, R. A. 2004b, ApJ, 604, 969

Jiménez-Reyes, S. J., Chaplin, W. J., Elsworth, Y., et al. 2007, ApJ, 654, 1135

Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy: Open Source Scientific Tools for Python, https://apisemanticscholar.org/CorpusID: 215874460

Karoff, C., Metcalfe, T. S., Santos, Â. R. G., et al. 2018, ApJ, 852, 46

Kiefer, R., Schad, A., Davies, G., & Roth, M. 2017, A&A, 598, A77

Komm, R. W., Gu, Y., Hill, F., Stark, P. B., & Fodor, I. K. 1999, ApJ, 519, 407

Kurtz, D. W., Saio, H., Takata, M., et al. 2014, MNRAS, 444, 102

Larson, T. P., & Schou, J. 2018, Sol. Phys., 293, 29

Li, G., Deheuvels, S., Ballot, J., & Lignières, F. 2022, Nature, 610, 43

Li, G., Deheuvels, S., & Ballot, J. 2024, A&A, 688, A184

Loi, S. T., & Papaloizou, J. C. B. 2018, MNRAS, 477, 5338

Marsden, S. C., Petit, P., Jeffers, S. V., et al. 2014, MNRAS, 444, 3517

Mathur, S., García, R. A., Morgenthaler, A., et al. 2013, A&A, 550, A32

Mathur, S., García, R. A., Ballot, J., et al. 2014, A&A, 562, A124

Mathur, S., García, R. A., Bugnet, L., et al. 2019, Front. Astron. Space Sci., 6, 46

McIntosh, S. W., Leamon, R. J., Egeland, R., et al. 2019, Sol. Phys., 294, 88

McQuillan, A., Mazeh, T., & Aigrain, S. 2014, ApJS, 211, 24

Mehta, T., Jain, K., Tripathy, S. C., et al. 2022, MNRAS, 515, 2415

Metcalfe, T. S. 2018, Phys. Today, 71, 70

Mursula, K., Usoskin, I. G., & Kovaltsov, G. A. 2001, Sol. Phys., 198, 51

Nigam, R., & Kosovichev, A. G. 1998, ApJ, 505, L51

Noraz, Q., Brun, A. S., Strugarek, A., & Depambour, G. 2022, A&A, 658, A144

Noraz, Q., Brun, A. S., Strugarek, A. 2024, A&A, 684, A156

Pallé, P. L., Régulo, C., & Roca Cortés, T. 1989, A&A, 224, 253

Reinhold, T., Cameron, R. H., & Gizon, L. 2017, A&A, 603, A52

Sakurai, K. 1979, Nature, 278, 146

Salabert, D., Chaplin, W. J., Elsworth, Y., New, R., & Verner, G. A. 2007, A&A, 463, 1181

Salabert, D., García, R. A., Pallé, P. L., & Jiménez-Reyes, S. J. 2009, A&A, 504, L1

Salabert, D., Régulo, C., Ballot, J., García, R. A., & Mathur, S. 2011a, A&A, 530, A127

Salabert, D., Ballot, J., & García, R. A. 2011b, A&A, 528, A25

Salabert, D., García, R. A., & Turck-Chièze, S. 2015, A&A, 578, A137

Salabert, D., Régulo, C., García, R. A., et al. 2016, A&A, 589, A118

Santos, Â. R. G., Campante, T. L., Chaplin, W. J., et al. 2018, ApJS, 237, 17

Santos, Â. R. G., García, R. A., Mathur, S., et al. 2019, ApJS, 244, 21

Santos, Â. R. G., Breton, S. N., Mathur, S., & García, R. A. 2021, ApJS, 255, 17

Santos, Â. R. G., Godoy-Rivera, D., Finley, A. J., et al. 2024, Front. Astron. Space Sci., 11, 1356379

Scherrer, P. H., Bogart, R. S., Bush, R. I., et al. 1995, Sol. Phys., 162, 129

Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, Sol. Phys., 275, 207

See, V., Roquette, J., Amard, L., & Matt, S. P. 2021, ApJ, 912, 127

Simoniello, R., Finsterle, W., Salabert, D., et al. 2012, A&A, 539, A135

Simoniello, R., Tripathy, S. C., Jain, K., & Hill, F. 2016, ApJ, 828, 41

Sonett, C. P. 1983, Nature, 306, 670

Tripathy, S. C., Jain, K., Hill, F., & Leibacher, J. W. 2010, ApJ, 711, L84

The pandas development team 2020, https://doi.org/10.5281/zenodo. 3509134

van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Comput. Sci. Eng., 13, 22

Vasil, G. M., Lecoanet, D., Augustson, K., et al. 2024, Nature, 629, 769

Wes McKinney 2010, in Proceedings of the 9th Python in Science Conference, 56

Wilson, O. C. 1978, ApJ, 226, 379
# Appendix A: Extraction of the p-mode parameters

We extracted the parameters of the individual p modes using three different methods that we describe below.

## A.1. GOLF: MCMC fitting

We used the Lorentzian-profile p-mode model as described in Eq. A.2 to perform an MCMC analysis of the GOLF power spectra with the apollinaire library (Breton et al. 2022b) which implements the ensemble sampler provided by the emcee library (Foreman-Mackey et al. 2013) for asteroseismic analyses. Given a set of mode parameters, $\theta$, the MCMC process provides a sampling of the posterior probability, $p(\theta|\mathbf{S}_x)$:

$$p(\theta|\mathbf{S}_x) = p(\mathbf{S}_x|\theta)p(\theta), \qquad (A.1)$$

where $p(\mathbf{S}_x|\theta)$ is the likelihood function and $p(\theta)$ the prior distribution of $\theta$. The marginal likelihood $p(\mathbf{S}_x)$ is omitted in the formula. In order to consider prior uniform distributions while preserving non-informative prior for every sampled parameter (Benomar et al. 2009), we choose to sample the logarithms $\ln \Gamma_{n,l}$ and $\ln H_{n,l}$ rather than sampling directly $\Gamma_{n,l}$ and $H_{n,l}$. The extracted mode parameters are taken as the median of the sampled distribution. The uncertainty we consider is the largest value between the absolute difference of the median value and the 16th centile value of the distribution and the absolute difference of the median value and the 84th centile value of the distribution. A detailed comparison of the ML and MCMC fitting results based on GOLF data can be found in Breton et al. (2022b).

## A.2. BiSON: ML fitting

The power spectrum of each time series was fitted assuming a $\chi^2$ with two degrees of freedom likelihood to estimate the mode parameters of the $l=0, 1, 2$, and 3 modes as described in Salabert et al. (2007). Each mode component of radial order $n$, angular degree $l$, and azimuthal order $m$ was parametrised with an asymmetric Lorentzian profile (Nigam & Kosovichev 1998), as

$$\mathcal{L}_{n,l,m}(\nu) = H_{n,l} \frac{(1 + b_{n,l}x_{n,l})^2 + b_{n,l}^2}{1 + x_{n,l}^2}, \qquad (A.2)$$

where $x_{n,l} = 2(\nu - \nu_{n,l})/\Gamma_{n,l}$, and $\nu_{n,l}$, $\Gamma_{n,l}$, and $H_{n,l}$ represent the mode frequency, linewidth, and height of the spectral density, respectively. The peak asymmetry is described by the parameter $b_{n,l}$. Because of their close proximity in frequency, modes are fitted in pairs (i.e. $l=2$ with 0 and $l=3$ with 1). While each mode parameter within a pair of modes is free, the peak asymmetry is set to be the same within pairs of modes. The $l=4$ and 5 modes were included in the fitted model when they are present within the fitting window, which is proportional to the mode linewidth (i.e. frequency dependent). An additive parameter $B$ is added to the fitted profile representing a constant background noise in the fitting window. Since SoHO and BiSON observe the Sun equatorward, only the $l+|m|$ even components are visible in Sun-as-a-star observations. The amplitude ratios between the $l=0, 1, 2$, and 3 modes and the $m$-height ratios of the $l=2$ and 3 multiplets correspond to those calculated in Salabert et al. (2011b). Finally, the mode parameters were extracted by maximizing the likelihood function, the power spectrum statistics being described by a $\chi^2$ distribution with two degrees of freedom. The natural logarithms of the mode height, linewidth, and background noise were varied resulting in normal distributions. The formal uncertainties

in each parameter were then derived from the inverse Hessian matrix.

## A.3. GONG: ML fitting

We used p-mode frequencies corresponding to the individual multiplets ($\nu_{n,l,m}$) obtained from the GONG time series, where $\nu$ is the frequency, $n$ is the radial order and $m$ is the azimuthal order, running from $-\ell$ to $+\ell$. The multiplets were estimated from the $m-\nu$ power spectra constructed from time series spanning 360 days (10 GONG months, where 1 GONG month corresponds to 36 days) with a spacing of 72 days between two consecutive data sets. The time series were processed through the standard GONG peak-fitting algorithm to compute the power spectra based on the multi-taper spectral analysis coupled with a fast Fourier transform (Komm et al. 1999). Finally, Lorentzian profiles were used to fit the peaks in the $m-\nu$ spectra using a minimisation scheme guided by an initial guess table (Anderson et al. 1990). In order to compare the mode frequencies for low-degree modes obtained from GOLF instrument, we averaged the multiplets over all $m$ values for a given $\ell$ and $n$ value to obtain m-averaged mode frequencies (i.e. quadruple modes were fitted for $m \pm 2$ including 0 and dipole using $m \pm 1$).

# Appendix B: Computed averaged frequency shifts

In Tables B.1, B.2, and B.3, the computed averaged frequency shifts $\langle\delta\nu\rangle$ for modes $\ell=0, 1$, and 2 are shown respectively for GOLF, BiSON, and GONG.

# Appendix C: Non-smoothed GOLF results showing the QBO

Figure C.1 shows the same smoothed frequency shifts and $F_{10.7}$ shown in Fig. 2 but including the non-filtered frequency shifts in dark blue. The behaviour of these frequency shifts is dominated by the presence of the QBO.

# Appendix D: Results for BiSON and GONG

The main differences between the instruments are summarised as follows (see Figs. D.1 and D.2):

- The GONG LFB of $\langle\delta\nu_{\ell=0}\rangle$ shows a clear QBO pattern that is not seen in the Sun-as-a-star instruments. Moreover, the filtered frequency $\langle\delta\nu_{\ell=0}\rangle$ follows perfectly $F_{10.7}$.

- The BiSON LFB of $\langle\delta\nu_{\ell=0}\rangle$ has the same behaviour as GOLF with an even faster rising phase of cycle 24 and a maximum earlier than in $F_{10.7}$.

- The BiSON and GONG LFBs of $\langle\delta\nu_{\ell=1}\rangle$ show similar or slightly higher amplitudes of the maximum of cycle 24 compared to cycle 23.

- The GONG LFB of $\langle\delta\nu_{\ell=2}\rangle$ has a very low signal-to-noise ratio. As a consequence, the scaled $F_{10.7}$ appears nearly flat.

- The BiSON and GONG MFBs of $\langle\delta\nu_{\ell=1}\rangle$ do not show any excess in cycle 24 near the maximum.

Although we note differences between the instruments, we favour GOLF results as the duty cycle of the subseries are in most cases above 95% while in GONG these vary between 78% to 92% and in BiSON these are often below 65%.
<table><caption>**Table B.1.** GOLF $\langle\delta\nu\rangle$ for modes $\ell = 0, 1$, and $2$ for the three frequency bands, LFB, MFB, and HFB.</caption><thead><tr><th rowspan="3">Time [a]</th><th colspan="9">GOLF</th></tr><tr><th colspan="3">&lang;δν&rang; LFB [μHz]</th><th colspan="3">&lang;δν&rang; MFB [μHz]</th><th colspan="3">&lang;δν&rang; HFB [μHz]</th></tr><tr><th>ℓ=0</th><th>ℓ=1</th><th>ℓ=2</th><th>ℓ=0</th><th>ℓ=1</th><th>ℓ=2</th><th>ℓ=0</th><th>ℓ=1</th><th>ℓ=2</th></tr></thead><tbody><tr><td>1996.770</td><td>0.014 ± 0.021</td><td>-0.017 ± 0.020</td><td>0.012 ± 0.030</td><td>...</td><td>...</td><td>...</td><td>-0.065 ± 0.037</td><td>0.001 ± 0.042</td><td>-0.027 ± 0.050</td></tr><tr><td>1997.019</td><td>-0.003 ± 0.020</td><td>-0.002 ± 0.021</td><td>0.014 ± 0.025</td><td>...</td><td>...</td><td>...</td><td>0.004 ± 0.038</td><td>-0.050 ± 0.041</td><td>-0.033 ± 0.049</td></tr><tr><td>1997.268</td><td>-0.020 ± 0.022</td><td>-0.017 ± 0.021</td><td>0.030 ± 0.027</td><td>...</td><td>...</td><td>...</td><td>0.020 ± 0.038</td><td>-0.029 ± 0.039</td><td>0.009 ± 0.051</td></tr><tr><td>1997.518</td><td>-0.006 ± 0.024</td><td>-0.001 ± 0.022</td><td>0.001 ± 0.027</td><td>...</td><td>...</td><td>...</td><td>0.062 ± 0.038</td><td>-0.018 ± 0.040</td><td>0.016 ± 0.052</td></tr><tr><td>1997.770</td><td>-0.018 ± 0.024</td><td>0.030 ± 0.026</td><td>-0.010 ± 0.028</td><td>...</td><td>...</td><td>...</td><td>0.073 ± 0.039</td><td>-0.001 ± 0.040</td><td>0.031 ± 0.054</td></tr><tr><td>1998.019</td><td>-0.007 ± 0.025</td><td>-0.002 ± 0.023</td><td>-0.039 ± 0.028</td><td>...</td><td>...</td><td>...</td><td>0.054 ± 0.041</td><td>0.087 ± 0.043</td><td>0.038 ± 0.052</td></tr><tr><td>1998.268</td><td>-0.003 ± 0.023</td><td>-0.020 ± 0.027</td><td>-0.074 ± 0.027</td><td>...</td><td>...</td><td>...</td><td>0.081 ± 0.041</td><td>0.124 ± 0.042</td><td>0.011 ± 0.049</td></tr><tr><td>1998.518</td><td>-0.014 ± 0.026</td><td>0.004 ± 0.026</td><td>-0.070 ± 0.027</td><td>...</td><td>...</td><td>...</td><td>0.125 ± 0.042</td><td>0.178 ± 0.044</td><td>0.018 ± 0.050</td></tr><tr><td>1998.770</td><td>-0.045 ± 0.027</td><td>-0.042 ± 0.023</td><td>0.036 ± 0.023</td><td>...</td><td>...</td><td>...</td><td>0.198 ± 0.042</td><td>0.240 ± 0.045</td><td>0.072 ± 0.048</td></tr><tr><td>1999.019</td><td>-0.027 ± 0.027</td><td>-0.027 ± 0.022</td><td>0.061 ± 0.021</td><td>...</td><td>...</td><td>...</td><td>0.270 ± 0.042</td><td>0.381 ± 0.047</td><td>0.057 ± 0.054</td></tr><tr><td>...</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** Time is given in fractional years and frequency shifts and errors are in μHz. The complete table with all significant digits is available online in a machine-readable format.

<table><caption>**Table B.2.** BiSON $\langle\delta\nu\rangle$ for modes $\ell = 0, 1$, and $2$ for the three frequency bands, LFB, MFB, and HFB.</caption><thead><tr><th rowspan="3">Time [a]</th><th colspan="9">BiSON</th></tr><tr><th colspan="3">&lang;δν&rang; LFB [μHz]</th><th colspan="3">&lang;δν&rang; MFB [μHz]</th><th colspan="3">&lang;δν&rang; HFB [μHz]</th></tr><tr><th>ℓ=0</th><th>ℓ=1</th><th>ℓ=2</th><th>ℓ=0</th><th>ℓ=1</th><th>ℓ=2</th><th>ℓ=0</th><th>ℓ=1</th><th>ℓ=2</th></tr></thead><tbody><tr><td>1996.770</td><td>-0.002 ± 0.035</td><td>-0.026 ± 0.034</td><td>-0.001 ± 0.027</td><td>...</td><td>...</td><td>...</td><td>-0.072 ± 0.048</td><td>0.027 ± 0.051</td><td>-0.060 ± 0.052</td></tr><tr><td>1997.019</td><td>0.008 ± 0.033</td><td>-0.025 ± 0.035</td><td>-0.005 ± 0.025</td><td>...</td><td>...</td><td>...</td><td>-0.017 ± 0.050</td><td>-0.056 ± 0.051</td><td>-0.066 ± 0.050</td></tr><tr><td>1997.268</td><td>-0.008 ± 0.034</td><td>-0.026 ± 0.033</td><td>0.013 ± 0.029</td><td>...</td><td>...</td><td>...</td><td>-0.003 ± 0.052</td><td>-0.046 ± 0.051</td><td>-0.023 ± 0.054</td></tr><tr><td>1997.518</td><td>-0.016 ± 0.036</td><td>-0.001 ± 0.034</td><td>-0.003 ± 0.026</td><td>...</td><td>...</td><td>...</td><td>0.073 ± 0.052</td><td>-0.023 ± 0.050</td><td>0.039 ± 0.056</td></tr><tr><td>1997.770</td><td>0.002 ± 0.038</td><td>0.031 ± 0.037</td><td>0.001 ± 0.028</td><td>...</td><td>...</td><td>...</td><td>0.087 ± 0.053</td><td>-0.028 ± 0.052</td><td>0.070 ± 0.056</td></tr><tr><td>1998.019</td><td>-0.014 ± 0.041</td><td>0.010 ± 0.036</td><td>-0.022 ± 0.027</td><td>...</td><td>...</td><td>...</td><td>0.112 ± 0.056</td><td>0.111 ± 0.051</td><td>0.066 ± 0.056</td></tr><tr><td>1998.268</td><td>0.025 ± 0.041</td><td>0.008 ± 0.038</td><td>-0.024 ± 0.028</td><td>...</td><td>...</td><td>...</td><td>0.189 ± 0.055</td><td>0.189 ± 0.050</td><td>0.118 ± 0.054</td></tr><tr><td>1998.518</td><td>0.039 ± 0.041</td><td>0.034 ± 0.037</td><td>-0.047 ± 0.031</td><td>...</td><td>...</td><td>...</td><td>0.255 ± 0.055</td><td>0.251 ± 0.051</td><td>0.135 ± 0.053</td></tr><tr><td>1998.770</td><td>-0.014 ± 0.041</td><td>0.036 ± 0.036</td><td>-0.015 ± 0.030</td><td>...</td><td>...</td><td>...</td><td>0.294 ± 0.053</td><td>0.309 ± 0.051</td><td>0.199 ± 0.052</td></tr><tr><td>1999.019</td><td>-0.022 ± 0.040</td><td>0.044 ± 0.035</td><td>0.034 ± 0.028</td><td>...</td><td>...</td><td>...</td><td>0.301 ± 0.051</td><td>0.324 ± 0.051</td><td>0.232 ± 0.054</td></tr><tr><td>...</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** Time is given in fractional years and frequency shifts and errors are in μHz. The complete table with all significant digits is available online in a machine-readable format.

<table><caption>**Table B.3.** GONG $\langle\delta\nu\rangle$ for modes $\ell = 0, 1$, and $2$ for the three frequency bands, LFB, MFB, and HFB.</caption><thead><tr><th rowspan="3">Time [a]</th><th colspan="9">GONG</th></tr><tr><th colspan="3">&lang;δν&rang; LFB [μHz]</th><th colspan="3">&lang;δν&rang; MFB [μHz]</th><th colspan="3">&lang;δν&rang; HFB [μHz]</th></tr><tr><th>ℓ=0</th><th>ℓ=1</th><th>ℓ=2</th><th>ℓ=0</th><th>ℓ=1</th><th>ℓ=2</th><th>ℓ=0</th><th>ℓ=1</th><th>ℓ=2</th></tr></thead><tbody><tr><td>1997.022</td><td>-0.007 ± 0.028</td><td>0.011 ± 0.068</td><td>0.030 ± 0.053</td><td>...</td><td>...</td><td>...</td><td>0.026 ± 0.040</td><td>-0.042 ± 0.066</td><td>-0.021 ± 0.075</td></tr><tr><td>1997.219</td><td>-0.030 ± 0.028</td><td>0.006 ± 0.066</td><td>0.026 ± 0.055</td><td>...</td><td>...</td><td>...</td><td>0.037 ± 0.041</td><td>-0.075 ± 0.070</td><td>0.004 ± 0.075</td></tr><tr><td>1997.416</td><td>-0.012 ± 0.028</td><td>0.060 ± 0.061</td><td>0.039 ± 0.051</td><td>...</td><td>...</td><td>...</td><td>0.064 ± 0.042</td><td>-0.070 ± 0.071</td><td>0.021 ± 0.075</td></tr><tr><td>1997.614</td><td>-0.008 ± 0.029</td><td>0.058 ± 0.062</td><td>0.046 ± 0.050</td><td>...</td><td>...</td><td>...</td><td>0.082 ± 0.042</td><td>-0.054 ± 0.074</td><td>0.039 ± 0.078</td></tr><tr><td>1997.811</td><td>0.013 ± 0.029</td><td>0.006 ± 0.063</td><td>0.038 ± 0.052</td><td>...</td><td>...</td><td>...</td><td>0.057 ± 0.043</td><td>-0.077 ± 0.075</td><td>0.088 ± 0.077</td></tr><tr><td>1998.008</td><td>0.039 ± 0.027</td><td>-0.022 ± 0.068</td><td>0.013 ± 0.054</td><td>...</td><td>...</td><td>...</td><td>0.044 ± 0.045</td><td>-0.019 ± 0.076</td><td>0.103 ± 0.078</td></tr><tr><td>1998.205</td><td>0.048 ± 0.030</td><td>-0.186 ± 0.055</td><td>-0.009 ± 0.054</td><td>...</td><td>...</td><td>...</td><td>0.104 ± 0.046</td><td>-0.004 ± 0.079</td><td>0.165 ± 0.079</td></tr><tr><td>1998.403</td><td>0.054 ± 0.028</td><td>-0.006 ± 0.069</td><td>0.004 ± 0.058</td><td>...</td><td>...</td><td>...</td><td>0.149 ± 0.045</td><td>0.076 ± 0.081</td><td>0.218 ± 0.079</td></tr><tr><td>1998.600</td><td>0.030 ± 0.031</td><td>0.042 ± 0.062</td><td>-0.002 ± 0.058</td><td>...</td><td>...</td><td>...</td><td>0.161 ± 0.046</td><td>0.101 ± 0.078</td><td>0.251 ± 0.078</td></tr><tr><td>1998.797</td><td>0.018 ± 0.028</td><td>0.066 ± 0.056</td><td>-0.002 ± 0.053</td><td>...</td><td>...</td><td>...</td><td>0.203 ± 0.046</td><td>0.179 ± 0.078</td><td>0.243 ± 0.077</td></tr><tr><td>...</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** Time is given in fractional years and frequency shifts and errors are in μHz. The complete table with all significant digits is available online in a machine-readable format.
<figure><img src="image_4.png" /><figcaption>**Fig. C.1.** Temporal evolution of the averaged but not smoothed GOLF frequency shifts (dark green lines). Errors are shown as light green regions. As a reference, the dark blue line represents the smoothed frequency shifts as shown in Fig. 2. Orange and vertical grey lines have the same meaning as in Fig. 2.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. D.1.** Temporal evolution of averaged but not smoothed BiSON frequency shifts (dark green lines) with errors represented by the light green regions. Same legend as in Fig. C.1.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. D.2.** Temporal evolution of averaged but not smoothed GONG frequency shifts (dark green lines) with errors represented by the light green regions. The signal-to-noise ratio of *l*=2 modes in the LFB is too low and the results are not reliable. Same legend as in Fig. C.1.</figcaption></figure>