Stellar magnetic activity and variability of oscillation parameters: An investigation of 24 solarlike stars observed by Kepler
^{1} KiepenheuerInstitut für Sonnenphysik, Schöneckstraße 6, 79104 Freiburg, Germany
email: kiefer@kis.unifreiburg.de
^{2} School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK
^{3} Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
Received: 9 March 2016
Accepted: 3 November 2016
Context. The Sun and solarlike stars undergo activity cycles for which the underlying mechanisms are not well understood. The oscillations of the Sun are known to vary with its activity cycle and these changes provide diagnostics on the conditions below the photosphere. Kepler has detected solarlike oscillations in hundreds of stars but as of yet, no widespread detection of signatures of magnetic activity cycles in the oscillation parameters of these stars have been reported.
Aims. We analysed the photometric short cadence Kepler time series of a set of 24 solarlike stars, which were observed for at least 960 d each, with the aim to find signatures of stellar magnetic activity in the oscillation parameters.
Methods. We analyse the temporal evolution of oscillation parameters by measuring mode frequency shifts, changes in the height of the pmode envelope, as well as granulation timescales.
Results. For 23 of the 24 investigated stars, we find significant frequency shifts in time. We present evidence for magnetic activity in six of these stars. We find that the amplitude of the frequency shifts decreases with stellar age and rotation period. For KIC 8006161 (the most prominent example), we find that frequency shifts are smallest for the lowest and largest for the highest pmode frequencies, as they are for the Sun.
Conclusions. These findings show that magnetic activity can be routinely observed in the oscillation parameters for solarlike stars, which opens up the possibility of placing the solar activity cycle in the context of other stars by asteroseismology.
Key words: asteroseismology / stars: activity / stars: oscillations / stars: solartype / stars: magnetic field / methods: data analysis
© ESO, 2017
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0),
which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Understanding the solar dynamo – exactly how and where it operates within the Sun – is a major open task in modern astrophysics. Helio and asteroseismology are outstanding tools to probe the interiors of the Sun and other stars. The vast amount of data that has become available from satellite missions like CoRoT (Baglin et al. 2006; Auvergne et al. 2009) and Kepler (Borucki et al. 2010; Koch et al. 2010) enables us to investigate stellar cycles for a wide range of stellar types. By putting the Sun and its magnetic activity cycle into the context of similar stars, we can widen our perspective on the basic processes at work in dynamos in different physical environments and may be able to decide between opposing dynamo theories.
Overview of the basic data of the investigated time series and their periodograms.
The better part of known activity cycles of mainsequence stars was detected from Mount Wilson Observatory (Wilson 1978; Baliunas et al. 1995). Their measurements of fluxes in the line cores of Ca II H and K, which can be used as proxies for stellar activity (Hall 2008, and references therein), suggest a relation between the age, rotation, and activity of stars. For a subset of solartype stars Baliunas et al. (1995) found that younger stars exhibit higher average levels of activity and rapid rotation rates, while older stars with a slower rotation usually have cycles with lower levels of magnetic activity. A relation of this kind has already been proposed by Skumanich (1972), who found that stellar Ca II emissions decay as the inverse square root of the age. Brandenburg et al. (1998) and Saar & Brandenburg (1999) made an empirical classification of stars from the Mount Wilson Observatory sample in a rotationactivity scheme. They found that most stars fall into one of three distinct branches when examining the ratio of cycle and rotation periods as a function of chromospheric activity. Karoff et al. (2013) used groundbased observations in combination with data from the Kepler satellite to test agerotationactivity relations.
Helioseismology found that the frequencies of solar acoustic oscillations (pmodes) are positively correlated with solar magnetic activity (e.g. Woodard & Noyes 1985; Libbrecht & Woodard 1990; JimenezReyes et al. 1998), whereas mode amplitudes are anticorrelated with it (e.g. Pallé et al. 1990; Komm et al. 2000; Broomhall et al. 2015). Activity related frequency shifts are larger for modes of higher frequency (JimenezReyes et al. 1998), because these modes have their maximum sensitivity to structural changes in the solar interior in shallower layers than modes of low frequency (Basu et al. 2012). This was recently used by Salabert et al. (2015) to study the change of subsurface solar activity as a function of time and for different depths.
In an analysis of photometric data of the CoRoT satellite, García et al. (2010) found evidence for an activity cycle with P_{cyc}> 120 d in the F5V star HD 49933. They found shifts in the frequencies of the star’s pmodes which showed a cyclelike behaviour and, in addition, changes in mode amplitudes which were anticorrelated to the frequency shifts. Later, Salabert et al. (2011) found the frequency shifts of HD 49933 to be dependent on mode frequency. The Kepler satellite produced a huge data set on solarlike stars with asteroseismic import, which shows great promise for many more detections of the kind García et al. made. Vida et al. (2014) found hints of cyclic activity in the Kepler data of nine fastrotating latetype stars. Recently, Salabert et al. (2016b) investigated Kepler and spectroscopic data from the young solar analog KIC 10644253. They found variability in the pmode frequencies of this star with a modulation of about 1.5 yr.
In this paper, we describe our efforts to find signatures of stellar magnetic activity in the oscillation parameters of solarlike stars, focussing on the shifts of pmode frequencies and the variation of the height of the pmode envelope. In Sect. 2 we describe our approach to measure these quantities and to estimate the errors on them. Section 3 is dedicated to a presentation of the results and a more detailed discussion of two stars from our sample. Section 4 gives a conclusion of our findings.
2. Methodology
2.1. Data and computation of periodograms
We analysed Kepler data from 24 stars, which were observed in the satellite’s short cadence mode. Only data from the Kepler data release 25 was used in this work. In this data release, the scrambling of the shortcadence collateral smear data, which was reported in the Global Erratum for Kepler Q0Q17 & K2 C0C5 ShortCadence Data, and the Dynablack calibration problem, which was reported in the Data Release 24 Notes Q0–Q17 Erratum, are corrected. These 24 stars were selected because they are known to exhibit solar like oscillations (e.g. Mathur et al. 2012; Appourchaux et al. 2012; Huber et al. 2013) and are known to be similar to the Sun regarding fundamental parameters like mass, radius, and effective temperature (cf. Tables A.1 and A.2 and references therein). No strict criteria were applied to define this special sample of stars. It is therefore just a first assessment of the variability of stellar pmode parameters in the Kepler short cadence sample. The full set of stars, which were observed in Kepler’s short cadence mode over extended periods of time, is considerably larger and is a worthwhile subject for future investigation.
The coverage of data for these stars ranges from 10 to 13 Kepler quarters of observation (960 to 1147 d). A list of the Kepler input catalogue (KIC) numbers of the stars and the available quarters of data is given in the first two columns of Table 1. Missing quarters are listed in the notes for Table 1. The third column gives the number of days the data for each star spans, counting from the first day of observation to the last. The light curves were cleaned of jumps, drifts, and outliers similar to the method described by García et al. (2011).
To investigate the temporal variation of pmode parameters, the time series were divided into consecutive, overlapping segments of equal length. To compromise between frequency resolution and number of independent segments, we chose a typical length of 150 d for the segments. From one segment to the next, the starting point was shifted by 50 d, resulting in a three time overlap for 150 d long segments. If the errors on the frequency shifts that were eventually measured allowed a decrease in segment length and hence frequency resolution of the periodograms, we reduced it to 100 d. The shift of the starting point remained at 50 d, resulting in a two time overlap. The lengths of the segments are given in the fourth column of Table 1.
Since the time series are slightly unevenly sampled and are affected by observational gaps, we used the LombScargle method (Lomb 1976; Scargle 1982) in the fast implementation of Press (2007) to calculate the periodograms (LSperiodogram). The resulting periodograms are scaled to comply with Parseval’s theorem.
The observational gaps lead to spectral leakage of mode power. To study the effect of the gaps in the Kepler time series on our analysis, we used the inpainting software of Pires et al. (2015) to fill gaps up to a length of 20 d (García et al. 2014b). For none of our targets did the high frequency noise level decrease significantly. Also, the results described in the subsequent sections, frequency shifts and mode heights, do not change in a way that would make necessary the inpainting of the time series. This is because the Kepler gap structure does not have a large effect on the periodograms of solarlike stars as long as the stars are no fast rotators (García et al. 2014b) and the periodograms are computed from time series segments with a length of only 100 or 150 d.
To get the same frequency axis for the LSperiodogram of each data segment of a given star, the LSperiodograms were interpolated onto the frequency grid of the segment with the lowest fill factor. This way, the LSperiodograms can later be correlated with each other. The minimum fill factor of segments was set to 50%. Thereby, we ensured a reasonable frequency resolution for the LSperiodograms of all the segments of a star. In Fig. 1, the 0.81 μHz (7 bins) boxcar smoothed segments between 3570−3610 μHz of the LSperiodograms of the first and the nineteenth segment of the time series of KIC 8006161 are shown. Even by visual inspection of the peak at ~ 3590 μHz, a shift towards higher frequency can be spotted for the periodogram of the nineteenth segment.
Fig. 1 Section of the LSperiodograms of the first and nineteenth segment of KIC 8006161 in black and red colour, respectively. The LSperiodograms are boxcar smoothed over 0.81 μHz. 

