Issue 
A&A
Volume 640, August 2020



Article Number  L10  
Number of page(s)  8  
Section  Letters to the Editor  
DOI  https://doi.org/10.1051/00046361/202038539  
Published online  04 August 2020 
Letter to the Editor
Predicting frequency changes of globalscale solar Rossby modes due to solar cycle changes in internal rotation
^{1}
MaxPlanckInstitut für Sonnensystemforschung, JustusvonLiebigWeg 3, 37077 Göttingen, Germany
email: crgoddard10@gmail.com, gizon@mps.mpg.de
^{2}
Institut für Astrophysik, GeorgAugustUniversität Göttingen, FriedrichHundPlatz 1, 37077 Göttingen, Germany
^{3}
Center for Space Science, NYUAD Institute, New York University Abu Dhabi, Abu Dhabi, UAE
Received:
30
May
2020
Accepted:
13
July
2020
Context. Largescale equatorial Rossby modes have been observed on the Sun over the last two solar cycles.
Aims. We investigate the impact of the timevarying zonal flows on the frequencies of Rossby modes.
Methods. A firstorder perturbation theory approach is used to obtain an expression for the expected shift in the mode frequencies due to perturbations in the internal rotation rate.
Results. Using the timevarying rotation from helioseismic inversions we predict the changes in Rossby mode frequencies with azimuthal orders from m = 1 to m = 15 over the last two solar cycles. The peaktopeak frequency change is less than 1 nHz for the m = 1 mode, grows with m, and reaches 25 nHz for m = 15.
Conclusions. Given the observational uncertainties on mode frequencies due to the finite mode lifetimes, we find that the predicted frequency shifts are near the limit of detectability.
Key words: Sun: rotation / waves / Sun: oscillations / Sun: interior / Sun: activity / hydrodynamics
© C. R. Goddard et al. 2020
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Open Access funding provided by Max Planck Society.
1. Introduction
Rossby waves are globalscale waves of vorticity which arise due to the conservation of the radial component of the absolute vorticity for rotating fluids. They are most commonly discussed in the context of planetary atmospheres and oceans (e.g. Vallis 2017). Detections of their stellar analogues, also known as r modes, have been proposed for Gamma Doradus stars (Van Reeth et al. 2016; Saio et al. 2018; Saio 2018). Recently Rossby waves have been unambiguously detected in the nearsurface flows of the Sun (Löptien et al. 2018); they are characterised by a robust dispersion relation which is close to that of classical Rossby waves in a uniformly rotating fluid (e.g. Saio 1982). Followup studies with helioseismic methods and data have confirmed the detection of modes with azimuthal orders m = 3–15 (Liang et al. 2019; Hanasoge & Mandal 2019; Hanson et al. 2020). The measured latitudinal eigenfunctions deviate from sectoral Legendre polynomials (Proxauf et al. 2020). Calculations in the equatorial βplane suggest that the effect of latitudinal differential rotation can account for the real part of the observed latitudinal eigenfunctions (Gizon et al. 2020). The observed Rossby waves can be described as equatorial modes trapped by latitudinal differential rotation.
The theory of r modes in the solar and stellar context has developed over the last few decades for various limiting cases (e.g. Provost et al. 1981; Wolff 1998; Zaqarashvili et al. 2007; Saio et al. 2018). Damiani et al. (2020) recently computed eigenmodes in slowly rotating polytropes and confirmed the r^{m} radial dependence of the displacement eigenfunctions proposed by Provost et al. (1981).
Globalscale Rossby waves are probes of the convection zone. In this paper we predict the effect of solar cycle variations in the zonal flows on their frequencies, since these frequencies are directly related to the solar rotation rate. This work is required to interpret future measurements of time variations, in particular to disentangle perturbations caused by solar cycle changes in rotation from perturbations caused by the evolving largescale magnetic field.
Helioseismology with data from the Michelson Doppler Imager on the Solar and Heliospheric Observatory (SOHO/MDI; Scherrer et al. 1995) and the Global Oscillation Network Group (GONG; Hill et al. 1996) allowed the internal rotation profile of the Sun to be measured in the convection zone. The latitudinal differential rotation seen at the surface extends throughout the convection zone (e.g. Schou et al. 1998). Additionally, data from the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory (SDO/HMI; Schou et al. 2012) confirmed these findings. All data sets revealed changes in the rotation rate that are associated with the solar cycle. These timevarying zonal flows or torsional oscillations were first noted at the surface (Howard & Labonte 1980) and have since been shown to extend throughout most of the convection zone (Vorontsov et al. 2002), although they weaken with depth (e.g. Howe 2009). The timevarying zonal flows consist of bands of slower and faster rotation that move towards the equator at low latitudes and towards the poles at high latitudes. These zonal flows vary significantly between the last two solar cycles, solar cycle 23 (SC23) and solar cycle 24 (SC24; Komm et al. 2014; Howe et al. 2018). Additionally, there is a difference in the average rotation between these two cycles. In SC24, the convection zone rotates more quickly at low latitudes, but more slowly at high latitudes (e.g. Basu & Antia 2019).
Section 2 describes the internal rotation data. In Sect. 3, we use a firstorder perturbation approach (Born approximation) to obtain a kernel for the sensitivity of Rossby waves to variations of an arbitrary rotation profile, and the predictions via a series of steady problems are presented in Sect. 4. Finally, the discussion and conclusions are given in Sect. 5.
2. Solar cycle variations of internal rotation
The rotation rate data used to prescribe the zonal flows are the 2D regularised leastsquares (RLS) global helioseismology inversions for solar rotation from MDI and HMI (Larson & Schou 2018). The data used here cover the period from 1996 June 06 to 2019 June 30, with a cadence of 72 days. The HMI data are used during the period of overlap between MDI and HMI. We work in sphericalpolar coordinates, where r is the radial coordinate, ϕ is longitude, and θ is the colatitude. A reference rotation profile, , is taken to be the unweighted mean rotation profile over the 23year time interval. Variations about this mean profile over time are thus
The few missing time steps are dealt with via linear interpolation. The rotation residual is then binned by taking the unweighted mean for each calendar year.
Figure 1 shows the binned rotation rate residual, δΩ(r, θ, t), at r/R_{⊙} = 0.90, 0.95, and 0.99 for the specified time window. The zonal flows and high latitude variations associated with the 11year solar cycle are seen alongside the asymmetry between SC23 and SC24. As mentioned previously, in the convection zone SC24 has higher rotation rates at low latitudes, but lower rotation rates at high latitudes. This asymmetry is not due to the inconsistencies between the MDI and HMI data, and is clearly seen in GONG data (Basu & Antia 2019).
Fig. 1. Rotation rate residuals, δΩ(r, θ, t)/2π, at r = 0.9, 0.95, and 0.99 R_{⊙} for 1996Jun.06 to 2019Jun.30. The data is binned for each calendar year, and the colour scale is saturated at ±3 nHz. The vertical dashed grey lines denote the positions of the maxima of solar cycles 23 and 24 (2002 and 2014). The vertical dashed brown line indicates the minimum between the two cycles (2009). 
3. Perturbation to Rossby mode frequencies due to variations in rotation
We utilise firstorder perturbation theory to obtain a kernel for the sensitivity of Rossby wave frequencies to small variations in the rotation profile. This method is commonly used to obtain sensitivity kernels for the effect of rotation on p modes (e.g. Gough & Thompson 1990; Howe 2009; Basu 2016).
3.1. Oscillation equation
We start from the linearised equation of motion for the displacement vector ξ, written in an inertial frame (LyndenBell & Ostriker 1967; Unno et al. 1989). Assuming a steady rotation rate Ω(r, θ) (the wave period is shorter than the timescale of the evolution of rotation) and a displacement vector proportional to exp(−iωt), we have
with
where
is the rotational velocity, r is the radial position vector, p is pressure, ψ is the gravitational potential, and d is the Lagrangian change operator (e.g. dp = p′+ξ ⋅ ∇p_{0}). Frictional and electromagnetic forces are not considered.
The following relation holds if we assume the azimuthal dependence of the displacement to be proportional to exp(imϕ), where m is the azimuthal order of a solar Rossby mode:
Using this relationship we can rewrite the operators that depend on u as
Substituting these into Eq. (2), and defining ℒ(ξ) = V(ξ)+P(ξ), gives
3.2. Firstorder oscillation equation
We now perturb Eq. (7) around an arbitrary reference rotation rate Ω_{0}(r, θ) using
The zeroorder eigenvalue problem for (ω_{0}, ξ_{0}) is
The firstorder eigenvalue problem for (δω, δξ) is
3.3. Firstorder frequency shift
We define an inner product of two functions as ⟨η, ξ⟩ = ∫_{V}η^{*} ⋅ ξρ dV, where V is the volume of the star, the asterisk denotes the complex conjugate, and ρ(r) is the radial density profile.
Dotting Eq. (10) with on the left, integrating over V, and utilising the selfadjointness of the operators (cf. LyndenBell & Ostriker 1967) gives
where is the unit vector along the rotation axis, . Using this expression the frequency shift due to variations in an arbitrary rotation profile may be calculated from the unperturbed (zeroorder) displacement eigenfunctions. Rearranging the terms that include δω, and defining 𝒮(ξ) to contain the centrifugal and δℒ terms,
we have
with
We note that the denominator is more complicated than in the standard pmode case with a nonrotating reference model. For the rest of the paper, we neglect the term 𝒮(ξ) in Eq. (13), which corresponds to neglecting the perturbation of the centrifugal acceleration and the perturbation to the pressure and gravity terms due to the change in rotation δΩ.
3.4. Rossby wave eigenfunction and resulting kernel
From here forwards we consider reference eigenfunctions ξ_{0} which correspond to incompressible global Rossby waves for the case of slow uniform rotation. In this case deviations of the solar structure (i.e. oblateness) from variations in rotation have a negligible impact on the eigenfrequencies (Damiani et al. 2020). The relevant zeroorder equation of motion, Eq. (9) without the centrifugal (third) term, is now
This is the inertial frame equivalent of the oscillation equation given in Provost et al. (1981) and Damiani et al. (2020), under the Cowling approximation.
For the case of slow uniform rotation, the eigenfunction of the Rossby mode with azimuthal wavenumber m is of the form
where ∇_{h} is the horizontal gradient, R(r) and L(θ) are functions to be specified for each mode, and L′(θ) = dL/dθ.
The frequency shift to the unperturbed eigenfrequency ω_{0} of mode m due to a change in rotation is given by
where the sensitivity kernels K are given by
For the displacement eigenfunction given in Eq. (17) and a realvalued L(θ), we have
and
Now we choose specific radial and latitudinal dependences for the eigenfunction given by Eq. (17), along with a corresponding reference eigenfrequency. Considering the case of uniform slow rotation
we take a sectoral spherical harmonic (ℓ = m) solution with no radial node:
This simple form for the radial eigenfunctions is suggested by the study of r modes in slowly (and uniformly) rotating polytropes in Damiani et al. (2020). Corrections may be needed to account for strong radial gradients in solar rotation and for the subadiabatic stratification below the convection zone; however, this is beyond the scope of this paper. We note that the kinetic energy density of most of the modes we consider here (m ≥ 4) is largely restricted to the convection zone.
For this case, the kernel simplifies to
This kernel does not depend on the uniform rotation rate Ω_{0} from the zeroorder model. Appendix A shows that this kernel predicts the correct answer for the case of a perturbation δΩ that is constant.
Figure 2 shows the kernel calculated with the smoothed solarlike density profile from Gizon et al. (2017). Figure 3 shows normalised latitudinal and radial cuts for azimuthal orders m = 3, 9, and 15. The frequency shifts are mostly sensitive to perturbations of the rotation rate in the low latitude convection zone for m ≥ 3. The negative higher latitude component of the kernel reduces in amplitude with increasing azimuthal order. The detected solar Rossby waves currently span m = 3–15, and hence the plots are restricted to m ≤ 15. Appendix B discusses the latitudinal dependence of the kernel in more detail.
Fig. 2. Kernels for the sensitivity of sectoral Rossby waves to variations in rotation rate for azimuthal orders m = 1, 3, 6, 9, 12, and 15 (from left to right, top to bottom). The different panels show the kernels in the meridional plane and the kernels are symmetric about the equator. The vertical axis is the rotation axis. The spatial grid is in units of R_{⊙}, and the grided region denotes the convection zone. 
Fig. 3. Radial (top) and latitudinal (bottom) dependences of the kernel for sectoral Rossby modes with azimuthal orders m = 3, 9, and 15, normalised by their maxima. The vertical dashed line in the top panel denotes the position of the base of the convection zone at r/R_{⊙} = 0.7. 
4. Results
4.1. Predicted frequency shifts over the last two solar cycles
Here we compute the Rossby wave frequency shifts δω for each mode m, which would be expected from the meansubtracted rotation residuals, δΩ(r, θ, t), defined in Eq. (1). The rotation residuals are binned by calendar year, as shown in Fig. 1. We do not consider the timeindependent frequency shift for each m caused by the difference between our zeroorder model (uniform rotation, Ω_{0}) and the timeaveraged differential rotation of the Sun () used to define δΩ(r, θ, t).
The frequencyshift time series for each azimuthal order were calculated from the binned rotation residuals and Eqs. (18) and (26) using trapezoidal rule integration. The perturbation to the rotation profile δΩ(r, θ, t) is set to zero in the radiative zone (r/R_{⊙} > 0.7), due to the inaccuracies in the inversions there and the lack of zonal flow signatures at these depths.
Figure 4 shows the obtained frequency shifts separately for low and high azimuthal orders, 1 ≤ m ≤ 7 and 8 ≤ m ≤ 15, respectively. For m = 1 no discernible trend is seen as this mode has little sensitivity to the strongest zonal flows in the upper convection zone. From m = 2 to m = 7 an apparent cycle with a period of > 22 years is observed: the frequency shift is negative for SC23 and positive for SC24. For m > 7 a doublepeaked trend is seen, similar to the 11year Sunspot cycle, but with the frequency during SC24 systematically higher than in SC23. This is a superposition of the 11year cycle of the zonal flows with the asymmetry of the two cycles (as seen for low m). Appendix C discusses the relationship between the frequency shift and the rotation residuals in more detail. Appendix D shows a comparison to frequency shifts computed from direct numerical calculations in the βplane (Gizon et al. 2020).
Fig. 4. Rossby mode frequencyshift time series due to variation in the solar rotation profile, plotted for m from 1 to 8 (top), and 9 to 15 (middle), calculated for each calendar year. The vertical grey dashed lines denote the positions of the maxima of solar cycles 23 and 24, to the nearest year. The vertical brown dashed line indicates the minimum between the two cycles. The smoothed sunspot number is plotted in the bottom panel. 
As one simple way to quantify these trends, we calculate the change in frequencies between two times for each m. The upper panel of Fig. 5 shows the difference in the frequencies between 2014 and 2002 (two solar maxima) and between 2014 and 2009 (maximum and minimum). The difference between maximum and minimum highlights the 11year component, and continues to increase for high m, reaching 23 nHz for m = 15. The difference between the cycle maxima highlights the asymmetry between the two cycles, which saturates for m > 8 at about 7 nHz.
Fig. 5. Top: predicted difference in the Rossby mode frequencies between the two solar maxima (grey curve, 2014 minus 2002) and between the maximum of cycle 24 and the preceding solar minimum (brown curve, 2014 minus 2009). The error bars are propagated errors from the scatter in the solar rotation measurements. Bottom: estimated errors in the Rossby mode frequencies due to the finite mode lifetimes as a function of observation time T. These errors can be used to determine whether the predicted frequency shifts from the top panel are detectable. 
4.2. Predicted measurement uncertainties due to measured finite mode lifetimes
Let us denote by σ the uncertainty on the measured frequency ω of a Rossby mode m due to the finite mode lifetimes. Under the assumption that Rossby modes are stochastically excited, the uncertainties σ scale as T^{−1/2}, where T is the observation time. Specifically, we have σ = (Γ/2T)^{1/2}, where Γ is the mode linewidth in frequency space (Libbrecht 1992).
We used the measurements of the mode linewidths from Liang et al. (2019) to calculate the measurement uncertainties σ for various values T. The estimate of the uncertainties in the r mode frequencies depends on the measured lifetimes. We note that the various methods of helioseismology and granulation tracking give essentially the same answers for the mode linewidths (Hanson et al. 2020). The bottom panel of Fig. 5 shows this for several values of m. We note the strong dependence on m. For some azimuthal orders and T = 3 years, for example, the uncertainties are smaller than the predicted solar cycle frequency shifts. Thus, the detection of the Rossby mode frequency shifts due to the timevarying zonal flows may be possible. However, as this is near the limit of detectability some averaging over several values of m may be required to reduce noise.
5. Conclusion
We used firstorder perturbation theory to obtain a kernel for the frequency shift of solar Rossby waves caused by perturbations in the rotation profile. We used these kernels together with solar torsional oscillation data to predict the time varying component of the frequency shift for solar Rossby waves. The predicted frequency shifts over the last two solar cycles have a strong dependence on the azimuthal order m. Modes with m ≤ 7 are mainly sensitive to the asymmetry between the two cycles, while those with m > 7 are mainly sensitive to the 11year cycle of the zonal flows. The predicted frequency difference between cycle maximum and minimum reaches 23 nHz for m = 15.
Several approximations and simplifications were made in the construction of the kernel, and in the chosen latitudinal and radial dependence for the sectoral Rossby waves. These are adequate for the problem at hand, but may require further development once the frequency shifts have been observed. For example, it might be useful to consider eigenfunctions from a differentially rotating reference model. In future work we intend to investigate potential solar cycle changes in the amplitudes and the lifetimes of the modes, as well as potential changes in the eigenfunctions due to the torsional oscillations or the magnetic field.
The predicted shift of the mode frequencies due to variation in the internal rotation over the last two solar cycles is near the limit of detection. A detection of this effect would help further constrain the physics of solar Rossby modes. Should there be much larger solar cycle variations in the mode frequencies, they will have to be attributed to the effect of the solar magnetic field in the convection zone.
Acknowledgments
We thank ZhiChao Liang and Chris S. Hanson for useful discussions. Author contributions: LG proposed the research. CRG performed the kernel calculations and the data analysis. ACB and LG contributed to the kernel calculations. DF performed the βplane calculations in Appendix D. CRG wrote the draft paper and all authors contributed to the final manuscript. This research is supported by ERC Synergy grant WHOLE SUN #810218. The data were processed at the German Data Center for SDO, funded by the German Aerospace Center (DLR). LG acknowledges partial support from the NYU Abu Dhabi Center for Space Science under grant G1502. The HMI and MDI helioseismology data products are used courtesy of the JSOC team. NumPy (van der Walt et al. 2011) and matplotlib (Hunter 2007) were used in this research.
References
 Basu, S. 2016, Liv. Rev. Sol. Phys., 13, 2 [Google Scholar]
 Basu, S., & Antia, H. M. 2019, ApJ, 883, 93 [Google Scholar]
 Damiani, C., Cameron, R. H., Birch, A. C., & Gizon, L. 2020, A&A, 637, A65 [EDP Sciences] [Google Scholar]
 Gizon, L., Barucq, H., Duruflé, M., et al. 2017, A&A, 600, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gizon, L., Fournier, D., & Albekioni, M. 2020, A&A, submitted [Google Scholar]
 Gough, D. O., & Thompson, M. J. 1990, MNRAS, 242, 25 [Google Scholar]
 Hanasoge, S., & Mandal, K. 2019, ApJ, 871, L32 [Google Scholar]
 Hanson, C. S., Gizon, L., & Liang, Z.C. 2020, A&A, 635, A109 [EDP Sciences] [Google Scholar]
 Hill, F., Stark, P. B., Stebbins, R. T., et al. 1996, Science, 272, 1292 [Google Scholar]
 Howard, R., & Labonte, B. J. 1980, ApJ, 239, L33 [Google Scholar]
 Howe, R. 2009, Liv. Rev. Sol. Phys., 6, 1 [Google Scholar]
 Howe, R., Hill, F., Komm, R., et al. 2018, ApJ, 862, L5 [Google Scholar]
 Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [Google Scholar]
 Komm, R., Howe, R., González Hernández, I., & Hill, F. 2014, Sol. Phys., 289, 3435 [Google Scholar]
 Larson, T. P., & Schou, J. 2018, Sol. Phys., 293, 29 [Google Scholar]
 Liang, Z.C., Gizon, L., Birch, A. C., & Duvall, T. L. 2019, A&A, 626, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Libbrecht, K. G. 1992, ApJ, 387, 712 [Google Scholar]
 Löptien, B., Gizon, L., Birch, A. C., et al. 2018, Nat. Astron., 2, 568 [Google Scholar]
 LyndenBell, D., & Ostriker, J. P. 1967, MNRAS, 136, 293 [Google Scholar]
 Provost, J., Berthomieu, G., & Rocca, A. 1981, A&A, 94, 126 [NASA ADS] [Google Scholar]
 Proxauf, B., Gizon, L., Löptien, B., et al. 2020, A&A, 634, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Saio, H. 1982, ApJ, 256, 717 [Google Scholar]
 Saio, H. 2018, ArXiv eprints [arXiv:1812.01253] [Google Scholar]
 Saio, H., Kurtz, D. W., Murphy, S. J., Antoci, V. L., & Lee, U. 2018, MNRAS, 474, 2774 [Google Scholar]
 Scherrer, P. H., Bogart, R. S., Bush, R. I., et al. 1995, Sol. Phys., 162, 129 [Google Scholar]
 Schou, J., Antia, H. M., Basu, S., et al. 1998, ApJ, 505, 390 [Google Scholar]
 Schou, J., Scherrer, P. H., Bush, R. I., et al. 2012, Sol. Phys., 275, 229 [Google Scholar]
 Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial Oscillations of Stars (Tokyo: University of Tokyo Press) [Google Scholar]
 Vallis, G. K. 2017, Atmospheric and Oceanic Fluid Dynamics: Fundamentals and LargeScale Circulation, 2nd edn. (Cambridge: Cambridge University Press) [Google Scholar]
 van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Comput. Sci. Eng., 13, 22 [Google Scholar]
 Van Reeth, T., Tkachenko, A., & Aerts, C. 2016, A&A, 593, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vorontsov, S. V., ChristensenDalsgaard, J., Schou, J., Strakhov, V. N., & Thompson, M. J. 2002, Science, 296, 101 [Google Scholar]
 Wolff, C. L. 1998, ApJ, 502, 961 [Google Scholar]
 Zaqarashvili, T. V., Oliver, R., Ballester, J. L., & Shergelashvili, B. M. 2007, A&A, 470, 815 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Case δΩ = const.
