arxiv: v1 [astro-ph.im] 3 Jul 2014

Transcription

1 Astronomy & Astrophysics manuscript no. Bravo_et_al_214 c ESO 218 February 2, 218 Wavelets: a powerful tool for studying rotation, activity, and pulsation in Kepler and CoRoT stellar light curves J. P. Bravo, S. Roque, R. Estrela, I. C. Leão, and J. R. De Medeiros Universidade Federal do Rio Grande do Norte, Departamento de Física, Natal, RN, Brazil jbravo@dfte.ufrn.br Received date / Accepted date arxiv: v1 [astro-ph.im] 3 Jul 214 ABSTRACT Aims. The wavelet transform has been used as a powerful tool for treating several problems in astrophysics. In this work, we show that the time frequency analysis of stellar light curves using the wavelet transform is a practical tool for identifying rotation, magnetic activity, and pulsation signatures. We present the wavelet spectral composition and multiscale variations of the time series for four classes of stars: targets dominated by magnetic activity, stars with transiting planets, those with binary transits, and pulsating stars. Methods. We applied the Morlet wavelet (6th order), which offers high time and frequency resolution. By applying the wavelet transform to the signal, we obtain the wavelet local and global power spectra. The first is interpreted as energy distribution of the signal in time frequency space, and the second is obtained by time integration of the local map. Results. Since the wavelet transform is a useful mathematical tool for nonstationary signals, this technique applied to Kepler and CoRoT light curves allows us to clearly identify particular signatures for different phenomena. In particular, patterns were identified for the temporal evolution of the rotation period and other periodicity due to active regions affecting these light curves. In addition, a beat-pattern signature in the local wavelet map of pulsating stars over the entire time span was also detected. Key words. methods: data analysis techniques: photometric stars: rotation stars: activity stars: oscillations (including pulsations) binaries: eclipsing 1. Introduction The CoRoT (Baglin et al. 29) and Kepler (Borucki et al. 21) space missions produced unique sets of light curves for about 3, stars, with excellent time sampling and unprecedented photometric precision. These data, in addition to the major scientific goals of the missions (asteroseismology and the search for exoplanets) open new perspectives for studying different stellar properties, including rotation, magnetic activity, and binarity. For extracting information from raw signals, several mathematical transformations can be applied, such as Laplace transform (Widder 1945), Z-transform (Jury 1964), Wigner distributions (Boashash 1988), and Fourier transform (Bochner & Chandrasekharan 1949), the last being the most widely used. The wavelet transform (Torrence & Compo 1998) is a more recent tool applied for treating a large number of phenomena in different areas, including geophysics, atmospheric turbulence, health (cardiology), and astrophysics. This transformation has a major advantage, since it allows analysis of frequency variations in time of a given signal. Analogous to sunspots and solar photospheric faculae, whose visibility is modulated by stellar rotation, stellar active regions consist of cool spots and bright faculae caused by the magnetic field of the star. Such starspots are well established as tracers of rotation, but their dynamic behavior may also be used to analyze other relevant phenomena, such as magnetic activity and cycles (e.g., Mathur et al. (214); García et al. (29); Willson & Mordvinov (1999)). The present work provides a time frequency analysis using the wavelet transform of different sort of stellar light curves from CoRoT and Kepler space missions in order to identify particular features associated to rotation, magnetic activity, and pulsation 1. This procedure allows us to obtain a distribution of the signal s energy, in time scale space, from which we can identify the temporal evolution of different phenomena affecting the light curves (such as active regions and possible beats related to pulsations or surface differential rotation). This paper is organized as follows. Section 2 discusses procedures and methods, with an analysis of artificial stationary/nonstationary signals and a description of the wavelet transform. Sect. 3 presents the primary results, including the main characteristics of the stars studied here, with conclusions in Sect Procedures and methods The light curve obtained from a star can be decomposed into a number of frequencies represented in the power, which allows us to determine the periodic components of data that may be related to the physical properties of the system. These properties may be, for instance, rotational modulation and several related dynamic phenomena on the stellar surface, pulsation, as well as planetary transits. From the application of Fourier transform to the signal, we can obtain its frequency amplitude representation. Nevertheless, since the stellar light curves present events not occurring periodically, such as growth and decay of an active region, our interest, in addition to obtaining the spectral composition, which offers an idea of rotation behavior and pul- 1 The processing of CoRoT and Kepler light curves is carried out by an I.C.L. routine called Coroect, with methods described in De Medeiros et al. (213), and wavelet analysis by a J.P.B. routine, with methods described in this paper, both using the Interactive Data Language (IDL). Article number, page 1 of 1