Open with DEXTER 
2.2. Measuring the frequency shifts
The first segment of each star was used as the reference for defining zero frequency shift. Due to this, the level of zero frequency shift can be at any level of a stellar activity cycle, e.g. its maximum. We computed the shift relative to the reference segment by calculation of the crosscorrelation of that frequency range of the LSperiodograms which contains the pmodes and fitting a Lorentzian to the crosscorrelation. We selected this frequency range by visually determining the lower and upper boundaries between which pmode peaks can be spotted in the smoothed LSperiodogram of the entire time series of each star (see last column of Table 1). Recently, a crosscorrelation technique for measuring mode frequency shifts was presented by Régulo et al. (2016). Their method to measure the frequency shifts is very similar to the approach presented here. However, our approach to estimate the error on the frequency shift differs in several key ways from that of Régulo et al. (2016).
In the top panel of Fig. 2, the part between ± 5 μHz of the crosscorrelation between the pmode region of two realisations (see Sect. 2.3) of the first segment of KIC 8006161 is depicted. The Lorentzian fit to the data is also shown. We have also indicated the zero frequency lag and the centre of the fitted Lorentzian. The bottom panel of Fig. 2 depicts the same as the top panel, but for the crosscorrelation between the first and the nineteenth segment of the time series from KIC 8006161. Here, the maximum of the crosscorrelation is clearly shifted towards positive values. Peaks of stochastically excited solar and stellar pmode oscillations appear to be of Lorentzian shape. The Lorentzian is invariant under convolution with a Lorentzian, hence the crosscorrelation of two of them can be fitted with a Lorentzian. In our analysis we used a symmetric Lorentzian profile. We found that the extra parameter introduced by a skewness does not decrease the error on the frequency shift or increase the stability of the fits.
Fig. 2 Top panel: crosscorrelation between the part of interest of the LSperiodograms of two realisations of the first segment of KIC 8006161 (black). The Lorentzian fit to the crosscorrelation is shown in red. The dotted black vertical line marks zero frequency lag, the dashed red line marks the centre of the fitted Lorentzian. Bottom panel: same as in the top panel, but for the crosscorrelation between a realisation of the first and one of the nineteenth segment. 

Open with DEXTER 
2.3. Resampling approach
Since the correlation structure of the estimated crosscorrelation function depends on the true unknown cross and autocorrelation function in a complicated manner, we use a resampling approach to estimate the error on the frequency shift (see Bartlett 1978 and e.g. Chaplin et al. 2007b). We generate a sample of new LSperiodograms, which follow the statistical properties of the original periodogram, and from which we compute the sample mean and sample error. In our resampling approach we neglect effects from uneven sampling, since it is small for Kepler data. The largest effects on our analysis are assumed to result from gaps in the data.
For each data segment i, we draw B = 200 complex random series for { ± ν_{k} } _{k = 1,...,N} from a zero mean normal distribution N(0,1) weighted by the estimated spectral power Ŝ_{i}(ν_{k}), such that and . For ν = 0, we set , reflecting a zero mean time series. An estimate Ŝ_{i} of the true unknown spectral density S_{i} is obtained by smoothing the LSperiodogram with a boxcar window. The width of the boxcar varies from star to star, ranging from 0.38–1.52 μHz (5–12 bins). We chose the width of the boxcar depending on the frequency resolution of the periodograms and the width of the pmodes of each individual star. The complex random series may be considered as realisations of Fourier transforms that would be obtained in the ideal case of evenly sampled data without observational gaps but with the same spectral density underlying the respective original data. We note that the real and the imaginary part of are by construction independent of each other and with respect to ν_{k}. In order to take into account the effect of gaps in the original data, we generate new time series on an evenly sampled grid { t_{n} } _{n = 1,··· ,N} by the inverse Fourier transform of . The grid in time approximates the one for the original data. Data points at the respective grid points with observation gaps are removed. These new time series are treated as described above to compute frequency shifts from the crosscorrelation function of LSperiodograms. We further note that these resampled time series have no physical meaning, since the phases between the real and imaginary parts of the Fourier transforms are random. However, the procedure described conserves the linear spectral properties of the original data.
Amplitudes of frequency shifts and their significance.
We emphasise that the LSperiodogram is computed with oversampling factor 1. In case of evenly sampled data without gaps, the LSperiodogram is then equivalent to the classical Fourier based periodogram and follows the same statistical distribution (Scargle 1982). An oversampling factor greater than 1 corresponds to an interpolation of the periodogram along frequency bins and gives correlated values at adjacent frequency bins. In this study, we aim to suppress such additional correlations since they lead to an unwanted complex correlation structure in the error analysis and result in a systematic underestimation of the error on the frequency shifts.
Fig. 3 Frequency shifts of solar pmodes from BiSON data with 1σ error bars (black). The seven month boxcar smoothed monthly sunspot number is shown in red. Source of the sunspot number: WDCSILSO, Royal Observatory of Belgium, Brussels. 

