# Transit timing variations in WASP-10b induced by stellar activity?

## Abstract

The hot-Jupiter WASP-10b was reported by Maciejewski et al. (2011a); Maciejewski et al. (2011b) to show transit timing variations (TTV) with an amplitude of minutes. These authors proposed that the observed TTVs were caused by a 0.1 perturbing companion with an orbital period of d, and hence, close to the outer 5:3 mean motion resonance with WASP-10b. To test this scenario, we present eight new transit light curves of WASP-10b obtained with the Faulkes Telescope North and the Liverpool Telescope. The new light curves, together with 22 previously published ones, were modelled with a Markov-Chain Monte-Carlo transit fitting code.

Transit depth differences reported for WASP-10b are thought to be due to star spot induced brightness modulation of the host star. Assuming the star is brighter at the activity minimum, we favour a small planetary radius. We find in agreement with Johnson et al. (2009) and Maciejewski et al. (2011b). Maciejewski et al. (2011a) and Husnoo et al. (2012) find no evidence for a significant eccentricity in this system. We present consistent system parameters for a circular orbit and refine the orbital ephemeris of WASP-10b.

Our homogeneously derived transit times do not support the previous claimed TTV signal, which was strongly dependent on 2 previously published transits that have been incorrectly normalised. Nevertheless, a linear ephemeris is not a statistically good fit to the transit times of WASP-10b. We show that the observed transit time variations are due to spot occultation features or systematics. We discuss and exemplify the effects of occultation spot features in the measured transit times and show that despite spot occultation during egress and ingress being difficult to distinguish in the transit light curves, they have a significant effect in the measured transit times. We conclude that if we account for spot features, the transit times of WASP-10 are consistent with a linear ephemeris with the exception of one transit (epoch 143) which is a partial transit. Therefore, there is currently no evidence for the existence of a companion to WASP-10b. Our results support the lack of TTVs of hot-Jupiters reported for the Kepler sample.

###### keywords:

stars: planetary systems – stars: individual (WASP-10) –stars: spots –techniques: photometric^{1}

^{2}

## 1 Introduction

The discovery of WASP-10b, a hot-Jupiter planet in a 3.09 day orbit, was reported by Christian et al. (2009). Its host star is a K5V type with K, [M/H] (Christian et al., 2009) and is relatively young, with an age of Myr (Maciejewski et al., 2011a). Using a high signal-to-noise transit light curve, Johnson et al. (2009) updated the stellar density and re-derived the stellar mass to and radius to . There has been some discussion about the planetary size, with Christian et al. (2009), Dittmann et al. (2010) and Krejčová et al. (2010) obtaining a planetary radius of while Johnson et al. (2009) derived a smaller radius of . Maciejewski et al. (2011b) argue that the discrepancy in previous results was due star spot induced brightness variability. In fact, WASP-10 has been shown to have stellar variability due to the rotation modulation of spots with a period of d (Smith et al., 2009) and semi-amplitude up to mmag. Star spots reduce the effective stellar disc area leading to an underestimation of the stellar radius and a larger transit depth. This variation of transit depth with spot coverage explains the differences in the measured planetary radius (e.g. Czesla et al. 2009). Using a transit taken close to maximum brightness that corresponds to a minimum spot coverage, Maciejewski et al. (2011b) obtained a maximum planetary radius of .

Interestingly, Maciejewski et al. (2011a) reported a periodic modulation of transit timing variations (TTVs) for WASP-10 with an amplitude of minutes. They suggested a 0.1 planet companion in the outer 5:3 mean motion resonant orbit as the cause of the measured TTVs. Later, Maciejewski et al. (2011b) confirmed the TTV solution and showed that only spots located close to the stellar limb would affect the transit times but they would also result in transit duration variations. Since no transit duration variations were found in their sample, the authors argued that spot occultation features did not significantly affect the measured transit times of WASP-10.

The possibility of finding additional planets by measuring their effect on the central transit times of known transiting planet systems was proposed by Holman & Murray (2005) and Agol et al. (2005). They showed that for systems near mean-motion resonances this method can be sensitive to planets less massive than the Earth. This has led many groups to search for TTVs from the ground. However, most of the searches resulted in only upper limits on the mass of possible companions (e. g. Miller-Ricci et al. 2008; Bean 2009; Gibson et al. 2009, 2010). The few initial reports of TTVs were unconfirmed later. For example, Díaz et al. (2008) found TTVs in OGLE-111b which were unconfirmed in later studies (Adams et al., 2010; Hoyer et al., 2011). At present, only three hot-Jupiter transiting systems show TTVs that were measured from ground-based data: WASP-3 (Maciejewski et al., 2010), WASP-10 (Maciejewski et al., 2011a), and WASP-5 (Fukui et al., 2011). Meanwhile, Kepler data has allowed the discovery of many TTV systems (e.g. Kepler-9 (Holman et al., 2010) and Kepler-11 (Lissauer et al., 2011) ). However, a recent TTV study of the Kepler sample revealed the lack of TTVs in hot-Jupiters, contrasting with the TTVs found for longer-period Jupiters or hot-Neptunes (Steffen et al., 2012). The authors suggested that this implied a different formation and/or evolution history for hot-Jupiters. Therefore, it is very important to confirm and understand the hot-Jupiter TTV detections that have been reported from ground-based data.

To confirm and possibly characterise WASP-10b’s companion we obtained new high precision transit light curves of WASP-10b. Six transit observations of WASP-10b were obtained with the RISE fast CCD camera mounted on the Liverpool Telescope and two more with the Merope camera on the Faulkes Telescope North. We combined the new observations with previously published 22 light curves of WASP-10 to homogeneously derive the transit times.

The new observations are described in Section 2. In Section 3, we discuss our transit model. We present the updated parameters of the system in Section 4 and in Section 5 discuss our results and show possible causes of spurious TTVs. Finally, we present our conclusions in Section 6.

## 2 Observations

We present further high precision transit observations of WASP-10b taken with Faulkes Telescope North (FTN) and with the Liverpool Telescope (LT).

Two full transits of WASP-10b were observed on the 2008-09-15 and 2008-09-19 with the Merope camera on the FTN using the SDSS-iâ² filter. The Merope camera consists of a e2V, pixel-CCD and has a field of view of 4.7’ 4.7’. Aperture photometry was performed with the DAOPHOT package (Stetson, 1987) within the IRAF4 environment. The differential photometry was performed relative to at least 5 comparison stars that were within the FTN field-of-view.

Another six transits were observed with the fast CCD camera RISE (Steele et al., 2008; Gibson et al., 2008) mounted on the 2.0m Liverpool Telescope on La Palma, Canary Islands. RISE has a wideband filter covering –nm which corresponds approximately to V+R. The pixel scale is 0.54 arcsec/pixel resulting in a 9.4’ 9.4’ field-of-view. Thanks to its frame transfer CCD, RISE has a deadtime of only 35 ms for exposures longer than 1 second. To decrease systematic noise due to poor guiding (Barros et al., 2011) the telescope was defocussed to -0.6mm and we used an exposure time of 9 seconds. However, in one of the observations (2010-11-10) we experimented with a higher defocussing of -1.0mm and longer exposure time of 34 seconds. The RISE data were reduced using the ultracam pipeline (Dhillon et al., 2007) following the same procedure as for WASP-21b (Barros et al., 2011). Details of the observations are given in Table 1.

The final new high precision transit light curves for WASP-10b are shown in Figure 1 along with the best-fit model described in Section 3.3. We overplot the model residuals and the estimated uncertainties, which are discussed in Section 3.2.

