Issue |
A&A
Volume 504, Number 2, September III 2009
|
|
---|---|---|
Page(s) | 561 - 564 | |
Section | Stellar atmospheres | |
DOI | https://doi.org/10.1051/0004-6361/200912329 | |
Published online | 15 July 2009 |
Transit mapping of a starspot on CoRoT-2
Probing a stellar surface with planetary transits
U. Wolter1 - J. H. M. M. Schmitt1 - K. F. Huber1 - S. Czesla1 - H. M. Müller1 - E. W. Guenther2 - A. P. Hatzes2
1 - Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg, Germany
2 -
Thüringer Landessternwarte Tautenburg, Sternwarte 5, 07778 Tautenburg, Germany
Received 14 April 2009 / Accepted 17 June 2009
Abstract
We analyze variations in the transit lightcurves of CoRoT-2b,
a massive hot Jupiter orbiting a highly active G star.
We use one transit lightcurve to eclipse-map a photospheric spot
occulted by the planet.
In this case study we determine
the size and longitude of the eclipsed portion of the starspot
and systematically study the corresponding uncertainties.
We determine a spot radius between
and
on the stellar surface and the spot longitude
with a precision of about
degree.
Given the well-known transit geometry of the CoRoT-2 system, this implies
a reliable detection of spots on latitudes typically covered by sunspots;
the size of the modelled spot is comparable
to large spot groups on the Sun.
We discuss the future potential of eclipse mapping by planetary transits for
the high-resolution analysis of stellar surface features.
Key words: stars: planetary systems - stars: activity - stars: late-type - stars: imaging - stars: individual: CoRoT-2
1 Introduction
The atmosphere of the Sun shows inhomogeneities down to the smallest scales currently accessible to solar observations, which are of the order of 50 km (e.g. Scharmer et al. 2002). Such intricate fine structure can also be expected in the atmospheres of other active stars. However, the best currently available stellar observations only resolve surface features down to the size of a few degrees on the surface, corresponding to several 10 000 km on a main sequence star.
The increasing number of known transiting extrasolar planets offers an outstanding opportunity to study surface inhomogeneities of their host stars with an unprecedented surface resolution. As an example of this new technique, we present the analysis of one transit lightcurve of the planetary system CoRoT-2, recently detected by the CoRoT satellite (Rouan et al. 1998). Our study indicates that under favourable conditions, the CoRoT lightcurves allow the study of e.g. starspots and/or faculae down to a sub-degree scale on the stellar surface.
Deformations of planetary transit lightcurves, attributable to dark spots, have been observed for several systems: HD 189 733, Pont et al. (2007) HD 209 458, Silva-Valio (2008) Tres-1, Rabus et al. (2009) and CoRoT-2, Lanza et al. (2009). Pont et al.'s study, based on HST data, for the first time showed that low-noise transit photometry of exoplanets can yield detailed information about surface features of their host star. We present a systematic analysis of the spot locations and extensions, including their uncertainties, that can be deduced from transit lightcurves.
CoRoT-2a (GSC 00465-01282) is an apparently young late G-dwarf star (Alonso et al. 2008, AL08 in the following);
it is unsually active and intrinsically variable among the
presently known planetary host stars.
Its rotation period of
d (Lanza et al. 2009)
amounts to less than three times the planetary orbit period
:
.
Densely sampled high-precision photometry in combination with
spectroscopic measurements
make the transit geometry of
CoRoT-2b exceptionally well known (AL08, Bouchy et al. 2008):
The planet is on a nearly circular orbit (
)
which is seen close to edge-on (
)
and oriented approximately perpendicular to the stellar rotation axis (deviating by
).
These attributes turn the planetary disk into an extremely well-defined
probe which periodically scans
a band on the stellar surface covering
in latitude. The orbital period of CoRoT-2b translates into
an orbital angular velocity of 0.002 deg/s.
As a result, close to the center of the stellar disk,
the planet moves
0.1 deg in longitude
during a 32 s time bin of the CoRoT data.
By order of magnitude, these values determine the attainable
surface resolution of the CoRoT observations of CoRoT-2.
![]() |
Figure 1: Section of the CoRoT lightcurve of CoRoT-2, purged of outlying values (see Sect. 2). The transits appear as recurrent dips overlying the modulation by the stellar rotation. Transits are numbered from the beginning of the CoRoT observations; our analysis focuses on transit No. 56. |
Open with DEXTER |
2 Observations and data analysis
The CoRoT (Auvergne et al. 2009) lightcurve of CoRoT-2 observed from 2007 May 16 to 2007 October 15 continuously samples 31 stellar rotations. After five days of observations the satellite's sampling cadence of CoRoT-2 was switched from 512 s to 32 s, resulting in a lightcurve that covers 79 planetary transits at this high time resolution.
Our analysis started out from the lightcurve data as delivered by the CoRoT pipeline. The pipeline sorts out defective data (e.g. due to crossings of the South Atlantic Anomaly, SAA) and performs a background subtraction. As a first step we added all three CoRoT photometry channels (red-green-blue) into a single lightcurve, because the individual channels are more affected by instrumental instabilities than their summed signal.
Similar to AL08, we then removed outlier points,
which show a pronounced non-normal distribution.
To this end, we
computed
the standard deviation
in several
narrow intervals,
yielding
.
We discarded all points deviating more than
from a boxcar-smoothed copy of the lightcurve,
thus rejecting nearly 2% of all points.
Finally, also following AL08, we corrected for the
contamination by a neighbouring
object;
a subinterval
of the resulting lightcurve is shown in Fig. 1.
The transits occur during different stellar rotation phases, i.e. at different levels and slopes of the lightcurve. In order to compare the transit lightcurves, it is convenient to normalize them to a common level. To this end, we carried out a linear interpolation of points adjacent to each transit and divided all transit lightcurves by their interpolating function. One lightcurve normalized in this way is shown in Fig. 2.
We note that this procedure introduces a systematic error that depends on the
spot coverage of the stellar disk visible during the transit.
If, e.g., the disk is covered by dark spots not occulted by the planet,
the transit depth will not yield the true ratio of radii
.
Instead, the stellar radius
will be underestimated relative to the planet, because
part of the stellar disk is dark and essentially invisible.
Given the typical amplitude of CoRoT-2's lightcurve of about 4% peak-to-peak,
this introduces a comparable uncertainty on the planetary radius
deduced from a single transit.
Consequently, this uncertainty also affects our spot size estimate
of the next section. We will discuss the influence of activity-induced lightcurve variations
on the determination of planetary parameters in a forthcoming paper.
![]() |
Figure 2:
Normalized lightcurve
during ``transit 56'', as a function of time from the transit center.
Gray symbols indicate the unbinned measurements, the red line shows them averaged in 224 s time bins with 1 |
Open with DEXTER |
3 Transit lightcurve modelling
To determine a spot configuration compatible with the deformations of a transit lightcurve, we selected one particular transit occuring close to JD 2454335.0 and referred to as ``transit 56''. Rendered in Fig. 2, it shows the most pronounced and isolated ``bump'' of the whole time series of CoRoT-2, suggesting a relatively narrow spot occulted close to the disk center. Additionally, the symmetry of the bump indicates a spot geometry that is largely symmetric in longitude.
3.1 The model
To model a transit lightcurve, we decompose the stellar surface into roughly square elements.
The integrated stellar flux is computed by summing the fluxes of all visible elements, respecting
their projected area
and the limb darkening (e.g. Wolter et al. 2005).
To model the planetary transit, all surface elements
occulted by the planetary disk at a given phase are removed
from the sum, thus treating the planet as a circular disk without any intrinsic emission. The surface resolution needs to be sufficient to
analyze the densely sampled movement of the planetary disk over the
stellar surface. We choose a surface grid with 750 elements at the equator,
yielding 178 868 elements in total
and a surface resolution of roughly
near the
equator. Also, this high resolution is required to compute lightcurve
models sufficiently free of jitter,
not exceeding a few 10-4 in our models.
Next, we adjusted the limb darkening to optimize the fit
to the transit ingress and egress as well as to a tentative lower envelope of
transit 56 and the surrounding transits.
We adopt a linear limb-darkening law with
;
we note, however, that this has little influence on our determined spot parameters,
since the spot occultation in transit 56 occurs close to the disk center.
We compute the position of the planetary disk on the stellar disk and the corresponding stellar rotation phase using the orbital parameters and planetary size given by AL08, as well as
a stellar rotation period of
d,
determined by a periodogram analysis and
compatible with the value of Lanza et al. (2009),
adopting the lightcurve maximum at JD 2 454 242.67 as the phase zero point.
To this end, we introduce Cartesian coordinates in units of the stellar radius
with their origin at the center of the star and
the observer located on the x-axis.
The stellar rotation axis is assumed to
coincide with the z-axis, while the axis of the planetary orbit is
tilted by
in the xz-plane.
Both the stellar rotation and the planetary orbit are adopted as right-handed
around the z-axis.
Using these coordinates,
and
describe the position of the center of the planetary disk projected onto the
stellar disk visible at the stellar rotation phase
.
The following relations apply:
where i and a describe the inclination and the semimajor axis of the planetary orbit, respectively.
where