Open with DEXTER 
Each realisation of the LSperiodogram of a given segment is crosscorrelated with a realisation of the LSperiodogram of the reference segment. For the reference segment, two samples of 200 realisations are drawn. The crosscorrelation analysis between these two sets yields the error bar on the zero level of the frequency shifts. The standard deviation of the centroid values of the Lorentzian profile of the 200 fits is used as the error for the shift of the pmodes, while their mean is used as the value of the shift. The value B = 200 was chosen, because it is a conservative number of realisations for which the variance of the estimator of the error is small compared to the variance of the distribution of the centroid values of the frequency shifts. Hence, the standard deviation of the 200 centroid values is a reliable estimate for the error on the measured frequency shifts.
2.4. Consistency check with BiSON data
For a consistency check of our method, we analysed data from the BiSON network (Davies et al. 2014). The time series encompasses the period from 1985–2014, covering about two and a half solar cycles. Figure 3 depicts the frequency shifts and the associated 1σ error bars as black data points. The monthly sunspot number is shown as a continuous red line (boxcar smoothed over seven months). It is scaled to match the amplitude of the variation of the pmode frequencies in this figure. In this case, we set the length of the segments to 200 d and the step length to 100 d in order to increase the clarity of the plot by reducing the number of data points. We selected the frequency range between 1800–3800 μHz for the crosscorrelation of the periodograms. The frequency shift closely follows the sunspot number, a proxy of solar magnetic activity.
For the amplitude of the solar pmode frequency shift we find a value of 0.62 ± 0.06 μHz. This is in agreement with previous investigations of the pmode frequency shifts obtained with the crosscorrelation method using BiSON data (Chaplin et al. 2007b).
2.5. Height of the mode envelope
As we are considering the temporal evolution of the frequencies of all pmodes in the periodogram, we compute a proxy for the amplitudes of the pmodes by measuring their excess over the background in the periodogram.
For this, we have fitted a background to each of the segments and recorded the summary statistics. We used a variant of “model F” as described in and recommended by Kallinger et al. (2014), which consists of two Harveylike profiles (Harvey 1985), a Gaussian pmode envelope, and an instrumental noise component. Following Kallinger et al., we restricted the exponents of the Harveylike profiles to a value of four.
Here our purpose was to study the impact of stellar activity. So, in order to isolate the activity signal of the pmodes from the low frequency activity signal due to surface modulation, we restricted the fitting range to frequencies above 400 μHz. This lower limit excludes regions of the periodogram that can also contain activity signals, but mitigates problems caused by a Kepler noise feature found at 374 μHz. As the surface activity signal is excluded, we did not include a powerlaw to account for surface activity in the model.
To fit the model, we used the affine invariant Markov Chain Monte Carlo code emcee (ForemanMackey et al. 2013). From the resulting Markov Chains we quote the summary statistics of the posterior probability distributions, which are normally distributed, as the median value and the 15.9th and 84.1th percentile to define the highest posterior density interval, that is, the credible interval. We used the standard Gaussian form for the mode envelope (3)where a is the height of the Gaussian, b is the frequency of maximum height, and c is the width of the Gaussian. In the subsequent sections, we will refer to a as the height of the pmode envelope, or simply height.
To exclude contamination of the measured mode heights with a temporally varying instrumental noise component, which might for example arise from the quarterly satellite rolls performed by Kepler during its nominal mission phase, we checked its temporal variation for all our targets. We found that there are no considerable jumps in the instrumental noise component at the quarter boundaries of the time series we used for the analyses.
3. Results
3.1. Frequency shifts and background parameters
In Table 2, the peaktovalley amplitude of the frequency shifts of all 24 stars are listed. The amplitude A was calculated as the difference between the minimum and maximum values of the frequency shifts. With the propagated error σ_{A}, we also calculated the significance of the observed shift amplitude A/σ_{A}, see Table 2. The observed frequency shifts in our sample range from ≈0.17 μHz for KIC 4914923 to ≈0.95 μHz for KIC 8006161. Only one star from our sample, KIC 9955598, does not exhibit a frequency shift amplitude which is larger than the error and is therefore compatible with zero for all the segments. In the top panels of Figs. B.1–B.23 the frequency shifts are plotted as a function of starting time of the segments, while in the bottom panels the temporal variation of the mode heights are depicted.
As described in Sect. 1, the temporal variation of pmode frequencies can be a proxy for stellar magnetic activity. For some stars, the frequency shifts follow a cyclelike pattern, e.g. KIC 8379927, KIC 8760414, KIC 10644253. Yet, the available time series are not sufficiently long to draw the conclusion that we found clear evidence for a complete cycle in any of the stars of this sample.
We calculated the Spearman rank correlation coefficients ρ between the variation of the mode frequency shifts, mode heights, the granulation timescales of the background Harvey profiles, and the noise level, as well as the corresponding pvalues p for independent segments of all 24 stars. The values are listed in Tables B.1–B.24. The correlation coefficients of the frequency shifts and the mode height variations and their pvalues are also given in Table 3.
For the Sun (e.g. Pallé et al. 1990; Komm et al. 2000; Broomhall et al. 2015) and one more star, the solarlike CoRoT star HD 49933 (García et al. 2010), the shifts in mode frequency and the variation of pmode amplitudes are anticorrelated over the course of their magnetic cycle. For six stars from our sample, the correlation coefficient between the frequency shifts and the variation of the height of the pmode envelope is less than −0.5 and the amplitude of the frequency shift is significant: KIC 6933899, KIC 8006161, KIC 8760414, KIC 9955598, KIC 10644253, KIC 11244118, and KIC 12258514. For these stars, the evidence for magnetic activity influencing the pmode parameters is strongest.
Correlation coefficients of frequency shifts and pmode envelope height variations and their pvalues
Fig. 4 Measured shift amplitudes as a function of effective temperature. Colour and symbol coding: black diamonds represent stars younger than 4 Gyr, red triangles are stars between 4–5 Gyr, orange circles between 5–6 Gyr, blue squares between 6–7 Gyr, and purple asterisks are stars older than 7 Gyr. 

Open with DEXTER 
Fig. 5 Measured shift amplitudes as a function of stellar age. 

Open with DEXTER 
Fig. 6 Measured shift amplitudes as a function of rotation period. 