Epoch | Date | telescope | exp time | number exposures | aperture radius |
---|---|---|---|---|---|

sec | arcsec | ||||

20 | 2008-09-15 | FTN | 125 | 115 | 6.1 |

21 | 2008-09-19 | FTN | 125 | 109 | 7.0 |

26 | 2008-10-04 | RISE | 5 | 1580 | 8.6 |

252 | 2010-09-03 | RISE | 9 | 1195 | 8.1 |

254b | 2010-09-09 | RISE | 9 | 1233 | 8.6 |

265 | 2010-10-13 | RISE | 9 | 1063 | 10.3 |

274 | 2010-11-10 | RISE | 34 | 339 | 11.4 |

285 | 2010-12-14 | RISE | 9 | 755 | 9.2 |

Previously published WASP-10b light curves were included in our analysis which were kindly made available by the authors. Details of these observations are given in Table 2. For completeness we also show the new best-fit model overploted on the previously published transit light curves in Figure 2, whose references for these data are given in Table 2.

Hereafter, we refer to each light curve by their epoch number relative to the ephemeris presented in Section 4.2. The epoch number is also given in the observation logs and light curve plots.

Epoch | Date | telescope | exp time | num exp | filter | ref |
---|---|---|---|---|---|---|

sec | arcsec | |||||

-103 | 2007-09-01 | Mercator | 30 | 170 | V | A |

-89 | 2007-10-15 | Tenagra | 125 | 125 | I | A |

-88 | 2007-10-18 | Tenagra | 125 | 91 | I | A |

0 | 2008-07-16 | 2.2m Hawaii | 50 | 194 | sloan z’ | B |

36 | 2008-11-04 | Slovak | 120 | R | D | |

56 | 2009-01-05 | 0.6m Torun | 50 | 254 | R | E |

123a | 2009-07-31 | 0.6m Jena | 70 | 144 | R | E |

123b | 2009-07-31 | 0.6m Rozhen | 120 | 116 | R | E |

124 | 2009-08-03 | 0.6m Rozhen | 90 | 122 | R | E |

132 | 2009-08-28 | 0.6m Rozhen | 90 | 99 | R | E |

133a | 2009-08-31 | Slovak | 165 | R | D | |

133b | 2009-08-31 | 0.6m Jena | 55 | 203 | R | E |

134 | 2009-09-03 | 0.6m Rozhen | 90 | 105 | R | E |

143 | 2009-10-01 | 2.0m Rozhen | 15 | 548 | V | E |

144 | 2009-10-04 | Slovak | 118 | R | D | |

145 | 2009-10-07 | Slovak | 114 | R | D | |

147 | 2009-10-14 | Kuiper | 10 | 579 | I | C |

158 | 2009-11-17 | Tenagra | 76 | 76 | R | E |

242 | 2010-08-03 | 2.2m Calar Alto | 40 | 165 | R | F |

243 | 2010-08-06 | 2.2m Calar Alto | 40 | 213 | R | F |

253 | 2010-09-06 | 2.2m Calar Alto | 50 | 164 | R | F |

254a | 2010-09-09 | 2.2m Calar Alto | 60 | 184 | R | F |

## 3 Data analysis

### 3.1 Re-analyse of the radial velocity observations

The discovery paper for WASP-10b (Christian et al., 2009) favoured an eccentric orbit for the planet and ). However, Maciejewski et al. (2011a) argued that eccentricity might have been overestimated due to stellar variability and favoured a circular orbit. Later, Husnoo et al. (2012) showed that the eccentricity detection was dependent on the first two radial velocity (RV) SOPHIE observations that could have a different zero point. Using a method of Bayesian model selection they find the eccentricity detection is not significant and obtained . Furthermore, the eccentricity is very correlated with the offset between the SOPHIE and FIES data. For all these reasons we assume a circular orbit. A better constraint on the eccentricity requires further RV observations. To estimate the stellar reflex velocity for a circular orbit we re-analysed the RV observations presented in Christian et al. (2009) using the new ephemeris given in Section 4.2.

Using an MCMC routine we fitted 5 parameters of a keplerian model for the host star reflex motion: orbital period, central transit time (T0), RV semi-amplitude (K), centre-of-mass velocity (), and an offset between SOPHIE and FIES data. We imposed a Gaussian prior to the orbital period and T0 using the new ephemeris given in Section 4.2 and obtained , , and offset. This solution has median rms of m/s which compares with the median RV uncertainty of m/s. Hence, the radial velocities have extra jitter which could be due to stellar activity.

Maciejewski et al. (2011a) argued that the stellar activity had a sinusoidal effect (semi-amplitude of m/s) in the measured RVs which should be accounted for when deriving the stellar reflex velocity due to WASP-10b. Following the procedure of Maciejewski et al. (2011a), we computed the Lomb-Scargle periodogram of the RV residuals after removing the circular model given above. This is shown in Figure 4 as well as the stellar rotation period of d determined by Smith et al. (2009) and the false alarm probability at (95.4%) confidence limit. We find no significant peak at the stellar rotation period, and hence we cannot confirm that the radial velocities are modulated by spot variability. Moreover, other host stars that show 1-2% photometric variability have smaller activity induced jitter in the measured RVs than the claimed m/s for WASP-10. For example, HD 189733b has a RV jitter semi-amplitude of m/s (Boisse et al., 2009) and CoRoT-7b has a RV jitter semi-amplitude smaller than m/s (Queloz et al., 2009). Therefore, we conclude that the extra jitter observed in the RVs of WASP-10 does not appear to be periodic. Including one jitter parameter for each data set in the radial velocity fit we find that the FIES data jitter is less than while the SOPHIE data jitter is . The latter is much higher than the expected jitter for SOPHIE supporting the result of Husnoo et al. (2012) about the possible offset of the two first SOPHIE points relative to the subsequent observations. Thus, we experimented discarding the first two SOPHIE points from the fit which resulted in a much lower estimated jitter for the SOPHIE data compatible with the expected values for this instrument. This was adopted as our final solution for the RVs with the remaining parameters being: FIES jitter , , , and offset. Our final solution for the stellar reflex velocity is in agreement with the value presented by Christian et al. (2009). The expected RV signal of the possible companion of m/s is below the median RV uncertainty, as also pointed out by Maciejewski et al. (2011a).

### 3.2 Photometric errors

An accurate estimate of the photometric errors is important to obtain reliable system parameters. We begin by scaling the errors of each light curve so that the of the best fitting model is 1.0. This is especially important in our case because we are combining data from different instruments which were reduced with different pipelines.

Exoplanet transit observations are often affected by correlated red noise which must be taken into account in error analyses. To estimate the time-correlated noise factor we followed the procedure of Winn et al. (2008). This consists of comparing the standard deviation of the residuals of the best fit solution for different time bins. For each light curve, we rescaled the errors by the estimated which is given in Table 3. This procedure does not distinguish between systematic noise and real stellar variability. All departures from the transit shape are considered as noise. For example, when we fit simultaneously all light curves (model 2), variations of the transit shape will also be considered as noise. Therefore, the values of for the simultaneous fit are higher than for the individual fits. If the shape variations are real this procedure will overestimate the errors. However, if they are indeed caused by systematic noise it will allow a more accurate error estimation.

Epoch | ind | simu |
---|---|---|

-103 | 1.0 | 1.0 |

-89 | 1.0 | 1.0 |

-88 | 1.0 | 1.0 |

0 | 1.0 | 1.0 |

20 | 1.0 | 1.3 |

21 | 1.5 | 1.7 |

26 | 1.5 | 2.0 |