2 sation modes, is to follow the time frequency behavior of those events and identify any specific signature to a particular stellar variability even if the light curve presents some kind of noise or singularities. For this, the time localization of the spectral components is necessary and the application of the wavelet transform to the signal will produce its time frequency representation (hereafter TFR). Amplitude Time (s) 2.1. Fourier transform The Fourier transform, named in honor of Joseph Fourier, is the extension of the Fourier series for nonperiodic functions; it decomposes any function into a sum of sinusoidal basis functions. Each of these functions is a complex exponential of a different frequency ν. The Fourier transform is defined as F(ν) = f (t) e 2πiνt dt. (1) The function F(ν) represents the amount of power inherent in f (t) at frequency ν, providing a frequency-based decomposition of the signal; that is, F(ν) tells us how much of each frequency exists in the signal, but offers no information on the existence of these frequencies over time. This information is not required when the signal is stationary. As an illustration, we simulated a stationary signal composed of 15-second sines (Fig. 1, top). This signal has different amplitudes and frequencies (1, 1, 5, and 2 Hz) at any given time instant. The power is obtained by applying the Fourier transform to this signal, as shown in Fig. 1 (bottom), where the four spectral components corresponding to frequencies 1, 1, 5, and 2 Hz are identified. Amplitude Time (s) Fig. 1. Artificial stationary signal composed of sines with different amplitudes and frequencies (1, 1, 5, and 2 Hz) (top panel) and its power (bottom panel). We also simulated a nonstationary signal with four different frequencies at three different time intervals, shown in Fig. 2 (top). In the first interval (up to 5.5 seconds (s)), the signal is composed of the four frequencies; in the second (from 5.5 s to 11 s), only two of the four frequencies are added (1 and 1 Hz); and in the last, the other two frequencies (2 and 5 Hz) make up the final part of the nonstationary signal. The corresponding power is shown in Fig. 2 (bottom). The two spectra in Figs. 1 and 2 are similar, in the sense that both show four spectral components at exactly the same frequencies (1, 1, 5, and 2 Hz). Nevertheless, although the first simulated signal contains these frequencies at all times, they are Fig. 2. Artificial nonstationary signal with four frequencies in three different time intervals (top panel) and its power (bottom panel). present in the second at different time intervals. This is one disadvantage of the Fourier transform, because it provides no information regarding the variation in these frequency components over time. However, globally it provides a more resolved power than the wavelet transform (Sect. 2.2), because it is useful for refining the period determination. On the other hand, the wavelet method allows a better interpretation of physical features (as shown below) prior to considering them for a period refinement. Gabor (1946) modified the Fourier transform, creating the so-called short-term Fourier transform (STFT) or Gabor transform. The mechanism consists of dividing the signal into small enough fixed-length segments, which can be assumed to be stationary. The function to be transformed (signal) is multiplied by a window function w, commonly a Gaussian function, and the resulting function is then processed with a Fourier transform to derive the TFR. Although the STFT has contributed significantly to the study of nonstationary signals, there was still a resolution problem to solve because the STFT does not show what frequency components exist at any given time. Indeed, we only know which frequency band exists at any given time interval (Hubbard 1996), which is a problem related to the width of the window function used. A wide window gives better frequency resolution but poor time resolution, while a narrow window has the opposite trade-off. This is interpreted as a limit on the simultaneous time frequency resolution one may achieve (Heisenberg uncertainty principle applied to time frequency information) The wavelet transform To overcome the resolution problem, the wavelet technique is a useful tool for analyzing nonstationary and nonperiodic signals, displaying characteristics that can vary in both time and frequency (or scale) (Burrus, Gopinath, & Guo 1998). The central idea of the wavelet is based on multiresolution analysis, from which the signal is analyzed at different frequencies with different resolutions showing details of the signal that characterize its different contributions. At low resolution, the details generally characterize large structures and, by increasing the resolution, we obtain more detailed information on the signal. In the 198s, Jean Morlet and Alex Grossman worked together on a mathematical function with two major characteris- Article number, page 2 of 1

3 Bravo et al.: wavelet procedure tics: having finite energy and subjected to dilation or compression (Grossmann & Morlet 1984). From a convolution between the wavelet and the signal, we can determine how much a section of the signal looks like the wavelet providing the TFR, also called wavelet local power or wavelet map (hereafter WVM). The analysis uses a function prototype, called mother wavelet ψ(t) that generates the other window functions. These functions are called daughter wavelets ψ a,b and are defined by translation and dilation (scale) of the mother wavelet ψ(t) as ψ a,b (t) = 1 ( ) t b ψ, a, b R, a, (2) a a where a and b are the scale and translation parameters, respectively. Scaling either dilates (large scales) or compresses (small scales) the signal. The constant number 1 a is an energy normalization factor so that the transformed signal will have the same energy at every scale (E signal = + f (t) 2 dt where f (t) is a continuous-time signal). There are two types of wavelet transform: the continuous wavelet transform (CWT) and the discrete wavelet transform (DWT) defined by CWT f (a, b) = f (t) ψ a,b dt = 1 ( ) t b f (t) ψ dt (3) a a and DWT f ( j, k) = a j 2 f (t) ψ ( a j t k b ) dt, (4) respectively, where j,k Z, a = a j, b = k b a j, and a >1 and b >1 are fixed (Foster 1996). The difference between Eqs. 3 and 4 is that the CWT operates on all possible scales and displacements, whereas the DWT uses a specific set of scales (or frequencies) and displacements (fixed values) (Daubechies 1992). In the present work, the CWT is used to achieve periods due to different phenomena, e.g., stellar rotation in some Kepler and CoRoT stellar light curves. The choice of the mother wavelet is imposed by the information that we want to emphasize in the signal. The most common continuous wavelets are the Morlet and the Mexican hat (Morettin 1999). The Morlet wavelet is a complex harmonic function contained within a Gaussian envelope as defined by Ψ(t) = e a[ν(t b)]2 e i2πν(t b) (5) where a and b are the scale and translation parameters, respectively, and ν is related to the order of the wavelet. The seconddegree exponential decay of the Gaussian function provides excellent spatial resolution, and its Fourier transform is a Gaussian distribution with very good frequency resolution. In this work, we use the sixth-order Morlet wavelet as the mother wavelet because of its good time localization and frequency resolution. To demonstrate the advantage of the wavelet technique, we applied this method to the artificial stationary and nonstationary signals of Figs. 1 and 2 respectively. Figure 3 shows the resulting WVMs (left panels) and their global wavelet spectra (GWS), which are obtained by time integration of the wavelet local power spectra (right panels). The horizontal and vertical axes correspond to the running time (in seconds) and logarithmic time scale (or the period in seconds (1/ν) where ν is the frequency in Hz), respectively. At the top, the WVM of the stationary signal shows the presence of four spectral components at all times, Period (s) Period (s) Artificial stationary signal Time (s) Artificial nonstationary signal >.96 s >.48 s > 9 s >.9 s >.96 s >.48 s > 9 s >.9 s Time (s). Fig. 3. Top panel: the wavelet map of the artificial stationary signal of Fig. 1. Bottom panel: the wavelet map of the artificial nonstationary signal of Fig. 2. The 6th-order Morlet wavelet was used. spectra are illustrated to the right. which was expected. At the bottom, in the case of the nonstationary signal, we can identify at which time each frequency is present or not. These spectral components are calculated via the GWS and illustrated by an arrow in Fig. 3. Therefore, from the WVM, we can identify which frequencies are predominant in the signal and at which instant they exist or not, resulting in improved accuracy on the timescale. Accordingly, using these WVMs via wavelet transform, we can identify several physical phenomena in stars from their light curves. For example, the technique allows us to determine the rotation period, to identify changes of active regions on the star due to growth or decay of spots and/or to differential rotation, as well as to analyze pulsation. 3. Results In this study, the wavelet method is applied to different Kepler and CoRoT public stellar light curves, including stars with planetary transit, binary systems, a variable star dominated by magnetic activity, and pulsating stars. Indeed, we present here the results of our analysis for a set of targets listed in different studies of variability reported in the literature, to compare our results with those produced by different procedures. First of all, we analyzed the CoRoT-2, which is a widely studied star, in order to better understand the surface phenomena behavior of this young, spotted yellow dwarf, and also Kepler-4, a Sun-like star presenting low changes in amplitude in its light curve allowing us to consider it as a quiet star. A Kepler apparently single star (KIC ) dominated by magnetic activity was Article number, page 3 of 1