Open with DEXTER 
There is no clear consensus in the literature, on whether the granulation and supergranulation properties are correlated with the solar activity cycle. Lefebvre et al. (2008) find no correlation between the activity level and granulation velocities or timescales in GOLF data, Muller et al. (2007) show that there is an inphase variation of the contrast of solar granulation with the activity cycle but no such correlation for granulation length scales, Meunier et al. (2008) find smaller supergranules at cycle maximum. Since granulation has different properties for different stellar parameters (e.g. Beeck et al. 2013), we expect different responses by the granulation to magnetic cycles for different stars. In the case of HD 49933, García et al. (2010) showed that the variation of the granulation timescale, the photon noise, and the mode amplitudes are uncorrelated.
The correlation coefficients we present in Tables B.1–B.24 do not present any systematic correlation between the oscillation and background parameters across the sample of stars. For the two stars, which we will discuss in detail in the subsequent sections, KIC 8006161 and KIC 10644253, both of which show evident signs of magnetic activity, we find very different correlations between the investigated parameters.
For KIC 10644253, which was shown to exhibit magnetic activity by Salabert et al. (2016b), we find that the frequency shifts are anticorrelated with the first granulation timescale with ρ = −0.71 and p = 0.07, but correlated with the second timescale with ρ = 0.86 and p = 0.01, cf. Table B.19. The measured mode heights show the opposite behaviour. It is noteworthy that for this evidentially active star (Mathur et al. 2014; Salabert et al. 2016b) the frequency shifts are correlated with changes in the high frequency noise with ρ = 0.50, while the mode heights are anticorrelated to these changes with ρ = −0.75.
In contrast to this, KIC 8006161, which exhibits the strongest frequency shifts of the investigated sample, shows virtually no significant correlations or anticorrelations between any of the background parameters, cf. Table B.11, where the only exception is a correlation between the two granulation timescales with ρ = 0.74 and p = 0.01.
For four stars, KIC 6106415, KIC 6116048, KIC 6603624 and KIC 7680114, the shifts in pmode frequency are strongly correlated with the variation in mode height, with correlation coefficients greater than 0.5. For magnetic activity, which is comparable to that observed on the Sun, this is not expected. As some of these stars (KIC 6116048, for example) show significant temporal variability in their pmode frequencies as well as their mode heights, this might either be due to a dynamo in these stars, which works differently to the solar case, or to another mechanism influencing the oscillation parameters (e.g. temporally varying flows, Roth & Stix 2003). Whether there is relationship between basic stellar parameters (e.g. mass, effective temperature, or thickness of the convection zone) and the correlation of seismic activity proxies as well as granulation parameters, will be investigated in a future study of an expanded sample of Kepler stars.
Stellar magnetic activity need not be necessarily cyclic in the solar sense and it might be that the variation of frequency shifts and mode heights is correlated in physical settings which are different from the solar reference case regarding, for example, differential rotation. The frequency shifts and mode heights of KIC 10018963 for example, see Fig. B.17, are moderately correlated with a correlation coefficient of 0.46. Both, frequencies and mode heights, exhibit significant temporal variations which might be caused by noncyclic sporadic stellar magnetic activity.
3.2. Results for the ensemble
With this systematic investigation of the frequency shifts for a larger sample of stars at hand, we try to find relationships between the measured shift amplitudes and fundamental stellar parameters. We assume here that the observed shifts are solely due to magnetic activity. The available values for stellar radius, mass, and age are listed in Table A.1. The spectral types, effective temperatures, rotation periods, and special features (e.g. binarity), are presented in Table A.2. The references for all values we used are given below these two tables.
In Fig. 4 the measured frequency shift amplitudes are plotted as a function of effective temperature with an indication of the age of the stars.
There are two proposed scaling models for the frequency shifts over a stellar activity cycle: Chaplin et al. (2007a) come to the conclusion that the frequency shifts are directly proportional to the strength of the activity cycle given by ΔR_{HK′}, which is given by the average fraction of the stellar luminosity that is emitted in the Ca II H and K line cores. The second model, presented by Metcalfe et al. (2007), assumes that the frequency shifts are proportional to D/I·ΔR_{HK′}, where D is the depth of the perturbations and I is the mode inertia. In Fig. 9 of Karoff et al. (2009), these two models are compared as a function of effective temperature and for different stellar ages.
The sample of stars presented here does not allow us to distinguish between the two models. The errors on the frequency shift amplitudes are too large and the sample of 24 stars is too small. However, there appears to be a slight tendency towards greater shift amplitudes for hotter stars. This would support the scaling model by Metcalfe et al. (2007). The outstanding exception from this is KIC 8006161, located at T_{eff} = 5258 K and a shift amplitude of 0.95 μHz. That KIC 8006161 shows outoftheordinary activity was also noticed by Karoff et al. (2013). An expansion of our study to all stars with Kepler short cadence data of sufficient length and a reduction of the uncertainty on the frequency shifts, for example, by peakbagging of the periodograms of all segments, will help to shed light on the question which of the two models is the correct one.
Figure 5 shows the frequency shift amplitudes as a function of the stellar ages, which are available in the literature, see Table A.1. As shown by, for example, Skumanich (1972) and recently by Salabert et al. (2016a), the strength of stellar activity is declining over time. Both scaling models for the cycle frequency shifts include the strength of the activity cycle via the quantity ΔR_{HK′}. It is the youngest star of the investigated sample, KIC 10644253, which shows the second greatest shift amplitude. Stars between 1.07–2.88 Gyr are missing from our sample. We excluded KIC 9025370 from this plot, as its age is only poorly determined. From Fig. 5, it can be seen that the amplitude of the frequency shifts is decreasing with stellar age. Again, KIC 8006161, at an age of 5.04 Gyr and a shift amplitude of 0.95 μHz, is the exception.
In Fig. 6 the shift amplitudes are plotted as a function of rotation period. Only 15 stars of the sample have measured rotation periods. As the rotation period is linked to the strength of stellar activity (e.g. Skumanich 1972; Wright et al. 2011), it is straightforward to assume that also the frequency shift amplitudes decrease with increasing rotation period. We note that once more KIC 8006161, found at P_{rot} = 29.79 days and a shift amplitude of 0.95 μHz, appears to deviate from the behaviour of the remaining sample. Indeed, even though the sample now consists of only 15 stars, there is a trend towards smaller shift amplitudes for longer rotation periods.
Fig. 7 Top panel: frequency shifts for KIC 8006161 for the frequency range between 2800–4400 μHz (top half) and height of the pmode envelope of KIC 8006161 (bottom half) as a function of time. Bottom panel: frequency shift as a function of height of the pmode envelope. 

Open with DEXTER 
Fig. 8 Frequency shifts of KIC 8006161 for three frequency ranges as a function of time: 2800–3300 μHz (left panel), 3300–3800 μHz (middle panel), and 3800–4400 μHz (right panel). 