36 | 1.4 | 1.4 |

56 | 1.0 | 1.0 |

123a | 1.1 | 1.1 |

123b | 1.1 | 1.1 |

124 | 1.1 | 1.0 |

132 | 1.0 | 1.0 |

133a | 1.4 | 1.7 |

133b | 1.7 | 1.8 |

134 | 1.3 | 1.3 |

143 | 1.5 | 1.8 |

144 | 1.4 | 1.6 |

145 | 1.9 | 2.3 |

147 | 2.4 | 4.3 |

158 | 1.0 | 1.0 |

242 | 1.5 | 1.8 |

243 | 1.8 | 2.3 |

252 | 2.0 | 2.0 |

253 | 1.7 | 2.0 |

254a | 1.0 | 1.2 |

254b | 3.0 | 3.0 |

265 | 2.7 | 3.0 |

274 | 1.8 | 2.2 |

285 | 1.3 | 1.5 |

### 3.3 Transit model and fitting procedure

To determine the planetary and orbital parameters, we fitted the light curves of WASP-10b with a transit model coupled with an Markov Chain Monte Carlo (MCMC) routine, a procedure similar to Barros et al. (2011) for WASP-21b. In our previous transit fitting analyses we fit all the light curves simultaneously. However, here we also present the individual light curve fits, mostly for comparison with previous results and call the latter model 1 and the former model 2. We used the Mandel & Agol (2002) transit model parametrised by the normalised separation of the planet, , ratio of planet radius to star radius, , orbital inclination, , and the transit epoch, , of each light curve. A circular orbit was adopted as discussed above. For each light curve, we also include two extra parameters ( and ) to account for a linear normalisation. We used the quadratic limb darkening coefficients (LDCs) deduced from the tables of Howarth (2011) for K, and [M/H] (Christian et al., 2009). Most of the light curves are of insufficient quality to allow for the fitting of the LDCs. To be consistent in the treatment of all the light curves we fixed the limb darkening coefficients to the theoretical values. Although this may lead to a slight underestimation of the errors, it was preferred to avoid spurious correlations. The values of the LDCs used for each filter are given in Table 4.

filter | linear coef | quadratic coef |
---|---|---|

sloan z’ | 0.391 | 0.210 |

I | 0.421 | 0.192 |

R | 0.602 | 0.133 |

rise | 0.612 | 0.146 |

V | 0.763 | 0.036 |

The MCMC chains were constructed following the Metropolis-Hastings algorithm. At each step of the chain each parameter is perturbed by a jump function, which is a random number sampled from a Gaussian distribution with mean zero and standard deviation equal to the respective parameter uncertainty. The new parameter set is accepted with probability: where is the difference in the of subsequent parameters sets. The jump functions are scaled by a common factor in order to ensure that the percentage of accepted steps lies between and . Further details about the fitting procedure can be found in Barros et al. (2011). For both models the initial 20% of each MCMC chain that corresponded to the burn-in phase were discarded and the remaining parts merged into a master chain. The best fit parameter was estimated as the median of its probability distribution and the limits as the value at which the integral of the distribution equals 0.341 from both sides of the median. To test convergence, the Gelman & Rubin (1992) statistic was computed for each fitted parameter and was found to be less than from unity for all the parameters in both models.

#### Model 1

The existence of an additional planet or moon affects the transit timings of the transiting planet, but can also affect the orbital parameters of the system. Hence, in previous analyses of the transit times of WASP-10, the authors opted for fitting each light curve individually (Dittmann et al., 2010; Maciejewski et al., 2011a; Maciejewski et al., 2011b). However, in some cases they were forced to fix the inclination of the system due to the poor quality of individual light curves. Here, we follow the same procedure primarily for comparison with previous methods. The results will also test possible variations of the light curve shape. Thus, we fitted each light curve individually using the MCMC method explained above with points per chain. Five chains were combined to obtain the final result. Instead of fixing the inclination we choose to impose a prior of the form:

(1) |

where we assume , taken from the simultaneous fit (model 2). This allows to account for the uncertainty of the inclination in the fit. We refer to this as model 1.

#### Model 2

For model 2 we assume that a possible additional planet has no significant effect on the transit shape during the time span of the observations. Therefore, we performed a simultaneous MCMC fit to all 30 light curves, fitting globally , , and and individually , and . We remind the reader that deviations from a constant shape are considered to be only due to systematic noise or intrinsic stellar variability and hence are treated as noise. Specifically, deviations of limb darkening parameters from the tabulated values are also considered as noise. However, since the brightness of WASP-10 has been shown to vary significantly due to star spots (Smith et al., 2009; Maciejewski et al., 2011a) a factor was included to account for out-of-transit flux variability between each light curve. This factor has been set arbitrarily to zero for our best quality light curve (Epoch 0, Johnson et al. (2009)) which is also the one taken in the reddest filter, and hence that should be less affected by stellar variability. Therefore, a total of parameters were fitted simultaneously. We computed seven MCMC chains each of points and different initial parameters.

## 4 Results

### 4.1 Model 1

The individual light curve fitted parameters are given in Table 5. Transit times were converted to Barycentric Dynamical Time (TDB) using the IDL codes kindly made available by Eastman et al. (2010). The weighted mean of the fitted parameters is: , , and degrees and the dispersion of inclination is equal to degrees as expected from the prior.

Comparing the values of for all light curves we conclude that individual deviations from the mean value of the normalised separation of the planet are smaller than except for 2 transits. These are epoch 20 which was observed in 2008-09-15 with FTN and epoch 147 (Dittmann et al., 2010). Both have higher than the mean by more than . Both of these transits show signs of correlated noise either due to systematics or stellar variability and will be discussed in section 5.2.

As expected, due to spot induced brightness variability of the host star, the parameter shows a wide spread. Two of the light curves (243, 254a) show a significantly higher () value of , while three light curves (0,133,158) show a significantly lower value. Therefore, we confirm transit depth variations reported by other authors and the need to include the extra parameter in the global fit. For completeness we also tested a model with free inclination and concluded that the inclination is constant within the uncertainties. Hence, we conclude that there are no significant shape variations apart from the transit depth, except for transits epoch 20 and 147.

Besides testing variations of transit shape, model 1 is also useful for direct comparison with previous transit times. This allows to better isolate the causes for the differences in the estimated transit times. The first three transits (-103,-89,-88) were previously combined together with the WASP data in a single ephemeris, and hence, they cannot be directly compared. Regarding the other previously published transit times, our results are within of previous results except for epochs 36 and 147 for which the difference is respectively and minutes. The disagreement can be explained by failure to account for the baseline function in the previous analysis of these transits. In our model a linear trend is fitted simultaneously with the transit model. When we do not fit a trend we recover the previous published results. T0 correlates with the baseline function because the latter affects the symmetry of the light curve. Therefore, including the baseline function is crucial to obtain reliable T0s

In general, we obtain slightly larger errors () than previous studies. For each case this is due to one or more of the following differences from the previous results: inclusion of a baseline function in the fit, inclusion of the red noise factor, and/or the treatment of the inclination. We argue that fixing or artificially constraining the inclination will cause underestimation of the transit time errors.

Despite the difference from previous results being small, the new transit times are in much better agreement with a linear ephemeris. The of a linear fit to the transit times is 91.95 for 28 degrees of freedom, i.e. which compares with found by Maciejewski et al. (2011a). However, the p-value is still low and hence the null hypothesis can be discarded at 5.7 . Four of the new transit times deviate more than from a linear ephemeris. These are epochs 20, 21, 143, and 147. Both epoch 21 and epoch 147 show trends in the light curve that can be affecting the transit times.