The frequency shift in any reference frame for the described sectoral Rossby modes caused by a change in the rotation rate is
The kernel given should reduce to this result in the case of uniform rotation (Ω_{0} = const) with a uniform perturbation (δΩ = const). This has been shown via numerical integration, once converged the resulting frequency shifts agree with Eq. (A.1) to machine precision.
The kernel may also be integrated directly to obtain the frequency shift for the case of uniform rotation and perturbation. Defining the frequency shift as
and taking the kernel given in Eq. (26), the radial component cancels from top and bottom giving
Utilising the identity , we obtain
which is the expected result.
Appendix B: Latitudinal form of the kernel
The figures in Sect. 3.4 show that the sign of the kernel switches at a given latitude, which decreases as the azimuthal order increases. This behaviour is due to the terms from the advection and coriolis force in the final kernel having opposite signs, and different latitudinal dependences. The remaining advective term in the kernel is proportional to m(sin θ)^{2m − 1} and the remaining coriolis term is proportional to −(m^{2} + 2)(cos θ)^{2}(sin θ)^{2m − 1}; therefore, they are equal and opposite at
Figure B.1 shows the latitudinal dependence of the coriolis and advective terms of the kernel, and their sum, for m = 3. The vertical dashed line shows the position as which the two terms are equal and opposite, summing to zero (θ_{s}). As m increases this position shifts towards the equator and the relative contribution of the coriolis term decreases.
Fig. B.1. Latitudinal dependences of the coriolis and advective terms that contribute to the kernel, and their sum, for m = 3. 
Appendix C: Explaining the behaviour of δω for high and low m values
The top three panels of Fig. C.1 show the result of multiplication of the rotation rate residual and the kernel, δΩ(r, θ) K(r, θ), with integration over the radial coordinate. These plots aid in understanding the connection between the rotation data and the calculated frequency shifts for the Rossby waves. We note that faster rotation than the mean profile will result in a positive frequency shift (except where the sensitivity kernel is negative).
Fig. C.1. Top three panels: rotation rate residual multiplied by the sensitivity kernel and integrated over the radial coordinate. Bottom panel: difference in the solar rotation rate in the convection zone between the maxima of cycles 24 and 23. 
For m = 3, we see that both high and low latitudes have a negative contribution to the frequency shift for SC23 and positive contribution for SC24. This is due to SC23 having lower rotation rates at low latitudes and higher rotation rates at high latitudes compared to SC24, in combination with the kernel being negative for higher latitudes (see Fig. 2). This is shown in the lower panel of Fig. C.1 where the difference between the rotation residual at the maxima of SC23 and SC24 is plotted. For m = 9 we see a combination of this, and a strong contribution at low latitudes from the zonal flows, which becomes even stronger for m = 15 as the kernel becomes more constrained to low latitudes near the surface (see Fig. 2). In all three plots the switch of the kernel from positive to negative is visible as the pair of horizontal lines of zero frequency at latitudes π/2 − θ_{s}.
Appendix D: Comparison with calculations in the βplane
Here we perform a comparison with recent direct numerical calculations of frequencies of viscous Rossby modes in the equatorial βplane (Gizon et al. 2020) to validate our firstorder frequency shifts. These calculations cannot inform us about the effects of the radial variations in rotation, but they can be performed for any latitudinal variations in the zonal flows. In addition, the βplane approximation is only valid for m ≫ 1.
To make a direct comparison, we define the zonal velocity in the βplane as a function of the coordinate y = R_{⊙}cosθ as
where Ω_{eq}(r) is the equatorial rotation rate as a function of radius. This flow is used as an input to the computations described in Gizon et al. (2020). In addition, we specify the viscosity through the Reynolds number Re = 300.
The βplane calculation is performed at the full 72day cadence of the rotation data, and the frequencyshift time series for each m is then mean subtracted to give a result equivalent to our firstorder calculations. These are then binned into calendar years.
Figure D.1 shows the resulting frequencyshift time series along with our firstorder results for m = 9, 12, and 15. There is good agreement for m = 12 and 15. The agreement is still reasonable for m = 9. For even smaller values of m the agreement becomes worse as the problem is less suited to the βplane approximation.
Fig. D.1. Frequencyshift time series from calculations in the βplane (orange) plotted along with the results from the kernel (blue) for m = 9, 12, and 15. 
All Figures
Fig. 1. Rotation rate residuals, δΩ(r, θ, t)/2π, at r = 0.9, 0.95, and 0.99 R_{⊙} for 1996Jun.06 to 2019Jun.30. The data is binned for each calendar year, and the colour scale is saturated at ±3 nHz. The vertical dashed grey lines denote the positions of the maxima of solar cycles 23 and 24 (2002 and 2014). The vertical dashed brown line indicates the minimum between the two cycles (2009). 