Open with DEXTER 
3.3. KIC 8006161
KIC 8006161 is a solarlike star with a mass of 1.04 ± 0.02 M⊙, a radius of 0.947 ± 0.007 R⊙, and an age of 5.04 ± 0.17 Gyr (model values from Metcalfe et al. 2014). It is a G8V star (MolendaŻakowicz et al. 2013) and its rotation period is 29.79 ± 3.09 d (García et al. 2014a). This star shows the most significant and greatest frequency shifts of our sample with A/σ_{A} = 4.23. The variation of the pmode frequencies in the range between 2800–4400 μHz is shown in the top half of the top panel of Fig. 7. From the zero level at the beginning of the time series, the mode frequencies increase by ≈0.95 μHz until the end of the time series. The error bars on the frequency shifts, which result from our resampling approach described in Sect. 2.3, are only of the order of ≈0.1 μHz. This small error is due to the narrow width of the pmodes in the LSperiodogram of KIC 8006161 and the broad frequency range where modes can be detected. The shifts reach a local maximum at around day 350, followed by a brief decrease. The shifts then increase again towards the end of the time series. This behaviour is reminiscent of the double maximum of the Sun’s activity cycle. Unfortunately, the available time series does not cover a full activity cycle of KIC 8006161. Therefore, we can only limit the cycle length to P_{cyc}> 1147 d.
The variation of the mode heights are depicted in the bottom half of the top panel of Fig. 7. We find a clear anticorrelation between the shift of mode frequencies and the variation of mode heights with a correlation value of ρ = −0.93 with a pvalue of 4 × 10^{5}. This linear anticorrelation is clearly seen in the bottom panel of Fig. 7, where the frequency shifts are plotted as a function of height of the mode envelope. There is a slight departure from the linear relationship towards higher values of height, which correspond to a state of lower magnetic activity. This behaviour, albeit less pronounced, is also observed for KIC 5184732, see Fig. C.1. As discussed before, the other parameters of the background fit show no correlations or anticorrelations, this can be compared with Table B.11, where the only exception is a correlation between the two granulation timescales with ρ = 0.74 and p = 0.01. We therefore conclude that the observed variations in mode frequencies and mode heights are due to magnetic activity on KIC 8006161.
Karoff et al. (2013) found that the S index of KIC 8006161, a measure for stellar chromospheric activity, is comparable to that of the Sun. They also found that, at the same time, KIC 8006161 exhibits an unusually low excess flux. This flux is associated with magnetic sources on the stellar surface. They argue that due to this odd combination of activity indices, KIC 8006161 might be in a minimum state of activity. Their measurements were taken between 2010 and 2012, which corresponds to Kepler Quarters 4 through 15, and therefore largely overlap with the time series we used for our analyses, which encompasses Q517.2^{1}. In the study of García et al. (2014a), KIC 8006161 was found to be more active than the Sun. García et al. used the standard deviation of the time series, a measure of photospheric activity, to determine the activity level.
Since this star is very similar to the Sun regarding its fundamental parameters, our findings for the activity related shift of its pmode frequencies, and especially the strong anticorrelation between mode frequencies and heights of the pmode envelope, suggest that the magnetic activity on KIC 8006161 is somewhat stronger than on the Sun, following the results presented by García et al. (2014a).
For the Sun, it was shown that the amplitude of frequency shifts over the solar cycle are larger for higher frequencies (e.g. Libbrecht & Woodard 1990; JiménezReyes et al. 2001; Chaplin et al. 2007b; Salabert et al. 2015). We looked for a similar behaviour for KIC 8006161 by dividing the frequency range where pmodes are visible in the periodogram into three smaller sections: a low frequencies range between 2800–3300 μHz, an intermediate range with frequencies 3300–3800 μHz, and high frequency range between 3800–4400 μHz. The results are shown in the three panels of Fig. 8. As in the solar case, the frequency shifts increase with increasing mode frequency: the modes in the low frequency range vary by only ≈0.6 μHz, while the modes in the intermediate range shift ≈1.1 μHz. The modes in the high frequency range even shift by as much as ≈1.2 μHz. The measuring of the frequency shifts via the crosscorrelation and the estimation of the error bars on them are sensitive to the mode widths and amplitudes. Due to the fact that mode widths increase with frequency and mode amplitudes follow a Gaussian envelope centred around the frequency of maximum power, the estimated error bars on the frequency shifts are smallest in the intermediate range and largest in the high frequency range (cf. Fig. 8, Table 2, and Régulo et al. 2016) The overall course of the variation of the frequencies is similar for all three frequency ranges. The characteristics of the variation for the entire frequency range are dominated by the variation of the modes in the intermediate frequency range, since these modes have the largest amplitudes and thus contribute strongest to the crosscorrelation. As for the entire frequency range, the frequency shifts of these subranges are strongly anticorrelated to the variation of the mode heights with correlation coefficients of –0.92 (low frequency range), –0.93 (intermediate frequency range), and –0.87 (high frequency range). This is only the fourth star for which a frequency dependence of the activity related shift of pmode frequencies has been detected after the Sun, HD 49933 (Salabert et al. 2011), and KIC 10644253 (Salabert et al. 2016b).
3.4. KIC 10644253
This star has a mass of 1.13 ± 0.05 M⊙, a radius of 1.108 ± 0.016 R⊙ (Mathur et al. 2014), and an age of 1.07 ± 0.25 Gyr (model values from Metcalfe et al. 2014). It is a G0V star (MolendaŻakowicz et al. 2013) with a surface rotation rate of 10.91 ± 0.87 d (García et al. 2014a).
The frequency shifts we found for this star follow a cyclelike pattern, see top half of the top panel of Fig. 9, and have the second greatest peaktovalley amplitude in our sample of stars with a value of ≈0.87 μHz. The heights of the mode envelope, see bottom half of the top panel of Fig. 9, are strongly anticorrelated to the frequency shifts with a correlation coefficient of ρ = −0.93 and a p value of 3 × 10^{3}. This linear anticorrelation is seen in the bottom panel of Fig. 9, where the frequency shifts of KIC 10644253 are plotted as a function of height of the mode envelope.
In their study of 22 F stars Mathur et al. investigated the long cadence Kepler time series of KIC 10644253 regarding signs of magnetic activity and activity cycles. In their magnetic proxy for this star, they found indications of two activity cycles during the first 16 quarters of Kepler data. They use the scale average variance, which is a projection of the wavelet power spectrum in a small region around the stellar rotation period onto the time axis, as a magnetic proxy. The overall magnetic activity of this star was among the highest from their sample of stars. Mathur et al. spotted the strongest peak in their cycle proxy of KIC 10644253 during Quarter 10, which corresponds to a time during day ≈400–450 for our data set. This was recently confirmed by Salabert et al. (2016b) who find a positive shift in pmode frequencies and a large stellar magnetic variability in this time period. Indeed, in our analysis the frequency shifts are at a local maximum and the mode heights are at a local minimum during this period. The local minima found in the frequency shifts before and after this maximum are also reflected in the magnetic proxy of Mathur et al. (2014) and the frequency shifts and stellar magnetic variability of Salabert et al. (2016b).
Interestingly, we find strong positive frequency shifts and small mode heights, which indicates strong magnetic activity, during the last 250 d of the time series. This is not reflected in the proxy of Mathur et al. (2014) or in the stellar magnetic variability of Salabert et al. (2016b). During the same time period, there is a discrepancy between the mode frequency shifts and the photometric activity proxy presented by Salabert et al. (2016b).
There are different values for the inclination angle of KIC 10644253 in the literature. Mathur et al. (2014) state i = 43.44 ± 14.48 °, whereas Salabert et al. (2016b) find values of i ≈ 11 ° derived from spectroscopic measurements^{2}, derived from the stellar radius, rotation period and the asteroseismic vsini, and derived with another value of the spectroscopic vsini. Therefore, the inclination of the rotation axis of KIC 10644253 to the line of sight seems to be rather low, with a weighted average value of i = 23 ± 6 °.
Due to this low inclination, it is conceivable that the regions of high activity are largely confined to the nearly outofsight hemisphere of the star during the last 250 d of the time series, where the discrepancy between the activity proxies are found. The global, low degree modes we use for our analyses are susceptible to magnetic activity throughout the star. Therefore, the frequency shifts and variation of the mode heights might reveal signatures of magnetic activity which are hidden to photometry based techniques such as that of Mathur et al. (2014).
Fig. 9 Top panel: frequency shifts for KIC 10644253 for the frequency range between 2450–3350 μHz (top half) and height of the pmode envelope of KIC 10644253 as a function of time (bottom half). Bottom panel: frequency shift as a function of height of the pmode envelope. 