4 also analyzed, in the interest of finding similar spot dynamics behavior as in CoRoT-2 with the transit removed 2. In addition, we applied the wavelet procedure to a Kepler eclipsing binary system, KIC , as well as to four pulsating variable stars, two targets observed by CoRoT (CoRoT and CoRoT ) and two by Kepler (KIC and KIC ). The light curves of the Kepler targets observed at a cadence of 3 minutes (with a mean total time span of 138 days) were reduced with the Pre-Search Data Conditioning (PDC) module of the Kepler data analysis pipeline, which tries to remove discontinuities, outliers, systematic trends, and other instrumental signatures (Twicken et al. 21). Those of the CoRoT stars are provided by the CoRoT public N2 data archive (Baudin et al. 26). TSI (W/m2) The Sun Because of its proximity, the Sun has become a standard model for studying stars. By analogy with the Sun and the solar magnetic cycles, active regions identified in other stars offer the possibility of studying stellar differential rotation, magnetic activity, dynamic of spots, and cycle variability. In this context, before dealing with our selected sample of stars, we briefly describe the results from the wavelet procedure applied to the total solar irradiance (TSI) time series from 1976 until 213, including cycles and the beginning of cycle 24, obtained from radiometers on different space platforms: HF on Nimbus7, ACRIM I on SMM, ACRIM II on UARS, and VIRGO on SOHO. The composite TSI time series, taken from the World Radiation Center, was expanded back to the minimum in 1976 using a model described in Fröhlich & Lean (24). Figure 4 shows the wavelet analysis for the Sun, where in each panel we present the time series at the top, its local map (modulus of the CWT and normalized to its maximum value) in the center, whose amplitudes are shown in terms of a color contour map, and the GWS (as the weighted average by time span) to the right. The WVM of the top panel exposes the 11-year cycle periodicity (or 384 days using our method), which is the most dominant feature of the even if in the GWS we see other periodicities and some subharmonics with lower power. Removing the long-term contributions, the intermediate- and short-term solar periodicities, as well as their changes over the entire time span, are clearly identified (bottom panel of Fig. 4). The dominant feature in the global is the 364-day periodicity, which is probably related to the 1.3-year periodicity at the base of the solar convection zone, as reported by Howe et al. (2), and also detected in sunspots areas and sunspots number time series studied using wavelet transforms by Krivova & Solanki (22), leading to an association of this period with an annual solar feature caused by magnetic fluxes generated deep inside the Sun. The other interesting features in the are the 158-day, 3-day, and 14-day periodicities. The 158-day flare occurrence period, called the Rieger-type period (Rieger et al. 1984), is stronger in cycle 21 but weaker for the subsequent cycles, and almost absent in the cycle 24. Krivova & Solanki (22) propose that Rieger-period is the third harmonic of the 1.3-year period. We also obtained the solar rotation period of 3 days, which is more evident in cycle 21 because of maximum activity but persists over the next three cycles. The identified 14- day periodicity seems to be a harmonic of the 3-day variation. Donnelly & Puga (199) demonstrate that such a cycle is also 2 The planetary and binary transits are removed using the I.C.L. routine based in (Drake 23) and (Tenenbaum et al. 21) methods. TSI (W/m2) Years Years. Fig. 4. Top panel: the composite TSI time series (at the top) and its wavelet local/global power spectra (in the center/to the right) considering long-term contributions. The most dominant feature of the is the 11-year cycle periodicity (P = 384 days). Bottom panel: the long-term contributions are removed. The intermediate- and short-term periodicities are clearly identified: 364 days (annual solar feature related to magnetic fluxes generated deep inside the Sun), 158 days (Rieger-type period), 3 days (solar rotation period) and 14 days (harmonic of the rotation period and associated to solar active regions). Contour levels are 9%, 8%, 7%,..., 2% and 1% of the map maximum. The contour levels are not plotted in the bottom panel for better viewing of periods. The 6th-order Morlet wavelet was used. associated to the active regions located opposite each other in solar longitude. The solar periodicities issued from the present analysis are in close agreement with those obtained by different authors, on the basis of different procedures for the treatment of the total solar irradiance (e.g., Willson (1997); Willson & Mordvinov (1999); Fröhlich & Lean (1998)) Stars with transiting planet One of the first planets detected with the CoRoT satellite, during its first long run in the galactic center direction (LRc1, time base 142 days), was CoRoT-Exo-2b, a hot Jupiter with a day orbit around a main-sequence G7V star. Because Article number, page 4 of 1