To test the model for the TTVs proposed by Maciejewski et al. (2011a) we fitted the sinusoidal model given by equation 1 of Maciejewski et al. (2011a) for frequencies: and , where P is the orbital period of WASP-10b. The preferred solution presented by Maciejewski et al. (2011a) for frequency had a p-value equal to which suggests that the model could be rejected at the level. However, we remind the reader that this sinusoidal model is only used to estimate the most probable TTV frequency to be used later in the 3-body simulations. Therefore, the p-value test cannot be directly applied. Nevertheless, we found that the sinusoidal model fit to the new transit times is worse, as it results in a much lower p-value for both frequencies and . Most importantly, for the new results we obtain 3.5 times smaller TTV semi-amplitudes: and .

In Figure 5 we show the Lomb-Scargle periodogram of the residual transit times of model 1 after subtracting a linear ephemeris. Notice that the TTV frequencies proposed by Maciejewski et al. (2011a) have disappeared which agrees with the much lower amplitudes found for the sinusoidal model fit.

Number | (degrees) | (hours) | T0 (days) | ||
---|---|---|---|---|---|

-103 | 12.02 0.21 | 0.1509 0.0034 | 88.67 0.12 | 2.198 0.043 | 4345.48630 0.00061 |

-89 | 12.02 0.17 | 0.1620 0.0030 | 88.65 0.12 | 2.219 0.033 | 4388.78602 0.00042 |

-88 | 11.13 1.06 | 0.1626 0.0046 | 88.67 0.12 | 2.411 0.245 | 4391.88193 0.00416 |

0 | 11.91 0.06 | 0.1575 0.0003 | 88.70 0.10 | 2.236 0.005 | 4664.03808 0.00005 |

20 | 12.35 0.13 | 0.1557 0.0011 | 88.64 0.12 | 2.143 0.017 | 4725.89175 0.00026 |

21 | 12.01 0.14 | 0.1570 0.0014 | 88.67 0.12 | 2.212 0.023 | 4728.98414 0.00038 |

26 | 11.15 0.41 | 0.1545 0.0021 | 88.72 0.11 | 2.392 0.097 | 4744.45206 0.00160 |

36 | 11.96 0.18 | 0.1641 0.0023 | 88.66 0.12 | 2.235 0.033 | 4775.37693 0.00051 |

56 | 11.93 0.15 | 0.1602 0.0017 | 88.65 0.12 | 2.233 0.027 | 4837.23105 0.00042 |

123 | 11.93 0.16 | 0.1572 0.0019 | 88.66 0.12 | 2.227 0.029 | 5044.44435 0.00044 |

123 | 11.90 0.16 | 0.1580 0.0026 | 88.68 0.12 | 2.238 0.031 | 5044.44407 0.00041 |

124 | 11.86 0.12 | 0.1559 0.0013 | 88.67 0.12 | 2.239 0.020 | 5047.53701 0.00033 |

132 | 11.67 0.51 | 0.1582 0.0057 | 88.66 0.12 | 2.281 0.116 | 5072.27832 0.00182 |

133 | 11.85 0.15 | 0.1598 0.0019 | 88.67 0.12 | 2.249 0.027 | 5075.37131 0.00041 |

133 | 12.11 0.21 | 0.1511 0.0024 | 88.67 0.12 | 2.181 0.038 | 5075.37145 0.00063 |

134 | 11.75 0.23 | 0.1623 0.0043 | 88.66 0.12 | 2.274 0.051 | 5078.46453 0.00080 |

143 | 11.76 0.14 | 0.1576 0.0024 | 88.68 0.11 | 2.265 0.028 | 5106.29690 0.00038 |

144 | 11.62 0.24 | 0.1629 0.0035 | 88.67 0.12 | 2.304 0.049 | 5109.39109 0.00080 |

145 | 11.50 0.26 | 0.1649 0.0034 | 88.67 0.12 | 2.334 0.055 | 5112.48544 0.00088 |

147 | 12.41 0.13 | 0.1601 0.0013 | 88.65 0.12 | 2.142 0.019 | 5118.66785 0.00030 |

158 | 12.00 0.15 | 0.1529 0.0018 | 88.66 0.12 | 2.206 0.026 | 5152.68920 0.00041 |

242 | 11.83 0.11 | 0.1608 0.0010 | 88.67 0.12 | 2.256 0.016 | 5412.47863 0.00024 |

243 | 11.81 0.09 | 0.1616 0.0009 | 88.68 0.12 | 2.264 0.013 | 5415.57161 0.00020 |

252 | 12.03 0.14 | 0.1550 0.0024 | 88.68 0.12 | 2.205 0.026 | 5443.40659 0.00035 |

253 | 11.80 0.09 | 0.1592 0.0009 | 88.68 0.11 | 2.260 0.013 | 5446.49881 0.00018 |

254 | 11.86 0.08 | 0.1610 0.0005 | 88.62 0.11 | 2.245 0.007 | 5449.59138 0.00007 |

254 | 11.86 0.12 | 0.1581 0.0014 | 88.66 0.12 | 2.245 0.020 | 5449.59120 0.00031 |

265 | 11.71 0.17 | 0.1614 0.0026 | 88.66 0.12 | 2.281 0.035 | 5483.61165 0.00054 |

274 | 11.92 0.10 | 0.1596 0.0008 | 88.67 0.12 | 2.236 0.014 | 5511.44560 0.00021 |

285 | 11.48 0.53 | 0.1573 0.0038 | 88.69 0.12 | 2.323 0.122 | 5545.46688 0.00205 |

### 4.2 Model 2

The individual transit times and the out-of-transit flux variability derived from a simultaneous fit of the 30 light curves of WASP-10 are given in Table 6.

Number | df | T0 | Time residuals |
---|---|---|---|

days | sec | ||

-103 | 0.047 0.065 | 4345.48665 0.00038 | -28.4 32.6 |

-89 | -0.067 0.030 | 4388.78596 0.00041 | 66.7 35.7 |

-88 | -0.040 0.057 | 4391.87917 0.00054 | 108.0 46.4 |

0 | 0.000 0.000 | 4664.03808 0.00005 | -2.4 4.6 |

20 | 0.030 0.032 | 4725.89179 0.00038 | -77.9 32.4 |

21 | 0.007 0.020 | 4728.98421 0.00043 | -104.6 37.3 |

26 | 0.065 0.024 | 4744.44923 0.00021 | 13.9 18.3 |

36 | -0.066 0.053 | 4775.37696 0.00054 | 51.0 46.7 |

56 | -0.029 0.021 | 4837.23102 0.00043 | 5.8 37.4 |

123a | -0.002 0.028 | 5044.44407 0.00031 | 20.8 26.9 |

123b | 0.004 0.024 | 5044.44434 0.00046 | 44.1 39.3 |

124 | 0.018 0.016 | 5047.53700 0.00030 | 38.0 26.3 |

132 | 0.034 0.044 | 5072.27905 0.00052 | 56.7 45.0 |

133a | 0.080 0.034 | 5075.37113 0.00073 | 0.9 62.7 |

133b | -0.022 0.025 | 5075.37133 0.00043 | 18.2 37.4 |

134 | -0.028 0.036 | 5078.46489 0.00062 | 90.3 53.5 |

143 | 0.023 0.020 | 5106.29655 0.00032 | -161.2 27.9 |

144 | -0.063 0.053 | 5109.39103 0.00095 | -10.2 82.0 |

145 | -0.110 0.069 | 5112.48544 0.00110 | 135.1 95.4 |