Open with DEXTER 
4. Summary and conclusion
We investigated the Kepler data of 24 solarlike stars for the temporal variation of pmode frequencies and height of the pmode envelope. For this, we split the stars’ time series into shorter segments. The shifts of pmode frequencies were measured by using a crosscorrelation analysis of the segments’ periodograms. The errors on the frequency shifts were estimated with a resampling approach. This approach will be further investigated and compared to similar approaches used in helio and asteroseismology in a future study. The temporal variation of the heights of the pmode envelope, the high frequency noise, and two granulation timescales were measured by fitting the segments’ periodograms.
We found significant variation of pmode frequencies (above 1σ), which can be signatures of stellar magnetic activity, during the observed period for 23 of the 24 stars of our sample. A cyclelike behaviour could be spotted in the frequency shifts of several stars (e.g. KIC 3632418, KIC 6933899, KIC 8760414, KIC 10644253). The correlation between the frequency shifts and the variation of mode heights was found to be strongly negative (below −0.5) for six stars (KIC 6933899, KIC 8006161, KIC 8760414, KIC 9955598, KIC 10644253, KIC 11244118, and KIC 12258514). Together, the significant shift in pmode frequencies and this anticorrelation strongly supports asteroseismic detection of solarlike magnetic activity on these stars. Due to the limited length of the time series, we refrain from assigning cycle periods to any of these stars. However, we can limit the cycle period of KIC 8006161 to P_{cyc}> 1147 d and state that more detailed and longer observations of the respective stars with evidence for cyclic activity could shed light on the period lengths of their magnetic activity cycles. Moreover, there are hints of sporadic activity on some of the investigated stars. The frequency shifts and mode heights of KIC 10018963 and KIC 11295426 exhibit significant temporal variations which might be caused by such noncyclic sporadic stellar magnetic activity. It is also noteworthy that the correlation of frequency shifts and mode heights is strongly negative, following the solar reference case, for only seven out of 24 stars, and even strongly positive for a couple of stars. Whether there is something to be learned from this about the connection of stellar magnetic activity or the dynamos at work in stars slightly different than the Sun to the temporal evolution of the mode frequency shifts and the mode heights is currently subject of further investigations.
The background parameters – specifically mode heights, granulation timescales, and high frequency noise level – and the frequency shifts do not show any systematic correlation across the sample of stars. We found very different correlations between the investigated parameters for the two stars, which we discussed in detail, KIC 8006161 and KIC 10644253, both of which show clear signs of magnetic activity. It will need an expanded sample of stars with a broader range of fundamental parameters, optimally all available Kepler stars with long enough short cadence coverage, to find the cause for this varying behaviour of the oscillation and background parameters.
The star KIC 10644253 exhibited shifts in pmode frequencies and changes of the mode heights which might be caused by two stellar cycles during the available data set, as was also recently discovered by Salabert et al. (2016b). When compared to the activity indices of Mathur et al. (2014) and Salabert et al. (2016b), we conclude that we found indications of stellar activity on the nearly outofsight hemisphere of this star during the last 250 d of Kepler observations.
For the solarlike star KIC 8006161 we were able to investigate the frequency dependence of the mode frequency shifts. We found that the shifts are larger for pmodes of higher frequencies, just like in the solar case (JimenezReyes et al. 1998), for HD 49933 (Salabert et al. 2011), and for KIC 10644253 (Salabert et al. 2016b). It is only the fourth star for which this frequency dependence of the activity related frequency shifts could be measured. KIC 8006161 is very similar to the Sun regarding its mass, age, and rotation period. It is therefore noteworthy that the frequency shifts are larger than for the Sun, even though the time series does clearly not cover a full cycle. This raises the question of what sets this star’s magnetic cycle apart from the Sun’s cycle. Further investigation of KIC 8006161 and its magnetic activity cycle is therefore worthwhile. When examining the frequency shifts as a function of mode height there is a departure from the negative linear relationship between these two quantities for higher values of mode height. This is also observed for KIC 5184732 and raises the question of whether this observed pattern is part of a hysteresis curve of the frequency shifts over a stellar cycle, similar as observed for the Sun (JimenezReyes et al. 1998).
We found that the frequency shift amplitude decreases with stellar age and rotation period. The results for our sample did not enable us to reject one of the two opposing scaling models for cycle frequency shift amplitudes by Chaplin et al. (2007a) and Metcalfe et al. (2007), while the data seems to show a slight tendency towards the scaling of Metcalfe et al. (2007).
Acknowledgments
The research leading to these results received funding from the European Research Council under the European Union’s Seventh Framework Program (FP/2007–2013)/ERC Grant Agreement No. 307117. G.R.D. acknowledges the support of the UK Science and Technology Facilities Council (STFC). This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France.
References
 Appourchaux, T., Chaplin, W. J., García, R. A., et al. 2012, A&A, 543, A54 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Baglin, A., Auvergne, M., Barge, P., et al. 2006, in The CoRoT Mission PreLaunch Status – Stellar Seismology and Planet Finding, eds. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, ESA SP, 1306, 33 (In the text)
 Baliunas, S. L., Donahue, R. A., Soon, W. H., et al. 1995, ApJ, 438, 269 [NASA ADS] [CrossRef] (In the text)
 Bartlett, M. S. 1978, An Introduction to Stochastic Processes: With Special Reference to Methods and Applications, 3rd edn. (Cambridge University Press) (In the text)
 Basu, S., Broomhall, A.M., Chaplin, W. J., & Elsworth, Y. 2012, ApJ, 758, 43 [NASA ADS] [CrossRef] (In the text)
 Beeck, B., Cameron, R. H., Reiners, A., & Schüssler, M. 2013, A&A, 558, A49 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [NASA ADS] [CrossRef] [PubMed] (In the text)
 Brandenburg, A., Saar, S. H., & Turpin, C. R. 1998, ApJ, 498, L51 [NASA ADS] [CrossRef] (In the text)
 Broomhall, A.M., Pugh, C., & Nakariakov, V. 2015, Adv. Space Res., 56, 2706 [NASA ADS] [CrossRef] (In the text)
 Bruntt, H., Basu, S., Smalley, B., et al. 2012, MNRAS, 423, 122 [NASA ADS] [CrossRef] (In the text)
 Chaplin, W. J., Elsworth, Y., Houdek, G., & New, R. 2007a, MNRAS, 377, 17 [NASA ADS] [CrossRef] (In the text)
 Chaplin, W. J., Elsworth, Y., Miller, B. A., Verner, G. A., & New, R. 2007b, ApJ, 659, 1749 [NASA ADS] [CrossRef] (In the text)
 Chaplin, W. J., Basu, S., Huber, D., et al. 2014, ApJS, 210, 1 [NASA ADS] [CrossRef] (In the text)
 Davies, G. R., Chaplin, W. J., Elsworth, Y., & Hale, S. J. 2014, MNRAS, 441, 3009 [NASA ADS] [CrossRef] (In the text)
 ForemanMackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [NASA ADS] [CrossRef] (In the text)
 García, R. A., Mathur, S., Salabert, D., et al. 2010, Science, 329, 1032 [NASA ADS] [CrossRef] [PubMed] (In the text)
 García, R. A., Hekker, S., Stello, D., et al. 2011, MNRAS, 414, L6 [NASA ADS] [CrossRef] (In the text)
 García, R. A., Ceillier, T., Salabert, D., et al. 2014a, A&A, 572, A34 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 García, R. A., Mathur, S., Pires, S., et al. 2014b, A&A, 568, A10 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Gilliland, R. L., Marcy, G. W., Rowe, J. F., et al. 2013, ApJ, 766, 40 [NASA ADS] [CrossRef] (In the text)
 Griffin, R. F. 2007, The Observatory, 127, 313 [NASA ADS] (In the text)
 Hall, J. C. 2008, Liv. Rev. Sol. Phys., 5, 2 (In the text)
 Harvey, J. 1985, in Future Missions in Solar, Heliospheric Space Plasma Physics, eds. E. Rolfe, & B. Battrick, ESA SP, 235 (In the text)
 Howell, S. B., Rowe, J. F., Bryson, S. T., et al. 2012, ApJ, 746, 123 [NASA ADS] [CrossRef] (In the text)
 Huber, D., Chaplin, W. J., ChristensenDalsgaard, J., et al. 2013, ApJ, 767, 127 [NASA ADS] [CrossRef] (In the text)
 Huber, D., Silva Aguirre, V., Matthews, J. M., et al. 2014, ApJS, 211, 2 [NASA ADS] [CrossRef] (In the text)
 JimenezReyes, S. J., Regulo, C., Palle, P. L., & Roca Cortes, T. 1998, A&A, 329, 1119 [NASA ADS] (In the text)
 JiménezReyes, S. J., Corbard, T., Pallé, P. L., Roca Cortés, T., & Tomczyk, S. 2001, A&A, 379, 622 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Kallinger, T., De Ridder, J., Hekker, S., et al. 2014, A&A, 570, A41 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Karoff, C., Metcalfe, T. S., Chaplin, W. J., et al. 2009, MNRAS, 399, 914 [NASA ADS] [CrossRef] (In the text)
 Karoff, C., Metcalfe, T. S., Chaplin, W. J., et al. 2013, MNRAS, 433, 3227 [NASA ADS] [CrossRef] (In the text)
 Koch, D. G., Borucki, W. J., Basri, G., et al. 2010, ApJ, 713, L79 [NASA ADS] [CrossRef] (In the text)
 Komm, R. W., Howe, R., & Hill, F. 2000, ApJ, 543, 472 [NASA ADS] [CrossRef] (In the text)
 Lefebvre, S., García, R. A., JiménezReyes, S. J., TurckChièze, S., & Mathur, S. 2008, A&A, 490, 1143 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Lépine, S., & Shara, M. M. 2005, AJ, 129, 1483 [NASA ADS] [CrossRef] (In the text)
 Libbrecht, K. G., & Woodard, M. F. 1990, Nature, 345, 779 [NASA ADS] [CrossRef] (In the text)
 Lomb, N. R. 1976, Ap&SS, 39, 447 [NASA ADS] [CrossRef] (In the text)
 Marcy, G. W., Isaacson, H., Howard, A. W., et al. 2014, ApJS, 210, 20 [NASA ADS] [CrossRef] (In the text)
 Mathur, S., Metcalfe, T. S., Woitaszek, M., et al. 2012, ApJ, 749, 152 [NASA ADS] [CrossRef] (In the text)
 Mathur, S., García, R. A., Ballot, J., et al. 2014, A&A, 562, A124 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 McQuillan, A., Mazeh, T., & Aigrain, S. 2013, ApJ, 775, L11 [NASA ADS] [CrossRef] (In the text)
 McQuillan, A., Mazeh, T., & Aigrain, S. 2014, ApJS, 211, 24 [NASA ADS] [CrossRef] [MathSciNet] (In the text)
 Metcalfe, T. S., Dziembowski, W. A., Judge, P. G., & Snow, M. 2007, MNRAS, 379, L16 [NASA ADS] [CrossRef] (In the text)
 Metcalfe, T. S., Creevey, O. L., Doğan, G., et al. 2014, ApJS, 214, 27 [NASA ADS] [CrossRef] (In the text)
 Meunier, N., Roudier, T., & Rieutord, M. 2008, A&A, 488, 1109 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 MolendaŻakowicz, J., Sousa, S. G., Frasca, A., et al. 2013, MNRAS, 434, 1422 [NASA ADS] [CrossRef] (In the text)
 Muller, R., Hanslmeier, A., & SaldañaMuñoz, M. 2007, A&A, 475, 717 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Pallé, P., Régulo, C., & Roca Cortés, T. 1990, in Progress of Seismology of the Sun and Stars, eds. Y. Osaki, & H. Shibahashi (Berlin Heidelberg: Springer), Lect. Not. Phys., 367, 129 (In the text)
 Pinsonneault, M. H., An, D., MolendaŻakowicz, J., et al. 2012, ApJS, 199, 30 [NASA ADS] [CrossRef] (In the text)
 Pires, S., Mathur, S., García, R. A., et al. 2015, A&A, 574, A18 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Press, W. H. 2007, Numerical recipes 3rd edn.: The art of scientific computing (Cambridge University Press) (In the text)
 Régulo, C., García, R. A., & Ballot, J. 2016, A&A, 589, A103 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Roth, M., & Stix, M. 2003, A&A, 405, 779 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Saar, S. H., & Brandenburg, A. 1999, ApJ, 524, 295 [NASA ADS] [CrossRef] (In the text)
 Salabert, D., Régulo, C., Ballot, J., García, R. A., & Mathur, S. 2011, A&A, 530, A127 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Salabert, D., García, R. A., & TurckChièze, S. 2015, A&A, 578, A137 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Salabert, D., Garcia, R. A., Beck, P. G., et al. 2016a, A&A, 596, A31 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Salabert, D., Régulo, C., García, R. A., et al. 2016b, A&A, 589, A118 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Scargle, J. D. 1982, ApJ, 263, 835 [NASA ADS] [CrossRef] (In the text)
 Silva Aguirre, V., Davies, G. R., Basu, S., et al. 2015, MNRAS, 452, 2127 [NASA ADS] [CrossRef] (In the text)
 Skumanich, A. 1972, ApJ, 171, 565 [NASA ADS] [CrossRef] (In the text)
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Vida, K., Oláh, K., & Szabó, R. 2014, MNRAS, 441, 2744 [NASA ADS] [CrossRef] (In the text)
 Wilson, O. C. 1978, ApJ, 226, 379 [NASA ADS] [CrossRef] (In the text)
 Woodard, M. F., & Noyes, R. W. 1985, Nature, 318, 449 [NASA ADS] [CrossRef] (In the text)
 Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48 [NASA ADS] [CrossRef] (In the text)