5 Bravo et al.: wavelet procedure its stellar mass, radius, and effective temperature are comparable to those of the Sun and because it is to the ZAMS (Bouchy et al. 28), which is possibly younger than.5 Gyr, CoRoT-2 (CoRoT , 2MASS ) has become a laboratory for our understanding of the magnetic activity behavior of the young Sun. The physical parameters of the star and the planet s characteristics were determined by Alonso et al. (28) and Bouchy et al. (28). Photometric analysis shows that modulation of the star was related to two active longitudes initially on opposite hemispheres, i.e., separated by 18. The first does not appreciably migrate, showing a rotation period of ±.24 days, and the other slowly migrates (retrograde migration) with a rotation period of d (Lanza et al. 29). Figure 5 shows the wavelet map with the associated global for this star. In the WVM of the top panel, the transits have a fine and deep droplet form. The periods related to the transit are indicated in the GWS by a red dashed line (.85 d and.56 d) in the top right hand panel of Fig. 5, which no longer appear in the GWS once the transit is removed (bottom right panel). Indeed, these are aliases of the planet orbital period. Because the periodicity associated to the magnetic activity is mixed with the energy contribution of the transits in the WVM, this hides the real orbital period and outstrips their aliases. The transits are removed because they can alter the periodogram when the orbital and the rotational periods (or their aliases) are synchronized or at least very close, preventing us from visualizing the persistence of the predominant periods. In this case we do not see any significant differences between both WVMs, but in some cases such as the binary system in section 3.4 which presents deeper transits, it is necessary to remove them. As seen in the WVM of the CoRoT-2 light curve, after removing the planet transit there is a clear signature showing the persistence over time of two semi-regular dune ranges (assemblages of color levels) representing the two predominant periods, which are calculated by a time integration of the local map and illustrated by a black dashed line in the GWS. Thus, we have 4.53 d as the rotation period and 2.27 d (approximately half of the main period) associated to spot emergence on opposite hemispheres of the star, most likely caused by differential rotation. This feature is therefore considered an indicator of rotational modulation related to the starspots. These periods are in accordance with the results obtained by Lanza s method (Lanza et al. 29) and also compared using the Lomb-Scargle method, which gives us a main period of d and a second of d, allowing us to adopt this type of signature in the WVM as a magnetic activity signature. Also calculated via Lomb-Scargle, there is another period of d, which could represent the variation in intensity of the second spot area or to be related to cyclic oscillation of the total spotted area, as reported by Lanza et al. (29). In contrast to the CoRoT-2 star, the second star with planetary transit analyzed in the present study, Kepler-4 (KIC , 2MASS ), a G-type star, is slightly brighter and is regarded in this work as a quiet star. Its planet Kepler-4b discovered in 21 has the size of Neptune and orbits its host star in 3.21 days (Borucki et al. 211). Figure 6 shows the wavelet map with the associated global for Kepler-4. We used here the Kepler Quarters 5-7 and Quarters 9-1, which yield 468-day time series once added. After removing the planetary transit (bottom panel), resulting in a cleaner wavelet map, we observe that the orbital period P orb = 3.1 d (obtained with our method) and its alias (red dashed line in the GWS of the top panel) no longer appear, thereby making three periodicities evident. Because the amplitude variations are very Fig. 5. Top panel: light curve of CoRoT-2 with transiting planet (at the top), its wavelet map (in the center) and global (to the right). Bottom panel: light curve of CoRoT-2 with transits removed (at the top), its wavelet map (in the center) and global (to the right). Contour levels are 9%, 8%, 7%,..., 2% and 1% of the map maximum. The 6th-order Morlet wavelet was used. small, we can hardly be certain which periodicity is related to rotation. A first guess is that rotation is associated with the day periodicity, and the others with their harmonics. In this case, no evident magnetic activity signature is identified, and the two dune ranges do not appear in the wavelet map for this quiet star Variable star dominated by magnetic activity Here we present an example of an apparently single active star, KIC (ra = 19 h 4 m 23.2 s, dec = , J2) in the search for a magnetic activity signature comparable to that of CoRoT-2. In fact, the light curve of this star shows significant variability features, which, in principle, could be associated to the pulsation or to rotational modulation caused by active regions. Its wavelet map with the associated global is displayed in Fig. 7, confirming the hypothesis that it is a fast rotator with many active regions, reflected by the semi-regular pattern observed in the light curve. Two main periods persist in the local map over the entire time span, both also evident in the Article number, page 5 of 1

