# Four winters of photometry with ASTEP South at Dome C, Antarctica

N. Crouzet, E. Chapellier, D. Mékarnia, A. Agabi, Y. Fanteï-Caujolle, L. Abe, J.-P. Rivet, F.-X. Schmider, F. Fressin, E. Bondoux, et al.

► To cite this version:

N. Crouzet, E. Chapellier, D. Mékarnia, A. Agabi, Y. Fanteï-Caujolle, et al. Four winters of photometry with ASTEP South at Dome C, Antarctica. Astronomy & Astrophysics - A&A, 2018, 619, pp.A116. 10.1051/0004-6361/201732565. hal-02305118

HAL Id: hal-02305118

https://hal.science/hal-02305118v1

Submitted on 6 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Four winters of photometry with ASTEP South at Dome C, Antarctica*

N. Crouzet<sup>1,2</sup>, E. Chapellier<sup>3</sup>, T. Guillot<sup>3</sup>, D. Mékarnia<sup>3</sup>, A. Agabi<sup>3</sup>, Y. Fanteï-Caujolle<sup>3</sup>, L. Abe<sup>3</sup>, J.-P. Rivet<sup>3</sup>
F.-X. Schmider<sup>3</sup>, F. Fressin<sup>4</sup>, E. Bondoux<sup>3</sup>, Z. Challita<sup>5</sup>, C. Pouzenc<sup>6</sup>, F. Valbousquet<sup>7</sup>, D. Bayliss<sup>8</sup>, S. Bonhomme<sup>3</sup>,
J.-B. Daban<sup>5</sup>, C. Gouvret<sup>3</sup>, and A. Blazit<sup>3</sup>

Received 29 December 2017 / Accepted 28 June 2018

## ABSTRACT

**Context.** Dome C in Antarctica is a promising site for photometric observations thanks to the continuous night during the Antarctic winter and favorable weather conditions.

**Aims.** We developed instruments to assess the quality of this site for photometry in the visible and to detect and characterize variable objects through the Antarctic Search for Transiting ExoPlanets (ASTEP) project.

**Methods.** Here, we present the full analysis of four winters of data collected with ASTEP South, a 10 cm refractor pointing continuously toward the celestial south pole. We improved the instrument over the years and developed specific data reduction methods.

**Results.** We achieved nearly continuous observations over the winters. We measure an average sky background of 20 mag arcsec⁻<fnref n="2" /> in the 579–642 nm bandpass. We built the lightcurves of 6000 stars and developed a model to infer the photometric quality of Dome C from the lightcurves themselves. The weather is photometric 67.1 ± 4.2% of the time and veiled 21.8 ± 2.0% of the time. The remaining time corresponds to poor quality data or winter storms. We analyzed the lightcurves of σ Oct and HD 184465 and find that the amplitude of their main frequency varies by a factor of 3.5 and 6.7 over the four years, respectively. We also identify 34 new variable stars and eight new eclipsing binaries with periods ranging from 0.17 to 81 days.

**Conclusion.** The phase coverage that we achieved with ASTEP South is exceptional for a ground-based instrument and the data quality enables the detection and study of variable objects. These results demonstrate the high quality of Dome C for photometry in the visible and for time series observations in general.

**Key words.** methods: observational – methods: data analysis – techniques: photometric – site testing – stars: variables: delta Scuti – binaries: eclipsing

## 1. Introduction

Dome C in Antarctica offers exceptional conditions for time-series observations thanks to the continuous night during the Antarctic winter and very favorable atmospheric conditions. Dome C is located on a summit of the high Antarctic plateau at an altitude of 3233 m, 1100 km away from the coast, at coordinates 75°06′S – 123°21′E. The site is equipped with the French-Italian Concordia station, which was built between 1999 and 2005 and runs all year long since 2005. The station is designed to host various scientific experiments, and can host about 80 people in summer and 15 people for winter-overs. Site testing for astronomy began in the early 2000s and revealed exceptional conditions: a clear sky most of the time including during

the winter, low wind speeds, a dry atmosphere, and low scintillation (Aristidi et al. 2003, 2005; Lawrence et al. 2004; Ashley et al. 2005; Geissler & Masciadri 2006; Kenyon et al. 2006). Experiments also found a free-atmosphere seeing above a boundary layer of about 30 m (Agabi et al. 2006; Trinquet et al. 2008; Aristidi et al. 2009). Night-time astronomical observations started in 2006. They revealed a high duty cycle for spectroscopy and showed that high quality time-series photometry was achievable from this site (Mosser & Aristidi 2007; Strassmeier et al. 2008; Briguglio et al. 2009). The continuous coverage during the winter is beneficial to the search for and study of variable stars and transiting exoplanets.

We developed the Antarctic Search for Transiting ExoPlanets (ASTEP) project to search for and characterize transiting exoplanets and to determine the quality of Dome C for photometry in the visible (Fressin et al. 2005). The project started in 2005 and two instruments were installed at Dome C. The main instrument, ASTEP 400, is a 40 cm Newton telescope designed and

<footnote n="1">¹ Instituto de Astrofísica de Canarias, C. Vía Láctea s/n, 38205 La Laguna, Tenerife, Spain
e-mail: ncrouzet@iac.es</footnote>

<footnote n="2">² Dept. de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain</footnote>

<footnote n="3">³ Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Laboratoire Lagrange, CS 34229, 06304 Nice Cedex 4, France</footnote>

<footnote n="4">⁴ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA</footnote>

<footnote n="5">⁵ Observatoire Midi-Pyrénées, 14 avenue Edouard Belin, 31400 Toulouse, France</footnote>

<footnote n="6">⁶ Concordia Station, Dome C, Antarctica</footnote>

<footnote n="7">⁷ Optique et Vision, 6 bis avenue de l'Esterel, 06160 Juan-les Pins, France</footnote>

<footnote n="8">⁸ Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK</footnote>

<footnote>* Lightcurves are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/619/A116</footnote>
built by our team to perform high precision photometry under the extreme conditions of the Antarctic winter (Daban et al. 2010). ASTEP 400 was installed in 2010 and observed nominally over four winters. The results are presented in Abe et al. (2013), Mékarnia et al. (2016), Crouzet et al. (2016), Chapellier et al. (2016) and the lightcurves of more than 300 000 stars are publicly available<fnref n="1" />. The 2017 and 2018 winter campaigns were dedicated to the photometric monitoring of β Pic during the transit of the Hill sphere of its planet β Pic b (Mékarnia et al. 2017).

A first instrument, ASTEP South, was installed in 2008 to acquire preliminary data for the project, to assess its feasibility, and to identify the challenges of photometric observations from Dome C. ASTEP South observed over four winters from 2008 to 2011. The results from a preliminary analysis of the first winter are reported in Crouzet et al. (2010). In this paper, we present the full analysis of the ASTEP South data.

## 2. Observations

ASTEP South consists of a 10 cm refractor, a front-illuminated 4096 × 4096 pixel CCD camera, and a simple mount in a thermally enclosed enclosure. The instrument is completely fixed and points continuously toward the celestial south pole during the Antarctic winter, which minimizes jitter noise and airmass variations during the observations. This setup leads to stars moving circularly on the CCD with a one-sidereal day period and to elongated point spread functions (PSF). The field of view is 3.88 × 3.88°<fnref n="2" />, corresponding to a pixel size of 3.41 arcsec on the sky. The exposure time is 29 s with 10 s overheads between each exposure. The PSF full width half maximum (FWHM) is nominally 2 pixels, not taking into account the elongation due to rotation. At the border of the field, the rotation rate of the stars on the CCD is 0.15 pixel s⁻<fnref n="1" />, corresponding to a PSF elongation of 4.33 pixel during one exposure. The observing bandpass is equivalent to a large R band (600–900 nm). The enclosure is thermalized to –20 °C and fans are mounted inside to homogenize the temperature. The enclosure is closed by a double glass window to minimize temperature fluctuations (the outside temperature varies between –50 and –80 °C during the winter at Dome C). Software programs were developed by our team to control the camera, run the acquisitions, transfer and save the data. This instrumental setup is presented in more detail in Crouzet et al. (2010). ASTEP South was installed at the Concordia station in January 2008 and observed during the Antarctic winters 2008, 2009, 2010, and 2011. Spurious reflections presumably due to moonlight were identified on images of the first campaign. We installed a 40 cm long baffle around the enclosure’s window in July 2009 to minimize parasitic light reaching the CCD. The data quality is evaluated daily by a software program that performs a basic data reduction, runs a preliminary analysis, and sends the results by email to the team in Nice. With no moving parts, a constant field of view, and dedicated software programs, this first ASTEP instrument observed in automatic mode during almost all the winters and required only limited maintenance.

The number of science frames collected by ASTEP South over the four observing seasons is given in Table 1. With a size of 32 MB per image, this yields 27.8 terabytes (TB) of science data. The amount of calibrated images is approximately the same because a bias frame was taken between each science frame. This yields a total of about 60 TB of data. The limited internet bandpass at Concordia prevents a full data transfer during the winter. Instead, the data are stored on external hard drives that

<table><caption>**Table 1.** Number of science images and amount of data acquired by ASTEP South.</caption><thead><tr><th>Observing season</th><th>Number of science frames</th><th>Total size (terabytes)</th></tr></thead><tbody><tr><td>2008</td><td>125 003</td><td>4.0</td></tr><tr><td>2009</td><td>223 756</td><td>7.1</td></tr><tr><td>2010</td><td>241 026</td><td>7.7</td></tr><tr><td>2011</td><td>279 295</td><td>8.9</td></tr><tr><td>Total</td><td>869 080</td><td>27.8</td></tr></tbody></table>

Notes. Calibration frames are not included here.

are shipped from Concordia to Nice at the end of the winter, and a backup copy is kept at Concordia.

## 3. Data reduction

### 3.1. Star catalog

The ASTEP South target star catalog is a combination of UCAC3, UCAC2, *Tycho-2*, and APASS. These catalogs were obtained from VizieR at the Centre de Données Astronomiques de Strasbourg<fnref n="2" /> and from the American Association of Variable Star Observers Photometric All-Sky Survey website<fnref n="3" />, via requests centered on the celestial south pole (RA = 0°, Dec = –90°) with a search radius of 2°. Our main catalog is UCAC3. Most of the bright stars are flagged as unreliable data; we use the *Tycho-2* entries for these stars. *Tycho-2* does not add any stars. We also noticed that a small fraction of stars in our images were missing in UCAC3 but were present in UCAC2; we use the UCAC2 entries for these stars. Finally, APASS is used to gather complementary information on the magnitudes of the targets. We keep all the magnitude and flag information available in these catalogs. For multiple entries, defined as having a separation of less than 7 arcsec (2 pixels in the ASTEP South images), only the brightest star is kept; the other ones may be catalog artifacts or stars that would not be resolved by ASTEP South. We use the UCAC3 *fmag* magnitude defined in the bandpass 579–642 nm as our reference magnitude, the *UCmag* defined in the same bandpass for stars present only in UCAC2, and the Sloan i' magnitude for stars present only in APASS.

Due to strong image vignetting (~25% at the edges) and to PSF elongation, the limit magnitude of stars that are visible in the ASTEP South images is lower at the border than at the center of the field. The limit flux $F_{lim}$ is well represented by a linear function after a certain distance from the center:

$$F_{lim} = \begin{cases} F_{lim,c} & \text{if } r \le r_0 \\ F_{lim,c}(1 - a(r - r_0)) & \text{if } r > r_0 \end{cases} \quad (1)$$

where $F_{lim,c}$ is the limit flux at the center of the image, $r$ is the radial distance from the center in pixels, $r_0 = 300$ pixels, and $a = 4.284 \times 10^{-4}$ pixel⁻<fnref n="1" />. This corresponds to a 1.5 magnitude difference between the center and the edges, which have a limit magnitude *fmag*<sub>lim</sub> of 16 and 14.5 respectively. We used this function to select only stars below the limit magnitude after calculating their nominal radial position in the images.

In total, the ASTEP South catalog comprises 5954 target stars including 5830 from UCAC3, 91 from *Tycho-2* (almost all the stars with *fmag* < 9.8), 29 from UCAC2, and four from

<footnote n="1">¹ https://astep-vo.oca.eu/</footnote>

<footnote n="2">² http://vizier.u-strasbg.fr/</footnote>

<footnote n="3">³ http://www.aavso.org/apass</footnote>
APASS. A visual inspection of these targets superimposed on a high quality ASTEP South image shows that this catalog is nearly complete for our purpose.

## 3.2. Image calibration

Calibration images were taken during the observations, with a sequence of one bias between every science image, and one dark every ten science images. Although very consuming in terms of disk space, the large number of bias was useful for a sanity check, for example to measure the bias level as a function of the camera's electronics temperature. We removed most of the oscillations shown in Fig. 10 of Crouzet et al. (2010) by improving the thermalization after the first campaign, but non-periodic temperature fluctuations remain. To correct for this, we computed the median value of each bias and interpolated them linearly to obtain the bias level for each science image. A “masterbias” and a “masterdark” were created for each day by stacking frames together with successive median filters. The dark current itself is negligible and was not detected. The median level of the darks is lower than that of the bias by about 50 analog to digital units (ADU) for an unexplained reason, and the median level of the science images matches that of the darks.