Appendix A: Tables of stellar parameters
Radius, mass, and age of the investigated stars.
Spectral type, effective temperature, and rotation period of the investigated stars.
Appendix B: Plots of frequency shifts and variations of the height of the pmode envelope, correlationcoefficients of the background model parameters and the measured frequency shifts
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 3632418.
Fig. B.1 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 3632418 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 3656476.
Fig. B.2 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 3656476 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 4914923.
Fig. B.3 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 4914923 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 5184732.
Fig. B.4 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 5184732 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 6106415.
Fig. B.5 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 6106415 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 6116048.
Fig. B.6 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 6116048 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 6603624.
Fig. B.7 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 6603624 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 6933899.
Fig. B.8 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 6933899 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 7680114.
Fig. B.9 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 7680114 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 7976303.
Fig. B.10 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 7976303 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 8006161.
Fig. B.11 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 8006161 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 8228742.
Fig. B.12 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 8228742 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 8379927.
Fig. B.13 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 8379927 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 8760414.
Fig. B.14 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 8760414 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 9025370.
Fig. B.15 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 9025370 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 9955598.
Fig. B.16 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 9955598 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 10018963.
Fig. B.17 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 10018963 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 10516096.
Fig. B.18 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 10516096 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 10644253.
Fig. B.19 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 10644253 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 10963065.
Fig. B.20 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 10963065 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 11244118.
Fig. B.21 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 11244118 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 11295426.
Fig. B.22 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 11295426 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 12009504.
Fig. B.23 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 12009504 as a function of time. 