6 Porb Fig. 7. Light curve (at the top) of Kepler star KIC , its local map (in the center), and its global wavelet (to the right). Contour levels are 9%, 8%, 7%,..., 2% and 1% of the map maximum. The 6th-order Morlet wavelet was used Fig. 6. Top panel: light curve of Kepler-4 with transiting planet (at the top), its wavelet map (in the center) and global (to the right). Bottom panel: light curve of Kepler-4 with transits removed (at the top), its wavelet map (in the center) and global (to the right). Contour levels are 9%, 8%, 7%,..., 2% and 1% of the map maximum. The 6th-order Morlet wavelet was used. GWS. The most significant, around 3.3 d, is related to the rotation, and the second, almost equal to half of the primary period, that is, 1.54 d, is associated to active regions that could be on opposite sides of the star, be growing and decaying, or be migrating, forming a double dip in the light curve (easily identified by visual inspection in some quarters as a local feature). One important aspect to underline here is that the period 1.54 d can be a potential indicator for two or more active regions contributing to the signal, with a period close to the main value of 3.3 d and suffering relative changes from one to another. From Lomb-Scargle in a prewhitening approach, Reinhold et al. (213) found two main periodicities for this case: P1 = 3.24 d and P2 = 3.57 d Binary system In this section we present the wavelet analysis for the Kepler binary system KIC (2MASS ), classified as an eclipsing binary by Prša et al. (211) and studied by Dimitrov et al. (212). Its wavelet map is displayed in Article number, page 6 of 1 Fig. 8 with the associated global. The orbital period Porb = d calculated via wavelet procedure and illustrated in the GWS (top panel) conforms with Porb = 18.6 d of Dimitrov et al. (212). Also, we observe some possible aliases (9.27 d and 4.63 d) of the transit period in both the WVM and the GWS (red dashed line). To search for stable periods, namely those that are persistent along the entire light curve, we removed the eclipses (bottom panel), whose depths are greater than the amplitudes of the rotational modulation contribution, distorting both the WVM and the GWS. We observe that the aliases are no longer present after the eclipses are removed. The regular changes in the light curve are represented by two semi-regular and continuous dune ranges over the 13-day window in the WVM, corresponding to the remaining 6.4-day and 3.2-day periodicities. The first is associated with rotational modulation caused by spots, in agreement with the rotation period computed by Dimitrov et al. (212), whereas the second is the second harmonic that may be caused by active regions located 18 apart on the stellar surface. As we have seen previously, this feature was also observed for the Sun and CoRoT-2 star. The light curve shows that some active regions emerge and fade during the entire coverage period with lower amplitude variation, which is characterized in the WVM by the dune ranges and their power index variations. Two other periodicities at d and d are also present in the WVMs, but it seems that both are caused by the recurrent gaps in the light curve Pulsating stars Changes in the luminosity of stars are also caused by fluctuations in stellar radius. This phenomenon is present in many intrinsic variable stars, such as RR Lyrae, Cepheids, and Delta Scuti, producing large and rather regular variations in amplitude in their light curves. Some variable stars, such as Gamma Doradus stars (γ Dor), are nonradial pulsators and have smaller pulsation amplitudes. Here we analyze two types of pulsating stars, CoRoT and KIC , which are typical examples of RR Lyrae stars, and CoRoT and KIC , presenting the typical behavior of γ Dor stars.

7 Bravo et al.: wavelet procedure PB Porb.9 P P/2 P/ Fig. 9. Light curve (at the top) of CoRoT star , an RRabtype Blazhko RR Lyrae, its local map (in the center) and global wavelet (to the right). The main periodicities detected are P =.56 d, associated to the radial fundamental pulsation period, P2 =.28 d and P = 8 d, the second and third harmonics, and PB = d, related 3 to the Blazhko period. Contour levels are 9%, 8%, 7%,..., 2% and 1% of the map maximum. The 6th-order Morlet wavelet was used Fig. 8. Top panel: light curve of KIC with binary transits (at the top) and its local/global wavelet map (in the center/to the right). Bottom panel: light curve of KIC with eclipses removed (at the top) and its local/global wavelet map (in the center/to the right). Contour levels are 9%, 8%, 7%,..., 2% and 1% of the map maximum. The 6th-order Morlet wavelet was used. The CoRoT star (ra = 18h 39m 3.8 s, dec = , J2), observed during the second long run in the galactic center direction (LRc2, time base 145 days), is a new RRab-type Blazhko RR Lyrae star (pulsation in the radial fundamental mode), analyzed by Guggenberger et al. (211). Their results are considered here as a comparison with our findings obtained via the wavelet procedure. The light curve and wavelet map with the associated global for the referred star are shown in Fig. 9. Its local map shows the long-term behavior of pulsation and its stability on low scales (high frequencies) of less than one day, represented by a track associated to the stronger power index. Some harmonics are illustrated by weaker power tracks. These periodicities are indicated in the global. Period P =.56 d (or frequency f = d 1 ) corresponds to the radial fundamental pulsation period, and the second and third harmonics are P2 =.28 d (or 2 f = d 1 ) and P3 = 8 d (or 3 f = d 1 ). Finally, PB = d (or fb =.3 d 1 ) is associated to the Blazhko modulation. We underline that the Blazhko effect is a variation in period and amplitude in RR Lyrae type variable stars (e.g., Szabó (214)). The CoRoT star exhibits clearly strong cycle-to-cycle changes in Blazhko modulation. In a continuous time span, 255 pulsations and more than 4 full Blazhko cycles were observed and investigated by Guggenberger et al. (211). These cycles are clearly observable in the local map displayed in Fig. 9 in color intensity and shape, forming a beat pattern with some circular and regular dune ranges. These periodicities are in accordance with those calculated by Guggenberger et al. (211) using Fourier analysis ( f = d 1, f1 = d 1, fb =.28 d 1, with fb the Blazhko frequency). Nevertheless when using our method, we do not find the additional period f1, considered as an independent mode by Guggenberger et al. (211), possibly owing to its very low amplitude. To assume this signature as typical of this type of pulsation, we also applied the wavelet analysis to the long-term light curve of the RR Lyrae KIC (2MASS ) or V36 Lyr (variable star designation in the GCVS Catalog (Samus et al. 29)) with a total time span of 1426 days. Figure 1 shows the corresponding WVM and GWS. The contour levels are not plotted here to avoid hiding the pulsation signature in the local map. The beat pattern of the previous RR Lyrae is still evident here; i.e., the dune ranges are circular and regular, comprising tracks associated to the primary period and the harmonics. This beat pattern characterizes the Blazhko cycles (27 full cycles) whose periodic variation could be associated to the 52.8-day periodicity (PB or fb =.19 d 1 ) in the GWS. The other periodicities are P =.54 d (or f = d 1 ) corresponding to the radial fundamental pulsation period, and the second and third harmonics P2 =.27 d (or 2 f = 3.74 d 1 ) and P 1 3 = 8 d (or 3 f = d ), respectively. For comparison, we find similar results to those obtained by Benko et al. (21) using Fourier analysis ( f = d 1 or P = d and PB = 51.4 d). The authors also find additional frequencies ( f1 = d 1 and f = d 1 ) that we do not obarticle number, page 7 of 1