Around noon local time, the Sun is high enough (below the horizon) to slightly increase the sky background (Crouzet et al. 2010). Although not detectable by naked eye, this illuminates enough the CCD to create “twilight” images, which we used to build flat-fields. We used images with median values between 8000 and 30 000 ADU (the pixel dynamic range is 65 535 ADU). Between a few tens and two hundred such images were recorded daily except in June where the Sun remains too low. Thanks to their motion around the CCD, stars could be removed by stacking images together and applying successive medians. We built one flat-field for each day with the following procedure: bright images were calibrated from bias and dark, corrected from hot pixels, normalized by their median in the 512×512 pixel central region, and divided into sets of 15 images. We interleaved the different sets instead of grouping consecutive images, so that images in the same set were optimally spaced in time to have stars as far away as possible between images. Images in each set were stacked together with a median filter. We obtained up to seven intermediate flat fields. We stacked them together with an outlier resistant mean to create one flat-field. We found that this method was significantly more efficient at removing stars than applying a unique median filter to all the images together. Finally, these daily flat-fields were stacked together with a median to create “masterflats” for each winter. This method yielded high quality masterflats, as shown in Fig. 1.

Random failures of the camera shutter began after two winters of operation at +5 °C, corresponding to ~393 000 open-close cycles. This issue became more severe in 2010 and 2011. We increased the temperature to +30 °C around the shutter for the 2011 campaign, which improved its reliability but did not solve the issue completely. Science images, darks, and bias were affected (Fig. 2). We developed a calibration method to correct for these effects. The masterbias and masterdark did not provide a satisfactory calibration, so we built a third type of calibration image. First, we cleaned the masterbias from residual light: we removed a low frequency component from each row and replaced it by a linear function between both sides (this linear trend also appeared on bias unaffected by shutter issues). Then, we adjusted the median value of the masterbias and of the corrected masterbias to match that of the masterdark; the medians were calculated in 200 pixel wide boxes located at the corners.

<figure><img src="image_2.png" /><figcaption>Fig. 1. Masterflat (left) and same image after subtracting the low frequency component (right) for 2009. The vignetting is ~25%. The color scale is arbitrary.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Fig. 2. From left to right: masterdark, masterbias, masterbias with shutter effect corrected, and mastermix for July 2nd, 2010. The color scale is arbitrary.</figcaption></figure>

Finally, we stacked these three images together with a median to create a “mastermix” (Fig. 2). This median can be seen as a conditional statement: each part of the mastermix corresponds to one of these images. We found that it was a good representation of the shutter effects seen on the science images, and provided a satisfactory calibration.

Our science images were calibrated as:

$$C = \frac{(I - M - b) \times g}{F} \qquad (2)$$

where C is the calibrated image, I is the raw image, M is the masterdark or the mastermix in case of shutter issues, b is the bias variations due to temperature fluctuations (taken around its median value for each day), g is the gain of the CCD, and F is the masterflat. We measured the gain during the characterization of the ASTEP cameras in the lab using the approach of Leach et al. (1980), Leach (1987). We took pairs of images at different illumination levels spanning the full dynamical range of the CCD. Then, the gain is given by $g = 2m/v$, where m is the mean of the image and v is the variance of the differential image (we neglected the read-out noise in this calculation). The $v(m)$ relation was well represented by a linear function below 31 360 ADU and by a third order polynomial above. This yielded a constant gain of 2 e-/ADU below 31 360 ADU and an increasing gain from 2 to 2.2 e-/ADU between 31 360 and 60 000 ADU. Finally, we corrected for hot pixels. They were present in the masterdarks but not in the masterbias. Pixels in the masterdarks that exceeded their value in the masterbias by more than 150 ADU were flagged for all the images of that day and were replaced by the average of the four adjacent pixels. An example of calibrated image is shown in Fig. 3.

## 3.3. Field matching algorithm

In each science image, we identified point sources using the FIND procedure in IDL, with a threshold intensity of five times the sky standard deviation. These point sources were compared
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Example of a calibrated image. The field of view is $3.88 \times 3.88^{\circ^2}$. The color scale is arbitrary.</figcaption></figure>

to the ASTEP South catalog using a home-made field matching algorithm. The star coordinates in the catalog (RA, Dec) are converted into image coordinates $(x, y)$. Then, we perform iterations on pairs of stars using the 20 brightest non saturated stars. For each pair in the ASTEP South image and each pair in the catalog, we calculate the geometrical transformation that matches one pair to the other, skipping iterations where both stars do not have approximately the same separation. Then, we apply this transformation to the 400 brightest stars and evaluate the matching by a $\chi^2$, where $\chi$ is the median distance between a star in the transformed image and the closest star in the catalog. The best matching is found very quickly and yields a typical precision of 0.4 pixel. Then, $\chi^2$ is minimized using the simplex method in order to better estimate the transformation parameters. This yields a precision of 0.2 pixel. Finally, we calculate and correct for distortions using a regression in $x$ and $y$ up to the third order, including crossed terms. The final precision on the star positions is 0.1 pixel.

### 3.4. Photometric algorithm

The ASTEP South PSFs are particularly challenging for high-precision photometry. We developed our own photometric algorithm to take into account the elongated PSFs, and their elongation directions and positions that vary in time as the stars rotate on the CCD. The aim was to maximize the measured flux while reducing the sky background, readout noise, hot pixels, cosmic rays, and contamination by surrounding stars compared to circular apertures.

We computed the apertures for each star in each image. The apertures consist of two half-circles joined by a rectangle which length and orientation depend on the position of the star relatively to the celestial south pole (Fig. 4). The aperture is quasi-circular for stars close to the south pole and is the most elongated at the edges of the field. To measure the stellar flux, we divide each pixel within the aperture into a grid of 11 × 11 subpixels, project the aperture onto this grid, and sum the flux at the subpixel level; this is equivalent to weighting each pixel by the

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Subimage around the star (left) and photometric aperture computed by the ASTEP South photometric algorithm (right) for stars of magnitude 8.03 (top left), 10.6 (top right), 11.04 (bottom left), and 12.08 (bottom right). The color scale is arbitrary.</figcaption></figure>

fraction of its area within the aperture. The sky background is calculated in an annulus of 11 pixels starting 5 pixels outside the star aperture, multiplied by the aperture area, and subtracted to the measured flux.

To calculate the optimal aperture size for each star, defined by the radius of the half-circles, we ran the photometric algorithm on a typical day with radii ranging from 0.5 to 10 pixels by steps of 0.05. The lightcurves were calibrated by a reference lightcurve, calculated from the median of the 30 non saturated brightest stars using an aperture radius of 5 pixels. For each star, we kept the aperture radius that yielded the smallest root mean square (RMS). In the case of neighboring stars, the resulting aperture was generally too large because it contained several stars, and the optimal aperture appeared as a local RMS minimum rather than as a global minimum. This was corrected by forcing the aperture radius to remain within ±2 pixels of the typical aperture for this star's magnitude. In the end, each star had its own optimized aperture size, which remained the same for all the photometric reduction.

After reducing all the days of a winter, we created a single lightcurve file for each star with the flux measurements and information on the aperture position, centroid position, and local sky background. We also built a reference lightcurve using 60 bright stars with low RMS that we identified from one typical day. We normalized these 60 lightcurves by their respective medians, and averaged them using an outlier resistant mean with a $3\sigma$ cut-off. The resulting reference lightcurve was also normalized by its median.

### 3.5. Correction of systematic effects

We divided the target stars' lightcurves by the reference lightcurve and binned them per 400 s intervals (6.7 min, approximately ten points). Then, we identified systematic variations and applied several corrections. The largest variations were due to the rotation of stars on the CCD and were on the order of a few percent, sometimes up to 10%. They formed a pattern with a one-sidereal day period, which was different for each star and was usually repeatable. However, several complications arose. First, we sometimes modified the pointing during the winter, for example after improving the thermal control or after other interventions on the instrument; this changed the one-day patterns. Second, each pattern was split into two or more patterns that are shifted in flux, corresponding to different states of the FWHM (see Sect. 4.2). Third, we noticed flat-field variations over the winter possibly due to mechanical deformations and temperature variations inside the box, and we built a flat-field for
each image calibration method (with or without shutter correction). In the end, we used between one and four flat-fields over the winter. These pointing, FWHM, and flat-field effects were mixed together and yielded variable one-sidereal day patterns. To recover them, we split the data into sets that were as consistent as possible according to these parameters; in practice, each set corresponds to a specific time interval during the winter.

Then, we recovered the patterns for each star and each set individually: we fold the star lightcurve over a given time interval at a one-sidereal day period, split it into 100 bins, calculate the median flux for each bin, interpolate them linearly, and divide the lightcurve by this function (Fig. 5). We applied this procedure twice. In some cases, because of poor weather, large FWHM variations, or small data sets, the patterns could not be recovered. The 2008 and 2011 data were divided into one and four sets respectively and were well corrected. The 2009 and 2010 data were much less stable and were divided into 12 and 16 sets respectively; some of them were completely discarded, whereas several large sets comprised most of the good data and were well corrected.

We converted fluxes into magnitudes and normalized them by their mean magnitude. From this point, we worked with the magnitude residuals. Then, we calibrated each star by a set of reference stars chosen as follows. We divided each lightcurve by all others and sorted the resulting lightcurves by their mean absolute deviation (MAD). The stars yielding the lowest MADs were used as reference stars and their lightcurves were averaged to build a reference lightcurve. We repeated this procedure for each star (so in the end each star had its own set of reference stars). We used three reference stars for the 600 brightest stars ($f_{\text{mag}} < 11.9$) and ten reference stars for fainter stars.

We corrected for remaining systematic variations. We gathered the lightcurves in a star, epoch array with stars sorted by ascending magnitudes. Systematic trends appeared along the epochs and were common to stars of similar magnitudes. We corrected for them on a row by row (epoch by epoch) basis by computing a median smoothing over 200 stars and subtracting it to the initial row. This method yielded a better correction than Systematic Removal (SYSREM, Tamuz et al. 2005) for these data at this stage of the pipeline, because it is insensitive to variable stars, which are often interpreted as systematic effects by SYSREM. Our method provided more uniform lightcurves overall. After these corrections, we computed the average standard deviation for each epoch $\sigma_e$ using stars of magnitudes between 12 and 13.5.

### 3.6. Data selection and final lightcurves

In the process, we flagged outliers departing by more than $4\sigma$ from each lightcurve, where $\sigma$ is the standard deviation of the lightcurve. We flagged entirely the epochs with more than 50% of outliers, as well as the epochs with less than 500 detected stars. Considering individual stars, we flagged entirely the days with more than 50% of outliers, because the remaining points were generally off the regular one-sidereal day pattern. We flagged the epochs with a median sky background greater than 500 ADU. We flagged the epochs with $\sigma_e > 2\%$ because they correspond to the tail of the $\sigma_e$ distribution (see Sect. 5). For 2009, a limit at 3% would be more appropriate considering all the data; however the distribution is bi-modal, one mode corresponding to a nominal FWHM (~2 pixels) and the other to a larger FWHM (>3 pixels). To keep only nominal data and for consistency with the other years, we kept the limit $\sigma_e < 2\%$. We also flagged the epochs at which more than 30% of the 1000 brightest stars had bad data (these epochs were often of poor

quality for all the stars), and those at which more than 10% of the 100 brightest stars had bad data (we found that this was a good proxy for data taken around full Moon). Finally, we flagged lonely epochs (epochs with less than 20% of good epochs within their 40 closest neighbors), which were often missed by other flagging constraints. The flagging information was used in the analysis of the Dome C photometric quality (Sect. 5) but the flagged epochs were discarded in the final lightcurves and in the search for variable objects (Sect. 7). The final lightcurves of the 5954 target stars are publicly available on the ASTEP lightcurve database website<fnref n="4" /> and at the CDS. Examples of lightcurves of variable objects are shown in Fig. A.1.

### 3.7. Lightcurve precision

We calculated the precision of the four-winter lightcurves over 6.7 minute bins. We computed the true lightcurve RMS ($\sigma$) and the point-to-point RMS ($\sigma_p$) which does not take into account correlated noise. We also computed theoretical estimates of the stellar photon noise, readout noise, sky background noise, scintillation noise, and their quadratic sum, first considering circular PSFs and circular apertures, then taking into account the PSF elongation and specific aperture of each star ($\sigma_t$). The scintillation noise is obtained from Eq. (1) of Young (1967; see also Dravins et al. 1998; Ryan & Sandler 1998). By analyzing turbulence profiles, Kenyon et al. (2006) suggest a potential gain of 3.6 at Dome C compared to Chilean sites. We did not apply this factor here: the scintillation noise remains well below other noise sources for ASTEP South, and using larger telescopes or studying a few bright stars would be necessary to confirm this gain. The lightcurve RMS as a function of stellar magnitude is shown in Fig. 6 and the average values for various magnitudes are reported in Table 2.

At the faint end, the point-to-point RMS is well explained by the theoretical RMS for elongated PSFs. At magnitudes 10 and 12, the point-to-point RMS is larger than the theoretical RMS by a factor 1.5 and 1.2 respectively, indicating that additional sources of uncorrelated noise is present and affect mostly bright stars. The true RMS is larger than the theoretical RMS by a factor 2.6, 1.8, 1.6, and 1.3 at magnitudes 10, 12, 14, and 15.5 respectively, indicating that correlated noise also remains in the lightcurves. For most stars, the true RMS of the four-winter lightcurves is within a factor of two of the theoretical estimates. This could be explained by the particular issues related to the challenging observing mode of ASTEP South. Overall, this shows that good photometry over very long periods is achievable at Dome C.

## 4. Instrumental parameters and sky background

Besides stellar fluxes, we extracted several parameters in order to evaluate the functioning of the instrument and the quality of Dome C for photometry. We present the results in this section, which complements the analysis of Crouzet et al. (2010) for the first winter.