147 | -0.026 0.029 | 5118.66795 0.00056 | -119.4 48.8 |

158 | 0.057 0.020 | 5152.68916 0.00044 | -16.6 37.7 |

242 | -0.043 0.014 | 5412.47858 0.00028 | -4.0 24.2 |

243 | -0.053 0.013 | 5415.57163 0.00026 | 23.4 22.4 |

252 | 0.009 0.021 | 5443.40631 0.00027 | 33.5 22.9 |

253 | -0.018 0.012 | 5446.49884 0.00021 | 16.1 18.4 |

254b | -0.010 0.021 | 5449.59123 0.00031 | -12.8 27.0 |

254a | -0.040 0.006 | 5449.59139 0.00009 | 0.7 8.1 |

265 | -0.025 0.026 | 5483.61125 0.00044 | -13.7 38.3 |

274 | -0.025 0.011 | 5511.44562 0.00025 | -29.9 21.8 |

285 | 0.037 0.025 | 5545.46534 0.00029 | -56.4 25.2 |

We conclude that three light curves 242, 243, 254a (Maciejewski et al. 2011b) have a significantly lower df (at ) than the chosen zero point. In our notation a negative df means the star is fainter which corresponds to a larger transit depth. Some light curves have a higher df, but this is non-significant. Therefore, we confirm that our chosen zero point is valid. However, this zero point might not correspond to the real maximum brightness of WASP-10, and hence the planet to star ratio could still be overestimated.

In Figure 6 we plot the difference between the our transit times derived with model 1 (black stars) and model 2 (red dots) and previously published transit times. We also show the ratio of transit time errors. The highest difference in times is found for transits 36, 132 and 147. For epochs 36 and 147 this is due to the baseline function not been accounted for in previous fits, as mentioned before. For epoch 132 the difference is within the errors and it is due to epoch 132 being a partial transit and hence very sensitive to the details of the model fit. The transit times derived from model 1 and model 2 are consistent within . The biggest differences in the results from the two models are in partial transits. In model 2, the global shape of the transit defined by the complete transits constrains the shape of partial transits, this allows a better constrain on the baseline function that correlates with the transit times.

From the bottom panel of Figure 6 it is also clear that our estimated errors are higher than previous ones. This is due to the inclusion of red noise in our analysis. For model 2 the uncertainties of T0s of partial transits are smaller than for model 1 because the shape of the transit in model 2 is almost fixed by the other light curves. In contrast, some T0 uncertainties of complete transits are higher because differences from the common shape are considered red-noise which is accounted for as explained in Section 3.2. This is the key difference between the models and implies that in this case model 2 will produce much more reliable results, as will be discussed later.

We compute a linear ephemeris using the new estimated transit times for model 2:

(2) |

where and . Time residuals (or TTVs) from the linear ephemeris are given in Table 6 and shown in Figure 7. The RISE data are show as stars, the new FTN data as squares, and the new results based on previously published light curves as circles.

The of the linear ephemeris is 86 for 28 degrees of freedom and the p-value is so we conclude that the linear fit is not a a good fit at confidence level. However, only epoch 143 deviates significantly (at ) from the linear ephemeris. This transit is the only one obtained with the 2.0m Rozhen telescope and is a partial transit. If we discard epoch 143 the decreases to 53 for 27 degrees of freedom and the p-value increases to 0.00176 and hence we can only reject the null hypotheses at the level.

For these new transit times derived with model 2 we also tested the sinusoidal model proposed by Maciejewski et al. (2011a) for the two proposed frequencies as done in the previous section for model 1 and found similar results. We find and while both p-values are less then and hence the model could be rejected at the level. We reiterate that the expected TTV variations are not strictly periodic and the p-value can not be directly applied. However, clearly the new results do not support the TTV solution claimed by Maciejewski et al. (2011a). We analyse a perturber model for the our derived transit times in Section 5.1.

The periodogram of the time residuals relative to a linear ephemeris are given in Figure 8. Comparing with the model 1 results, it is clear that the peaks around the star’s rotation period have decreased, implying that model 2 results are less affected by stellar variability.

### 4.3 System parameters

Although the main objective of this study is not deriving updated parameters for WASP-10, for consistency with the simulations presented in Section 5.1 we derive the system parameters from the geometric fitted parameters, given in the top of Table 7. The same procedure as Barros et al. (2012) is followed.

We estimate the stellar mass from models by interpolating the Yonsei-Yale stellar evolution tracks by Demarque et al. (2004) using the effective temperature K and metallicity dex from Christian et al. (2009) and our new derived stellar density . We obtain a sightly larger stellar density than derived by Johnson et al. (2009) ( ) which results in a slightly lower mass , but consistent with previous results.

Combining the stellar mass, stellar reflex velocity estimated in Section 3.1, and the geometric parameters derived with model 2, we derive the system parameters which are given on Table 7. The parameters were derived using the equations of Seager & Mallén-Ornelas (2003) and Kipping (2010). Specifically, the refers to the total transit duration and is the duration of the sky-projected centre of the planet overlapping the stellar surface and is defined by equation 15 of Kipping (2010).

Parameter | Value |
---|---|

Normalised separation | |

Planet/star radius ratio | |

Transit duration [hours] | |

Transit duration [hours] | |

Orbital inclination [degrees] | |

Impact parameter [] | |

Stellar density [] | |

Orbital semimajor axis [AU] | |

Stellar mass [] | |

Stellar radius [] | |

Planet mass [] | |

Planet radius [] | |

Planet density [] |

## 5 Discussion

### 5.1 Additional perturber

From the analysis of the TTVs, Maciejewski et al. (2011a) postulated the presence of a perturber with a mass of on a circular orbit with a period of days. This orbital period places the system close a 5:3 mean-motion resonance. The existence of the planetary companion was motivated by a large reduced chi-square associated to a linear fit of the mid-transit times.

In the present study, additional transit observations and new estimations of the error bars lead to for a linear ephemeris. This value is lower than the one obtained by Maciejewski et al. (2011a), but remains significantly larger than one, suggesting the existence of an additional signal. As in Maciejewski et al. (2011a), we make the hypothesis that this signal is produced by a planetary companion and try to constrain its mass and orbital properties. For that purpose, we performed numerical 3-body simulations with 400 values of the companion orbital period ranging uniformly between 1.5 to 10 , 60 values of the initial longitude (at the time of the first transit Epoch=-103) uniformly distributed over 360 degrees, and 100 masses between 0.01 to 5 Jupiter masses with logarithmic steps. In all simulations, the initial conditions correspond to circular coplanar orbits. The best fit to the timing residuals is very isolated in the parameter space and gives days, , and deg with . The TTV signal produced by this set of parameters is superimposed to the timing residuals in Fig. 9.

Since our best solution is completely different from that of
Maciejewski et al. (2011a), we analyse its robustness to the
mid-transit time uncertainties. We generate 100 sets of artificial timing residuals
with the same 27 Epochs as the observations. The value of each point
is then taken from a normal distribution whose mean corresponds to our
best solution and width equals the uncertainty of that Epoch.
To reduce the cpu time, we forgo the 3-body simulations in favour of
the perturbation method ptmeth (Nesvorný
& Morbidelli, 2008).
This is a powerful analytical method about faster than numerical
simulations which provides accurate results outside of mean-motion
resonances (Nesvorný
& Morbidelli, 2008). It is also the
technique used by Maciejewski et al. (2011a) to obtain a good
estimate of the location of the best parameters before performing
numerical simulations. The best solutions derived for each of the 100
artificial data sets are plotted in Fig. 10. Each point
correspond to a different data set. The most remarkable result is that
these best solutions, with reduced chi-square^{3}