8 P B P P P /2 P 3 P P 1 P 5 P / Fig. 1. Light curve (at the top) of Kepler RR Lyrae star KIC or V36 Lyr, its local map (in the center) and global wavelet (to the right). The main periodicities detected are P =.54 d, associated to the radial fundamental pulsation period, P =.27 d and 2 P = 8 d, the second and third harmonics, and P 3 B = 52.8 d, related to the Blazhko period. The 6th-order Morlet wavelet was used P tain using our method, possibly due to their small amplitudes. Clearly the wavelet pattern and signatures observed for CoRoT star , a well defined RRab-type Blazhko RR Lyrae type, as discussed in the previous paragraph, are also observed for KIC The CoRoT star , observed during the first scientific run anti-center pointing which lasted about 6 d (IRa1), is a 12.4 magnitude eclipsing binary (ra = 6 h 48 m 54.3 s, dec = , J2), which is considered to be a γ Dor pulsator that shows modulated oscillations and narrow eclipses. For this star, an orbital period of d (F orb = 139 d 1 ) is reported by Maceroni et al. (21, 213). These authors also present the frequencies obtained from a Fourier analysis and detected a nearly equidistant frequency spacing of about.5 d 1. Figure 11 depicts the wavelet analysis of CoRoT , considering the binary transits (top panel) and the one with the eclipses removed (bottom panel). There are no significant differences between both maps, with only a few variations in scale intensity caused by the eclipses. However, the pulsation frequency f 1 = 1.22 d 1 (P 1 =.82 d) is still significant in both WVMs. As shown by Maceroni et al. (21, 213), the primary star pulsates with typical γ Dor frequencies, a result compatible with our wavelet analysis. The main periodicities illustrated in the GWS of the bottom panel are P 1 =.82 d (or f 1 = 1.22 d 1 ), associated to the pulsation period with the highest amplitude, and P 2 = d ( f 2 =.5 d 1, corresponding to.5f orb ) related to the beat pattern. Also, P 3 = 4.34 d ( f 3 =.23 d 1 ) and P 4 = 2.17 d ( f 4 =.46 d 1 ) remain after removing the eclipses, leading to the conclusion that they must be related to the beat pattern, although they are not exactly equal to the harmonics of the P 2. Finally P 5 =.41 d ( f 5 = 2.44 d 1 ) is one of the harmonics of the pulsation period. All these periodicities are in accordance with those obtained by Maceroni et al. (21). The 1.38-day periodicity in the GWS of the top panel could be associated to the orbital period because it no longer appears once eclipses are removed, whereas another period of 1.31 days with low power appears in both local maps, assuming that it is also Fig. 11. Top panel: light curve of CoRoT (a γ Doradus pulsator) with binary transits (at the top) and its local/global wavelet map (in the center/to the right). Bottom panel: light curve of CoRoT with eclipses removed (at the top) and its local/global wavelet map (in the center/to the right). The main periodicities detected and illustrated in the GWS are P 1 =.82 d, corresponding to the nonradial fundamental pulsation period, P 2 = d, related to the beat pattern, P 3 = 4.34 d and P 4 = 2.17 d, also related to the beat pattern (although they are not exactly equal to the harmonics of the P 2 ), and P 5 =.41 d, an overtone of P 1. Contour levels are 9%, 8%, 7%,..., 2% and 1% of the map maximum. The 6th-order Morlet wavelet was used. related to the beat pattern of low amplitude. The pulsation signature that we observe here presents a semiregularity of dune ranges in the WVM, putting in doubt that the modulation variations are caused by pulsation or by rotation accompanied by the presence of spots. Indeed, the observed semiregularity could be the result of the short-term light curve, with a coverage time limited to 57 days, which seems to be very short for identifying evident pulsation signatures by just looking at the local map. Finally, we applied the wavelet procedure to the longterm light curve of the Kepler star KIC (2MASS ), classified as a γ Dor star by Tkachenko et al. (213). From the wavelet analysis, illustrated by the corresponding WVM and GWS given in Fig. 12, we observe one P 3 P 4 P 1 P 5 Article number, page 8 of 1