### 4.1. Pointing variations

We investigated the mechanical stability of ASTEP South by measuring the motion of the celestial south pole on the CCD (Fig. 7). Within one day, the south pole is stable within one pixel: the standard deviation from its average position is typically 0.6

<footnote n="4">⁴ https://astep-vo.oca.eu/</footnote>
<figure><img src="image_6.png" /><figcaption>Fig. 5. Correction of the one-sidereal day variations. *From left to right*: winters 2008 to 2011. One star is given as an example for each winter; the stars' IDs are indicated. The lightcurves are phase-folded at a one sidereal day period. *First row*: initial lightcurve. *Second row*: split into different data sets. Each set is represented by a color and the recovered pattern is shown as a line of the same color. *Third row*: lightcurve after division by the patterns. *Fourth row*: identification of outliers, indicated in red. *Fifth row*: lightcurve with outliers removed. *Sixth row*: final lightcurve after applying all the photometric correction and data selection steps.</figcaption></figure>

and 0.5 pixel in *x* and *y*, respectively. The long term motion over the winters is given in Table 3 and averages to 0.18 and 0.09 pixel per day in *x* and *y*, respectively. This long term motion is oriented close to the north-south axis which is the vertical axis of the instrument; this is particularly evident in 2008 and 2010. This motion may correspond to a bending

of the instrument under its own weight over time, modulated by thermal fluctuations inside the box affecting its deformation coefficients. A motion of the ice below the instrument is also possible. Overall, this motion is very small and would be completely negligible for any instrument equipped with a guiding system operating at Dome C.
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** RMS of the four winter lightcurves as a function of stellar R magnitude. The theoretical stellar photon noise, readout noise, sky background noise, scintillation noise, and their quadratic sum are shown as long-dashed, dot-dashed, dotted, dashed, and plain lines respectively, in the case of circular PSFs and apertures. The theoretical noise for each star taking into account the PSF and aperture elongation is shown as black dots. The true and point-to-point RMS are shown as red and blue dots, respectively.</figcaption></figure>

<table><caption>**Table 2.** RMS of the four-winter lightcurves over 6.7 minute bins for different magnitudes, in percent.</caption><thead><tr><th><i>R</i> mag</th><th>10</th><th>12</th><th>14</th><th>15.5</th></tr></thead><tbody><tr><td>σ<sub>t</sub></td><td>0.16</td><td>0.45</td><td>1.54</td><td>4.08</td></tr><tr><td>σ<sub>p</sub></td><td>0.23</td><td>0.55</td><td>1.79</td><td>4.14</td></tr><tr><td>σ</td><td>0.41</td><td>0.81</td><td>2.41</td><td>5.47</td></tr></tbody></table>

**Notes.** σ<sub>t</sub> is the theoretical RMS, σ<sub>p</sub> is the point-to-point RMS, and σ is the true RMS.

<table><caption>**Table 3.** Motion of the south pole on the CCD in the *x* and *y* directions for each winter (Δ*x* and Δ*y*) calculated over the longest time interval without repointing the instrument (Δ*t*).</caption><thead><tr><th></th><th>Δx<br>(pixel unit)</th><th>Δy<br>(pixel unit)</th><th>Δt<br>(day)</th></tr></thead><tbody><tr><td>2008</td><td>25</td><td>10</td><td>128</td></tr><tr><td>2009</td><td>13</td><td>8</td><td>86</td></tr><tr><td>2010</td><td>25</td><td>11</td><td>115</td></tr><tr><td>2011</td><td>33</td><td>16</td><td>211</td></tr></tbody></table>

**Notes.** The pixel scale is 3.41 arcsec per pixel.

## 4.2. Point spread function variations

We measured the PSF variations using the 20 brightest stars in the 1000×1000 pixel central region of the image, where the elongation is negligible. We fitted the PSFs with a Gaussian function to derive the FWHM. This calculation was done by the day-to-day analysis software program and was used as an instrument control parameter. The nominal FWHM for ASTEP South is 2 pixels, and its variations over the four winters are shown in Fig. 8. Thermal stability is critical to a stable FWHM. For example, during the first part of the 2009 winter, the FWHM suffered extreme variations because of a fan failure. The recorded temperature inside the enclosure (around the probe) varied between −22 and −18 °C but the actual variations and gradients were probably much stronger, inducing thermal deformations of the instrument.

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Position of the celestial south pole on the CCD during the winters 2008, 2009, 2010, and 2011. For each winter, only the longest period without repointing the instrument is shown. Increasing time is represented by black, purple, blue, green, yellow, orange, and red. The orientation of the CCD (north, south, east, west) is indicated.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** PSF FWHM during the winters 2008, 2009, 2010, and 2011. Short-term variations (timescale < 1 day) are due to seeing variations, and long-term variations (timescale > 1 day) are due to poor thermal stability causing mechanical deformations of the instrument.</figcaption></figure>

After restarting the fans, the recorded temperature remained stable within 0.2 °C of the setting point and the FWHM remained stable between 2 and 2.5 pixels. Long-term variations (timescale > 1 day) in 2008 and 2010 are also correlated with thermal variations. Because any intervention on the instrument during the Antarctic winter is extremely delicate, the thermal control was not modified during these winters. Under stable thermal conditions, the FWHM recorded by ASTEP South is directly correlated with the atmospheric seeing, which yields short-term variations (timescale < 1 day). This has been shown for the 2008 winter by Crouzet et al. (2010) and is also evident during the second half of the 2009 winter and during the 2011 winter, where increases of the FWHM correspond to poor seeing conditions.

As for any photometric instrument, periods affected by large FWHM variations result in poor photometry. In particular,
<figure><img src="image_10.png" /><figcaption>Fig. 9. Sky background and related parameters for the winters 2008 to 2011 and for all the winters combined, from left to right. From top to bottom: sky background as a function of time, sky background distribution, sky background as a function of Sun altitude, epoch RMS σ_e (see text) as a function of sky background, and σ_e as a function of Sun altitude. The sky background histograms were calculated in flux units and are displayed in magnitude units. For the last three rows, we computed an outlier resistant mean using 50 bins across the abscissa range (red lines).</figcaption></figure>

FWHM variations are correlated with the one-sidereal day patterns, with different patterns for FWHM below and above a limit of either 2.5 or 3 pixels. Above this value, the pattern can be highly variable from day to day and was not always recovered, in particular in 2009 and 2010. These data were generally eliminated from the lightcurves.

## 4.3. Sky background

The sky background was computed in each image using the IDL SKY procedure, corrected for the zenith angle, converted from photo-electrons to instrumental magnitude, and calibrated to a fmag magnitude (579–642 nm bandpass). We performed this latter calibration by fitting a linear function between the stars’ instrumental magnitudes and their UCAC3 fmag magnitudes, independently for each winter. The sky background and

its correlations with other parameters are shown in Fig. 9. The first row shows the sky background as a function of time. The Sun increases the sky background every day around noon and the Moon cycle is clearly visible, as already reported in Crouzet et al. (2010). The second row displays the histograms of sky values. They are composed of three regions: a broad distribution centered around the nominal dark sky value (~20 mag arcsec⁻²), a region between 19.5 and 18 mag arcsec⁻² that corresponds to the illumination around full Moon, and a long tail for values brighter than 18 mag arcsec⁻² that corresponds to the illumination by the Sun around noon. To measure the average dark sky background at Dome C, we combined the distributions of the four winters and fitted the main peak with a Gaussian function. We did this calculation in flux units and converted the mean and 1σ values back to magnitudes. We find a dark sky background at Dome C of 20.0±0.3 mag arcsec⁻² in the 579–642 nm
bandpass. The uncertainty is partly due to real variations of the dark sky background: in 2011, the sky is darker around mid-winter (~20.5 mag arcsec⁻²) and brighter closer to the equinoxes (~19.5 mag arcsec⁻²). This trend is less clear for the other win- ters due to less consistent and less complete coverage. For com- parison, the dark sky background in mag arcsec⁻² in the R band is 20.4 at Mauna Kea<fnref n="5" />, 20.9 at Paranal (Patat 2008), 21.0 at La Palma<fnref n="6" />, and 21.2 at Cerro Tololo and La Silla (Kriisciunas et al. 2007). Thus, the sky background at Dome C is comparable to that of the best observing sites.

The sky background as a function of the Sun altitude $h_{\text{Sun}}$ is shown in the third row of Fig. 9. A limit $h_{\text{Sun}} < -15^{\circ}$ can be used to define dark sky conditions, whereas the sky back- ground increases quickly for $h_{\text{Sun}} > -13^{\circ}$. Finally, we calculated the lightcurve dispersion at each epoch $\sigma_e$ using the stars in the magnitude interval 12–13.5, which represents the photometric quality at that epoch (see Sect. 5). We show $\sigma_e$ as a function of the sky background and Sun altitude in the fourth and fifth rows, respectively. We find that $\sigma_e$ starts increasing noticeably when the sky background becomes brighter than 18 mag arcsec⁻². We find that $\sigma_e$ remains at its average level while $h_{\text{Sun}} < -12^{\circ}$ and increases steeply for $h_{\text{Sun}} > -10^{\circ}$. These numbers give the con- ditions under which good quality photometry can be achieved at Dome C as seen by our 10 cm diameter instrument.

## 5. Photometric quality of Dome C

The long, continuous, and consistent database provided by ASTEP South is well suited to measure the weather quality dur- ing the winter at Dome C. The clear sky fraction for the winter 2008 was measured by Crouzet et al. (2010) using the number of detected stars in the field of view. Here, we use the quality of the lightcurves themselves to derive the photometric quality. In the following calculations, we used only the subset of stars with R magnitudes between 12 and 13.5 (1090 control stars). The final lightcurves are in the form of magnitude residuals. We calculated the standard deviation of the residuals at each epoch $\sigma_e$ using the control stars. This parameter provides a direct measurement of the data quality at each epoch: large dispersions in the residuals correspond to low quality epochs. As shown in Fig. 12, the $\sigma_e$ distribution is composed of a bulk and a tail containing the high and low quality data, respectively. To differentiate the bulk and the tail, we built a model for $\sigma_e$ and compared it to the mea- sured distribution. In the following, $\sigma_{e,m}$ and $\sigma_{e,t}$ refer to the measured and theoretical $\sigma_e$ respectively. We calculated $\sigma_{e,t}$ in several steps as described next.

### 5.1. Theoretical noise

First, we calculated the theoretical noise that was expected for each star at each epoch taking into account the Poisson noise, read-out noise, and sky background noise. We did not consider scintillation noise because it is negligible in the magnitude range 12–13.5 (see Fig. 6). Each star had its own photometric aperture which remained fixed in size for the four years of data. The frac- tion of the PSF included in the aperture depended on its FWHM, which varied in time. We calculated this fraction for each star

at each epoch by integrating numerically the PSF over the pho- tometric aperture. The PSF was taken as a two-dimensional Gaussian function elongated in one direction, its FWHM was that measured by ASTEP South, and the elongation was calcu- lated as a function of the distance of the star from the celestial south pole. The photometric aperture was the same as used in the reduction pipeline.

We calculated the nominal flux for each star as the median of the flux measured under excellent conditions, which we defined as non-flagged epochs falling within the 20% largest number of stars and the 20% smallest FWHM. Then, we calculated the the- oretical flux for each star at each epoch by scaling the nominal flux to the fraction of PSF within the aperture, and used it to derive theoretical estimates of the stellar Poisson noise, read-out noise, and sky background noise, relatively to the stellar signal. In these calculations and in the following, we included a 10% uncertainty on the measured FWHM and sky background values to account for measurement errors.

As shown in Sect. 3.7, these theoretical noises do not explain fully the measured lightcurves’ RMS. Thus, we added unidenti- fied noise sources and split their contributions into two compo- nents $\sigma_w$ and $\sigma_r$, defined as:

$$ \sigma_w^2 = \sigma_p^2 - \sigma_t^2 \quad \text{and} \quad \sigma_r^2 = \sigma^2 - \sigma_p^2 \qquad (3) $$

where $\sigma_t$ is the theoretical RMS, $\sigma_p$ is the point-to-point RMS, and $\sigma$ is the true RMS. $\sigma_w$ corresponds to a white (uncorre- lated) noise component and $\sigma_r$ as a red (correlated) noise com- ponent. We obtained $\sigma_w$ and $\sigma_r$ for each star directly from their lightcurve (see Fig. 6); these lightcurves were cleaned from bad data points including a cut-off in $\sigma_e$ so they correspond to good weather conditions. We expressed the noise $\sigma_{s,e}$ for star s at epoch e as:

$$ \sigma_{s,e}^2 = \sigma_{\text{phot}, s, e}^2 + \sigma_{\text{read}, s, e}^2 + \sigma_{\text{sky}, s, e}^2 + \sigma_{w, s}^2 + \sigma_{r, s}^2 \qquad (4) $$

where $\sigma_{\text{phot}, s, e}$, $\sigma_{\text{read}, s, e}$, and $\sigma_{\text{sky}, s, e}$ are the stellar photon noise, read-out noise, and sky background noise, and $\sigma_{w, s}$ and $\sigma_{r, s}$ are the additional white and red noise contributions.