It can be noted that the most massive companions at a given orbital period, in Fig. 10, are those with the lowest eccentricity. This is expected since eccentricities increase the orbital interactions and hence the TTVs. This property has already been used on different systems to set an upper boundary on the mass of a perturber as a function of its semi-major axis (e.g., Steffen & Agol, 2005) . For instance, in the case of WASP-10, from the -body simulations we obtain a 3 sigma upper mass equal to at the 5.28 day period claimed by Maciejewski et al. (2011a). This is not stringent enough to rule out the solution of Maciejewski et al. (2011a), and neither does it confirm it since it lies at the 3 sigma limit and not at the best solution. Moreover, it should be stressed that any perturber with mass beyond would have been detected by radial velocities. They can thus be discarded.

This analysis illustrates the difficulties of interpreting the perturbations produced by a non-transiting companion in very low signal to noise data (see, e.g., Boué et al., 2012). The detection and the characterisation of non-transiting planets do require strong TTV signals as in the cases of Kepler-19 and KOI-872 (Ballard et al., 2011; Nesvorný et al., 2012).

### 5.2 The effect of spots in measured TTVs

Spot occultation features have been shown to affect measured transit times (Miller-Ricci et al., 2008; Alonso et al., 2009; Sanchis-Ojeda et al., 2011). The transit light curves of WASP-10b show features that might be due to the occultation of stellar spots by the planet. Maciejewski et al. (2011b) shows that the features in epochs 242, 243, 253 are compatible with spots along the transit cord with an area of approximately 25% -37% that of planetary disc. However, the same shapes features can also be caused by systematic noise (Barros et al., 2011). In this section, we will test the effect of spot features in the light curve shape and consequently in the fitted parameters, particularly in the measured transit times. Since we are analysing data from different instruments, filters, and signal-to-noise ratio we will not attempt to test if the ”spot features” are due to real spots occultations or systematic noise. Instead, we assume that independently of its causes, a ”spot feature” will have similar effect on the shape of the light curve.

We modelled the spot features as inverted transits, a procedure similar to Kundurthy et al. (2011) and Maciejewski et al. (2011b). We used the transit model described above, but for the spot features we assume no limb darkening and zero planet-spot impact parameter. Values for and were fixed to those of Table 7 and an extra factor was included to account for the spot-to-star contrast. This is a very simplified spot model, but reproduces the observed features quite well. To obtain typical values for the spot parameters for our simulations we used the spot seen in light curve epoch 265 as a template. According to our tests of light curve systematics this appears to be a real spot.

To simulate the effects of a spot in a transit light curve we added the model spot to the model light curve of WASP-10b at different in-transit phases. The resulting light curve was then fitted with a transit model not including a spot to test how the spot affects the fitted parameters. Errors were not included in the simulated light curves. We also do not account for geometrical foreshortening (Sanchis-Ojeda et al., 2011), i.e. the same spot amplitude was used for each phase. The derived parameters as a function of central spot time relative to the mid-transit time are shown in Figure 11, where evidently a spot similar to the one in light curve epoch 265 is able to create fake transit time variations up to 1.2 minutes. Clearly, the spots will also significantly affect the derived shape parameters: , , and the inclination. For the best signal-to-noise light curves, the spot effect will be more than of the errors quoted in Table 5.

Noteworthily, the effect of the spot is higher when the spot is in the limb of the star, although spots at longitudes up to 30 degrees from the limb will significantly affect the measured transit times. Furthermore, spot features during ingress and egress are harder to detect in the light curve both analytically and by eye. This is illustrated in Figure 12 which shows an example of a light curve where the planet occults a spot during egress. The combined light curve is well fitted by a transit model with a shorter duration and earlier egress time which results in an earlier T0. Moreover, in individual fits this light curve would appear to have low red noise.

### 5.3 Spot modelling

As seen in the previous section, the light curves that show higher TTVs due to spots also show a larger (or shorter duration). Hence, a high value is an indication that the light curve could be affected by spots. Moreover, in this case the red noise might be underestimated for individual fits. To check if the TTVs presented in Table 6 are affected by spots we tested the light curves that have a high red noise and also the ones with high .

From the results in Table 5, as mentioned before, we conclude that two light curves have higher than the mean value. These are epochs 20 and 147. Table 3 shows that five light curves have red noise higher than 2 for model 2, namely epochs 143, 145, 147, 254b, and 265. These 6 light curves were fitted with a transit model including a spot model as described in the previous section. The impact parameter was fixed to 0 and we assume no limb darkening, the contrast, phase, and duration of the spot were fitted. When we include a spot feature we find variations on the fitted transit times up to 2.2 minutes (for epoch 20). We extended the analysis to all the light curves, but found T0 variations only in another light curve; epoch 21. The difference between the derived T0 when including a spot model and the T0 derived with model 2 are given in Table 8 as well as the corresponding red noise for the one spot model parametrised by . For curiosity, we also give the spot fitted parameters in Table 8, however the fits are highly degenerate and dependent on the starting parameters. We confirm that for all the light curves, the time correlated features can be well described by a spot, given that the estimate of the red noise significantly decreases. However, in some cases such as epoch 147 residual red noise remains in the light curve after accounting for one spot. Although the fits are highly degenerate, we confirm that only spots close to the limb of the star have a significant effect on the measured transit times. Although the spot induced TTV can be large, the new transit times agree with model 2 results. However, this would not be the case if we had not included red noise. Furthermore, the existence of a spot explains the significant TTVs found for model 1 which does not penalise shape variations. Hence, in this case model 1 underestimates the errors in the mid-transit times.

In Figure 13, we show the timing residuals for the one spot model. The of a linear ephemeris is 76 for 28 degrees of freedom and the p-value is which would lead to a rejection of the model. Epoch 143 is still inconsistent () with a linear ephemeris since this discrepancy cannot be explained with a spot model. The ingress of this transit is clearly earlier than predicted by the linear ephemeris which cannot be caused by a spot. If we exclude epoch 143 the of the linear ephemeris reduces to 43 for 27 degrees of freedom (= 2.56) and the p-value is 0.028 and hence a linear ephemeris is no longer rejected. Therefore, we conclude that if we account for “spot features” the transit times of WASP-10 are consistent with a linear ephemeris.

Figure 14 shows the Lomb-Scargle periodogram of the residuals from a linear ephemeris. Notice that the peaks close to the rotation period of the star have decreased. The main difference between this figure and the results for model 2 is due to the “spot features” during egress found for epochs 20, 21 and 147. Given that these “spot features” add power at the rotation period of the star, this suggests that these spots are real rather than systematic noise.

However, we caution the reader to an over interpretation of these results because our model is very simplified, the fits are degenerate, and we cannot confirm that all the features present in the light curve are real spots. A complete spot model of WASP-10 is beyond the scope of this paper, our objective being to estimate the mid-transit timing errors that can be caused by spots.

In summary, stellar spots or correlated noise can significantly bias the derived parameters of the light curves especially if the spot occurs during ingress or egress. However, the errors due to spots can be accounted for if we include red noise due to shape variations as done for model 2.

Epoch | delta T0 | heigh | |||
---|---|---|---|---|---|

(sec) | (min) | (mmag) | (min) | ||

20 | -135 | 1. | 50.72 | 7.98 | 31.5 |

21 | -78 | 1. | 36.97 | 3.13 | 62.58 |

143 | -6 | 1.4 | 1.47 | 1.61 | 46.16 |

145 | 5 | 1.3 | -10.50 | 8.22 | 101.94 |