9 Bravo et al.: wavelet procedure Fig. 12. Light curve (at the top) of Kepler γ Dor KIC , its local map (in the center), and global wavelet (to the right). The main periodicities detected are P =.95 d, associated to the nonradial fundamental pulsation period and d, related to the beat pattern. The 6th-order Morlet wavelet was used. evident track showing a regular continuity of dune ranges associated to pulsation modes. We do not find the same dunes as the previous analyzed RR Lyrae cases because of the difference in the type of pulsations, but the observed regularity points to a clear pulsation pattern. The predominant periods are P =.95 d, the fundamental pulsation period and d, possibly associated to beat pattern, as shown in the GWS. The contour levels are not plotted here to avoid hiding the evidence of γ Dor pulsation signature. 4. Conclusions In the present work we have carried out a time frequency analysis using the wavelet transform of different sorts of stellar light curves obtained in the scope of the CoRoT and Kepler space missions. The procedure was applied to CoRoT-2 and Kepler-4, two well-known main-sequence stars with planetary transits, to a Kepler apparently single main-sequence star, KIC , which is dominated by magnetic activity, to a Kepler eclipsing binary system KIC , as well as to four pulsating variable stars, two RR Lyrae, CoRoT and KIC , and two γ Dor, CoRoT and KIC The procedure allowed us to obtain the distribution of the signal s energy, in time scale space, from which it was possible to identify the temporal evolution of different phenomena affecting the light curves (such as active regions and possible beats related to pulsations or surface differential rotation). Wavelet analysis in the solar case gives us a first idea about what to expect if we apply this method to other active stars. As dominant features, reported well by different studies (Howe et al. 2; Krivova & Solanki 22; Donnelly & Puga 199; Willson & Mordvinov 1999), we identified the 364-day periodicity probably related to an annual solar feature caused by magnetic fluxes generated deep inside the Sun, the solar rotation period, the Rieger-type period, and the 14-day periodicity associated to active regions located opposite each other in solar longitude. The 11-year cycle is also detected when considering the long-term contributions in the local map P From the wavelet analysis of the CoRoT-2 light curve, in addition to the orbital period, we identified in particular a main period corresponding to the rotation period and another that is nearly half of the primary period, which is associated to active regions at different longitudes evolving over time. In fact, instead of being considered as an harmonic of the rotation period, this second periodicity could be a presumed effect of active regions moving to the opposite side of the star, most probably due to differential rotation as in solar case (Donnelly & Puga 199). In the wavelet maps, we also distinguish two semi-regular and continuous dune ranges over the entire time span, which is a strong indicator of the dynamic of starspots yielding to assume these features as the more typical rotation and magnetic activity signature. In addition to these periodicities, a third period is sometimes evident (for example in the case of CoRoT-2) and related to long-term cycles of the stellar activity. However, we can notice that in some cases such a period can be hidden by other contributions in the local map owing to gaps in the light curve (as seen in Fig. 8). In contrast to CoRoT-2, with clear signatures of an active star, the wavelet analysis for Kepler-4 shows no evident signatures of rotation and magnetic activity, reflecting its quiet magnetic activity behavior, responsible for the observed light curve low amplitude variation. Then, we analyzed KIC to compare the identified wavelet signatures with those for CoRoT-2 without transits. In addition to the rotation periodicity, our analysis reveals also the presence of two or more active regions, pointing to a clear dynamic of starspots as in CoRoT-2. In addition to the confirmation of the orbital period already report in the literature, the wavelet analysis for the Kepler eclipsing binary KIC has also revealed that different periodicity signatures, including rotation, are better defined after removing the transits or eclipses. By comparing the wavelet analysis of the light curves with transits and with removed transits, for stars with planetary and binary transits, it is clear that the presence of periodicity signatures in the light curves are featured more after the transit has been removed, especially when their depths are greater than the amplitude of the rotational modulation. In the case of the pulsating stars CoRoT , V36 Lyr, CoRoT , and KIC , there are solid similarities between their wavelet maps, both RR Lyrae and γ Dor stars clearly showing the pulsation period with its harmonics and a beat pattern illustrated by continuous and regular dune ranges. The pulsation pattern is different between the two types of pulsating stars. The beat pattern for the RR Lyrae stars is represented by more circular and regular dune ranges, whereas the γ Dor stars are characterized by more compact but very regular dune ranges (or tracks). We also note that for a short total time span, the γ Dor pulsation signature could be confused by rotational modulation. Finally, to establish the observed regularity of dunes as a pulsation pattern in the referred pulsating stars, we extend our wavelet analysis to the following additional pulsating stars: KIC RR Lyrae star, KIC , KIC , KIC , KIC , KIC 18943, KIC γ Doradus stars, KIC RR-δ Scuti star, and KIC Cepheid star. All these are studies of their pulsating nature reported in the literature by other authors. The resulting wavelet maps confirm that the regularity of dunes in the maps is a major trace of a pulsation pattern. In summary, this study has shown that the wavelet technique offers a detailed interpretation of stellar light curves, giving additional information on different physical phenomena present in the signal. Semi-regular patterns represent changes of active regions due to growth or decay of spots and/or to differential rotation, whereas regular patterns indicate events that are more stable Article number, page 9 of 1