Open with DEXTER 
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 12258514
Fig. B.24 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 12258514 as a function of time. 

Open with DEXTER 
Appendix C: Plot of frequency shifts as a function of height of the pmode envelope for KIC 5184732
Fig. C.1 Frequency shifts of KIC 5184732 as a function of height of the mode envelope. The correlation between the presented quantities is negative with a tendency towards no correlation for larger heights. 

Open with DEXTER 
All Tables
Overview of the basic data of the investigated time series and their periodograms.
Correlation coefficients of frequency shifts and pmode envelope height variations and their pvalues
Spectral type, effective temperature, and rotation period of the investigated stars.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 3632418.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 3656476.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 4914923.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 5184732.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 6106415.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 6116048.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 6603624.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 6933899.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 7680114.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 7976303.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 8006161.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 8228742.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 8379927.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 8760414.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 9025370.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 9955598.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 10018963.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 10516096.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 10644253.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 10963065.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 11244118.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 11295426.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 12009504.
Correlation coefficients between parameters of the background model and the frequency shifts for KIC 12258514
All Figures
Fig. 1 Section of the LSperiodograms of the first and nineteenth segment of KIC 8006161 in black and red colour, respectively. The LSperiodograms are boxcar smoothed over 0.81 μHz. 

Open with DEXTER  
In the text 
Fig. 2 Top panel: crosscorrelation between the part of interest of the LSperiodograms of two realisations of the first segment of KIC 8006161 (black). The Lorentzian fit to the crosscorrelation is shown in red. The dotted black vertical line marks zero frequency lag, the dashed red line marks the centre of the fitted Lorentzian. Bottom panel: same as in the top panel, but for the crosscorrelation between a realisation of the first and one of the nineteenth segment. 

Open with DEXTER  
In the text 
Fig. 3 Frequency shifts of solar pmodes from BiSON data with 1σ error bars (black). The seven month boxcar smoothed monthly sunspot number is shown in red. Source of the sunspot number: WDCSILSO, Royal Observatory of Belgium, Brussels. 

Open with DEXTER  
In the text 
Fig. 4 Measured shift amplitudes as a function of effective temperature. Colour and symbol coding: black diamonds represent stars younger than 4 Gyr, red triangles are stars between 4–5 Gyr, orange circles between 5–6 Gyr, blue squares between 6–7 Gyr, and purple asterisks are stars older than 7 Gyr. 

Open with DEXTER  
In the text 
Fig. 5 Measured shift amplitudes as a function of stellar age. 

Open with DEXTER  
In the text 
Fig. 6 Measured shift amplitudes as a function of rotation period. 

Open with DEXTER  
In the text 
Fig. 7 Top panel: frequency shifts for KIC 8006161 for the frequency range between 2800–4400 μHz (top half) and height of the pmode envelope of KIC 8006161 (bottom half) as a function of time. Bottom panel: frequency shift as a function of height of the pmode envelope. 

Open with DEXTER  
In the text 
Fig. 8 Frequency shifts of KIC 8006161 for three frequency ranges as a function of time: 2800–3300 μHz (left panel), 3300–3800 μHz (middle panel), and 3800–4400 μHz (right panel). 

Open with DEXTER  
In the text 
Fig. 9 Top panel: frequency shifts for KIC 10644253 for the frequency range between 2450–3350 μHz (top half) and height of the pmode envelope of KIC 10644253 as a function of time (bottom half). Bottom panel: frequency shift as a function of height of the pmode envelope. 

Open with DEXTER  
In the text 
Fig. B.1 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 3632418 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.2 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 3656476 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.3 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 4914923 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.4 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 5184732 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.5 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 6106415 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.6 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 6116048 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.7 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 6603624 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.8 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 6933899 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.9 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 7680114 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.10 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 7976303 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.11 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 8006161 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.12 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 8228742 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.13 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 8379927 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.14 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 8760414 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.15 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 9025370 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.16 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 9955598 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.17 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 10018963 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.18 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 10516096 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.19 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 10644253 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.20 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 10963065 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.21 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 11244118 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.22 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 11295426 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.23 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 12009504 as a function of time. 

Open with DEXTER  
In the text 
Fig. B.24 Frequency shifts (top half) and height of the pmode envelope (bottom half) for KIC 12258514 as a function of time. 

Open with DEXTER  
In the text 
Fig. C.1 Frequency shifts of KIC 5184732 as a function of height of the mode envelope. The correlation between the presented quantities is negative with a tendency towards no correlation for larger heights. 

Open with DEXTER  
In the text 