147 | -126 | 2. | 52.38 | 5.93 | 52.60 |

243 | 7 | 1.1 | 39.73 | 3.03 | -10.28 |

254b | 17 | 1. | -26.48 | 1.97 | 55.69 |

265 | 2 | 1. | 32.35 | 3.99 | 28.71 |

### 5.4 Effects of partial transits

Several authors have noted that the uncertainties in the T0s derived from partial transit are underestimated (e.g Gibson et al., 2009; Barros et al., 2012). Given that the one transit that shows significant transit timing variations is a partial transit (epoch 143), in this section we illustrate how the derived transit times are affected in partial transits. We assume here that the light curves are fitted simultaneously (model 2), and hence the shape of the transits is set by the complete transits, i.e. we fixed and in our test fits. To test the TTVs induced by partial transits, we used some of our complete light curves, from which we subtract model 2 estimated T0 ( Table 6). The light curves were then gradually cut at a time -Tp or Tp which corresponds to cutting the beginning or end of the observation respectively. Resulting partial light curves were fitted to derived the induced T0 variation. Red noise was estimated for each fake light curve and included in the errors. In Figure 15 we show our test applied to the light curves epoch 0, 145, and 147 as examples. For other light curves the derived TTVs are enclosed in the examples presented. We conclude that light curves that do not have red noise such as epoch 0 show partial transit induced TTVs that are consistent with zero. However, for light curves with high red noise the induced TTVs can reach several minutes and the derived errors are indeed underestimated. The errors decrease with Tp due to red noise underestimation. For example, the TTVs derived from fake partials of epoch 147 are significant at even when we include red noise. This is caused by the correlation between red noise and the normalisation parameters and in partial light curves. Notice that TTVs are induced even when Tp is larger than half of the transit duration, i.e. the transit needs to include some out of transit data. Therefore this could explain the 2.7 minutes TTV found for epoch 143, though we remind that this is also the only transit taken with 2.0m Rozhen telescope.

If we discard all partial transits, the transit residuals of model 2 are consistent with a linear ephemeris, with a for 17 degrees of freedom ( ) and a p-value . Model 1 results are still inconsistent with a linear ephemeris at the level since transits show significant transit time variations due to spots that are not accounted for in this model as was shown in the previous section. The one spot model timing residuals of the complete transits are also consistent with a linear ephemeris: for 17 degrees of freedom () and a p-value . The increase of the p-value when we discard partial transits indicates that errors in transit times from partial transits are indeed underestimated, as discussed above.

In summary, it is possible to obtain a good estimation of the mid-transit times of partial transits if they do not have red noise. However, in the presence of red noise the errors of the mid-transit times of partial transits can be significantly underestimated due to our inability to recognise red noise in partial light curves. Unfortunately, complete transits are sometimes unattainable due to weather and technical issues. For example, all scheduled RISE observations were intended to cover the complete transit. We conclude that transit times derived from partial transits have to be regarded with extreme care.

## 6 Conclusion

We have presented eight new transit observations of WASP-10b with RISE mounted on the Liverpool telescope and Merope mounted on the Faulkes Telescope North. The new observations were combined with 22 previously published transit light curves to derive the transit times homogeneously. System parameters were re-derived assuming a circular orbit and assuming that the star is less active at maximum brightness. However, due to bright spots some stars, like the Sun, are actually brighter at the activity maximum. No bright spot is seen in the transit cord of WASP-10, therefore, we assume that the observed stellar modulation is due to dark spots.

We show that the significant transit timing variations for WASP-10b reported by Maciejewski et al. (2011a); Maciejewski et al. (2011b) were due to the normalisation parameters being neglected in some of previous fits to the data. We find a time difference of 120 and 49 seconds for epochs 36 and 147, respectively when we account for a baseline function in the transit fit. It is crucial to account for the baseline function because it correlates strongly with the estimated transit times.

We present two models to derive the transit times. Model 1 where the transits are fitted individually and model 2 where the transits are fitted simultaneously. Comparison between the two models lead us to suspect that the transits times were affected by spot features or systematic noise. We show that spot features or systematic noise during ingress and egress have a significant effect on the measured transit times leading to spurious TTV detections. Since they also cause transit duration variations they are treated as red noise in model 2 which explains the differences found between the two models. Nevertheless, the transit times derived with model 2 show higher dispersion than expected for a linear ephemeris.

We tested if the transit times could be explained by a perturber model. As expected the addition of a planetary companion in the system reduces the rms of the residuals. However, there is not more evidence for the TTV signal claimed by Maciejewski et al. (2011a); Maciejewski et al. (2011b) neither in the periodogram nor in the n-body analysis of the TTVs. In fact, we found that other configurations with perturber masses ranging from 0.01 to 5 and period between 1.5 and the orbital period of WASP-10b fit the data equivalently well. Therefore, we conclude that the transit times of WASP-10b have insufficient accuracy and precision to validate the presence of any planetary companion.

An alternative model for the dispersion seen in the transit times of WASP-10b is stellar activity induced TTVs. We tested this hypothesis by including a spot feature in the transit model fit. Our simple model shows that transits 20, 21 and 147 are consistent with having a spot feature during egress that results in transit time variations of -135, -78 and -126 seconds respectively. This explains the dispersion found for the transit times derived with model 2. We suggest that this spot induce transit variations are due to real spots rather than systematics because they add power to the residuals at periods close to the orbital period of the host star. Moreover, other spots are seen during transit, in particular epochs 242, 243 and 253 which were studied in detail by Maciejewski et al. (2011b) that concluded that they were consistent with real spots. When we account for these spot features, the derived transit times are consistent with a linear ephemeris with the exception of epoch 143, which we show it could be due to this transit being partial.

We also show that partial transits may produce fake TTVs up to 3 minutes, and hence T0s derived from partial transits should be regarded with extreme care. When we discard partial transits the dispersion of the transit times relative to a linear ephemeris decreases significantly both for model 2 results and for the one spot model results.

## 7 Acknowledgements

We thank Damian Christian, John Johnson, Jason Dittmann, Tereza Krejcova, and Gracjan Maciejewski for kindly sharing their light curves on WASP-10 and for interesting discussions. We thank the anonymous referee for contributing to the improvement of this manuscript. We thank Tom Marsh for the use of the ULTRACAM pipeline. SCCB thanks John Southworth and Rodrigo Díaz for interesting discussions. The RISE instrument mounted at the Liverpool Telescope was designed and built with resources made available from Queen’s University Belfast, Liverpool John Moores University and the University of Manchester. The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias with financial support from the UK Science and Technology Facilities Council.

### Footnotes

- pagerange: Transit timing variations in WASP-10b induced by stellar activity?–References
- pubyear: 2002
- The values of the reduced chi-square obtained with ptmeth are statistically lower than 1.77 found with the 3-body simulations. This is expected since the model of the perturbation method contains two additional parameters (eccentricity and longitude of periastron) in comparison to the slowest -body integration method. Moreover, the mock data sets fitted by ptmeth are not affected by stellar activity.

### References