10 in time, like pulsations. This method has an advantage in relation to the Fourier technique (Lomb-Scargle used), because in addition to identifying transits or eclipses, it is possible to identify the signature of the dynamic of different star characteristics associated to the observed light curves. Acknowledgements. The CoRoT space mission was developed and is operated by the French space agency CNES, with the participation of ESA s RSSD and Science Programs, Austria, Belgium, Brazil, Germany, and Spain. This paper includes data collected by the Kepler mission. Funding for the Kepler mission is provided by the NASA Science Mission directorate. Some of the data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS Support for MAST for non-hst data is provided by the NASA Office of Space Science via grant NNX9AF8G and by other grants and contracts. Research activities at the Federal University of Rio Grande do Norte are supported by continuous grants from the national funding agencies CNPq, and FAPERN and by the INCT-INEspaço. J.P.B. acknowledges a CAPES graduate fellowship. S.R. and R.E. acknowledge CNPq undergraduate fellowships. I.C.L. acknowledges a CNPq/PNPD fellowship. Tenenbaum, P., Bryson, S. T., Chandrasekaran, H., et al. 21, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 774, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series Tkachenko, A., Aerts, C., Yakushechkin, A., et al. 213, A&A, 556, A52 Torrence, C. & Compo, G. P. 1998, Bulletin of the American Meteorological Society, 79, 61 Twicken, J. D., Chandrasekaran, H., Jenkins, J. M., et al. 21, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 774, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series Widder, D. V. 1945, The American Mathematical Monthly, 52, 419 Willson, R. C. 1997, Science, 277, 1963 Willson, R. C. & Mordvinov, A. V. 1999, Geophys. Res. Lett., 26, 3613 References Alonso, R., Auvergne, M., Baglin, A., et al. 28, A&A, 482, L21 Baglin, A., Auvergne, M., Barge, P., et al. 29, in IAU Symposium, Vol. 253, IAU Symposium, ed. F. Pont, D. Sasselov, & M. J. Holman, Baudin, F., Baglin, A., Orcesi, J.-L., et al. 26, in ESA Special Publication, Vol. 136, ESA Special Publication, ed. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, 145 Benkő, J. M., Kolenberg, K., Szabó, R., et al. 21, MNRAS, 49, 1585 Boashash, B. 1988, IEEE Transactions on Acoustics, Speech, and Signal Processing, 36, 1518 Bochner, S. & Chandrasekharan, K. 1949, Fourier Transforms (Princeton University Press) Borucki, W. J., Koch, D., Basri, G., et al. 21, Science, 327, 977 Borucki, W. J., Koch, D. G., Basri, G., et al. 211, ApJ, 736, 19 Bouchy, F., Queloz, D., Deleuil, M., et al. 28, A&A, 482, L25 Burrus, C. S., Gopinath, R., & Guo, H. 1998, Introduction to Wavelets and Wavelet Transforms (Prentice-Hall) Daubechies, I., ed. 1992, Ten lectures on wavelets De Medeiros, J. R., Ferreira Lopes, C. E., Leão, I. C., et al. 213, A&A, 555, A63 Dimitrov, D., Kyurkchieva, D., & Radeva, V. 212, Bulgarian Astronomical Journal, 18, 3 Donnelly, R. F. & Puga, L. C. 199, Sol. Phys., 13, 369 Drake, A. J. 23, ApJ, 589, 12 Foster, G. 1996, The Astronomical Journal, 112, 179 Fröhlich, C. & Lean, J. 1998, Geophys. Res. Lett., 25, 4377 Fröhlich, C. & Lean, J. 24, A&A Rev., 12, 273 Gabor, D. 1946, J. Inst. Electr. Engineering, 93, 429 García, R. A., Régulo, C., Samadi, R., et al. 29, A&A, 56, 41 Grossmann, A. & Morlet, J. 1984, SIAM Journal Mathematical Analysis, 15, 723 Guggenberger, E., Kolenberg, K., Chapellier, E., et al. 211, MNRAS, 415, 1577 Howe, R., Christensen-Dalsgaard, J., Hill, F., et al. 2, Science, 287, 2456 Hubbard, B. 1996, The world according to wavelets. The Story of a Mathematical Technique in the Making (A. K. Peters Ltd, Wellesley, Massachusetts) Jury, E. 1964, Theory and application of the z-transform method (Robert E. Krieger Publishing) Krivova, N. A. & Solanki, S. K. 22, A&A, 394, 71 Lanza, A. F., Pagano, I., Leto, G., et al. 29, A&A, 493, 193 Maceroni, C., Cardini, D., Damiani, C., et al. 21, ArXiv e-prints Maceroni, C., Montalbán, J., Gandolfi, D., Pavlovski, K., & Rainer, M. 213, A&A, 552, A6 Mathur, S., García, R. A., Ballot, J., et al. 214, A&A, 562, A124 Morettin, P. A. 1999, Ondas e Ondeletas: da análise de Fourier à análise de ondeletas (EdUSP) Prša, A., Batalha, N., Slawson, R. W., et al. 211, AJ, 141, 83 Reinhold, T., Reiners, A., & Basri, G. 213, ArXiv e-prints Rieger, E., Kanbach, G., Reppin, C., et al. 1984, Nature, 312, 623 Samus, N. N., Durlevich, O. V., & et al. 29, VizieR Online Data Catalog, 1, 225 Szabó, R. 214, in IAU Symposium, Vol. 31, IAU Symposium, ed. J. A. Guzik, W. J. Chaplin, G. Handler, & A. Pigulski, Article number, page 1 of 1