![]() |
(3) |
with






3.2 Transit mapping of a starspot
Our initial tests showed that the given transit lightcurves do not significantly constrain the spot shape and contrast. Hence, we limit our models to circular spots or circle segments. Furthermore, we tentatively adopt two values for the spot-to-photosphere contrast: As a ``dark spot'' scenario we choose a spot flux of 30% of the photospheric flux. Based on Bouchy et al. (2008)'s effective temperature of CoRoT-2a of 5625 K and Planck's law at a wavelength of 6000 Å, this corresponds to a spot approximately 1200 K cooler than the photosphere. This is motivated by spot temperatures found for other highly active G and K stars (e.g. Strassmeier & Rice 1998; O'Neal et al. 2004). As a ``bright spot'' we adopt a value of 75% of the photospheric flux. This roughly represents the average contrast of large spot groups on the Sun at visible wavelengths (Walton et al. 2003; Albregtsen et al. 1984). Our two spot contrast scenarios are comparable to those adopted for CoRoT-2a by Lanza et al. (2009).
Defining the spot contrast and shape in this way, three
free parameters describe a given spot, namely the spot radius r as well as
its central longitude
and colatitude
.
However, as illustrated by Fig. 3,
the longitude is closely confined by the transit lightcurve,
irrespective of the spot contrast and colatitude.
This is due to the well defined maximum of the transit bump.
Thus only two undetermined parameters remain to optimize
the fit of transit 56's lightcurve: r and
.
We calculated model lightcurves for grids in these two parameters;
the resulting goodness-of-fit
is shown in Fig. 4.
![]() |
Figure 3:
Fit quality of our transit lightcurve models
as a function of spot longitude, shown for the models BC (blue stars), BN (red plus symbols) and DS (black curve) of Table 1.
Our adopted limiting value of
|
Open with DEXTER |
![]() |
Figure 4:
Fit quality |
Open with DEXTER |
To calculate
we rebinned the lightcurve
into 224 s bins, and estimated the resulting errors
assuming a Gaussian error propagation:
.
Here, M=3 is the number of free model parameters, namely the spot longitude, colatitude and radius.
N=8 gives the number of time bins used to constrain the modelled spot;
see the interval limits indicated in Fig. 2.
Reducing
in this way is non-unique because of the
significant correlation of the spot radius and colatitude,
illustrated by the slanted contours of Fig. 4.
Tentatively we choose a limit of
for models judged as compatible with the observations.
As an example, Fig. 2 shows a solution yielding
.
3.3 Results
As the -contours in the upper panel of Fig. 4 show,
for the ``bright spot'' scenario, the central spot colatitude
is confined to
,
i.e. inside the stellar surface belt transited
by the planet. The spot radii are slightly smaller than, or comparable
to, the size of the planetary disk:
.
This scenario is illustrated by Fig. 5 and the exemplary solutions BN, BC and BS
(``bright north, central and south '') of Table 1.
For the ``dark spot'' scenario, on the other hand, only spot centers away from the center
of the transit belt yield proper fits to the transit lightcurve:
and
.
This is illustrated by the lower panel of Fig. 4
which also shows that the resulting spot radii are smaller than
in the ``bright spot'' case (
).
Spots with centers outside the transit band also yield feasible solutions. An example is illustrated by the red arc in Fig. 5, showing solution DEQ of Table 1. As illustrated by the figure, concerning area and longitude extent they do not differ significantly from solutions with centers inside the band. We do not discuss them further since their radii do not describe the transited extension and area of the spot appropriately.
![]() |
Figure 5: A starspot on CoRoT-2a occulted by the planet during transit 56, as reconstructed by transit mapping. The black and gray circles represent the planetary disk and the spot for our ``bright spot'' scenario, i.e. adopting a spot flux of 75% relative to the photosphere. The purple arcs illustrate the northern- and southernmost solutions for this spot flux. The red arc illustrates a ``dark spot'' solution (DEQ in Table 1), see text for discussion. |
Open with DEXTER |
Table 1: Parameters of characteristic spot solutions discussed in the text; longitudes and colatitudes are given for the spot center.
4 Discussion
The shape of the transit lightcurves of CoRoT-2 exhibits highly significant variations between different transits. This indicates the ubiquitous presence of starspots in the surface regions transited by the planet. Furthermore, as illustrated by Fig. 1 and the analysis of Lanza et al. (2009), the overall lightcurve of CoRoT-2 continuously changes in amplitude during the complete CoRoT observations. This shows that the stellar surface of CoRoT-2a persistently evolves on timescales shorter than its rotation period. Such a fast spot evolution is interesting physically and makes CoRoT-2a a favourable case for the study of stellar activity and, potentially, a landmark system.
We concentrate our analysis on a single planetary transit whose lightcurve shows a pronounced and isolated bump close to the transit center. Assuming this bump is caused by a circular starspot, we determine parameter ranges for this spot that reproduce the observed transit lightcurve.
While, similar to Pont et al.'s analysis of HD189733, the spot contrast is only weakly constrained, the spot longitude and radius are closely confined by the transit lightcurve. The spot thus reconstructed on CoRoT-2b is comparable in extent with large spot groups on the Sun which cover up to about 1% of the solar surface (Baumann & Solanki 2005; Norman 2005).
Given the nearly
inclination of CoRoT-2b's orbit,
the transit-covered belt on its surface lies close to the equator.
Our analysis proves that CoRoT-2a, like the Sun, exhibits spots in this region.
Such well-constrained latitude measurements
near the equator are difficult or impossible with other surface reconstruction
methods like Doppler imaging. Using long-term transit observations,
this may allow one to study activity cycles analogous to the solar butterfly diagram.
Also, concerning possible indications of differential rotation on CoRoT-2a
(Lanza et al. 2009), the spots in the transit-covered belt
could supply additional information.
Solar umbrae have diameters up to about 10 Mm, corresponding to one degree in heliographic coordinates; penumbrae reach approximately twice this size. Our study indicates that a surface resolution of potentially better than one degree can be achieved for host stars of eclipsing planets when applying transit mapping to low-noise and fast-sampled lightcurves. Thus, such observations offer a new opportunity to study ``solar-like'' surface structures on other stars; they may for example permit one to measure the umbra/penumbra contrasts of their spots.
Acknowledgements
U.W. acknowledges financial support from DLR, project 50 OR 0105.
References
- Albregtsen, F., Joras, P. B., & Maltby, P. 1984, SoPh, 90, 17 [NASA ADS] (In the text)
- Alonso, R., Auvergne, M., Baglin, A., et al. 2008, A&A, 482, L21 [NASA ADS] [CrossRef] [EDP Sciences] (AL08) (In the text)
- Auvergne, M. P., Bodin, L., Boisnard, J. T., et al. 2009, A&A, accepted (In the text)
- Baumann, I., & Solanki, S. K. 2005, A&A, 443, 1061 [NASA ADS] [CrossRef] [EDP Sciences]
- Bouchy, F., Queloz, D., Deleuil, M., et al. 2008, A&A, 482, L25 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Charbonneau, D., Winn, J. N., & Latham, D. W. 2006, ApJ, 636, 445 [NASA ADS] [CrossRef] (In the text)
- Lanza, A. F., Pagano, I., Leto, G., et al. 2009, A&A, 493, 193 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Norman, P. 2005, JBAA, 115, 194 [NASA ADS]
- O'Neal, D., Neff, J., Saar, S., & Cuntz, M. 2004, AJ, 128, 1802 [NASA ADS] [CrossRef] (In the text)
- Pont, F., Gilliland, R. L., Moutou, C., et al. 2007, A&A, 476, 1347 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Rabus, M., Alonso, R., & Belmonte, J. A. 2009, A&A, 494, 391 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Rouan, D., Baglin, A., Copet, A., et al. 1998, in Earth, Moons and Planets, 81, 79 (In the text)
- Scharmer, G. B., & Gudiksen, B.V. 2002, Nature, 420, 151 [NASA ADS] [CrossRef] (In the text)
- Strassmeier, K., & Rice, J. 1998, A&A, 330, 685 [NASA ADS] (In the text)
- Silva-Valio, A. 2008, ApJ, 683, L179 [NASA ADS] [CrossRef] (In the text)
- Walton, S. R., Preminger, D. A., & Chapman, G. A. 2003, SoPh, 213, 301 [NASA ADS] (In the text)
- Winn, J. N., Holman, M. J., & Roussanova, A. 2007, 2007, ApJ, 657, 1098 [NASA ADS] [CrossRef]
- Wolter, U., Schmitt, J. H. M. M., & van Wyk, F. 2005, A&A, 435, 261 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
All Tables
Table 1: Parameters of characteristic spot solutions discussed in the text; longitudes and colatitudes are given for the spot center.
All Figures
![]() |
Figure 1: Section of the CoRoT lightcurve of CoRoT-2, purged of outlying values (see Sect. 2). The transits appear as recurrent dips overlying the modulation by the stellar rotation. Transits are numbered from the beginning of the CoRoT observations; our analysis focuses on transit No. 56. |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Normalized lightcurve
during ``transit 56'', as a function of time from the transit center.
Gray symbols indicate the unbinned measurements, the red line shows them averaged in 224 s time bins with 1 |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Fit quality of our transit lightcurve models
as a function of spot longitude, shown for the models BC (blue stars), BN (red plus symbols) and DS (black curve) of Table 1.
Our adopted limiting value of
|
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Fit quality |
Open with DEXTER | |
In the text |
![]() |
Figure 5: A starspot on CoRoT-2a occulted by the planet during transit 56, as reconstructed by transit mapping. The black and gray circles represent the planetary disk and the spot for our ``bright spot'' scenario, i.e. adopting a spot flux of 75% relative to the photosphere. The purple arcs illustrate the northern- and southernmost solutions for this spot flux. The red arc illustrates a ``dark spot'' solution (DEQ in Table 1), see text for discussion. |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.