In the text 
Fig. 2. Kernels for the sensitivity of sectoral Rossby waves to variations in rotation rate for azimuthal orders m = 1, 3, 6, 9, 12, and 15 (from left to right, top to bottom). The different panels show the kernels in the meridional plane and the kernels are symmetric about the equator. The vertical axis is the rotation axis. The spatial grid is in units of R_{⊙}, and the grided region denotes the convection zone. 

In the text 
Fig. 3. Radial (top) and latitudinal (bottom) dependences of the kernel for sectoral Rossby modes with azimuthal orders m = 3, 9, and 15, normalised by their maxima. The vertical dashed line in the top panel denotes the position of the base of the convection zone at r/R_{⊙} = 0.7. 

In the text 
Fig. 4. Rossby mode frequencyshift time series due to variation in the solar rotation profile, plotted for m from 1 to 8 (top), and 9 to 15 (middle), calculated for each calendar year. The vertical grey dashed lines denote the positions of the maxima of solar cycles 23 and 24, to the nearest year. The vertical brown dashed line indicates the minimum between the two cycles. The smoothed sunspot number is plotted in the bottom panel. 

In the text 
Fig. 5. Top: predicted difference in the Rossby mode frequencies between the two solar maxima (grey curve, 2014 minus 2002) and between the maximum of cycle 24 and the preceding solar minimum (brown curve, 2014 minus 2009). The error bars are propagated errors from the scatter in the solar rotation measurements. Bottom: estimated errors in the Rossby mode frequencies due to the finite mode lifetimes as a function of observation time T. These errors can be used to determine whether the predicted frequency shifts from the top panel are detectable. 

In the text 
Fig. B.1. Latitudinal dependences of the coriolis and advective terms that contribute to the kernel, and their sum, for m = 3. 

In the text 
Fig. C.1. Top three panels: rotation rate residual multiplied by the sensitivity kernel and integrated over the radial coordinate. Bottom panel: difference in the solar rotation rate in the convection zone between the maxima of cycles 24 and 23. 

In the text 
Fig. D.1. Frequencyshift time series from calculations in the βplane (orange) plotted along with the results from the kernel (blue) for m = 9, 12, and 15. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.