Finally, we computed the theoretical noise for epoch e, $\sigma_{e,t}$, as the median of $\sigma_{s,e}$ over the control stars. As shown in Fig. 10, the resulting $\sigma_{e,t}$ distribution is well centered on the $\sigma_{e,m}$ distri- bution but is very narrow. This is because the correlation of $\sigma_{s,e}$ between stars (the correlative nature of $\sigma_{r,s}$) was not accounted for: $\sigma_{r,s}$ is only an amplitude. Increasing the random errors on the FWHM and sky background measurements broadened only slightly the $\sigma_{e,t}$ distribution. In fact, more complex correlated noise was present and was not captured yet.

### 5.2. Addition of a blind correlative component

A complete model of the variations affecting the ASTEP South lightcurves was beyond the scope of this work. Instead, we used an empirical approach and included a correlative component in the $\sigma_r$ parameter. We modeled the correlations in the fol- lowing way. We divided the lightcurves into segments of about 30 minute duration. We associated each segment to a correlated noise value randomly drawn from a normal distribution of mean zero and standard deviation $\beta$, noted $N(\beta)$. These values were the same for all stars, meaning that the lightcurve residuals var- ied in the same way for all stars at a given epoch. For star s, we defined the red noise distribution $\mathcal{R}$ along the segments by:

$$ \mathcal{R} = \sigma_{r,s} (1 + N(\beta)) \qquad (5) $$

<footnote n="5">5 http://www.gemini.edu/sciops/telescopes-and-sites/observing-condition-constraints/optical-sky-background</footnote>

<footnote n="6">6 http://www.ing.iac.es/Astronomy/observing/conditions/skybr/skybr.html/#vsol</footnote>
<figure><img src="image_11.png" /><figcaption>Fig. 10. Construction of a theoretical model for the distribution of σ<sub>e</sub>. The blue line shows a model taking into account the theoretical noises plus the white and red noise values, as defined in Eq. (4). The brown line shows a model after introducing a correlative component between stars, as defined in Eq. (6). The red line shows a model after introducing correlations with the FWHM and the sky background, as defined in Eq. (10). The black line is the measured distribution σ<sub>e,m</sub>. Only the 2011 winter is shown as an example.</figcaption></figure>

The red noise value σ<sub>R,s</sub> for star s was taken as the realization of R on the segments and remained constant on each segment. The theoretical noise σ<sub>s,e</sub> for star “s” at epoch “e” was then defined as:

$$
\sigma_{s,e}^2 = \sigma_{\text{phot},s,e}^2 + \sigma_{\text{read},s,e}^2 + \sigma_{\text{sky},s,e}^2 + \sigma_{\text{w},s}^2 + \sigma_{\mathcal{R},s,e}^2 \quad (6)
$$

where σ<sub>R,s,e</sub> was then correlated between stars. The theoretical
noise for epoch e, σ<sub>e,t</sub>, was taken as the median of σ<sub>s,e</sub> over the
control stars.

The amount of correlation encoded in the β parameter was
unknown; the larger it was, the broader was the σ<sub>e,t</sub> distribution.
We adjusted β by matching the left side of the σ<sub>e,t</sub> histogram
to that of σ<sub>e,m</sub>. Indeed, the broadening on the left side of the
σ<sub>e,m</sub> histogram was due to noise that was present even in good
weather conditions, whereas the broadening on the right side was
also due to bad weather, high sky background, etc. The resulting
σ<sub>e,t</sub> distribution reproduced well the bulk of the σ<sub>e,m</sub> distribu-
tion, as shown in Fig. 10.

In this model, the sizes of the segments did not matter as
long as the stars were affected in the same way. The β parame-
ter represented the amount of correlated noise at a given time
in the star domain, whereas the correlated noise is generally
measured on a lightcurve, that is for a given star in the time
domain (e.g., Gillon et al. 2006). The value of β depended on
the number of correlated stars: if less stars were correlated,
β had to be larger to reproduce the σ<sub>e,m</sub> distribution; differ-
ent groups of stars may also have had different correlations.
Algorithms such as SYSREM (Tamuz et al. 2005), trend fil-
tering algorithm (TFA, Kovács et al. 2005) or others using a
principal component analysis aim at identifying and removing
such correlations (although the use of SYSREM on the ASTEP
South data did not improve the lightcurves compared to our
own procedures). Our approach was equivalent to convoluting
the σ<sub>e,t</sub> distribution from Eq. (4) with a Gaussian function, but
related this broadening to an amount of correlated noise that
could be included in the noise budget. A more general approach
based on the covariance matrix could also have been used
(Pont et al. 2006).

5.3. Correlations with the FWHM and sky background

The tail of the σ<sub>e,m</sub> histogram contained the bad weather data,
but also included other noise sources which effects were not
symmetrical with respect to the mode of the histogram and
which were not included in our model. We investigated poten-
tial correlations of σ<sub>e,m</sub> with the various recorded parameters
for data taken under good weather conditions, and found clear
correlations with the PSF FWHM and the sky background.

In spite of the correction by reference stars, FWHM varia-
tions affected the photometry more than predicted in Sect. 5.1.
To derive the correlation between σ<sub>e,m</sub> and the FWHM, we
divided the FWHM measurements into ten intervals, and calcu-
lated the median FWHM and the median σ<sub>e,m</sub> in each interval.
In this process, we removed the σ<sub>e,m</sub> outliers using a ten itera-
tion sigma-clipping procedure with a threshold of 1.5 times the
standard deviation, to ensure that the correlation was extracted
from good weather data only. The extracted correlation function
was denoted *f* and is shown for the winters 2010 and 2011 in
Fig. 11. We interpolated *f* linearly to each measured FWHM
(after including a 10% uncertainty on the measurements) and
calculated the associated noise σ<sub>F,e</sub> by:

$$
\sigma_{F,e}^2 = f^2 - f_0^2 \quad (7)
$$

where $f_0$ was the minimum of $f$. This ensured that $\sigma_{F,e}$ only accounted for additional noise caused by FWHM variations. We used the same $\sigma_{F,e}$ value for all stars at a given epoch.

We used the same method for the sky background noise.
We calculated the correlation function *g* using 22 intervals, with
smaller sizes at lower sky values. *g* is shown for the winters 2010
and 2011 in Fig. 11. We interpolated *g* at each measured sky
value (after including a 10% uncertainty on the measurements)
and calculated the associated noise *σ*<sub>S,e</sub> by:

$$
\sigma_{S,e}^2 = g^2 - g_0^2 \tag{8}
$$

where g₀ was the minimum of g. This empirical model was close
to the theoretical estimate from Poisson noise but with some dif-
ferences at low sky values, possibly related to image calibration.

We also had to modify the theoretical relative photon noise,
read-out noise, and sky noise terms in Eqs. (4) and (6) because
part of their contributions was then contained in σ<sub>F,e</sub> and σ<sub>S,e</sub>.
To account for this, we calculated the stellar fluxes for a unique
value of FWHM (that of f<sub>0</sub>), take a unique sky value (that of
g<sub>0</sub>), and recomputed the theoretical noise estimates. These cal-
culations were still different for each star because each had its
own photometric aperture, but they were equal for all epochs.
The FWHM and sky background variations were also responsi-
ble for a certain amount of red noise, so the contribution of σ<sub>R,s,e</sub>
had to be reduced. Some amount may have been considered as
white noise (larger Poisson noise for larger FWHM or sky back-
ground), but we reduced only σ<sub>R,s,e</sub> and not σ<sub>w</sub> (the final result
was equivalent). We defined a new red noise distribution R' as:

$$
\mathcal{R}' = \alpha \mathcal{R} \tag{9}
$$

with α < 1 and ℛ defined in Eq. (5). The corresponding red noise
value was noted σ<sub>R',s,e</sub>. Our final model was:

$$
\sigma_{s,e}^2 = \sigma_{\text{phot},s}^2 + \sigma_{\text{read},s}^2 + \sigma_{\text{sky},s}^2 + \sigma_{\text{w},s}^2 + \sigma_{\text{F},e}^2 + \sigma_{\text{S},e}^2 + \sigma_{\mathcal{R}',s,e}^2 \quad (10)
$$

As before, the theoretical noise for epoch e, σ<sub>e,t</sub>, was taken as the median of σ<sub>s,e</sub> over the control stars. We derived α by matching the left sides of the σ<sub>e,t</sub> and σ<sub>e,m</sub> histograms. The final σ<sub>e,t</sub> distributions for each winter are shown in Fig. 12. We also defined
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Correlations of $\sigma_e$ with the PSF FWHM (top) and the sky background (middle and bottom) for the winters 2010 (left) and 2011 (right). Measurements are indicated as black dots. The correlations extracted from good weather data are shown as red lines. The sky background correlation plots are shown in two different scales.</figcaption></figure>

the ratio $r$ of the joint contribution of $\sigma_{F,e}$ and $\sigma_{S,e}$ in the red noise budget compared to the unknown sources $\sigma_{R'}$ as:

$$r = \sqrt{1/\alpha^2 - 1} \qquad (11)$$

The two parameters of our model, $\beta$ and $\alpha$, and the ratio $r$ are reported for each winter in Table 4. For the four winters, $\sigma_{F,e}$ and $\sigma_{S,e}$ have a larger contribution than $\sigma_{R',s,e}$ ($r > 1$), so they are the main sources of correlated noise. For 2011, $r$ is lower than for the other winters, which is consistent with the stable FWHM during that winter. Following this approach, we could in principle have split entirely the red noise into its various components. However, no obvious correlations were found between $\sigma_{e,m}$ and other recorded parameters that were independent of the FWHM and the sky background, so we kept $\sigma_{R',s,e}$ as a global representation of the remaining components.

## 5.4. Weather statistics

We attributed to veiled weather the values of $\sigma_{e,m}$ in excess of our model. We defined the “photometric” and “veiled” weather fractions as the areas of $\sigma_{e,m}$ below and above $\sigma_{e,t}$ respectively, as illustrated in Fig. 12. We repeated this procedure 100 times with different random distributions for the definition of the correlation segments, for $\mathcal{R}$, and for the uncertainties on the FWHM and sky background measurements. We took the mean and standard deviation of the results to obtain the photometric and veiled weather statistics and their uncertainties. In addition, small fractions of the lightcurves were of reasonable quality but were not well corrected by the initial one-sidereal day period variation removal. They were kept but flagged early in the pipeline, and did not have a $\sigma_e$ value. Whether these data were affected by

<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Normalized distribution of $\sigma_e$ for the winters 2008, 2009, 2010, and 2011 (black lines). Periods of photometric and veiled weather as defined by our model are indicated in red and orange respectively.</figcaption></figure>

<table><caption>**Table 4.** Model parameters $\beta$ and $\alpha$ for the epoch noise distribution $\sigma_{e,t}$ and ratio of the red noise contributions $r$.</caption><thead><tr><th>Winter</th><th>β</th><th>α</th><th>r</th></tr></thead><tbody><tr><td>2008</td><td>0.22 ± 0.09</td><td>0.50 ± 0.01</td><td>1.73</td></tr><tr><td>2009</td><td>0.17 ± 0.07</td><td>0.38 ± 0.01</td><td>2.44</td></tr><tr><td>2010</td><td>0.39 ± 0.09</td><td>0.37 ± 0.01</td><td>2.49</td></tr><tr><td>2011</td><td>0.21 ± 0.02</td><td>0.65 ± 0.01</td><td>1.18</td></tr></tbody></table>

**Notes.** The values and uncertainties were obtained after repeating the procedure 100 times.

instrumental effects or poor weather conditions could not be inferred by our method. We classified these periods manually by interpreting the variations using the recorded parameters, the daily reports of the winter-over crew, and knowledge of the instrument behavior. We adopted a relative uncertainty of 10% for the corresponding fractions. These periods tend to increase the veiled weather fraction, and some of them correspond to bad weather. Furthermore, bad data were also eliminated during the photometric calibration process, and the instrument was stopped during “white-outs” (winter storms lasting one to a few days). We included these periods and adopted a relative uncertainty of 10% for the corresponding fractions. Finally, in the first part of 2009 (JD < 2455032), the instrument and lightcurves suffered extreme variations and our model did not reproduce the measured distribution accurately, so we derived the weather statistics only from the second part of the winter.

After these adjustments, we computed the final weather statistics for the winters at Dome C, which are reported in Table 5. The cumulative distributions of photometric weather, veiled weather, and bad weather plus white-out fractions are shown in Fig. 13. We also calculated the average weather statistics for the winter at Dome C by taking the mean and standard deviation over the four winters. On average, we find a photometric weather fraction of 67.1 ± 4.2%, a veiled weather fraction of 21.8 ± 2.0%, a bad weather fraction of 3.3 ± 2.0%, and a white-out fraction of 7.9 ± 6.0%. The photometric weather fraction is comparable to that of the best astronomical sites at temperate latitudes (see Crouzet et al. 2010, for a more detailed comparison). These statistics are consistent with those obtained from a different method based on the number of detected stars using the
<table><caption>**Table 5.** Weather statistics for the winter at Dome C as observed with ASTEP South.</caption><thead><tr><th></th><th>2008</th><th>2009</th><th>2010</th><th>2011</th><th>Average</th></tr></thead><tbody><tr><td colspan="6">Excluding white-outs:</td></tr><tr><td>Photometric</td><td>75.6 ± 2.2</td><td>69.6 ± 1.2</td><td>75.3 ± 2.7</td><td>70.9 ± 1.2</td><td>72.9 ± 3.1</td></tr><tr><td>Veiled</td><td>23.1 ± 2.2</td><td>25.8 ± 1.1</td><td>22.3 ± 1.7</td><td>23.3 ± 1.1</td><td>23.7 ± 1.5</td></tr><tr><td>Bad</td><td>1.3 ± 0.1</td><td>4.5 ± 0.5</td><td>2.3 ± 0.2</td><td>5.8 ± 0.6</td><td>3.5 ± 2.0</td></tr><tr><td colspan="6">Including white-outs:</td></tr><tr><td>Photometric</td><td>63.0 ± 2.2</td><td>64.7 ± 1.2</td><td>72.5 ± 2.6</td><td>68.0 ± 1.2</td><td>67.1 ± 4.2</td></tr><tr><td>Veiled</td><td>19.2 ± 1.8</td><td>24.0 ± 1.1</td><td>21.5 ± 1.6</td><td>22.4 ± 1.1</td><td>21.8 ± 2.0</td></tr><tr><td>Bad</td><td>1.1 ± 0.1</td><td>4.2 ± 0.4</td><td>2.2 ± 0.2</td><td>5.5 ± 0.6</td><td>3.3 ± 2.0</td></tr><tr><td>White-out</td><td>16.6 ± 1.7</td><td>7.1 ± 0.7</td><td>3.8 ± 0.4</td><td>4.1 ± 0.4</td><td>7.9 ± 6.0</td></tr></tbody></table>

**Notes.** All numbers are in percent.

<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Cumulative normalized distribution of $\sigma_e$ for the winters 2008, 2009, 2010, and 2011 (black lines). Periods of photometric and veiled weather are indicated in red and orange respectively, and include adjustments for periods without a $\sigma_e$ value (see text). Periods of bad weather and white-out are indicated in white.</figcaption></figure>

2008 winter (Crouzet et al. 2010). In particular, both methods are in excellent agreement for the 2008 winter.

## 5.5. Validation

We validated the results of our method by verifying the epoch quality with respect to several parameters. To infer the quality of individual epochs, we calculated the ratio $\sigma_e, m/\sigma_e, t$ and we considered that an epoch is photometric or veiled if this ratio was below or above a certain limit, respectively. This limit was calculated for each winter to match the weather statistics. Figure 14 shows $\sigma_e$ as a function of the FHWM and of the sky background for 2010 and 2011 (same as Fig. 11 now showing the quality of each epoch). As expected, the bulk of data corresponds to photometric weather, and the variations of $\sigma_e$ with the FHWM and sky background are well taken into account in differentiating the photometric and veiled weather. We also show the number of detected stars as a function of these parameters for the 2011 winter (Fig. 15). This number is a good proxy for the weather quality and was used by Crouzet et al. (2010) to derive the weather statistics at Dome C for the 2008 winter, before extracting the lightcurves. As expected, the number of detected stars is correlated with $\sigma_e$. Interestingly, the correlation has two tails,

<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Correlations of $\sigma_e$ with the PSF FWHM (top) and the sky background (middle and bottom) for the winters 2010 (left) and 2011 (right). Photometric and veiled epochs as inferred from our model are indicated in red and orange, respectively. Bad epochs are indicated in blue. The sky background plots are shown in two different scales.</figcaption></figure>

one that extends the bulk of photometric weather and one that corresponds to veiled weather. This is also the case for 2010 but is less obvious for 2008 and 2009. This reflects the sky background distribution: the photometric (lower) tail corresponds to the correlation between large sky backgrounds and $\sigma_e$ whereas the veiled (upper) tail corresponds to low sky backgrounds for which the number of stars was affected by poor weather conditions. In these plots, the photometric, veiled, and bad epochs are distributed as expected. In particular, the bulk of photometric data inferred from our model matches very well the bulk of data with a nominal number of detected stars, which confirms that they correspond to excellent weather. A few epochs have almost nominal parameters (high number of stars, low sky background, nominal FWHM) but were classified as bad. These epochs were identified as outliers before calculating $\sigma_e$ (see Sect. 3.6) and were often adjacent to or surrounded by bad weather periods;
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Correlations of the number of detected stars with $\sigma_e$ (top left), the FWHM (top right), and the sky background (bottom left and right, in two different scales) for the 2011 winter. Photometric and veiled epochs as inferred from our model are indicated in red and orange, respectively. Bad epochs are indicated in blue; most of them do not have a calculated $\sigma_e$.</figcaption></figure>

they were not suitable for photometry. Overall, these results validate our model and show that the approach used here and that of Crouzet et al. (2010) are consistent.

# 6. Duty cycle

The duty cycle of ASTEP South for the winters 2008–2011 is represented in Fig. 16: the limit due to the Sun elevation, the fraction of time where ASTEP South was operating, and the photometric, veiled, and bad weather fractions including whiteouts are shown for each day of the winter. In total, we acquired 11415.1 h of data with ASTEP South including 6699.8 h under photometric conditions and nominal functioning of the instrument. Considering only times with a Sun elevation lower than −9°, these numbers are 9451.2 h and 6308.5 h respectively. The detail for each year is reported in Table 6. The increase of the operating time from the 2008 to the 2011 winter is evident and reflects the improvements of the system (hardware and software) made over the years. The first half of the 2008 winter was hampered by numerous technical issues while setting up the instrument. A long white out period occurred at the end of August and lasted for about ten days. In 2009, the beginning of the winter was dedicated to setting up the instrument and improving the thermalization; data acquired during this time were not used in the lightcurves. In 2010, we lost two weeks of data taken in June, either in the data transfer between computers and external disks or while shipping the disks from Concordia to Nice. In 2011, the installation was fast and the system acquired high quality data nearly continuously for the whole winter, with almost no human interventions. This demonstrates that a simple instrument can operate almost continuously for several winters at Dome C. In addition, as noted in Crouzet et al. (2010), periods of bad weather are generally concentrated in white-out episodes lasting from one to a few days; between these episodes, high quality photometric observations can be achieved for extended periods of time. Overall, the weather statistics and duty cycle at Dome C demonstrate the high quality of this site for continuous photometric observations during the winter.

# 7. Study of variable stars and eclipsing binaries

## 7.1. Search for periodic signals

We searched for periodic signals in the lightcurves using the box least square (BLS) algorithm (Kovács et al. 2002). This program searches for square shapes representing transit-like events, and also detects stellar variations. We performed the search over a wide period range ($0.4 < P < 100$ days) which we divided into 15 intervals with sizes following approximately a logarithmic scale. We built a frequency comb for each interval using a constant spacing in $\Delta f/f$, a number of frequencies set by the expected duty cycle, and we ran BLS in each interval separately. During the search, we applied locally a refined comb with nine times more frequencies, first around the forty best frequencies, then around the best frequency. This adaptive comb yielded a more precise determination of the best frequency. We ran BLS twice. After the first run, the histogram of periods showed peaks that were centered around integer or half-integer values, indicating harmonics of the “day-night” cycle (the sky background increases around noon even during the winter, see Sect. 4.3). Another peak around 29 days was due to the Moon cycle. We excluded periods around these peaks for the second BLS run. Finally, we combined the 15 intervals into three period ranges: 0.4–1 day, 1–10 days, and 10–100 days. For each lightcurve, we compared the results obtained for the intervals within each range and kept only the period that yielded the largest power in the BLS spectrum. Every lightcurve ended up with a best period, but most of them did not contain a valid signal. We ranked them by decreasing signal to noise ratios (S/Ns), where the signal was calculated from the amplitude and duration of the variations and the noise was obtained from the standard deviation of the lightcurve residuals. Then, we inspected the lightcurves visually to identify variable stars, eclipsing binaries, and transiting exoplanet candidates. Examples of lightcurves of variable objects are shown in Fig. A.1.

## 7.2. Analysis of variable objects

We detected periodic signals in the lightcurves of 211 stars. A similar search was performed by the Chinese Small Telescope ARray (CSTAR, Yuan et al. 2008; Zhou et al. 2010) from Dome A, Antarctica, using data acquired during the 2008, 2009, and 2010 winters (Wang et al. 2011, 2013, 2015; Oelkers et al. 2015; Yang et al. 2015). We cross-matched our objects with the variable stars reported by CSTAR and by the GCVS (Samus’ et al. 2017), VSX (Watson et al. 2006), and ASAS (Pojmanski 2002; Richards et al. 2012) catalogs. The results are given in Table 7. We identified 96 objects that were not reported and we analyzed their lightcurves with the Period04 software program (Lenz & Breger 2005). First, we analyzed the 2011 data alone because they were more abundant and of better quality, then we analyzed the four winters together. In the Fourier spectra, we found systematic peaks around $1.003$ and $0.034 \text{ d}^{-1}$ caused by the day-night and Moon cycles, respectively (these frequencies were not completely avoided in our BLS search probably due to the frequency refinement procedure). We eliminated stars showing variations only at these frequencies. We also discarded stars with frequencies that were common to many stars, indicating non-astrophysical origin, and stars with S/Ns that were too low to extract accurate periods and amplitudes. This left 42 objects with confirmed variability. We analyzed them year by year to distinguish variable stars and eclipsing binaries: the amplitude of variations due to stellar activity varies from year to year but remains the same for eclipsing binaries. We also analyzed the
<figure><img src="image_17.png" /><figcaption>**Fig. 16.** Daily observing time fraction for ASTEP South for the winters 2008, 2009, 2010, and 2011. The light blue and dark blue regions indicate the fraction of time with a Sun elevation lower than −9° and −13° (9° and 13° below the horizon), respectively. Periods of photometric and veiled weather as observed by ASTEP South are indicated in red and orange, respectively. Periods of bad weather including white-outs are indicated in white. Periods in yellow indicate that data were recorded but not used in the photometric analysis, either because they were taken during the instrument setup at the beginning of the winter or were lost during the data transfer.</figcaption></figure>

<table><caption>**Table 6.** Amount of data acquired by ASTEP South and subsample of data taken under photometric conditions and nominal functioning of the instrument, with or without the constraint of a Sun elevation lower than −9°.</caption><thead><tr><th rowspan="2">Winter</th><th rowspan="2">All data</th><th rowspan="2">Photometric data</th><th>All data</th><th>Photometric data</th><th rowspan="2">h<sub>Sun</sub> &lt; −9°</th></tr><tr><th>h<sub>Sun</sub> &lt; −9°</th><th>h<sub>Sun</sub> &lt; −9°</th></tr></thead><tbody><tr><td>2008</td><td>1640.5</td><td>887.0</td><td>1376.3</td><td>864.8</td><td>3030.6</td></tr><tr><td>2009</td><td>2977.1</td><td>1692.8</td><td>2453.4</td><td>1575.7</td><td>3030.9</td></tr><tr><td>2010</td><td>3155.5</td><td>1829.9</td><td>2627.7</td><td>1749.2</td><td>3031.0</td></tr><tr><td>2011</td><td>3642.0</td><td>2290.1</td><td>2993.8</td><td>2118.8</td><td>3031.0</td></tr></tbody><tfoot><tr><td>Total</td><td>11415.1</td><td>6699.8</td><td>9451.2</td><td>6308.5</td><td>12123.5</td></tr></tfoot></table>

**Notes.** The amount of time with a Sun elevation lower than −9° for the whole year at Dome C is indicated for comparison. All numbers are in hours.

lightcurves of five known objects: we investigated in detail $\sigma$ Oct and HD 184465, we improved the parameters of EN Oct, and we re-classified two CSTAR variables as eclipsing binaries.

In this study, we computed the frequency and amplitude uncertainties using the formulae proposed by Montgomery & Odonoghue (1999) providing $3\sigma$ estimates. These uncertainties correspond to ideal cases (constant noise, regular spacing, no gaps) and underestimate the true uncertainties. Based on our experience, we used the $3\sigma$ values for our uncertainties.

Finally, we counted the number of previously reported variables in the ASTEP South field of view that were missed by our search (that were not in our list of 211 potential variables). We missed 103 objects including 90 from CSTAR and 13 from GCVS, VSX, or ASAS. Some were not in our initial target list, some were too faint, and some were simply not identified as variable. CSTAR is composed of four Schmidt-Cassegrain telescopes, each with a pupil entrance aperture of 145 mm and a field of view of 4.5° in diameter (20$^{\circ}$2 in total) equipped with different filters (g, r, i, and no filter). The CSTAR target list contains around 18 000 objects. Although technical issues generally prevented the acquisition of data by the four telescopes simultaneously, each CSTAR telescope already collects about twice more flux than ASTEP South. This set up can explain that many variables were detected by CSTAR and missed by ASTEP South. The photometric data reduction pipelines are also different, and Oelkers et al. (2015) use differential images. In addition, our search for periodic signals was primarily oriented toward the detection of transiting exoplanet candidates and we used only the BLS algorithm to search for variability. Other algorithms or metrics that were used for CSTAR such as the Welch-Stetson variability index (Stetson 1996), the analysis of variance method (Schwarzenberg-Czerny 1996), or the Lomb-Scargle method (Lomb 1976; Scargle 1982) may

reveal more variables in the ASTEP South lightcurves. Besides, we restrained the search range to periods larger than 0.4 days, and most $\delta$ Scuti were not identified by their main frequency: some were found from their harmonics and some may have been missed. In contrast, the longer coverage by ASTEP South compared to CSTAR may explain that some variables were detected only by ASTEP South.

Our goal is not to detect all the variables in the ASTEP South field of view, which has already been studied extensively by CSTAR. Instead, we complement the CSTAR search and focus in the following on a few interesting objects that highlight the long photometric coverage that we achieved with ASTEP South.

### 7.2.1. $\sigma$ Oct

$\sigma$ Oct, also known as HR 7228 and HD 177482, is the brightest star in the field ($V = 5.42$). It is a $\delta$ Scuti reported first by McInally & Austin (1978) and re-observed only once by Coates et al. (1981). Its properties are poorly known as it had been observed only four nights in total before the ASTEP South and CSTAR experiments. Its known pulsation period is on the order of 0.1 day with an amplitude of 15 mmag. The ASTEP South lightcurve of $\sigma$ Oct is shown in Fig. 17; we note that this star was saturating the CCD so the measured amplitudes should be taken with caution. We analyzed the data with the Period04 software program and detected 21 frequencies with S/Ns larger than four (Fig. 18, Table B.1). Seven-teen of them are between 8.023 and 11.756 $d^{-1}$. The dominant peak is at 10.493 $d^{-1}$ and corresponds to $\delta$ Scuti-type oscillations with an average amplitude of 3.73 mmag. Interestingly, four other frequencies are between 0.612 and 2.848 $d^{-1}$ with amplitudes below 0.7 mmag, and could correspond to $\gamma$
<table><caption>**Table 7.** Summary of our search for variable objects in the ASTEP South lightcurves.</caption><tr><td>Number of target stars</td><td>5954</td></tr><tr><td>Objects with detected variability</td><td>211</td></tr><tr><td>New objects</td><td>96</td></tr><tr><td>Confirmed variability</td><td>42</td></tr><tr><td>New variable stars</td><td>34</td></tr><tr><td>New eclipsing binaries</td><td>8</td></tr><tr><td>Unconfirmed variability</td><td>54</td></tr><tr><td>Known objects</td><td>115</td></tr><tr><td>From CSTAR</td><td>95</td></tr><tr><td>From GCVS, VSX, and ASAS</td><td>20</td></tr><tr><td>Objects missed</td><td>103</td></tr><tr><td>From CSTAR</td><td>90</td></tr><tr><td>From GCVS, VSX, and ASAS</td><td>13</td></tr></table>

<figure><img src="image_18.png" /><figcaption>**Fig. 17.** Lightcurve of σ Oct as observed with ASTEP South over four winters, phase-folded at a period of 0.0953 day. The black points correspond to a time sampling of 10 min. The red line is a binning with 100 points over the phase range. The period and V magnitude are indicated on the plot.</figcaption></figure>

Dor-type pulsations. Thus, σ Oct is apparently an hybrid δ Scuti – γ Dor with very low amplitudes for the γ Dor modes. We analyzed the ten main δ Scuti frequencies year by year (Fig. 19, Table B.2). Four of them show amplitude variations. In particular, the amplitude of the dominant frequency decreases from 9.74 to 2.76 mmag (a factor of 3.5) between 2008 and 2011, and this decrease is not compensated by an equivalent increase of the amplitudes of other frequencies. Thus, the cause is not an energy transfer between oscillation modes but is an actual diminution of the total amplitude of the star's oscillations. A similar behavior was found in 603 δ Scuti out of 983 using four year lightcurves from the *Kepler* mission and no explanation is known to date (Bowman et al. 2016).

### 7.2.2. HD 184465

HD 184465 is a star of V magnitude 8.82 and spectral type A1 V, located in the δ Scuti instability strip. Its ASTEP South lightcurve and Fourier spectrum are shown in Figs. 20 and 21. We detected ten δ Scuti-type frequencies with S/Ns greater than four and we analyzed them year by year (Fig. 22, Table B.3). The 9.671 d⁻¹ frequency largely dominates the Fourier spectrum,

<figure><img src="image_19.png" /><figcaption>**Fig. 18.** Fourier spectrum of σ Oct obtained from its four winter lightcurve.</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>**Fig. 19.** Amplitude variations of the ten main δ Scuti frequencies of σ Oct from 2008 to 2011, with one point per winter. The 3σ uncertainties are shown for the main frequency and are the same for all frequencies.</figcaption></figure>

and its amplitude increases from 0.9 to 6.02 mmag (a factor of 6.7) between 2008 and 2011. This increase is not compensated by an equivalent decrease of the amplitudes of other frequencies. The amplitudes of other frequencies are generally too small to detect significant variations between years. We analyzed the 2011 data by 30-day segments (Table B.4, Fig. 22). The amplitudes of the 9.671 d⁻¹ and 6.465 d⁻¹ frequencies vary at similar rates in 2011 compared to the four winters together. Interestingly, the amplitude of the 8.691 d⁻¹ frequency increases in 2011 from 0.7 to 2.7 mmag (a factor of 3.8). Overall, this star does not show energy transfer between oscillations modes, but shows actual variations of the amplitude of its oscillations on several timescales.

### 7.2.3. New variable stars and eclipsing binaries

We identified 34 new variable stars (Table C.1). Their amplitudes are typically between 1 and 5 mmag and several have periods longer than ten days. More than half are active stars with variations induced by their rotation. In many cases, the S/N is too low to identify the precise type of variable. We also identified eight new eclipsing binaries and reclassified two CSTAR variables as eclipsing binaries (Table C.2). Six of them have orbital periods between 23 and 81 days, which are much
<figure><img src="image_21.png" /><figcaption>**Fig. 20.** Lightcurve of HD 184465 as observed with ASTEP South over four winters, phase-folded at a period of 0.1034 day. The black points correspond to a time sampling of 10 min. The red line is a binning with 100 points over the phase range. The period and V magnitude are indicated on the plot.</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Fig. 21.** Fourier spectrum of HD 184465 obtained from its four winter lightcurve.</figcaption></figure>

longer than those usually detected from the ground. In addition, we refined the period and eclipse depth of the known eclipsing binary EN Oct compared to previous studies (Knigge & Bauernfeind 1967; Otero 2004). We also detected variable stars and eclipsing binaries already reported by CSTAR and other catalogs, so we do not report them here. Figure A.1 shows examples of lightcurves of variable stars and eclipsing binaries obtained with ASTEP South. Despite challenges in the photometric reduction due to the unusual observing mode of ASTEP South, these lightcurves are of good quality, and the large amount of data enables the detection of variations over a wide range of magnitudes and periods (in Fig. A.1, the V magnitudes range from 6.52 to 15.04 and the periods range from 0.12 to 40.53 days). Overall, these detections demonstrate the benefit of the long photometric coverage that can be achieved during the winter at Dome C.

### 7.2.4. Inventory of variable objects

We gathered the new variable stars and eclipsing binaries detected by ASTEP South with those reported by CSTAR,

<figure><img src="image_23.png" /><figcaption>**Fig. 22.** Top: amplitude variations of the ten δ Scuti frequencies of HD 184465 with signal to noise ratios greater than four from 2008 to 2011, with one point per winter. Bottom: amplitude variations of three frequencies with two time intervals for each of the 2008, 2009, and 2010 winters and six 30-day intervals for the 2011 winter. The $3\sigma$ uncertainties are shown for the main frequency and are the same for all frequencies.</figcaption></figure>

GCVS, VSX, and ASAS to make an inventory of variable objects in the ASTEP South field of view. Their classification is given in Table 8. Further details on these objects can be found in Wang et al. (2011, 2013, 2015), Oelkers et al. (2015), Yang et al. (2015), Samus’ et al. (2017), Watson et al. (2006), Pojmanski (2002), Richards et al. (2012) and in Tables C.1 and C.2 of this paper.

### 7.2.5. A 75.6 day period G–M eclipsing binary

Eclipsing binaries are the most common by-products of transiting exoplanet searches. In our search for transiting objects, we detected an eclipsing binary with a 75.6 day orbital period and a 13.6 h primary eclipse (star ID: UCAC3 2-3056, fmag: 12.76; see Table C.2). These period and duration are much longer than those usually detected from the ground, which shows the benefit of the exceptional phase coverage that we achieved from Dome C. This object was not detected by Period04 because of its too low S/N but it was among the transit candidates that we identified from the BLS search and visual inspection. We conducted follow-up observations in radial velocity and photometry, which confirmed the orbital period and the presence of the eclipses. This system is a G9–M2.5 eclipsing
<table><caption>**Table 8.** Inventory of variable objects in the ASTEP South field of view (Dec < −88°) including those detected by ASTEP South and those reported by CSTAR, GCVS, VSX, and ASAS.</caption><thead><tr><th>Total number of variable objects</th><th>260</th></tr></thead><tbody><tr><td>Variable stars</td><td>209</td></tr><tr><td>VRot</td><td>96</td></tr><tr><td>Ceph</td><td>1</td></tr><tr><td>Delta Scuti</td><td>19</td></tr><tr><td>Gamma Dor</td><td>10</td></tr><tr><td>RR Lyrae</td><td>4</td></tr><tr><td>Mira</td><td>2</td></tr><tr><td>Multiperiodic</td><td>7</td></tr><tr><td>Irregular</td><td>38</td></tr><tr><td>Unclassified</td><td>32</td></tr><tr><td>Eclipsing binaries</td><td>47</td></tr><tr><td>EA</td><td>19</td></tr><tr><td>EB</td><td>10</td></tr><tr><td>EW</td><td>18</td></tr><tr><td>Transit-like</td><td>4</td></tr></tbody></table>

**Notes.** “VRot” are rotational variables and include the ACV, BY, SP, ELL classes and the “VRot” and “VRot?” of Table C.1. “Unclassified” include the “Var” of Table C.1. “Transit-like” objects are those reported by CSTAR; the transit candidates identified by ASTEP South are not reported here and are under investigation. This table does not include the 54 objects initially suspected as variable with ASTEP South that were not confirmed by a finer analysis (see Sect. 7.2).

binary with an eccentric orbit (e = 0.26) and a relatively large semi-major axis (a = 0.42 au). Thus, tidal coupling between both stars is weak and the M dwarf should behave as a single star. This object will provide an ideal benchmark to constrain the evolution models of low mass stars, and will be presented in a forthcoming paper (Crouzet et al., in prep.).

## 8. Conclusion

We presented the full analysis of four winters of observations collected with ASTEP South, a 10 cm refractor installed at Dome C, Antarctica. The instrument ran continuously for most of the winters. The functioning of the instrument and the lightcurve quality improved over the years. A main limitation was shut-ter malfunctions that started after two winters and required spe-cific temperature adjustments and image calibration methods. Another factor that affected the data quality was imperfect ther-malisation, in particular during the first half of the 2009 winter, caused by deficient air circulation leading to variable temperature and gradients inside the box. This resulted in unstable PSFs and poor quality lightcurves for these periods. By acquiring knowl-edge of the instrument’s behavior over the years and improving the experiment (hardware and software), we were able to observe longer and longer, from half of the winter in 2008 to the full win-ter with almost no intervention in 2011. As previously shown for the 2008 winter, we confirmed that the instrument was very stable over the ice and that PSF stability was achieved when the ther-malisation was nominal. We also measured the sky background magnitude and its variations over the winters.

We developed a specific data reduction pipeline to take into account the motion and elongated shape of the PSFs on the CCD. We built the lightcurves of nearly 6000 stars over four winters, which are available publicly. The experimental design was particularly sensitive to flat-field errors, which we

mitigated by removing the one-sidereal day period variations in the lightcurves. Most of the data are of good quality. The final lightcurve RMS is standard for a 10 cm instrument, and the large quantity of data yields a drastic improvement of the precision when folding the lightcurves at given periods and binning the data points. We measured the photometric quality of Dome C from the lightcurves themselves using a method that we devel-oped specifically for this work, which may be applied to other experiments.

We detected periodic signal in the lightcurves of around two hundred objects. Half of them were already reported by the CSTAR experiment. We discovered new variable stars and eclipsing binaries, and conducted a detailed study of a few objects that was enabled by the long and continuous coverage of our observations. These results demonstrate the high quality of Dome C for photometry in the visible and the technical fea-sibility of running astronomical instruments for several years at this site.

**Acknowledgements.** The ASTEP project was funded by the Agence Nationale de la Recherche (ANR), the Institut National des Sciences de l'Univers (INSU), the Programme National de Planétologie (PNP), and the Plan Pluri-Formation OPERA between the Observatoire de la Côte d'Azur and the Université de Nice-Sophia Antipolis. The logistics at Concordia is handled by the French Institut Paul-Emile Victor (IPEV) and the Italian Programma Nazionale di Ricerche in Antartide (PNRA). This work made use of the SIMBAD and VizieR databases operated at CDS, Strasbourg, France, the NASA's Astrophysics Data System bibliographic services, and the NASA's IDL astronomy library. Software: Period04 (Lenz & Breger 2005).

## References

* Abe, L., Gonçalves, I., Agabi, A., et al. 2013, A&A, 553, A49

* Agabi, A., Aristidi, E., Azouit, M., et al. 2006, PASP, 118, 344

* Aristidi, E., Agabi, A., Vernin, J., et al. 2003, A&A, 406, L19

* Aristidi, E., Agabi, K., Azouit, M., et al. 2005, A&A, 430, 739

* Aristidi, E., Fossat, E., Agabi, A., et al. 2009, A&A, 499, 955

* Ashley, M. C. B., Lawrence, J. S., Storey, J. W. V., & Tokovinin, A. 2005, in EAS Pub. Ser., eds. M. Giard, F. Casoli, & F. Paletou, 14, 19

* Bowman, D. M., Kurtz, D. W., Breger, M., Murphy, S. J., & Holdsworth, D. L. 2016, MNRAS, 460, 1970

* Briguglio, R., Tosti, G., Strassmeier, K. G., et al. 2009, Mem. Soc. Astron. It., 80, 147

* Chapellier, E., Mékarnia, D., Abe, L., et al. 2016, ApJS, 226, 21

* Coates, D. W., Halprin, L., Moon, T. T., & Thompson, K. 1981, IBVS, 2047

* Crouzet, N., Guillot, T., Agabi, A., et al. 2010, A&A, 511, A36

* Crouzet, N., Mékarnia, D., Guillot, T., et al. 2016, in Ground-Based and Airborne Instrumentation for Astronomy VI, Proc. SPIE, 9908, 99080U

* Daban, J.-B., Gouvret, C., Guillot, T., et al. 2010, in Ground-based and Airborne Telescopes III, Proc. SPIE, 7733, 77334T

* Dravins, D., Lindegren, L., Mezey, E., & Young, A. T. 1998, PASP, 110, 610

* Fressin, F., Guillot, T., Bouchy, F., et al. 2005, in EAS Pub. Ser., eds. M. Giard, F. Casoli, & F. Paletou, 14, 309

* Geissler, K., & Masciadri, E. 2006, Proc. SPIE, 6267, 626719

* Gillon, M., Pont, F., Moutou, C., et al. 2006, A&A, 459, 249

* Kenyon, S. L., Lawrence, J. S., Ashley, M. C. B., et al. 2006, PASP, 118, 924

* Knigge, R., & Bauernfeind, H. 1967, Bamberg Veroeffentlichungen der Remeis-Sternwarte, 6, 45

* Kovács, G., Zucker, S., & Mazeh, T. 2002, A&A, 391, 369

* Kovács, G., Bakos, G., & Noyes, R. W. 2005, MNRAS, 356, 557

* Krisciunas, K., Semler, D. R., Richards, J., et al. 2007, PASP, 119, 687

* Lawrence, J. S., Ashley, M. C. B., Tokovinin, A., & Travouillon, T. 2004, Nature, 431, 278

* Leach, R. W. 1987, Opt. Eng., 26, 1061

* Leach, R. W., Schild, R. E., Gursky, H., et al. 1980, PASP, 92, 233

* Lenz, P., & Breger, M. 2005, Commun. Asteroseismol., 146, 53

* Lomb, N. R. 1976, Ap&SS, 39, 447

* McInally, C. J., & Austin, R. R. D. 1978, MNRAS, 184, 885

* Mékarnia, D., Guillot, T., Rivet, J.-P., et al. 2016, MNRAS, 463, 45

* Mékarnia, D., Chapellier, E., Guillot, T., et al. 2017, A&A, 608, L6

* Montgomery, M. H., & Odonoghue, D. 1999, Delta Scuti Star Newsletter, 13, 28

* Mosser, B., & Aristidi, E. 2007, PASP, 119, 127

* Oelkers, R. J., Macri, L. M., Wang, L., et al. 2015, AJ, 149, 50
Otero, S. A. 2004, IBVS, 5532

Patat, F. 2008, A&A, 481, 575

Pojmanski, G. 2002, Acta Astron., 52, 397

Pont, F., Zucker, S., & Queloz, D. 2006, MNRAS, 373, 231

Richards, J. W., Starr, D. L., Miller, A. A., et al. 2012, ApJS, 203, 32

Ryan, P., & Sandler, D. 1998, PASP, 110, 1235

Samus', N. N., Kazarovets, E. V., Durlevich, O. V., Kireeva, N. N., & Pastukhova, E. N. 2017, Astron. Rep., 61, 80

Scargle, J. D. 1982, ApJ, 263, 835

Schwarzenberg-Czerny, A. 1996, ApJ, 460, L107

Stetson, P. B. 1996, PASP, 108, 851

Strassmeier, K. G., Briguglio, R., Granzer, T., et al. 2008, A&A, 490, 287

Tamuz, O., Mazeh, T., & Zucker, S. 2005, MNRAS, 356, 1466

Trinquet, H., Agabi, A., Vernin, J., et al. 2008, PASP, 120, 203

Wang, L., Macri, L. M., Krisciunas, K., et al. 2011, AJ, 142, 155

Wang, L., Macri, L. M., Wang, L., et al. 2013, AJ, 146, 139

Wang, S., Zhang, H., Zhou, X., et al. 2015, ApJS, 218, 20

Watson, C. L., Henden, A. A., & Price, A. 2006, SAS Annual Symp., 25, 47

Yang, M., Zhang, H., Wang, S., et al. 2015, ApJS, 217, 28

Young, A. T. 1967, AJ, 72, 747

Yuan, X., Cui, X., Liu, G., et al. 2008, in Ground-based and Airborne Telescopes II, Proc. SPIE, 7012, 70124G

Zhou, X., Wu, Z.-Y., Jiang, Z.-J., et al. 2010, Res. Astron. Astrophys., 10, 279
Appendix A: Lightcurve examples

<figure><img src="image_24.png" /><figcaption>**Fig. A.1.** Examples of phase-folded lightcurves of variable stars and eclipsing binaries obtained with ASTEP South, sorted by magnitude. The black points correspond to a time sampling of 10 min. A binning with 100 points over the phase range is shown as a red line, except for eclipsing binaries with narrow eclipses. The star ID, the period in days, and the V magnitude (UCAC3 fmag) are indicated.</figcaption></figure>
## Appendix B: Analysis of δ Scuti pulsations of σ Oct and HD 184465

<table><caption>**Table B.1.** Twenty-one main frequencies detected in the Fourier spectrum of σ Oct with their amplitudes and phases.</caption><thead><tr><th>Frequency (d<sup>-1</sup>)</th><th>Amplitude (mmag)</th><th>Phase (rad)</th></tr></thead><tbody><tr><td>10.492693 (13)</td><td>3.73 (11)</td><td>0.449 (30)</td></tr><tr><td>9.719016 (16)</td><td>3.11 (11)</td><td>2.972 (36)</td></tr><tr><td>10.741796 (18)</td><td>2.77 (11)</td><td>2.953 (40)</td></tr><tr><td>8.799398 (19)</td><td>2.61 (11)</td><td>4.475 (42)</td></tr><tr><td>10.251893 (20)</td><td>2.47 (11)</td><td>4.517 (45)</td></tr><tr><td>11.429052 (22)</td><td>2.27 (11)</td><td>2.818 (49)</td></tr><tr><td>9.361590 (28)</td><td>1.80 (11)</td><td>2.375 (61)</td></tr><tr><td>9.771399 (30)</td><td>1.65 (11)</td><td>1.183 (67)</td></tr><tr><td>10.440784 (37)</td><td>1.37 (11)</td><td>6.276 (81)</td></tr><tr><td>10.058734 (52)</td><td>0.96 (11)</td><td>1.24 (12)</td></tr><tr><td>11.756385 (69)</td><td>0.73 (11)</td><td>4.06 (15)</td></tr><tr><td>9.139954 (75)</td><td>0.67 (11)</td><td>1.88 (17)</td></tr><tr><td>0.960644 (81)</td><td>0.62 (11)</td><td>3.39 (18)</td></tr><tr><td>0.611525 (82)</td><td>0.61 (11)</td><td>3.31 (18)</td></tr><tr><td>0.898446 (85)</td><td>0.59 (11)</td><td>4.89 (19)</td></tr><tr><td>11.189760 (93)</td><td>0.54 (11)</td><td>5.06 (21)</td></tr><tr><td>2.848108 (99)</td><td>0.50 (11)</td><td>2.09 (22)</td></tr><tr><td>11.26188 (10)</td><td>0.49 (11)</td><td>4.08 (23)</td></tr><tr><td>8.23509 (11)</td><td>0.47 (11)</td><td>5.52 (24)</td></tr><tr><td>8.02296 (14)</td><td>0.35 (11)</td><td>2.15 (32)</td></tr><tr><td>11.02859 (16)</td><td>0.32 (11)</td><td>5.41 (35)</td></tr></tbody></table>

**Notes.** The reference for the phase is BJD 2454600. Numbers in parentheses indicate the 3σ uncertainties on the last two digits.

<table><caption>**Table B.2.** Amplitude of the ten main δ Scuti frequencies found in the Fourier spectrum of σ Oct for each winter and for the four winters together.</caption><thead><tr><th>Frequency (d<sup>-1</sup>)</th><th>2008</th><th>2009</th><th>2010</th><th>2011</th><th>All</th></tr></thead><tbody><tr><td>10.492693 (13)</td><td>9.69 (32)</td><td>5.42 (26)</td><td>3.01 (22)</td><td>0.99 (16)</td><td>3.73 (11)</td></tr><tr><td>9.719016 (16)</td><td>2.58 (32)</td><td>2.78 (26)</td><td>3.07 (22)</td><td>3.32 (16)</td><td>3.11 (11)</td></tr><tr><td>10.741796 (18)</td><td>2.84 (32)</td><td>2.74 (26)</td><td>2.67 (22)</td><td>2.79 (16)</td><td>2.77 (11)</td></tr><tr><td>8.799398 (19)</td><td>2.91 (32)</td><td>2.56 (26)</td><td>2.56 (22)</td><td>2.52 (16)</td><td>2.61 (11)</td></tr><tr><td>10.251893 (20)</td><td>2.29 (32)</td><td>2.14 (26)</td><td>2.27 (22)</td><td>2.58 (16)</td><td>2.47 (11)</td></tr><tr><td>11.429052 (22)</td><td>2.39 (32)</td><td>2.24 (26)</td><td>2.12 (22)</td><td>2.53 (16)</td><td>2.27 (11)</td></tr><tr><td>9.361590 (28)</td><td>2.99 (32)</td><td>1.67 (26)</td><td>1.51 (22)</td><td>1.55 (16)</td><td>1.80 (11)</td></tr><tr><td>9.771399 (30)</td><td>1.56 (32)</td><td>1.55 (26)</td><td>1.82 (22)</td><td>1.75 (16)</td><td>1.65 (11)</td></tr><tr><td>10.440784 (37)</td><td>1.21 (32)</td><td>1.24 (26)</td><td>1.36 (22)</td><td>1.72 (16)</td><td>1.37 (11)</td></tr><tr><td>10.058734 (52)</td><td>1.38 (32)</td><td>1.31 (26)</td><td>0.93 (22)</td><td>0.77 (16)</td><td>0.96 (11)</td></tr></tbody></table>

**Notes.** Numbers in parentheses indicate the 3σ uncertainties on the last two digits. The amplitudes are in mmag.

<table><caption>**Table B.3.** Amplitude of the ten δ Scuti frequencies with signal to noise ratios greater than four found in the Fourier spectrum of HD 184465 for each winter and for the four winters together.</caption><thead><tr><th>Frequency (d<sup>-1</sup>)</th><th>2008</th><th>2009</th><th>2010</th><th>2011</th><th>All</th></tr></thead><tbody><tr><td>9.671458 (12)</td><td>0.86 (32)</td><td>2.24 (24)</td><td>3.60 (19)</td><td>5.98 (14)</td><td>3.94 (10)</td></tr><tr><td>6.465094 (56)</td><td>1.33 (32)</td><td>0.90 (24)</td><td>0.93 (19)</td><td>0.67 (14)</td><td>0.85 (10)</td></tr><tr><td>8.691272 (64)</td><td>0.26 (32)</td><td>0.34 (24)</td><td>0.28 (19)</td><td>1.41 (14)</td><td>0.74 (10)</td></tr><tr><td>4.934047 (69)</td><td>0.79 (32)</td><td>0.57 (24)</td><td>0.72 (19)</td><td>0.69 (14)</td><td>0.68 (10)</td></tr><tr><td>3.785687 (87)</td><td>0.36 (32)</td><td>0.38 (24)</td><td>0.63 (19)</td><td>0.61 (14)</td><td>0.55 (10)</td></tr><tr><td>14.20371 (12)</td><td>0.07 (32)</td><td>0.06 (24)</td><td>0.19 (19)</td><td>0.74 (14)</td><td>0.40 (10)</td></tr><tr><td>4.35554 (13)</td><td>0.24 (32)</td><td>0.40 (24)</td><td>0.36 (19)</td><td>0.53 (14)</td><td>0.37 (10)</td></tr><tr><td>7.70589 (18)</td><td>0.26 (32)</td><td>0.38 (24)</td><td>0.38 (19)</td><td>0.18 (14)</td><td>0.26 (10)</td></tr><tr><td>10.03363 (25)</td><td>0.39 (32)</td><td>0.29 (24)</td><td>0.19 (19)</td><td>0.32 (14)</td><td>0.19 (10)</td></tr><tr><td>11.82126 (26)</td><td>0.19 (32)</td><td>0.27 (24)</td><td>0.10 (19)</td><td>0.22 (14)</td><td>0.19 (10)</td></tr></tbody></table>

**Notes.** Numbers in parentheses indicate the 3σ uncertainties on the last two digits. The amplitudes are in mmag.

<table><caption>**Table B.4.** Amplitude variations of three δ Scuti frequencies of HD 184465 during the four winters, with two time intervals for each of the 2008, 2009, and 2010 winters and six 30-day intervals for the 2011 winter.</caption><thead><tr><th rowspan="2">JD - 2450000</th><th colspan="3">Frequency (d<sup>-1</sup>)</th><th rowspan="2"></th></tr><tr><th>9.671</th><th>8.692</th><th>6.465</th></tr></thead><tbody><tr><td>4652</td><td>0.84</td><td>0.32</td><td>1.30</td><td>(36)</td></tr><tr><td>4702</td><td>1.04</td><td>0.41</td><td>1.40</td><td>(52)</td></tr><tr><td>5003</td><td>2.23</td><td>0.39</td><td>0.79</td><td>(35)</td></tr><tr><td>5074</td><td>2.28</td><td>0.38</td><td>0.85</td><td>(31)</td></tr><tr><td>5329</td><td>3.27</td><td>0.35</td><td>0.81</td><td>(32)</td></tr><tr><td>5419</td><td>3.91</td><td>0.53</td><td>0.99</td><td>(25)</td></tr><tr><td>5652</td><td>5.59</td><td>0.76</td><td>0.75</td><td>(39)</td></tr><tr><td>5682</td><td>5.76</td><td>0.65</td><td>0.83</td><td>(40)</td></tr><tr><td>5712</td><td>5.59</td><td>1.37</td><td>0.77</td><td>(26)</td></tr><tr><td>5742</td><td>6.12</td><td>1.24</td><td>0.49</td><td>(29)</td></tr><tr><td>5772</td><td>6.02</td><td>1.88</td><td>0.89</td><td>(33)</td></tr><tr><td>5802</td><td>6.10</td><td>2.74</td><td>0.61</td><td>(51)</td></tr></tbody></table>

**Notes.** Dates are the central time of each interval. Numbers in parentheses indicate the 3σ uncertainties on the last two digits for the amplitudes (they are the same for all amplitudes at a given time interval). The amplitudes are in mmag.
Appendix C: New variables stars and eclipsing binaries identified with ASTEP South

<table><caption>Table C.1. New variable stars identified with ASTEP South.</caption><thead><tr><th>ID<br>(UCAC3)</th><th>RA<br>(J2000)</th><th>Dec<br>(J2000)</th><th>fmag<br>(mag)</th><th>Frequency<br>(d<sup>-1</sup>)</th><th>Period<br>(d)</th><th>Amplitude<br>(mmag)</th><th>Type</th></tr></thead><tbody><tr><td>004-005139</td><td>13:31:52.604</td><td>-88:22:17.75</td><td>10.23</td><td>0.14252 (3)</td><td>7.02</td><td>2.3 (0.1)</td><td>Var</td></tr><tr><td>003-002049</td><td>07:54:08.982</td><td>-88:31:04.01</td><td>10.91</td><td>5.83351 (6)</td><td>0.17</td><td>1.1 (0.1)</td><td>δ Scuti</td></tr><tr><td>003-004378</td><td>15:31:11.257</td><td>-88:31:17.40</td><td>10.93</td><td>0.10227 (3)</td><td>9.78</td><td>3.7 (0.2)</td><td>VRot?</td></tr><tr><td>004-006257</td><td>16:13:29.743</td><td>-88:12:25.00</td><td>11.82</td><td>0.08001 (2)</td><td>12.50</td><td>4.5 (0.2)</td><td>Var</td></tr><tr><td>002-000843</td><td>05:12:39.072</td><td>-89:10:02.73</td><td>12.18</td><td>0.13586 (3)</td><td>7.36</td><td>3.0 (0.5)</td><td>VRot</td></tr><tr><td>002-000214</td><td>01:15:41.760</td><td>-89:00:19.72</td><td>12.21</td><td>0.02587 (2)</td><td>38.65</td><td>6.1 (0.2)</td><td>VRot</td></tr><tr><td>002-000853</td><td>05:17:21.064</td><td>-89:21:19.78</td><td>12.25</td><td>1.05540 (6)</td><td>0.95</td><td>1.6 (0.2)</td><td>VRot?</td></tr><tr><td>004-000689</td><td>01:57:52.063</td><td>-88:13:15.56</td><td>12.55</td><td>0.12161 (4)</td><td>8.22</td><td>3.8 (0.3)</td><td>VRot?</td></tr><tr><td>003-002304</td><td>08:39:17.699</td><td>-88:44:20.12</td><td>12.59</td><td>0.08060 (2)</td><td>12.41</td><td>5.3 (0.3)</td><td>Var</td></tr><tr><td>003-004226</td><td>15:05:46.232</td><td>-88:39:23.33</td><td>12.67</td><td>0.02679 (3)</td><td>37.33</td><td>3.9 (0.3)</td><td>VRot</td></tr><tr><td>003-005314</td><td>19:03:58.030</td><td>-88:59:04.65</td><td>12.69</td><td>0.02492 (3)</td><td>40.13</td><td>3.4 (0.2)</td><td>VRot?</td></tr><tr><td>001-000287</td><td>05:43:45.773</td><td>-89:33:57.57</td><td>12.69</td><td>0.10116 (3)</td><td>9.89</td><td>4.5 (0.2)</td><td>Var</td></tr><tr><td>001-000176</td><td>03:35:19.894</td><td>-89:41:08.84</td><td>12.80</td><td>0.04475 (3)</td><td>22.35</td><td>4.0 (0.2)</td><td>Var</td></tr><tr><td>003-001672</td><td>06:30:27.586</td><td>-88:59:11.08</td><td>12.83</td><td>0.04459 (2)</td><td>22.43</td><td>5.6 (0.2)</td><td>Var</td></tr><tr><td>003-003813</td><td>13:42:59.863</td><td>-88:53:36.27</td><td>12.96</td><td>0.07838 (2)</td><td>12.76</td><td>6.7 (0.3)</td><td>VRot</td></tr><tr><td>003-000054</td><td>00:16:06.515</td><td>-88:34:08.37</td><td>12.98</td><td>0.01286 (1)</td><td>77.76</td><td>2.1 (0.5)</td><td>VRot</td></tr><tr><td>002-000558</td><td>03:32:17.485</td><td>-89:15:06.15</td><td>13.51</td><td>0.09741 (4)</td><td>10.27</td><td>4.6 (0.4)</td><td>VRot?</td></tr><tr><td>003-001584</td><td>06:12:58.011</td><td>-88:35:13.08</td><td>13.52</td><td>0.12314 (4)</td><td>8.12</td><td>6.5 (0.5)</td><td>VRot?</td></tr><tr><td>003-000486</td><td>01:59:49.979</td><td>-88:54:05.91</td><td>13.64</td><td>0.02878 (2)</td><td>34.75</td><td>13.3 (0.6)</td><td>VRot</td></tr><tr><td>002-002489</td><td>14:58:37.642</td><td>-89:07:37.09</td><td>13.81</td><td>0.04184 (2)</td><td>23.90</td><td>18.5 (0.6)</td><td>VRot</td></tr><tr><td>003-004679</td><td>16:35:40.811</td><td>-88:43:46.12</td><td>13.83</td><td>0.02809 (5)</td><td>35.60</td><td>88.9 (1.0)</td><td>VRot</td></tr><tr><td>003-004666</td><td>16:32:18.171</td><td>-88:36:33.72</td><td>14.21</td><td>0.02950 (3)</td><td>33.90</td><td>18.4 (1.1)</td><td>VRot</td></tr><tr><td>004-004408</td><td>11:55:09.767</td><td>-88:23:58.50</td><td>14.32</td><td>0.23194 (4)</td><td>4.31</td><td>7.7 (0.7)</td><td>VRot?</td></tr><tr><td>004-004603</td><td>12:17:46.280</td><td>-88:05:49.99</td><td>14.44</td><td>1.03020 (5)</td><td>0.97</td><td>9.6 (1.1)</td><td>Var</td></tr><tr><td>004-007265</td><td>19:04:41.581</td><td>-88:29:17.88</td><td>14.46</td><td>2.11961 (4)</td><td>0.47</td><td>5.4 (0.5)</td><td>Var</td></tr><tr><td>003-003146</td><td>11:28:48.149</td><td>-88:50:41.89</td><td>14.48</td><td>0.42148 (4)</td><td>2.37</td><td>11.5 (1.0)</td><td>Var</td></tr><tr><td>002-001119</td><td>06:54:02.979</td><td>-89:05:10.41</td><td>14.54</td><td>5.89217 (4)</td><td>0.17</td><td>10.8 (1.0)</td><td>Var</td></tr><tr><td>003-006325</td><td>23:21:41.025</td><td>-88:50:24.29</td><td>14.91</td><td>0.64068 (2)</td><td>1.56</td><td>35.6 (1.4)</td><td>Var</td></tr><tr><td>001-000349</td><td>06:48:03.486</td><td>-89:32:05.52</td><td>14.99</td><td>5.64772 (6)</td><td>0.18</td><td>97.0 (1.3)</td><td>δ Scuti</td></tr><tr><td>001-000856</td><td>15:36:30.014</td><td>-89:47:17.82</td><td>15.05</td><td>0.80258 (3)</td><td>1.25</td><td>18.9 (1.2)</td><td>VRot</td></tr><tr><td>002-001114</td><td>06:52:18.556</td><td>-89:01:33.13</td><td>15.09</td><td>0.12765 (2)</td><td>7.83</td><td>14.9 (0.6)</td><td>Var</td></tr><tr><td>002-003252</td><td>19:53:39.345</td><td>-89:01:49.63</td><td>15.20</td><td>1.82823 (2)</td><td>0.55</td><td>43.0 (1.4)</td><td>VRot</td></tr><tr><td>002-001915</td><td>11:32:03.119</td><td>-89:04:57.12</td><td>15.21</td><td>0.95488 (2)</td><td>1.05</td><td>28.6 (1.4)</td><td>Var</td></tr><tr><td>002-001216</td><td>07:28:33.245</td><td>-89:03:48.27</td><td>15.69</td><td>0.77189 (3)</td><td>1.30</td><td>14.1 (0.9)</td><td>Var</td></tr></tbody></table>

Notes. The ID, coordinates, and magnitudes fmag are extracted from the UCAC3 catalog. For each object, we report the frequency, the corresponding period, the semi-amplitude of the variations, and the type of variable. Numbers in parentheses indicate the $3\sigma$ uncertainties on the last digit for the frequencies, and the $3\sigma$ uncertainties for the amplitudes. “Var” indicates that the type of variable cannot be determined and “?” indicates a tentative classification.
<table><caption>**Table C.2.** New eclipsing binaries identified with ASTEP South and known objects with improved classification and parameters.</caption><thead><tr><th>ID<br>(UCAC3)</th><th>RA<br>(J2000)</th><th>Dec<br>(J2000)</th><th>fmag<br>(mag)</th><th>Period<br>(d)</th><th>T<sub>0</sub><br>(BJD)</th><th>Depth<br>(mmag)</th><th>Type</th></tr></thead><tbody><tr><td>003-001765<sup>a</sup></td><td>06:53:22.343</td><td>-88:34:04.28</td><td>10.04</td><td>0.4574110 (48)</td><td>4627.4151 (34)</td><td>2.6 (0.4)</td><td>EB</td></tr><tr><td>002-002671</td><td>16:04:20.122</td><td>-89:00:40.74</td><td>10.99</td><td>13.4892 (25)</td><td>4627.70 (13)</td><td>5.7 (0.4)</td><td>EB</td></tr><tr><td>002-003763</td><td>23:23:51.240</td><td>-89:25:17.62</td><td>11.25</td><td>23.3853 (53)</td><td>4636.35 (18)</td><td>8.6 (0.4)</td><td>EB</td></tr><tr><td>004-006956</td><td>18:03:31.003</td><td>-88:16:36.79</td><td>11.86</td><td>36.263 (24)</td><td>4634.18 (23)</td><td>5.9 (0.4)</td><td>EB</td></tr><tr><td>004-000272<sup>b</sup></td><td>00:45:26.547</td><td>-88:25:21.03</td><td>11.98</td><td>34.559 (19)</td><td>4627.69 (19)</td><td>15.0 (0.7)</td><td>EB</td></tr><tr><td>001-000823<sup>c</sup></td><td>14:50:14.729</td><td>-89:46:58.96</td><td>12.18</td><td>0.28891080 (20)</td><td>4627.2346 (02)</td><td>479 (8)</td><td>EW</td></tr><tr><td>001-000701</td><td>12:41:48.571</td><td>-89:42:30.30</td><td>12.69</td><td>54.133 (15)</td><td>4668.37 (20)</td><td>38 (2)</td><td>EB</td></tr><tr><td>002-003056<sup>d</sup></td><td>18:36:30.712</td><td>-89:11:02.84</td><td>12.76</td><td>75.57626 (20)</td><td>4629.9895 (60)</td><td>30 (2)</td><td>EA</td></tr><tr><td>001-001090</td><td>20:37:01.418</td><td>-89:44:34.88</td><td>12.82</td><td>23.7052 (70)</td><td>4632.38 (10)</td><td>10 (1)</td><td>EW</td></tr><tr><td>001-000120</td><td>02:26:38.190</td><td>-89:31:31.72</td><td>14.45</td><td>81.051070 (20)</td><td>5058.9223 (24)</td><td>335 (5)</td><td>EA</td></tr><tr><td>001-000682</td><td>12:27:06.862</td><td>-89:48:58.91</td><td>14.59</td><td>1.048938 (30)</td><td>4627.2738 (97)</td><td>28.0 (0.8)</td><td>EB</td></tr></tbody></table>

**Notes.** The ID, coordinates, and magnitudes fmag are extracted from the UCAC3 catalog. For each object, we report the period, the mid-eclipse date of the first detected eclipse, the depth of the primary eclipse (or of the secondary eclipse in the potential case of eccentric binaries with no primary eclipses), and the type of eclipsing binary. Numbers in parentheses indicate the $3\sigma$ uncertainties on the last two digits for the periods and ephemerides, and the $3\sigma$ uncertainties for the depths. <sup>(a)</sup>This object was reported as a rotating ellipsoidal variable by Wang et al. (2015). <sup>(b)</sup>This object was reported as a variable of unknown type by Wang et al. (2011) and was not identified by Wang et al. (2013). <sup>(c)</sup>This object is EN Oct: we improved its period and eclipse depth compared to previous studies by Knigge & Bauernfeind (1967), Otero (2004). <sup>(d)</sup>This object is discussed in Sect. 7.2.5.