Misleading variations in estimated rotational frequency splittings of solar p modes: Consequences for helio- and asteroseismology
The aim of this paper is to investigate whether there are any 11-yr or quasi-biennial solar cycle-related variations in solar rotational splitting frequencies of low-degree solar p modes. Although no 11-yr signals were observed, variations on a shorter timescale were apparent. We show that the variations arose from complications/artifacts associated with the realization noise in the data and the process by which the data were analyzed. More specifically, the realization noise was observed to have a larger effect on the rotational splittings than accounted for by the formal uncertainties. When used to infer the rotation profile of the Sun these variations are not important. The outer regions of the solar interior can be constrained using higher-degree modes. While the variations in the low- splittings do make large differences to the inferred rotation rate of the core, the core rotation rate is so poorly constrained, even by low- modes, that the different inferred rotation profiles still agree within their respective uncertainties. By contrast, in asteroseismology, only low- modes are visible and so higher- modes cannot be used to constrain the rotation profile of stars. Furthermore, we usually only have one data set from which to measure the observed low- splitting. In such circumstances the inferred internal rotation rate of a main sequence star could differ significantly from estimates of the surface rotation rate, hence leading to spurious conclusions. Therefore, extreme care must be taken when using only the splittings of low- modes to draw conclusions about the average internal rotation rate of a star.
The frequencies of the Sun’s acoustic (p-mode) oscillations vary throughout the solar cycle with the frequencies of the most prominent modes being at their highest when solar activity is at its maximum [?]. Mid-term signals can also be observed in p-mode frequencies [?]. Solar cycle variations are also observed in other p-mode parameters such as powers and lifetimes [?]. With the advent of long, continuous, and high-quality asteroseismic observations of solar-like stars [?] by the Convection Rotation and Planetary Transits [?] space mission and Kepler [?] it is possible to study changes in p-mode properties [?] and to measure rotational splittings [?].
Solar rotation splits p modes into azimuthal orders . The difference in frequency between azimuthal orders of a mode can be used to infer the Sun’s internal rotation profile. Variations in the splittings of intermediate- modes have been used to observe the torsional oscillation [?]. No noticeable 11-yr solar cycle changes in the p-mode splittings of low- data have been observed to date [?]. [?] used the rotational splittings of medium- p modes to infer the evolution of the rotation rate at the base of the convection zone through the solar cycle. [?] found a periodicity in the data, however, the signal seemed to disappear after 2000 [?]. Furthermore, the result was not confirmed by other investigations [?]. [?] observed mid-term variations in the rotational splittings of low- modes. Here we examine whether solar cycle-related changes, potentially including mid-term variations, can be observed in the most up-to-date low- p-mode splittings. We also discuss the consequences of observed variations in the rotational splitting on solar rotation profile inversions.
Acoustic p modes are not the only oscillations that propagate through the interior of a star. Gravity (g) modes can propagate only in regions of stable stratification and so are trapped within the central regions of a star, beneath the convection zone. Mixed modes show characteristics of both p and g modes, carrying information on both the outer envelope and the central core. Mixed modes have been detected in stars that have evolved off the main sequence and the splittings of mixed modes have been used to infer the rotation rate of stellar cores [?]. However, here we consider only the inferred average internal rotation rate of main sequence stars, like the Sun, that do not have detectible g modes and mixed modes. We must, therefore, rely on the splittings of low- p modes only and so any variations in the splittings could be important.
The structure of this paper is as follows: In Section Section 2 we describe how the frequency splittings were obtained and discuss the presence of any periodicities. In Section 3 we compare the results obtained from the real data with those obtained from simulated data. In Section Section 4 the consequences of the variations in the splittings for inversions of the internal rotation profile of the Sun and other stars are discussed. The main results of the paper are summarized in Section Section 5.
2Extraction of rotational splitting frequencies
We have analyzed the p-mode rotational splittings observed by the Birmingham Solar-Oscillations Network [?] during the last two solar cycles in their entirety i.e. from 1986 April 14 to 2010 April 7. The Global Oscillations at Low Frequencies [?] instrument on board the ESA/NASA Solar and Heliospheric Observatory (SOHO) spacecraft has been collecting data since 1996 and so we have been able to analyze the velocity data (following the methods described in [?]) covering almost the entirety of solar cycle 23, i.e., from 1996 April 11 to 2010 April 7.
The precision with which p-mode parameters, such as the rotational splittings, can be determined is directly related to the length of data set under consideration. Consequently, p-mode frequencies are often determined from data sets whose lengths are of the order of years. However, a compromise must be made here regarding the appropriate length of data set for study so that changes as the solar cycle evolves can be resolved. Therefore the observations made by GOLF and BiSON were divided into 182.5-day-long independent subsets. After 1996 April 11, when both sets of data were available, we ensured that the start times of the subsets from each observational program were the same.
Estimates of the mode splittings were extracted from each subset by fitting a modified Lorentzian model to the data using a standard likelihood maximization method, where the rotational splitting was a free parameter in the model.
The unweighted mean mode splitting for a particular subset was then determined by averaging estimates in the frequency range . The lower limit of this frequency range (i.e., ) was the lowest frequency for which the signal-to-noise value allowed good fits to be obtained in all the subsets. Furthermore, below this frequency the mode frequencies experience almost no solar cycle shift [?]. Above it becomes difficult to accurately determine the mode splittings as the widths of the modes become larger than the splitting itself and so the mode components overlap in frequency, eventually becoming a single peak. The mean splittings were calculated for the and modes separately.
Two different fitting codes have been used to extract the mode frequencies [?], both giving similar results. For GOLF the mean absolute difference in the estimated average splittings was for and for , and the difference in the estimated splittings from the two codes was always within . For BiSON, the mean absolute difference in the estimated splittings was for and for . Differences in the way the fitting procedures dealt with the BiSON window function meant that before 1992, when the BiSON duty cycle was often below 50 per cent, the splittings obtained from the BiSON data were separated by up to . However, after 1992 the maximum separation of the BiSON splittings obtained by the two different fitting methods was . For clarity, we only show the results of one method, which is described in [?].
Many authors have examined the uncertainties associated with estimating splittings from Sun-as-a-star observations and the biases associated with the estimated splittings are now well known [?]. We follow the advice of [?] in order to minimize these biases. For example, [?] found that one of the main sources of bias in estimated splittings arises if the amplitude ratios of the components inside a multiplet are not set correctly. We therefore fix these ratios at previously determined and accepted values for BiSON and GOLF [?]. Furthermore, the amplitude of the majority of known systematic errors become appreciable above approximately , which is above the upper limit on the frequency range we examine here.
2.1Variation with time of the rotational splitting frequencies
Let us now consider whether the rotational splitting measured with the methodology outlined above show any significant variations with time. Figure 1 shows the mean determined splittings as a function of time, , for GOLF and BiSON. There is no obvious evidence for an 11-yr solar cycle effect on the splittings. This is consistent with previous results [?]. However, there are discernible mid-term variations, which is also in agreement with the results of [?]. There is significant correlation between the variations in the splittings observed in the two data sets (see Table 1), such that there is less than a 0.05 per cent probability of the correlations occurring by chance.
There is no correlation between the and mode splittings for either the GOLF or the BiSON data. The splittings of the modes show more variation than the mode splittings. This is not surprising and indicates that, in Sun-as-a-star data, the rotational splittings estimated from modes are more stable than the splittings estimated using modes [?]. However, as we now show, this difference is accounted for by the formal uncertainties associated with the estimated splittings. We define to be the temporal mean of and . Fig. Figure 2, which shows (where is the formal uncertainty associated with ), demonstrates that the scatter in the plots is similar for and 2, indicating that the relative variation in the splittings is comparable. Table ? gives the standard deviations of the error normalized residuals, . The standard deviations of the normalized residuals are close to unity for both GOLF and BiSON data, which is what we would expect to observe for a totally random process with appropriate error bars.
2.2Periodograms of the observed splittings
To further investigate whether the observed variations are significant, we computed periodograms of the observed splittings. These periodograms are plotted in Fig. ?. In calculating the periodograms the data have been oversampled by a factor of 10. The cut-off frequency at the low end of the periodograms, below which no information can be obtained, is and for GOLF and BiSON respectively. Also plotted in Fig. ? are the 1 per cent false alarm significance levels [?], which were determined using Monte Carlo simulations. 200,000 noise data sets were simulated, using a normal distribution random number generator, to mimic those plotted in Fig. Figure 1. The standard deviation of each point on the data set was taken as the uncertainty associated with the rotational splittings plotted in Fig. Figure 1. The simulated data sets were then used to create periodograms and the amplitudes observed in the simulated periodograms were used to define the 1 per cent false alarm significance levels. The 1 per cent significance level is somewhat arbitrary, but was chosen a priori so that the number of expected false detections is less than unity [?].
There are no significant peaks (at a 1 per cent level) in the periodogram of the GOLF data indicating that there are no significant periodicities in the GOLF rotational splittings. However, there are three significant peaks in the BiSON periodogram and one significant periodicity in the BiSON splittings. The significant peak in the BiSON periodogram (at ) does not correspond to any of the significant peaks in the BiSON periodogram.
Table 3 contains the probabilities of those peaks in the BiSON periodogram where the probability of observing the peaks by change is less than 1 per cent. Although two of the peaks are only marginally significant the other peak (at ) and the peak are well below the significance level. The uncertainties associated with the splittings would have to be increased by 10 per cent for the modes and 17 per cent for the modes before these peaks are no longer significant.
3Comparison with solarFLAG data
To help determine whether the signals observed in the BiSON data are really solar in origin (e.g. perhaps associated with the solar cycle) or simply an artifact of the data analysis procedure we examined data that were simulated for the solar Fitting at Low Angular degree Group [?]. The data were simulated in the time domain and were designed to mimic Sun-as-a-star observations. The simulated oscillations were stochastically excited and damped with lifetimes analogous to those of real solar oscillations. Five sets of solarFLAG data were examined. The input mode frequencies were the same in each simulated data set and were constant in time. The mode splittings were fixed at . However, the data simulated the stochastic nature of the mode excitation giving access to different mode and noise realizations. The simulated data were approximately 9 yr in length but were split into subsets of 182.5 d for analysis.
The splittings obtained from a typical solarFLAG data set (FLAGa) are plotted in Figure 4. Almost as much variation is observed in the splittings obtained from the solarFLAG data as was observed in the real data. Table ? shows the mean observed splittings and the maximum absolute deviation from the mean for both the real and simulated data. The maximum absolute deviations of the BiSON data are slightly larger than all of the simulated solarFLAG data sets, while the GOLF data are in reasonable agreement with the solarFLAG data. Since the only source of variation in the solarFLAG data is the realization noise these results imply that the changes seen in the real data are more than likely due to realization noise.
3.1Impact of the BiSON window function
The discrepancy between the maximum absolute deviations observed in the BiSON results and the GOLF and solarFLAG results can be explained by the impact of the window functions of the respective data sets. The duty cycles of the BiSON subsets range from 22 to 88 per cent, with the average duty cycle being 69 per cent. Normally the duty cycle of the GOLF data is above 90 per cent and the mean duty cycle is 95 per cent. However, the lowest observed duty cycle of the GOLF subsets was 42 per cent, which coincides with when control of SOHO was temporarily lost, known as the “SOHO vacation”. BiSON window functions were applied to the GOLF data such that the 182.5 d BiSON subset from which the window function was extracted matched in time the GOLF subset to which it was applied. Table 4 shows that this increases the maximum absolute deviations of the GOLF data. Furthermore, although not shown here, a marginally significant peak (at a 1 per cent significance level) was observed at in the periodogram of the splittings.
Various BiSON window functions were also imposed on the solarFLAG data, with duty cycles covering the range observed here. In addition to the splittings observed in the 100 per cent duty cycle solarFLAG data, Figure 4 also shows the splittings obtained from one set of solarFLAG data when a 182.5-d-long BiSON window function, with a duty cycle of 88 per cent, was imposed upon each 182.5-d-long subset. This is just one example of the window functions that were imposed. The same BiSON window function was imposed on each FLAG subset so that the average duty cycle could be varied systematically. When some window functions were imposed the maximum absolute deviations were observed to increase to values similar to those observed in the BiSON data. However, in other cases the maximum absolute deviation was observed to decrease. Table 5 shows that, although altering the window function changes the maximum absolute deviation, the obtained deviation is not dependent on the value of the duty cycle. For example, the rotational splittings observed in FLAGa varied more when the duty cycle was 88 per cent than when the duty cycle was 47 per cent.
This is consistent with the surmise that the observed variations are due to realization noise and that the window function could affect the variations in the rotational splittings both constructively and destructively via its interaction with the noise. Consider a period of time when the realization noise shifts the obtained splitting away from the true value. If the window function blanks this period of time (i.e. there is no data available) then the estimated rotational splitting will be close to the true value and the observed variation in the rotational splittings will decrease. However, the window function could equally favour times when the realization noise shifts the observed rotational splitting away from the true value. The variation in the rotational splittings will then increase.
This hypothesis is strengthened further by analysis of the periodograms of the solarFLAG data. Fig. ? shows periodograms of the splittings obtained from the solarFLAG data sets when various BiSON window functions were imposed. We start by considering the top-left panel of Figure 5, which shows periodograms of the splittings obtained from the FLAGc data when BiSON window functions were imposed such that the duty cycle of the data was 59 per cent, 68 per cent and 76 per cent respectively. Even though the splittings included in the simulated FLAG data were constant, significant peaks (at a 1 per cent false alarm level) are observed in the periodograms. Significant peaks are observed when the duty cycle was both 59 per cent and 76 per cent but not 68 per cent. Similarly in the top middle panel of Fig. Figure 5 significant peaks (at a 1 per cent false alarm level) are observed in periodograms of the splittings obtained from the FLAGd data when the duty cycle is 68 per cent but not 47 per cent or 100 per cent. This indicates that the likelihood of observing a significant peak is not simply a function of the duty cycle. The bottom left and middle panels show that the same is true for the splittings: The bottom left panel shows that significant peaks are observed in the splittings obtained from the FLAGe data when the duty cycle was 39 per cent and 100 per cent but not 59 per cent. The bottom middle panel shows that significant peaks are observed in the splittings obtained from the FLAGb data when the duty cycle was 76 per cent but not 47 per cent or 100 per cent. This shows that more important than the actual duty cycle itself is the manner in which the window function interacts with the realisation noise. The top and bottom right-hand panels of Fig. ? show just how different the periodograms can look simply by changing the window function of the data.
|GOLF (BiSON window function)|
|Duty cycle (%)|
We note here that Tables Table 4 and Table 5 indicate that the mean rotational splittings are consistently higher than the input value of , because the rotationally split components overlap in frequency giving rise to a positive bias [?].
The results of the solar-FLAG simulations show that significant quasi-periodic variations can be observed in the splittings even when no actual, underlying variations in the splittings are present. The BiSON window function can act to both enhance or reduce the significance of the variation.
4Consequences for rotation profile inversions
Let us examine the consequences of these apparent variations in the splittings for both helioseismic rotation inversion profiles and asteroseismic inferences of the mean internal rotation rate. Each rotational splitting is a spatially weighted average of certain volumes of a star’s internal rotation rate and various inversion techniques have been developed to infer the Sun’s internal rotation profile from the observed splittings. We now consider the consequences of the observed variations in the splittings on inferences that can be made about the rotation profiles of the Sun and other stars.
4.1Inversions of the Sun’s internal rotation profile
We have seen that, although the mean BiSON splitting is , the maximum and minimum observed splittings are and respectively. We have already shown that the observed variation in the splittings is just an artifact associated with the realization noise in the data and the process by which the data were analyzed. However, one may erroneously ascribe the splitting variations to a physical change in the rotation rate of a region of the solar interior. We have split the solar interior into 4 regions; the solar core , the radiative zone , the convection zone , and the near-surface shear layer . We have then determined by how much the rotation rate of that layer would need to change, relative to an average 2D Regularized Least Squares (RLS) inversion [?] of the solar rotation profile, to explain the observed differences in the rotational splittings. The inversions [?] were made using 15 yrs of Global Oscillations Network Group [?] data. This analysis is similar to that performed for the Sun by [?] and more details can be found in Appendix ?.
Using the above mentioned maximum and minimum measured splittings we find that the maximum departure from the mean BiSON splitting was . Table 6 shows the change in rotation rate, , required to change the observed splitting of and 2, modes by (results for other modes can be found in Appendix ?). The results indicate that the changes required in the rotation rate of near-surface regions to explain a difference in the splitting of are not unreasonably large. However, observations of higher- modes rule out variations in the near-surface regions even of this magnitude (see Appendix ?). Information from higher- modes dominate and constrain inversions of the rotation in outer regions of the solar interior at better than a 0.5 per cent level. Although variations with time in the higher- splittings are detectable (due to the torsional oscillation) these variations are less than [?] and so they are approximately 100 times smaller than the deviations observed here. Large changes in the core rotation rate are required to explain the observed change in the rotational splittings. Although the required changes are so large as to be unlikely
|Shell||Inner radius of|
|surface shear layer||0.95||0.49||0.49|
4.2The rotation profile of a solar-type star
Asteroseismic data often consist of just a few months of observations in which only the low- modes are visible. Thus asteroseismic data are often similar to the 182.5 d Sun-as-a-star subsets analyzed here. We note here that this analysis is relevant for main sequence stars for which no mixed modes or g modes can be detected and used to aid inferences of the internal rotation rate. Therefore it is difficult to make inversions of the internal rotation profile of a star. However, lightcurves of a star, such as those observed by CoRoT and Kepler, can show rotational modulation due to the presence of magnetic stellar activity, such as spots and active regions. The modulation of the light curve can, therefore, be used to determine the surface rotation rate of a star [?], which can then be compared to the mean internal rotation rate determined from the oscillations. The sidereal rotation rate of the Sun, as estimated from the rotation of sunspots at the surface, is approximately 27 d. The mean rotation rate of the solar interior, determined from low-l rotational frequency splittings, is 28 d. However, if we were to use the extreme values of the splittings obtained from a single 182.5 d subset, similar to those available in asteroseismic studies, we would find the mean rotation rate to be either or . One may, therefore, erroneously conclude that the average rotation rate of the stellar interior is different to the surface rotation rate.
As with previous studies, [?], we find no 11-yr solar cycle variation in the rotational splittings. Discernible quasi-periodic mid-term variations are present in the splittings. However, we show that these variations are due to realization noise.
Many authors have discussed the biases associated with estimating splittings from Sun-as-a-star data [?]. We have shown that, despite taking the advice of [?], the realization noise has a larger effect on rotational splittings than accounted for by formal uncertainties, implying that the uncertainties on the fitted splittings are underestimated. The effect of the realization noise produces an erroneous apparent mid-term signal in the observed splittings that could, potentially, look like the quasi-biennial solar-cycle related signal that is observed in the mode frequencies [?]. Although we have shown that the mid-term signal is an artifact associated with the realization noise and the data analysis process the periodicities in the splittings may be mistakenly interpreted as a physical variation in the rotation rate. Such variations are unlikely because they would require either very large changes in the core rotation or changes in the near surface that, although smaller, are ruled out through observations of higher- modes.
It is not surprising that the splittings observed by BiSON and GOLF are highly correlated, since they both are observing the same star. In fact we expect the splittings to be correlated whether these variations are of solar origin or whether, as we have shown, they are due to realization noise. [?] found that there was a strong degree of correlation between numerous mode parameters, such as mode linewidths and power densities, observed in GOLF and BiSON data. They concluded that this showed that the estimated mode parameters were dominated by the same mode realization noise (the signature of stochastic excitation). [?] observed that, after the solar-activity dependence had been removed from the frequencies, significant fluctuations in the frequencies were highly correlated between BiSON, GONG, and MDI data. [?] explained these correlations in terms of the stochastic nature of the mode excitation, which could be interpreted as a form of realization noise. The high level of correlation observed here between the GOLF and BiSON splittings is also due to the mode realization noise. This supports the conclusion that the mode realization noise is responsible for the observed variation in the splittings, which again demonstrates how careful one needs to be, particularly when data have correlated noise sources.
Although the spurious variations do not have a significant effect on inferred rotation profiles of the Sun they could be important for any asteroseismic inferences of the average internal rotation rate of main sequence stars for which g modes and mixed modes have not been detected. Firstly because, unlike for the Sun, we must rely on the splittings estimated from low- modes only. Secondly, because we often have only one relatively short data set, of the order of a few months, from which to extract the splittings. An additional consideration is that the accuracy of the determined splitting is dependent on the size of the splittings relative to the width of the modes i.e. whether the modes overlap in frequency [?]. We have shown that care must be taken in using average splittings (extracted from data sets of a few months duration) to infer the mean internal rotation rate (which, for example, might be compared in asteroseismic studies to the surface rotation rate inferred directly from rotational modulation of the light curve, e.g. [?]). Artificial data, such as those available through AsteroFLAG [?], should be used to test for potential biases in the results. Care should be taken when interpreting the observed splittings of low- modes in both astero- and helioseismic data.
AInversion process: Forward calculations
A rotation profile of the Sun was simulated, based on the average of 2D inversion profiles computed over 15 yrs using a Regularized Least Squares [?] inversion method and Global Oscillations Network Group data [?]. The GONG instrument make spatially resolved observations of the Sun and is thus able to observe much higher- modes than those examined here. The GONG observations were used to build a latitudinally dependent rotation profile of the Sun between 0.5 and 1.0 solar radii . Below the rotation was taken to be constant and equal to the mean inferred rotation rate across all latitudes at . We varied the rotation rate in four different regions of the solar interior; the solar core , the radiative zone , the convection zone , and the near-surface shear layer . The splittings were calculated when the rotation rate in each region was varied between . The and 2 splittings, , were then determined [?].
To first order, the magnitudes of the splittings vary linearly with the rotation rate in a given region. That linear factor depends on the position of the specific region in the solar interior and so can be used to infer the change in rotation rate in a specific region that is required to produce a change in splitting of . These values are plotted in Figure 6 and are listed in Table 6 for , and , (corresponding to modes at 2561 and respectively). The results in Figure 6 and Table Table 6 indicate by how much the rotation rate in a particular region would have to deviate from the values obtained from the 15-yr 2D RLS inversions to explain the observed results. The required deviation in rotation rate change is similar for the radiative zone, the convection zone and the surface shear layer and is approximately constant across the range of examined here. However, the change in rotation rate required in the core is substantially larger and depends on .
As expected from the rotation kernels, the observed splittings are much more sensitive to surface regions than to the core. The changes in Table 6 can be compared to the mean rotation rate at the surface, which is approximately . The rotation rate in the core would have to be approximately 5 times faster than the mean solar surface rotation rate to account for an increase in the splitting of , while the rotation rate in the other regions of the solar interior would have to change by more than 100 per cent of the mean solar value. In other words, even if the shell responsible for the change in splitting is close to the surface, a large, but not inconceivable, change in rotation rate is required to explain a difference in the observed splitting of .
BVariations in the inferred solar rotation profile
We have combined the low-degree splittings observed in the BiSON and GOLF data with a set of intermediate-degree splittings made from the average of three 72 d Michelson Doppler Imager [?] sets observed between 1996 and 1997. Each set of splittings was then inverted using 2D RLS methods [?]. Figure 7 shows residuals in the inferred rotation profile at the equator once a mean profile has been subtracted. In regions where no information is available (i.e. the core) the 2D RLS inversion method extrapolates a solution. In these regions the averaging kernels, the functions that describe how the estimate of the solution corresponds to a spatial average of the underlying true solution, are not well localized. The rotation profiles plotted in Figure 7 were generated using two sets of splittings measured by BiSON and GOLF. These were chosen because they correspond to epochs when the observed splitting deviates noticeably from the mean. In each case the rotation profiles differ from the mean, and from each other, in the deep interior but these differences are less than the formal uncertainties for the inversions. Therefore, the spurious changes in the observed splittings do not have a significant effect on the inferred rotation profile of the Sun. This is not entirely unexpected because even though the (model-dependent) spatial weighting functions (kernels) of low- modes do probe the solar core, they are still more sensitive to the surface regions. Indeed, the sound speed in the deep interior is much larger than near the surface, so the time the modes dwell in the core is small relative to the outer regions.
We thank the referee for their insightful comments. This paper utilizes data collected by the Birmingham Solar-Oscillations Network (BiSON), which is funded by the UK Science Technology and Facilities Council (STFC). We thank the members of the BiSON team, colleagues at our host institutes, and all others, past and present, who have been associated with BiSON. 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. A.M.B., W.J.C., and Y.E. acknowledge the financial support of STFC. D.S. acknowledges the support from the Spanish National Research Plan (grant PNAyA2007-62650) and from CNES. R.A.G. thanks the support of the CNES/GOLF grant at the CEA/Saclay. R.H. thanks the National Solar Observatory for computing resources. We thank R. New for helpful comments and discussions. NCAR is supported by the National Science Foundation.
- The rotation rate in the core is not well constrained through p modes [?]and so, in this paper, we assumed a rate similar to the rest of the radiative zone. However, we note that g modes would provide tighter constraints on the core rotation rate [?].