- Adams E. R., López-Morales M., Elliot J. L., Seager S., Osip D. J., 2010, ApJ, 714, 13
- Agol E., Steffen J., Sari R., Clarkson W., 2005, MNRAS, 359, 567
- Alonso R., Aigrain S., Pont F., Mazeh T., CoRoT Exoplanet Science Team 2009, in IAU Symposium Vol. 253 of IAU Symposium, Searching for the secondary eclipse of CoRoT-Exo-2b and its transit timing variations. pp 91–96
- Ballard S., Fabrycky D., Fressin F., Charbonneau D., Desert J.-M., Torres G., Marcy G., Burke C. J., Isaacson H., Henze C., Steffen J. H., Ciardi D. R., et al. 2011, ApJ, 743, 200
- Barros S. C. C., Pollacco D. L., Gibson N. P., Howarth I. D., Keenan F. P., Simpson E. K., Skillen I., Steele I. A., 2011, MNRAS, 416, 2593
- Barros S. C. C., Pollacco D. L., Gibson N. P., Keenan F. P., Skillen I., Steele I. A., 2012, MNRAS, 419, 1248
- Bean J. L., 2009, A&A, 506, 369
- Boisse I., Moutou C., Vidal-Madjar A., Bouchy F., Pont F., Hébrard G., Bonfils X., Croll B., Delfosse X., Desort M., Forveille T., Lagrange A.-M., et al. 2009, A&A, 495, 959
- Boué G., Oshagh M., Montalto M., Santos N. C., 2012, MNRAS, 422, L57
- Christian D. J., Gibson N. P., Simpson E. K., Street R. A., Skillen I., Pollacco D., Collier Cameron A., Joshi Y. C., Keenan F. P., Stempels H. C., et al 2009, MNRAS, 392, 1585
- Czesla S., Huber K. F., Wolter U., Schröter S., Schmitt J. H. M. M., 2009, A&A, 505, 1277
- Demarque P., Woo J., Kim Y., Yi S. K., 2004, ApJS, 155, 667
- Dhillon V. S., Marsh T. R., Stevenson M. J., Atkinson D. C., Kerry P., Peacocke P. T., Vick A. J. A., Beard S. M., Ives D. J., Lunney D. W., McLay S. A., Tierney C. J., Kelly J., Littlefair S. P., Nicholson R., Pashley R., Harlaftis E. T., O’Brien K., 2007, MNRAS, 378, 825
- Díaz R. F., Rojo P., Melita M., Hoyer S., Minniti D., Mauas P. J. D., Ruíz M. T., 2008, ApJL, 682, L49
- Dittmann J. A., Close L. M., Scuderi L. J., Morris M. D., 2010, ApJ, 717, 235
- Eastman J., Siverd R., Gaudi B. S., 2010, PASP, 122, 935
- Fukui A., Narita N., Tristram P. J., Sumi T., Abe F., Itow Y., Sullivan D. J., Bond I. A., Hirano T., Tamura M., Bennett D. P., Furusawa K., et al. 2011, PASJ, 63, 287
- Gelman A., Rubin D., 1992, Statistical Science, 7, 457
- Gibson N. P., Pollacco D., Simpson E. K., Barros S., Joshi Y. C., Todd I., Keenan F. P., Skillen I., Benn C., Christian D., Hrudková M., Steele I. A., 2009, ApJ, 700, 1078
- Gibson N. P., Pollacco D., Simpson E. K., Joshi Y. C., Todd I., Benn C., Christian D., Hrudková M., Keenan F. P., Meaburn J., Skillen I., Steele I. A., 2008, A&A, 492, 603
- Gibson N. P., Pollacco D. L., Barros S., Benn C., Christian D., Hrudková M., Joshi Y. C., Keenan F. P., Simpson E. K., Skillen I., Steele I. A., Todd I., 2010, MNRAS, 401, 1917
- Holman M. J., Fabrycky D. C., Ragozzine D., Ford E. B., Steffen J. H., Welsh W. F., Lissauer J. J., Latham D. W., Marcy G. W., et al. 2010, Science, 330, 51
- Holman M. J., Murray N. W., 2005, Science, 307, 1288
- Howarth I. D., 2011, MNRAS, 413, 1515
- Hoyer S., Rojo P., López-Morales M., Díaz R. F., Chambers J., Minniti D., 2011, ApJ, 733, 53
- Husnoo N., Pont F., Mazeh T., Fabrycky D., Hébrard G., Bouchy F., Shporer A., 2012, MNRAS, 422, 3151
- Johnson J. A., Winn J. N., Cabrera N. E., Carter J. A., 2009, ApJL, 692, L100
- Kipping D. M., 2010, MNRAS, 407, 301
- Krejčová T., Budaj J., Krushevska V., 2010, Contributions of the Astronomical Observatory Skalnate Pleso, 40, 77
- Kundurthy P., Agol E., Becker A. C., Barnes R., Williams B., Mukadam A., 2011, ApJ, 731, 123
- Lissauer J. J., Fabrycky D. C., Ford E. B., Borucki W. J., Fressin F., Marcy G. W., Orosz J. A., Rowe J. F., Torres G., Welsh W. F., Batalha N. M., et al. 2011, Nature, 470, 53
- Maciejewski G., Dimitrov D., Neuhäuser R., Niedzielski A., Raetz S., Ginski C., Adam C., Marka C., Moualla M., Mugrauer M., 2010, MNRAS, 407, 2625
- Maciejewski G., Dimitrov D., Neuhäuser R., Tetzlaff N., Niedzielski A., Raetz S., Chen W. P., Walter F., Marka C., Baar S., Krejcová T., Budaj J., Krushevska V., Tachihara K., Takahashi H., Mugrauer M., 2011a, MNRAS, 411, 1204
- Maciejewski G., Raetz S., Nettelmann N., Seeliger M., Adam C., Nowak G., Neuhäuser R., 2011b, A&A, 535, A7
- Mandel K., Agol E., 2002, ApJL, 580, L171
- Miller-Ricci E., Rowe J. F., Sasselov D., Matthews J. M., Kuschnig R., Croll B., Guenther D. B., Moffat A. F. J., Rucinski S. M., Walker G. A. H., Weiss W. W., 2008, ApJ, 682, 593
- Nesvorný D., Kipping D. M., Buchhave L. A., Bakos G. Á., Hartman J., Schmitt A. R., 2012, Science, 336, 1133
- Nesvorný D., Morbidelli A., 2008, ApJ, 688, 636
- Queloz D., Bouchy F., Moutou C., Hatzes A., Hébrard G., Alonso R., Auvergne M., Baglin A., Barbieri M., Barge P., Benz W., Bordé P., et al. 2009, A&A, 506, 303
- Sanchis-Ojeda R., Winn J. N., Holman M. J., Carter J. A., Osip D. J., Fuentes C. I., 2011, ApJ, 733, 127
- Seager S., Mallén-Ornelas G., 2003, ApJ, 585, 1038
- Smith A. M. S., Hebb L., Collier Cameron A., Anderson D. R., Lister T. A., Hellier C., Pollacco D., Queloz D., Skillen I., West R. G., 2009, MNRAS, 398, 1827
- Steele I. A., Bates S. D., Gibson N., Keenan F., Meaburn J., Mottram C. J., Pollacco D., Todd I., 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series Vol. 7014 of Edited by McLean, Ian S.; Casali, Mark, RISE: a fast-readout imager for exoplanet transit timing. pp 70146J–70146J
- Steffen J. H., Agol E., 2005, MNRAS, 364, L96
- Steffen J. H., Ragozzine D., Fabrycky D. C., Carter J. A., Ford E. B., Holman M. J., Rowe J. F., Welsh W. F., Borucki W. J., Boss A. P., Ciardi D. R., Quinn S. N., 2012, Proceedings of the National Academy of Science, 109, 7982
- Stetson P. B., 1987, PASP, 99, 191
- Winn J. N., Holman M. J., Torres G., McCullough P., Johns-Krull C., Latham D. W., Shporer A., Mazeh T., Garcia-Melendo E., Foote C., Esquerdo G., Everett M., 2008, ApJ, 683, 1076