Free Access
Issue
A&A
Volume 564, April 2014
Article Number A58
Number of page(s) 18
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/201323041
Published online 04 April 2014

© ESO, 2014

1. Introduction

Observing planetary transits gives access to a set of properties of the star+planet system (e.g. Winn 2009). The relative size of the planet with respect to the star, given by the depth of the transit light-curve, is one of the most straightforward system observables. Spectrally resolved observations have shown for gaseous giant exoplanets that transit the brightest stars that the transit depth varies with wavelength (the first to report this were Charbonneau et al. 2002; Vidal-Madjar et al. 2003, 2004) which reveals the signature of the planetary atmosphere. At certain wavelengths, the atmosphere indeed absorbs the stellar light and becomes opaque, which makes the planet look larger in front of the star and produces a deeper transit. The atmosphere of smaller terrestrial exoplanets close to or in the habitable zone (HZ) is in principle also accessible by this method (early works by Schneider 1992, 1994, and more recent models by Ehrenreich et al. 2006, 2012; Kaltenegger & Traub 2009; García-Muñoz et al. 2012; Snellen et al. 2013; Bétrémieux & Kaltenegger 2013). Recent simulations show that the next generation of instruments in the infrared or visible domains, that is the James Webb Space Telescope (JWST) and the European Extremely Large Telescope (E-ELT), should be able to offer a first access to characterizing terrestrial planets in the HZ and detecting biosignatures (here atmospheric O2 and O3) and molecules linked to habitability (Rauer et al. 2011; Hedelt et al. 2013; Snellen et al. 2013; Rodler & López-Morales 2014).

In this context, it is of prime interest to observe Earth as a transiting exoplanet. Vidal-Madjar et al. (2010) observed the August 2008 partial lunar eclipse with the high-resolution SOPHIE échelle spectrograph at the 193 cm Haute-Provence telescope (Perruchot et al. 2008; Bouchy et al. 2009) and showed that the Earth atmosphere thickness versus wavelength can be derived from such observations. Earth observed from the Moon during a lunar eclipse indeed transits in front of the Sun and gives the possibility to study the Earth atmosphere as during a transit. These first observations also allowed one to identify molecular species in the atmosphere, that is water vapour (H2O), a molecule linked to Earth habitability, biogenic oxygen (O2), and ozone (O3). Ozone and oxygen are considered as reliable biosignatures for Earth-like planets (Owen 1980; Léger et al. 1993; Seager et al. 2013) except for terrestrial planets outside the liquid-water habitable zone (Schindler & Kasting 2000). These species were also detected by Pallé et al. (2009) from umbra spectra taken during the same eclipse. However, Vidal-Madjar et al. (2010) pointed out that the Pallé et al. spectrum is not the transmission spectrum of Earth observed as a distant transiting planet. The umbra signal indeed results from photons that are highly scattered and refracted by the deep atmosphere. These photons strongly deviate from the star-planet direction and cannot be collected during the observation of a transit. From penumbra observations, Vidal-Madjar et al. (2010) additionally detected the thin high-altitude sodium Na I layer. The larger diameter of Earth in the blue due to increasing Rayleigh scattering at shorter wavelengths is also well visible, although it is partially blended with residuals from spectrograph orders correction.

In this paper, we analyze new observations obtained during the lunar eclipse on Dec. 21, 2010 with the ESO spectrographs HARPS (from 3800 to 6900 Å, Mayor et al. 2003) and UVES (from 3200 to 10 400 Å, Dekker et al. 2000).

We present first the principle of the observation and different reduction methods: the original Vidal-Madjar et al. (2010) data reduction method (method 1), which makes use of a pair of penumbra spectra, and a second method (method 2), which requires only one single penumbra spectrum to recover the atmosphere thickness profile (Sect. 2). We also present a third alternative method to derive the thickness profile, which makes use of penumbra and umbra spectra. This third approach is biased by the refraction in the Earth atmosphere (Sect. 3, method 3), but can be corrected a posteriori with the result of method 1. An observation log and more details about the spectrographs are given in Sect. 4. We give details on the data reduction in Sects. 5 and 6 for HARPS and UVES, respectively. The altitude profiles versus wavelength obtained with the three different methods are detailed in Sect. 7, as well as simulations of the E-ELT observing a transiting Earth twin at 10 pc.

2. Principle of the observation (method 1 and 2)

The principle of Vidal-Madjar et al. (2010) method is the following: for an observer on the Moon in the deep penumbra, the Sun appears above the Earth limb as a small crescent. The light in the penumbra is a mix of direct solar light that did not interact with Earth with light that passed through the (mostly upper) Earth atmosphere. The nature of this mixed light in the penumbra is similar to what is observed during a transit of an exoplanet, here Earth in front of the Sun. Nevertheless, as seen from the Moon, Earth is larger than the Sun, and the Earth radius versus wavelength cannot be reconstructed from the transit depth et all wavelengths, as done with transiting exoplanets. The method consists of observing two different positions of the Sun above the Earth limb (during ingress or egress) – that is, of recording Moon spectra at two different depths in the penumbra (method 1). For these two different positions on the Moon, the solar crescent above the Earth limb as well as the arc of the Earth atmosphere intersecting the solar disk, are different. We have shown (Vidal-Madjar et al. 2010, Eq. (14) therein) that h, the thickness profile of the atmosphere, or the altitude under which the atmosphere can be considered as opaque, can be written h(λ)=(EA(λ)EA(λ0)×FA(λ0)FA(λ)EB(λ)EB(λ0)×FB(λ0)FB(λ))\begin{eqnarray} h (\lambda)& =& \left( \frac{E_A (\lambda)}{E_A (\lambda_0)} \times \frac{F_A (\lambda_0)}{F_A(\lambda)} - \frac{E_B (\lambda)}{E_B(\lambda_0)} \times \frac{F_B (\lambda_0)}{F_B (\lambda)} \right) \nonumber\\ && \div \left( \frac{L_B}{S_B} - \frac{L_A}{S_A} \right), \label{ratioAB} \end{eqnarray}(1)where E and F are the fluxes during the eclipse and outside the penumbra (almost full-Moon just before or after the eclipse), respectively. Subscripts A and B stand for the two different positions in the penumbra. The value λ0 is an arbitrary wavelength for which we define h(λ0) = 0, and the symbols S and L represent geometrical parameters equal to the surface of the Sun crescent above the Earth limb and the length of the arc of the Earth atmosphere intersecting the solar disk for h(λ0), respectively (see Fig. 4 in Vidal-Madjar et al. 2010). Equation (1) shows that h(λ) can thus be written as a function of differences between the measured fluxes divided by a factor of geometrical values. The numerical values for S are recovered from the measured fluxes too, E(λ0) and F(λ0), following Eq. (11) of Vidal-Madjar et al. (2010), E(λ0)=F(λ0)×SS,\begin{eqnarray} E (\lambda_0) = F (\lambda_0) \times \frac{S}{S_\odot}, \label{S} \end{eqnarray}(2)where S is the surface of the solar disk. At last, once the surfaces S are known, the lengths L can be calculated quite simply from geometrical considerations knowing the apparent diameters of Earth and Sun, so h(λ) can finally be calculated.

Equation (1) above is obtained from the differences between Eqs. (12) and (13) in Vidal-Madjar et al. (2010), each being written for one position on the Moon. These equations can be rewritten, for position A for example, h(λ)=(1EA(λ)EA(λ0)×FA(λ0)FA(λ))×SALA·\begin{eqnarray} h (\lambda)= \left( 1 - \frac{E_A (\lambda)}{E_A (\lambda_0)} \times \frac{F_A (\lambda_0)}{F_A (\lambda)} \right) \times \frac{S_A}{L_A} \cdot \label{ratioA} \end{eqnarray}(3)We note that h(λ) is fully calculable from Eq. (3) above, thus with only one penumbra spectrum (method 2). Nevertheless, because E and F are recorded at different times in the night, the barycentric Earth radial velocity (BERV) leads to a different Doppler shift between the solar and telluric lines in E and F spectra. Since we are interested in the signature of the Earth atmosphere, all spectra are re-aligned on the telluric lines, which slightly shifts the solar lines between E and F. Moreover, the solar lines integrated over the solar crescent seen from the Moon during the partial eclipse are also Doppler-shifted with respect to the lines of the full (uneclipsed) solar disk. The solar lines thus do not cancel out correctly in the ratio EA(λ)EA(λ0)\hbox{$\frac{E_A (\lambda)}{E_A (\lambda_0)}$}/FA(λ)FA(λ0)\hbox{$\frac{F_A(\lambda)}{F_A (\lambda_0)}$}, and the obtained h(λ) profile is noisier, although the overall shape is correct. With the difference in the left part of Eq. (1), the solar lines cancel out much better and the h(λ) profile is cleaner. This is discussed in more details in Sect. 7.2.

3. An alternative method that approximates h(λ) (method 3)

This section describes another method, different from that of Vidal-Madjar et al. (2010), to derive the thickness profile of Earth atmosphere. In a lunar eclipse, the radius of Earth is given by the radius of the umbra. The umbra radius can be calculated geometrically for a solid sphere that casts its shadow to the Moon, but in practice, the Earth atmosphere that scatters light along Earth limb, lead to a blurred edge between umbra and penumbra on the Moon surface. The radiance profile across the umbra edge has been measured in the past by several observers (especially by Dubois from 1950 to 1980, Danjon 1925; Combaz et al. 1950; Reaves & Walker 1952; Walker & Reaves 1957). The irradiance of the eclipsed Moon as well as the umbra edge have also been modelized (Link 1969; Vollmer & Gedzelman 2008; García-Muñoz & Pallé 2011) by an analysis of the light propagation from the Sun to the Moon through the Earth atmosphere. The observed magnitude versus position with respect to the umbra center looks approximately like a hyperbolic function (or sigmoid). The observations mentioned above indicate that magnitude typically drops by about 8 mag from penumbra to umbra, the steepest region being typically 2 arcmin wide (Dubois & Poumeyrol 1954). Unfortunately, the spatial resolution of published radial photometry across the umbra edge is only about 1 arcmin, and not sufficient to reveal a dependance of the umbra radius with wavelength. To fix the order of magnitude, with an Earth umbra radius of 43 arcmin, a variation of 10 km in the Earth radius (of 6371 km) translates into a change of the umbra radius of 4 arcsec.

With ESO’s HARPS and UVES spectrographs, we recorded spectra from deep penumbra to umbra over 7 arcmin with HARPS and 19 arcmin with UVES. Because the telescope tracked neither the star nor the Moon but the Earth umbra projected to the Moon, we were able to integrate the signal at fixed distances from the umbra edge, while the Moon itself slowly shifted in front of the spectrographs during the exposures. To compare the flux in these spectra taken at different positions in the umbra or in the penumbra, we calculate their corresponding magnitude −2.5  log 10(E/F) relative to the full-Moon spectra recorded before the eclipse. Then we plot for each wavelength the magnitudes versus their distance from the umbra center, as Dubois did. But while Dubois observed at three or four different wavelengths through broadband blue, green, yellow, and red Wratten filters, HARPS and UVES allow to do this now at a much higher spectral resolution. The points follow the same overall shape as Dubois’ profiles (Fig. 1). We fit a hyperbolic function through these points ζ(r)=a0+a1tanh[a3(ra2)],\begin{eqnarray} \zeta(r) = a_0+a_1\tanh[ a_3(r-a_2)], \label{eq:model} \end{eqnarray}(4)where ζ is the model giving the magnitude in the radial direction through the umbra edge, r the distance to the umbra center, and ai the parameters defining the hyperbolic function. The function ζ(r) varies between a0 − a1 and a0 + a1. We may arbitrarily define the abscissa a2 of the inflection point (at ordinate a0) of this function as a measurement of the umbra radius, and thus of the Earth radius at the observing wavelength.

Nevertheless, Fig. 1 shows that the inflection point corresponds to a zone about 3 mag (15×) fainter than shallower penumbra. But in a real transit, only the beams of light that interact with the top of the atmosphere, and thus only weakly absorbed, leave a signature in the transit. If a beam goes deeper to the atmosphere, it is significantly dimmed and scattered, and finally absent from the transit signature of the planet. Therefore the informations from the 3 mag fainter light we obtain at the inflection point are not those we obtain during a real transit. The inflection point is thus not the best definition of the edge of the Earth’s umbra on the Moon, at least in the context of our Earth-transit analysis. It is more relevant to define a point slightly more in the brighter penumbra for the edge of Earth’s umbra, as shown by the vertical dotted line in Fig. 1, and to build up the Earth radius versus wavelength from this point instead of from the inflection point. More details are given in Sect. 7.3.

thumbnail Fig. 1

Measured magnitude across the umbra edge with HARPS at 5500 Å. Crosses and diamonds stand for A and B fiber, respectively. The horizontal error bars are 2σ = 20 arcsec long and represent the typical pointing error of the telescope. Vertically, the 2σ error bars are the magnitude errors evaluated from the S/N of the spectra at that wavelength and the evaluated error in the local Moon albedo. The solid line is the model from Eq. (4) fitted through the data. The vertical dotted line marks our definition of the umbra edge. See Sect. 7.3 for details.

We understand that this method makes use of information that is not available in an exoplanet transit: indeed it uses information from spectra of the umbra to build the function ζ, but, as mentioned in the introduction, the photons from the umbra result from highly angularly scattered or refracted beams, which are not available in the observation of an exoplanet transit. Moreover, the umbra photons, especially the red ones, have been refracted by the deepest Earth atmospheric layers, and the Earth shadow is thus shrunk by refraction at these wavelengths, leading to a smaller Earth radius in the red. In other words, the gradient in the altitude profiles measured from blue to red wavelengths will be steeper due to refraction. But we show in Sect. 7.3 that the altitude profiles h(λ) obtained by this method can be calibrated (scaled) to the profiles obtained with method 1.

4. Observing log

HARPS provides spectra with a spectral resolution of 115 000 and covers the visible range (Table 1, Mayor et al. 2003). UVES was used in dichroic-1 and dichroic-2 modes (DIC1 and DIC2). In each of these modes, both BLUE and RED arms of the spectrograph are used simultaneously, and with standard settings, two exposures (one in each mode) give a full [3100–10 400] Å spectrum (Table 1, Dekker et al. 2000). The slit width was set to 0.3 arcsec, which gives a spectral resolution of 120 000.

Table 1

HARPS and UVES spectral ranges.

Tables 2 and 3 give information on the observations made from La Silla with HARPS. The position with respect to the umbra edge was measured with Guide-91. Tables 4 and 5 give the same information but for UVES at Paranal. With UVES, the telescope was not positioned as expected on the Moon during the eclipse and the expected positions are indicated in parentheses in Table 5.

Table 2

Series of n Moon reference spectra taken with HARPS before the eclipse.

Table 3

HARPS exposures and fiber positions in the penumbra (negative values) and umbra (positive values) with respect to the umbra edge given by Guide-9.

Table 4

Series of three Moon reference spectra taken with UVES before the eclipse.

Table 5

UVES exposures or series of n exposures, and measured slit positions in the penumbra (negative values) and umbra (positive values) with respect to the umbra edge given by Guide-9.

As an example, we show in Fig. 2 the Sun seen from where HARPS fiber-A was positioned on the Moon at 07:06:35 UT, at the beginning of the 300 s exposure (Table 3). The Sun is setting in the Pacific Ocean about 2000 km southeast of Japan, at coordinates ≈ 156° E and 23° N. The daytime cloud fraction for 2010 Dec. 21 given by Aqua/MODIS2 for these coordinates is 0.7.

thumbnail Fig. 2

Earth and Sun seen from the Moon deep penumbra where HARPS fiber A was positioned at 07:06:35 UT, at the beginning of the 300 s exposure. Image with Guide-9.

5. HARPS data processing

The data processing makes use of so-called s1d reduced data output files directly from the HARPS data reduction software (DRS). The HARPS instrumental response is calibrated with a high signal-to-noise ratio (S/N) G2V star spectrum.

5.1. Correction for the Moon reflectivity

To accurately maintain a given position of the spectrograph with respect to the umbra edge during the exposures, we needed to track the umbra, not the Moon itself. The Moon was thus moving in front of the spectrograph pair of fibers at a speed of about 0.55 arcsec/s. The shift of the Moon during a given exposure contributed to smooth out any relative spectral reflectance differences between the lunar terrains that drifted in front of the spectrograph. But different terrains with different albedos were seen by each fiber during the different exposures. To correct each spectrum for the flux variations induced by these albedo variations, we measure the mean lunar reflectivity along the path of the slit across the Moon surface during the exposures a posteriori. The path is reconstructed from the date and telescope coordinates available in the metadata. We build a picture of the Moon as seen during the eclipse from La Silla at a resolution of 1 arcsec/pixel: the reconstructed scene is an orthographic projection of a Clementine albedo map3 at 750 nm (Eliason et al. 1999) and takes into account the correct libration angles, the correct phase angle (close of zero degree), and a Lambert-Lunar limb-darkening law (Mc Ewen 1996). The mean lunar reflectivity is then measured with Audela-2.04 along the path of each fiber across the lunar surface as shown in Fig. 3.

thumbnail Fig. 3

Example of a path of the HARPS fibers across the reconstructed Moon image, during the 300 s exposure in the penumbra. The darkening from the penumbra and umbra are not superimposed on the image here, since we are only interested in recovering the moon reflectivity. The A and B fibers remained at 232 and 125 arcsec from the umbra edge, while the fibers path across the Moon was 164 arcsec long (length of the yellow lines). The mean reflectivity is measured along the path, taking into account a 10 arcsec pointing error in both RA and Dec directions. The terrain reflectivity fluctuates by σ = 21% along the path for Fiber A and σ = 15% for Fiber B in this highly contrasted lunar region. The error on the mean reflectivity due to pointing uncertainty is estimated to be about σ = 4%. The figure shows that while Fiber A went through the dark region of Mare Nectaris, Fiber B remained on a region measured 39% brighter between Catharina and Cyrillus craters, pointing out the relevance of correcting the flux for lunar local reflectivity.

5.2. Correction for the signature of the atmosphere above the telescope

To correct our eclipse spectra for the signatures of the atmosphere above the telescope, we need a reference transmission spectra for the atmosphere T(λ). To derive T(λ), we proceed as in Vidal-Madjar et al. (2010) by calculating the ratio of two full-Moon spectra, F1 and F2, taken before the eclipse at different airmass AM1 and AM2 < AM1. The ratio can be scaled to AM = 1 following T(λ)=[F1AM1(λ)F2AM2(λ)](1/(AM1AM2))·\begin{eqnarray} T(\lambda) = {\Bigg[ \frac{F_1^{AM_1} (\lambda)} {F_2^{AM_2} (\lambda)} \Bigg]}^{(1/(AM_1-AM_2))}\cdot \label{eq:am} \end{eqnarray}(5)We calculate T(λ) with AM1 = 2.2 and AM2 = 1.7 (Table 2). Since we are interested only in the accurate measurement of the telluric lines in T(λ), the two spectra in Eq. (5) are corrected for the BERV to properly align the telluric lines of both spectra (the solar lines become slightly misaligned). Figure 4 shows the obtained average reference transmission T(λ) that we now use in the following data analysis of HARPS data to re-evaluate all eclipse spectra as if they were collected from outside the atmosphere.

thumbnail Fig. 4

La Silla atmosphere transmission function T(λ) for AM = 1 as evaluated from full-Moon spectra taken with HARPS just before the eclipse (Table 2). The spectrum is normalized at 6866 Å. To verify this experimental T(λ) spectrum, a transmission fit based on the model of Hayes & Latham (1975) for the Rayleigh and aerosol components normalized at 6866 Å is displayed (blue dashed line). We also add an ozone transmission spectrum to the model (light-blue solid line) that correctly fits the observed wide Chappuis band. The best least-squares fit for Rayleigh and ozone over La Silla that night is obtained for an airmass of AM = 0.77 for ozone and Rayleigh, and an optical depth of 0.181 at 5320 Å for the aerosols, between 4000 and 6866 Å. The fit for ozone is made only where O2 and H2O absorption bands are absent. The O2γ and B absorption bands are clearly seen, around 6300 and 6900 Å respectively, the other deep bands are due to H2O. The small gap in HARPS data is visible between 5300 and 5340 Å (Table 1).

6. UVES data processing

The UVES instrumental response is calibrated with the reference star HD 80170 (K5 III-IV, Teff = 4537 K, log   (g) = 2.1, Randich et al. 1999).

6.1. Correction for the Moon reflectivity

To measure the mean lunar reflectivity along the path of the UVES slit across the Moon surface, we build a picture of the Moon as seen during the eclipse from Paranal, as we did for the HARPS observations. An example of the slit path is shown in Fig. 5.

thumbnail Fig. 5

Example of a path of the UVES slit across the reconstructed Moon image during the 600 s exposure in the umbra. The slit, 10 arcsec long and perpendicular to the path, remains at 936 arcsec from the umbra edge, while the slit path over the Moon was 330 arcsec long (length of the yellow line). The mean reflectivity is measured along the path, taking into account a 10 arcsec pointing error in both RA and Dec directions. The terrain reflectivity fluctuates by σ = 13% along the path in this bright lunar region, and the error on the mean reflectivity due to pointing uncertainty is estimated to be about σ = 4%.

6.2. Correction for the signature of the atmosphere above the telescope

We can not proceed exactly as with the HARPS data because full-Moon spectra show that the water-vapour lines were stronger while the Moon was higher in the sky, while it should have been lower for the Moon at a smaller airmass. Figure 6 indeed shows that the amount of water vapour above the telescope increased during the first part of the night, between the two sets of full-Moon observations. We describe below how we consequently corrected our UVES data.

thumbnail Fig. 6

Part of the spectrum showing the ratio of a high airmass over a low airmass UVES full-Moon spectra. The oxygen B and A bands at 6900 and 7600 Å  respectively, are in absorption as expected, but the water vapour bands at 7200 and 8200 Å appear in emission, showing that the amount of water vapour along the line of sight was higher when the airmass was lower, indicating that the amount of water vapour above Paranal increased during the first part of the night, before the eclipse. To calculate this ratio, the spectra are aligned by cross-correlation to match the radial velocity of the telluric lines, so the solar lines do not properly cancel out in the ratio, and appear with a wing in emission and the other in absorption.

Table 6

Measured PWV before the eclipse and at the end the observation, and reduced PWV for airmass =1.

We estimate the precipitable water vapour (PWV) in our spectra by using the spectral-fitting method described in Querel et al. (2011). This method compares the measured spectrum with a reference spectrum while varying the reference water-vapour content. The reference spectrum is generated using an atmospheric model configured with the local environmental parameters at the time of measurement (ambient temperature, pressure, elevation, etc.). A comparison is then made over several small spectral intervals (2–5 nm) where the only significant features are caused by absorption by water vapour. A least-squares fit is performed between the measured data and reference spectra for a range of PWV values. The reference spectrum with the smallest fit residual represents the water-vapour content (or PWV) most similar to that of the measured spectrum. Since this absorption is along the line-of-sight and assuming an isotropic distribution of water vapour above the telescope, the absorption must then be corrected back to zenith by dividing the best-fit PWV by the airmass of the measurement.

Table 6 indicates that three hours before the eclipse started the PWV measured on the full-Moon spectrum was about 6.4 mm for an airmass of 2. The PWV for airmass =1 was therefore 3.2 mm. Later in the night, and just before the penumbra ingress, we find that PWV increased to 8.00 mm at airmass =1.56, indicating the increase of moisture above the telescope. The reduced PWV for airmass =1 is thus 8.00/1.56 = 5.1 mm before the eclipse started. Immediately after the eclipse observations, we observed HD 80170 and measured a reduced PWV of 5.71 mm. We therefore assume that during the eclipse, the PWV at zenith remains within the [5.1; 5.7] mm range and we chose the mean value 5.4 mm as the reference PWV at airmass =1 during the eclipse. PWV values of 5 mm are relatively high for Paranal, and we understand that assuming such a high water vapour content is isotropic and constant above the telescope for three hours (the duration of our observations) may be a strong assumption. Nevertheless, it is the best we can do and the whole process of data reduction leading to h(λ) assumes the stability of the PWV above Paranal during the observations (this is obviously also true for the HARPS observations at La Silla).

To correct our eclipse spectra for the signatures of the atmosphere above the telescope, we need distinct reference transmission spectra, one for water and and one for oxygen. In the following, we describe how these reference spectra are obtained. These reference transmission spectra do not need to include Rayleigh and ozone signatures because they are already corrected for by UVES reduction pipeline.

The UVES transmission water-vapour spectrum TH2O(λ) for airmass =1 (PWV = 5.4 mm) is obtained by the ratio of the full-Moon spectra F for which we know the value of PWV (Table 6), TH2O(λ)=(F(PWV=8.00)F(PWV=6.43))5.4/(8.006.43)·\begin{eqnarray} T_{\rm H_2O}(\lambda) = \bigg( \frac{F(PWV=8.00)}{F(PWV=6.43)} \bigg) ^{5.4/(8.00-6.43)}\cdot \label{ratio_H2O} \end{eqnarray}(6)To derive a spectrum for pure water vapour (i.e. without the lines from other elements or molecules), TH2O(λ) is calculated only when water-vapour absorption is present (T ≤ 99%) in the water-vapour synthetic spectrum computed to measure the PWV. To avoid TH2O(λ) to become too noisy in the deepest absorption bands of water vapour, the spectrum is not calculated when water-vapour transmission becomes T ≤ 2.5%.

The transmission spectrum for O2 for airmass =1 is obtained by taking into account the variation of airmass AM between the two F spectra (Table 6), as done by Vidal-Madjar et al. (2010): TO2(λ)=(F(AM=2.01)F(AM=1.56))1/(2.011.56)·\begin{eqnarray} T_{\rm O_2}(\lambda) = \bigg( \frac{F(AM=2.01)}{F(AM=1.56)} \bigg) ^{1/(2.01{-}1.56)}\cdot \label{ratio_O2} \end{eqnarray}(7)And as for TH2O, TO2 is calculated only when O2 absorption is present (99.9% ≥ T ≥ 2%) in a AM = 1 MODTRAN5 synthetic spectrum.

Once TH2O and TO2 are known, we can correct our spectra simply by the division s=S(TH2O×TO2)AM,\begin{eqnarray} s = \frac{S}{(T_{\rm H_2O}\times T_{\rm O_2})^{AM}} \label{corr_AM}, \end{eqnarray}(8)where s is the corrected spectrum and S is a full-Moon or eclipse spectrum acquired at airmass AM. The final T(λ) spectrum for Paranal is shown in Fig. 7. The continuum shows neither the Rayleigh nor the ozone signature in this spectrum (on the opposite of what is obtained with HARPS, Fig. 4) because the UVES reduction pipeline automatically removed from all spectra tabulated values for the Rayleigh and ozone signatures at the corresponding airmass of the observation. Of course, the additional Rayleigh and ozone signatures present in the eclipse spectra and produced by the light path through the limb of the Earth atmosphere have not been removed by the UVES pipeline.

thumbnail Fig. 7

Paranal atmosphere transmission function T(λ) for AM = 1 as evaluated from full-Moon spectra taken with UVES before the eclipse, following the reduction process described in Sect. 6.2. The UVES reduction pipeline automatically removes from all individual spectra tabulated values for the Rayleigh, aerosol and ozone signatures at the corresponding airmass of the observation. These signatures therefore do not appear here as for HARPS (Fig. 4). A weak residual curvature of the continuum (5% maximum at 7800 Å, with no resemblance to a residual signature of ozone or Rayleigh) was fitted with a second-order function. The spectrum shown here is normalized to this fit to have a flat continuum at 1 (horizontal line). The O2γ, A, and B bands are clearly seen (in red), the other bands are due to H2O (in blue).

To calculate TH2O and TO2, the full-Moon spectra are aligned with a cross-correlation to match the radial velocity of the telluric lines.

7. Results and discussion

7.1. Results from Eq. (1) (method1)

From the HARPS data, we build five pairs of penumbra spectra suitable for Eq. (1), [98; 204], [125; 232], [66; 125], [66; 204], and [66; 232] where the numbers in brackets are the positions described in Table 3. The possible pair [66; 98] is not used: it produces a significantly noisier altitude profile probably because the positions in the penumbra are too close and the difference between the spectra is too small.

The penumbra spectra are affected by the wavelength-dependent solar limb darkening LD(λ) (Hestroffer & Magnan 1998). Equation (1) makes use of geometrical parameters S and L for two different solar crescents visible above Earth. The solar LD(λ) makes the crescents fainter than the mean solar brightness. The size of the solar crescents estimated from the penumbra spectra fluxes is thus underestimated if the penumbra fluxes are not corrected for the solar LD. To correct for this effect, we first evaluate the size of the solar crescents (parameters S and L) from the flux ratio (penumbra/full-Moon) at λ = 5798.8 Å which allows us to evaluate an average LD that affects each solar crescent following Hestroffer & Magnan (1998). For these crescent sizes and the LD obtained at λ = 5798.8 Å, a wavelength-dependant limb darkening LD(λ) is then re-evaluated from the low spectral resolution model of Hestroffer & Magnan (1998). It covers the spectral ranges [3033–3570] and [4160–10 990] Å, with a spectral resolution of a few tens (we do a linear interpolation between 3570 and 4160 Å). The model is then used to calculate the LD-corrected penumbra spectra EA(λ)/LD(λ) and EB(λ)/LD(λ) from which new S and L are evaluated for λ0 with Eq. (2). Equation (1) is then evaluated again to derive the final h(λ) profile.

thumbnail Fig. 8

Altitude profile with HARPS as obtained from Eq. (1) (method 1). A 10 Å bin is superimposed in light-blue over the full-resolution profile. The profile shows the large bump due to ozone Chappuis band around 6000 Å, the oxygen γ- and B-band at 6300 at 6900 Å  respectively, and also the water-vapour band around 6500 Å. The expected increased of h(λ) in the blue due to Rayleigh is visible as a negative slope between 4000 and 5000 Å. Oxygen and water-vapour absorption spectra at arbitrary scale are shown.

The mean altitude profile from the five pairs of fibers obtained from Eq. (1) is shown in Fig. 8. The profile is shifted upward by 30 km to have an altitude of 23.8 km in the [4520–4540] Å domain, as predicted by Ehrenreich et al. (2006). The profile shows the large bump due to ozone Chappuis band centered around 5900 Å and peaking here at about 38 km, and the oxygen γ and B bands at 6300 and 6900 Å, respectively. Water vapour also appears around 6500 Å at an altitude above 25 km. The increase of h(λ) in the blue due to Rayleigh scattering is visible as a negative slope between 4000 and 5000 Å (Ehrenreich et al. 2006; Kaltenegger & Traub 2009). The blue part of the profile shows many residual solar absorption lines that bias the binned profile towards lower altitudes. Numerous solar lines are visible, for example Hβ at 4861 Å and Hα at 6563 Å. The Fig. 9 is a zoom on the Na I doublet. The 2 Å-wide photospheric wings of the solar Na lines are still visible. The presence of all these solar features indicates that the solar signature is not properly removed from the h(λ) profile. This could result from the E/F ratios in Eq. (1). The penumbra spectra E indeed contain the spectrum from the solar crescent, that is, the edge of the Sun where solar absorption lines are deeper than in full-Moon spectra F, which correspond to a full disk of the Sun illuminating the Moon. We therefore tried to adjust a power p > 1 on F to compensate for the solar residuals with ratios E/(Fp) in Eq. (1) instead of E/F. But we found p to be very close to 1 and could not detect a significant difference between the average depth of the absorption lines in E and F.

thumbnail Fig. 9

Zoom in the Na I doublet region of the altitude profile from Fig. 8 with HARPS and method 1. The 2 Å-wide photospheric wings of the Na lines are still present, showing that the solar signature is not properly removed from this h(λ) profile. A residual solar line of Ni I is visible at 5892.9 Å (Wallace et al. 2011).

Vidal-Madjar et al. (2010) observe the same type of solar residuals in the SOPHIE data (residual wings of the Na I doublet) and the explanation invoked is the Ring effect. The depth of solar Fraunhofer lines in scattered light is lower than in direct sunlight. This effect is known as the Ring effect (Grainer & Ring 1962) and is mainly induced by rotational Raman scattering (Vountas et al. 1998). The line filling-in for forward scattering (as experienced by the photons that interact with the Earth atmosphere and are scattered towards the eclipsed Moon) is about 2.1% and 3.8% of the continuum for single and multiple scattering, respectively, for the Ca II K line (Vountas et al. 1998). Kattawar et al. (1981) report a filling-in of 2.5% at 6301.5 Å for a Fe-I line. Brinkmann (1968) calculates a filling-in of about 1.7 in the center to 1.9% in the wings for a hypothetical line at 4000 Å, underlining the continuous and smooth effect across the line profile, which is also measured by Pallamraju et al. (2000). Brinkmann (1968) reports observations showing a decrease of the filling-in with wavelength by a factor of 2 between 4383 and 6563 Å, and more recent measurements (Pallamraju et al. 2000, Fig. 3 therein) shows that the Ring effect decreases approximately linearly by a factor of 3 from 4285 to 6560 Å for a given line strength (normalized depth to continuum × width).

Several methods exist to correct for the Ring effect (Wagner et al. 2001). The simplest is that of Noxon et al. (1979), which consists of subtracting a small constant, close to 2%, from the spectrum. This was the solution applied in Vidal-Madjar et al. (2010) where a constant is empirically adjusted to remove the residual solar Na I lines (the constant was 4% and 0.6% of the continuum in the Na I region, for the deepest and the shallowest positions in the penumbra, respectively). We also correct this Ring effect in our data.

We observed that removing a Ring value xa from the shallowest spectrum in the penumbra and xb( > xa) to the deepest spectrum gives almost the same result as removing xb − xa to the deepest only (all spectra were corrected from the instrumental response, so Noxon’s approach applies correctly). The Ring signature is thus, as expected, stronger in the deeper spectrum (xb > xa), where collected photons have been more scattered in Earth atmosphere than photons collected in the shallower penumbra. In our empirical correction for the Ring effect, we therefore only apply a correction on the spectrum in the deepest position in the penumbra. We consider empirical Ring corrections following: i) Noxon et al. (1979); ii) Noxon et al. (1979) except in deep telluric bands (T < 0.5) where the Ring effect is negligible because Raman-scattered photons are strongly absorbed (Sioris & Evans 2000); iii) a λ-2 power law to mimic the chromatic dependance observed by Pallamraju et al. (2000) or Brinkmann (1968); iv) a λ-2 power law subtracted except in deep telluric bands (T < 0.5); and v) a subtraction of the reciprocal of the spectrum (Johnston & McKenzie 1989) that mimics the Ring spectrum. Here again, the obtained Ring spectrum is not removed from deep telluric lines (T < 0.5).

thumbnail Fig. 10

Altitude profile with HARPS and method 1 without Ring correction (upper panel) and different empirical Ring corrections. The amount of correction is 1% of the continuum at 5880 Å: i) Noxon et al. (1979); ii) Noxon et al. (1979) except in deep telluric bands (T < 0.5) where the Ring effect is negligible because Raman-scattered photons are strongly absorbed (Sioris & Evans 2000); iii) a λ-2 power law to mimic the chromatic dependance observed by Pallamraju et al. (2000) or Brinkmann (1968); iv) a λ-2 power law subtracted except in deep telluric bands (T < 0.5); and v) a subtraction of the reciprocal of the spectrum (Johnston & McKenzie 1989), which mimics the Ring spectrum. Here again, the obtained Ring spectrum is not removed from deep telluric lines (T < 0.5). The Na I doublet peaks at 5890 Å for all types of Ring corrections tested.

thumbnail Fig. 11

Same as Fig. 10 but zoomed-in on the Na I doublet region of the altitude profile without Ring correction (upper panel) and different empirical Ring corrections.

We observe that removing 1% of the continuum at 5880 Å to the deepest spectrum of each pair reduces the solar residuals over the full altitude profile (Fig. 10) and slightly reduces the Na wings (Fig. 11). The correction also leads to the rise of the Na lines to altitude values of 80 km, a value in agreement with the literature (Plane 2003; Moussaoui et al. 2010). Nevertheless, the width of the terrestrial Na I peaks are observed to be 0.4 Å. This is larger than the 0.2 Å observed for the other telluric lines, but might be induced by the presence of the strong Na I line in the solar spectrum. The extraction process of h(λ) is then more disturbed in the case of the Na I lines, easily inducing a change in the width, particularly because the solar lines cannot be perfectly positioned in all spectra since these have been Doppler-shifted to have the telluric lines aligned. An additional effect is the limb-darkening effect that changes the shape of the solar Na I lines emitted by different size crescent above the Earth limb. This is also why the broad Na I stellar wings are not perfectly corrected in the presented extraction procedure, which only includes a low spectral resolution limb-darkening model. Finally, we also note that only the Na I lines are detected in Fig. 10, not the strong Mg I triplet at 5170 Å nor Hα, for example, produce a peak in the altitude profile. We therefore conclude that the terrestrial Na layer is detected.

Interestingly, García-Muñoz et al. (2012) point out that the lunar eclipse of 16 August 2008 occurred only four days after the maximum of the Perseid meteor shower. It is indeed known that meteor showers can temporarily increase the amount of sodium in the atmosphere by a factor up of 10 (Plane 2003; Moussaoui et al. 2010). This corroborates the detection of the sodium layer by Vidal-Madjar et al. (2010) in August 2008. Nevertheless, the Geminid meteor shower maximum occurs on December 13. It is an active meteor shower, with a particle density similar to that of the Perseids6, that might have increased the sodium absorption during the eclipse of 21 December 2010. Nevertheless, our observations just show that the Na layer is indeed present, as during our first observations in August 2008. Since we have only one measurement, it is unfortunately not possible to detect any change of the atmospheric Na I content, either in 2008 or in 2010, and thus to confirm a link with any meteor shower.

thumbnail Fig. 12

Final altitude profile obtained with HARPS, method 1, and with a Noxon correction for the Ring effect, expect in the deep telluric lines (Sioris & Evans 2000). A 10 Å bin is superimposed in light-blue over the full-resolution profile; 30 and 100 Å bins are shown too, shifted upwards for clarity. The Na I peaks are visible near 5900 Å. Oxygen and water-vapour absorption spectra at arbitrary scale are shown. The altitude profile has been shifted to have an altitude of 23.8 km in the [4520–4540] Å domain, as predicted by Ehrenreich et al. (2006) (dots and dashes in green – dashes to improve the visibility of the curve around 6000 Å).

thumbnail Fig. 13

Zooms on the HARPS final altitude profile from method 1 (Fig. 12) with a Noxon correction for the Ring effect, except in the deep telluric lines (Sioris & Evans 2000). Oxygen and water vapour absorption spectra at arbitrary scale are shown. Oxygen shows up very clearly. Water vapour is barely visible around 6500 Å above 25 km and weak lines are visible in the higher Na I D1 D2 region. Note that our Noxon correction does not completely remove the wings of the Na I solar lines. The wings disappear with a stronger correction, but the two Na I peaks then reach an altitude well above 100 km which is not realistic.

It can be seen from Fig. 10 that the Ring correction is not homogeneous across the full spectral range: the Hα line remains undercorrected, while lines in the bluer part of the h(λ) profile are overcorrected. We therefore adopt the simplest Noxon approach, a subtraction of a constant, except in deep telluric lines (T < 0.5). The final result and the model from Ehrenreich et al. (2006) are shown in Fig. 12. The model is flatter than the observed profile, which peaks at about 35 km near 6000 Å, about 8 km above the model. This difference might be a refraction effect that amplifies the observed profile, as with method 3, although this effect is weaker here because we use only penumbra spectra. If the profile is calculated only with the two pairs [98; 204] and [125; 232] (Table 3), leaving out the spectrum from the deepest penumbra at only 66 arcsec from the umbra, the observed profile is indeed flatter (but noisier): it still peaks at 35 km, higher that the altitude expected from the model, but the lowest altitude at 6900 Å becomes 17 km instead of 12. We note that our observations of the 2008 eclipse agrees better with the model (Vidal-Madjar et al. 2010, Fig. 9 therein) although SOPHIE fibers were between 100 and 250 arcsec from the umbra edge, as for HARPS. The additional bins in Fig. 12 show that the B-band is barely detectable in the 30 Å bin, while the broad O3 Chappuis band and Rayleigh slope in the blue remain visible in the 100 Å bin. Figure 13 shows zooms in the Na I lines and the O2  γ and B band regions. Oxygen shows up very clearly, with peaks reaching h = 35–40 km. Some water-vapour peaks are visible but weak, which is not surprising because of the high altitude at which they are detected and also because of the lower strength of the water-vapour absorption in the HARPS wavelength range compared with λ > 7000 Å. We detect neither Li (6708 Å) nor Ca II (Hunten 1967; Plane 2003).

thumbnail Fig. 14

Final altitude profile obtained with UVES and method 1. No Ring correction is applied. A 10 Å bin is superimposed in light-blue on the full-resolution profile. Reference oxygen and water-vapour absorption spectra obtained with UVES are shown at arbitrary scale. The altitude profile has been shifted to have an altitude of 23.8 km in the [4520–4540] Å domain, as predicted by Ehrenreich et al. (2006) (green dashed line, no water vapour in the model).

With UVES, the coordinates of the telescope indicate that the telescope unfortunately did not move in the penumbra between the exposures in order to record spectra at 90 and 180 arcsec from the umbra edge as expected, but instead stayed 200 arcsec from the umbra (Table 5). The different recorded spectra are very similar and the result from Eq. (1) is noisy, especially in the blue. Residuals from spectral orders are also visible above 7000 Å. There is nevertheless a spectrum taken in UVES DIC2 configuration, while the telescope was thought to be in the umbra, but was too close to the Moon edge and left the moon about 20 s after the exposure started. The exposure was then stopped, but the telescope coordinates indicate that the telescope was in fact in the penumbra. The exposure time for that unexpected exposure is only about 21 and 38% of what was required for the BLUE and RED exposures, respectively. The altitude profile obtained with UVES is shown in Fig. 14. The overall shape agrees with the model of Ehrenreich et al. (2006), but the altitude observed in the UV is far higher than that of the model. Moreover, we had to shift the DIC1 BLUE profile by –110 km to connect it to the DIC2 BLUE profile. The DIC2 RED profile has also to be shifted by 35 km to be connected to the altitude obtained for the DIC1 RED profile. The oxygen bands are barely visible and the water-vapour bands appears in absorption. We did not consider any Ring correction for this noisy profile. Note also that the model of Ehrenreich et al. (2006) shown in the figures of this paper does not include water vapour.

7.2. Results from Eq. (3) (method 2)

The HARPS spectra are now used with Eq. (3) to calculate h from only one penumbra spectrum. The eclipse spectra are corrected for limb-darkening and Ring effect as for method 1. When the profile is calculated with the spectra in Eq. (3) aligned on the telluric lines, the solar lines become consequently slightly misaligned, which generates noise in the profile (Fig. 15 upper panel). We therefore build a composite profile where all lines – telluric and solar – are aligned. The composite uses two different profiles, one aligned on the telluric lines and a second aligned on the solar line. In the first profile, every misaligned solar line is recognized with the help of a synthetic G2 spectrum and replaced by its aligned residual taken from the second profile. Some solar lines are missing in the G2 spectrum in the most crowded telluric bands, therefore some residual noise due to uncorrected solar lines is still present in the composite inside these telluric bands. The procedure still slightly cleans the h profile (Fig. 15 lower panel). The overall shape of the profile is similar to the profile obtained with Eq. (1) shown in Fig. 12. Ozone appears about 6 km above the model. When zooming-in (Fig. 16), only the oxygen B and γ bands are correctly detected with peaks reaching h = 35–40 km. Water-vapour bands at 6500 and 5890 Å are not visible, nor the atmospheric Na I D-lines, which are maybe hidden in the strong solar residual lines. Note that the profile shown here is for the deeper position in the penumbra. The other positions lead to a noisier result, although the overall shape is correct.

thumbnail Fig. 15

Altitude profile with HARPS and Eq. (3) (method 2). Only the spectrum of the deepest penumbra, taken at 66 arcsec from the umbra, is used here (Table 3). Shallower penumbra spectra give similar but noisier profiles. A 10 Å bin is superimposed in light-blue on the full-resolution profiles. The upper panel shows the profile obtained with the HARPS spectra aligned on the telluric lines. Solar lines are consequently slightly misaligned and generate noise in the profile. The lower panel shows a composite profile showing the profile calculated with the spectra aligned on the telluric lines, expect when a solar line is present, where the profile from spectra aligned on the solar lines is used. Oxygen and water-vapour absorption spectra at arbitrary scale are shown. The green dashed line is the model of Ehrenreich et al. (2006).

thumbnail Fig. 16

Zooms on the HARPS altitude profile from method 2 shown in Fig. 15. Oxygen and water-vapour absorption spectra at arbitrary scale are shown to facilitate the identification of the spectral lines. Oxygen is correctly detected, while water vapour is not detected. Atmospheric Na I D-lines are not detected either, maybe they are hidden by the strong solar residual lines.

thumbnail Fig. 17

Altitude profile with UVES and method 2. A 10 Å bin is superimposed in light-blue on the full-resolution profile. The upper panel shows the profile obtained with the UVES spectra aligned on the telluric lines, with solar lines consequently slightly misaligned and thus generating noise in h. The lower panel shows a composite profile showing the profile aligned on the solar lines when no telluric line is present, and the profile aligned on the telluric lines when a telluric line is present. Oxygen and water-vapour absorption spectra at arbitrary scale are shown to facilitate the identification of the spectral lines. The model of Ehrenreich et al. (2006) is shown in green.

thumbnail Fig. 18

Zooms on the oxygen A and B-band regions of the UVES profil shown in Fig. 17 (method 2). An oxygen absorption spectrum at arbitrary scale is shown.

Equation (3) is also applied on UVES spectra to calculate h. As for HARPS, we compute a composite profile. The composite is the profile obtained with the spectra aligned on the solar lines in which the oxygen and water-vapour lines are replaced by the oxygen and water-vapour lines picked up in the profile obtained with the spectra aligned on the telluric lines. This procedure cleans the profile (Fig. 17 lower panel), although compromising the detection of atmospheric species different from oxygen and water. This may explain, at least partially, why the Na I D-lines do not appear clearly. The oxygen γ, B, and A bands are visible up to 80 km (zooms in Fig. 18), but water-vapour bands are not detected. We subtract a Ring effect from the UVES eclipse spectra of 0.5% of the continuum at 5880 Å. When we subtract a Ring effect of 1% as we did for HARPS, the profile does not change, except from the oxygen A and B bands which peak at 150 km.

Qualitatively speaking, the overall shape of the profile from 3000 to 7000 Å is in acceptable agreement with the model, with the ozone Chappuis bump at 6000 Å and a Rayleigh increase of h in the blue and near UV. The DIC2 BLUE altitude is set to h = 23.8 km in the [4520–4540] Å domain, as predicted by Ehrenreich et al. (2006). The other parts of the profile need to be shifted by about ± 6 km to be properly continuously connected to the DIC2 BLUE profile.

thumbnail Fig. 19

Altitude profiles obtained with HARPS and method 3 for different positions across the umbra edge, given by parameter c in Eq. (10) and shown in Fig. 20. The vertical axis is the measured altitude, set to 23.8 km at at [4520–4540] Å for the profile calculated for c = 0. Other profiles are shifted vertically for clarity by steps of 40 km. The profiles in this figure are not built at HARPS native resolution, but at a smaller 0.6 Å/pixel resolution. The light-blue lines are 10 Å bins.

7.3. Results from the umbra edge radiance profile (method 3)

With method 3, we make use of all available umbra and penumbra spectra. The penumbra spectra are corrected from the limb-darkening in the same way as for method 1 and 2. For each sampled wavelength in the HARPS spectra, we fit a sigmoid curve (Eq. (4)) through the eight available measured points (position, magnitude) in the umbra and deep penumbra region (Table 3). An example of the obtained fit is shown in Fig. 1. We now need to define the radius of the umbra edge re from which we build the Earth atmosphere thickness h(λ) simply given by h(λ)=re(λ)×rru+cte,\begin{eqnarray} h(\lambda) = r_{\rm e}(\lambda) \times \frac{r_\oplus } {r_{\rm u}} + cte, \label{eq:measured_h} \end{eqnarray}(9)where re is expressed in arcmin, ru = 42.7 arcmin is the umbra radius from the ephemeris, and r = 6371 km is the Earth radius. The constant cte is set to give h = 23.8 km at [4520–4540] Å.

thumbnail Fig. 20

Positions across the umbra edge used to build the altitude profiles shown in Fig. 19. The positions are defined by parameter c (Eq. (10)) and ranging from c =  −0.6 to 0.9. The definition for the umbra edge is taken for c = 0.6 (vertical solid line), at 1 arcmin in the penumbra side with respect to the inflection point, here fixed at abscissa =0. The umbra edge is only about one magnitude fainter than the penumbra at 4 arcmin from the inflection point.

As explained in Sect. 3, the inflection point abscissa a2 of the fit is not the best estimate of re (in our context of Earth transit analysis, at least). To explore the impact of re on h(λ), we write re such that ζ(re)=a0+c×a1.\begin{eqnarray} \zeta(r_{\rm e})= a_0+c\times a_1. \label{eq:re} \end{eqnarray}(10)For the inflection point, the coefficient c in front of half-amplitude a1 of the tanh function is c = 0 and ζ(a2) = a0. Positive values of c lead to positions on the penumbra side of the inflection point. Figure 19 shows a set of altitude profiles built for different positions across the umbra edge, with −0.6 < c < 0.9 as shown in Fig. 20. The profiles versus c change substantially: for c =  −0.6, the profile is built from data on the umbra side of the inflection point, and the bump from the ozone Chappuis band is almost invisible. The mean negative slope of the profile also underlines the effect of refraction that shrinks the Earth umbra at redder wavelengths. This effect decreases when c increases to positive values, with the profiles built from data on the penumbra side of the inflection point. The broad Chappuis band is detected. We finally define the umbra radius such as ζ(re) = a0 + 0.6   × a1 (c = 0.6). For this value, the noise in the altitude profile is close to minimum, and the bias from shallow umbra (that is, the mean negative slope of h(λ)) is reduced. With this definition, Fig. 20 puts the umbra edge at about 1 arcmin in the penumbra side with respect to the inflection point, and about one magnitude fainter than the penumbra at 4 arcmin from the inflection point. This choice is consistent with our reasoning developed Sect. 3.

thumbnail Fig. 21

Altitude profile obtained with HARPS and method 3. The left vertical axis is the measured altitude, set to 23.8 km at at [4520–4540] Å. The right vertical axis is the corrected altitude, if one considers the lowest region of h close to altitude 12 km as in Fig. 12 and again h = 23.8 km at [4520–4540] Å. The horizontal dotted line represents the altitude 23.8 km.The light-blue line is a 10 Å bin superimposed on the full-resolution profile. Oxygen and water-vapour absorption spectra at arbitrary scale are shown.

thumbnail Fig. 22

Zooms on the HARPS altitude profile from method 3 shown in Fig. 21. The left vertical axis is the measured altitude, set to 23.8 km at [4520–4540] Å. The right vertical axis is the corrected altitude, if one considers the lowest region of h close to altitude 12 km as in Fig. 12 and again h = 23.8 km at [4520–4540] Å. Oxygen and water-vapour absorption spectra at arbitrary scale are shown. Water vapour is not detected above 16 km, suggesting that the amount of water vapour is too low to absorb enough to be visible at this altitude. The sodium layer Na I is not detected.

thumbnail Fig. 23

Measured magnitude across the umbra edge with UVES at 10 000 Å. The horizontal error bars are 20 arcsec long and represent the typical pointing error of the telescope. Vertically, the 2σ error bars are the magnitude errors given by the UVES pipeline and the evaluated error in the local Moon albedo. The solid line is the model from Eq. (4) fitted through the data with parameter a3 set to the mean value of a3 found for the HARPS data. The vertical dotted line marks the edge of the umbra following our definition.

The result with HARPS and method 3 is shown in Figs. 21 and 22. Wide and narrow features of the Earth atmosphere are visible (Rayleigh, ozone Chappuis band, oxygen γ and B-bands). Water vapour is not detected, suggesting it is not present above 12 km. Compared with the profile obtained with method 1 (Fig. 12), the profile with method 3 shows a stronger negative slope from blue to red, leading to an altitude range about 3.3 times larger than with method 1. We attribute this effect to the fact that method 3 makes use of umbra spectra. The light in the umbra is reddened by Rayleigh scattering in the deep Earth atmospheric layers and refracted towards the umbra center. Therefore the umbra edge in red light is more refracted towards the umbra center than at blue wavelengths. The umbra radius (Earth radius on the Moon) is thus smaller in the red than in the blue, due to refraction and Rayleigh scattering. This effect amplifies the measured altitude range across the full visible spectrum. To correct for this effect, the right vertical axis of Fig. 21 is calibrated with the method 1 result, by setting the altitude for two different points, that is, the lowest region of h at altitude 12 km and h = 23.8 km at [4520–4540] Å.

thumbnail Fig. 24

Altitude profile obtained with UVES and method 3. The left vertical axis is the measured altitude, set to 23.8 km at [4520–4540] Å. The right vertical axis is the corrected altitude, if one considers the lowest region of h close to altitude 0 km and again h = 23.8 km at [4520–4540] Å. The horizontal dotted line represents the altitude 23.8 km. A 10 Å bin is superimposed in light-blue on the full-resolution profile, and a 100 Å bin in red is also shown, shifted upwards for clarity. The profile for λ ≤  ≈ 4300 Å is not shown, because the very weak umbra signal induces noise in the altitude profile which becomes noisy and h unexpectedly decreases. Oxygen and water-vapour absorption spectra at arbitrary scale are shown.

We process the UVES data in the same way. Nevertheless, since we have less observations points with UVES than with HARPS, it is not possible to keep four free parameters for the fitting function, the sigmoid ζ. We have to fix the a3 parameter to the mean value obtained with HARPS. The parameter a3 describes the steepness of the fit near the inflection point and is quite stable for HARPS at visible wavelengths: in the [6500–6900] Å range, its average value is a3=9.88×10-3\hbox{$\overline{a_3}=9.88\times10^{-3}$} with an RMS of σ = 2.2 × 10-4. We fix a3 to that mean value for UVES. An example of the fit through the UVES observations at 10 000 Å is shown in Fig. 23. Figure 24 shows the altitude profile obtained with UVES, showing the same features as for HARPS, with clear oxygen A and B bands, and water-vapour bands up to 10 000 Å. As with the HARPS profile (Fig. 21), the profile obtained with UVES displays a strong global negative slope. Here the profile is 4.5 times larger than with method 1. The right axis in Fig. 24 is calibrated considering that the lowest region of h is close to altitude 0 km and again that h = 23.8 km at [4520–4540] Å. The factor of 4.5 is larger than with HARPS (3.3) maybe because the slit positions with UVES are deeper in the umbra than with HARPS (see Tables 5 and 3).

In the blue (λ ≤    ≈ 4300 Å), the signal in the umbra spectra is weak and the altitude profile becomes noisy and unexpectedly decreases (this part of the profile is not shown in Fig. 24). The weak signal in the deepest lines of oxygen A-band and water vapour (9300–9500 Å) in the umbra spectra is also responsible for spikes above 100 km at these wavelengths. To cancel out most of these spikes, we set a flux threshold for the umbra spectra of 1% of the continuum below which the algorithm does not evaluate h. At last, to have h = 23.8 km at [4520–4540] Å, we subtract 37 km from the DIC2 BLUE profile. To have all parts of the profile connected in Fig. 24, we shift the DIC1 RED profile by –23 km and the DIC2 RED profile by +25 km.

The Figs. 25 and 26 show zooms in the UVES full profile shown in Fig. 24. In Fig. 25, the three oxygen bands are clearly detected. Note that the core of the oxygen line for the A-band could not be calculated because the signal in the line core is too weak. The sodium layer is not detected, and no water vapour appears above 19 km near the sodium lines. At λ > 7000 Å in Fig. 26, water vapour appears very well, even at an altitude higher than 20 km, because the water-vapour absorption bands become stronger in these longer wavelengths. Water vapour indeed presents a vertical distribution in the Earth atmosphere and is thus also present at higher altitudes but with lower volume densities. The altitude at which it is detected (where an optical thickness of about 1 is reached for a grazing line of sight) is of course higher in the stronger bands.

It is worth noting that the residuals of the Na I lines are very similar apart from the noise for HARPS and UVES (Figs. 22 and 25 respectively). This is expected because HARPS and UVES observed roughly the same solar crescents and thus experienced the same Doppler shift of the solar lines of the crescents that induced these residuals in the altitude profile. Note also that no correction of the Ring effect is implemented in method 3, either for HARPS or UVES. All spectra are indeed used here, all at different depths in the penumbra and umbra. A simple empirical approach cannot reasonably be considered. A physical model of the Ring effect is much needed, taking into account the path length through the Earth atmosphere for the different positions observed in the umbra and penumbra.

7.4. Possibilities with a future giant telescope

We analyze how our results can be extrapolated to what a future ground-based giant telescope might be able to detect in the visible transmission spectrum of an Earth-like exoplanet observed in transit. We use a full model of the future E-ELT, a 39 m telescope. The focal instrument has the same throughput, pixel sampling, and detector characteristics as the EPICS instruments. The model takes into account the photon noise from the star, the thermal and readout detector noises, the zodiacal noise, and the absorption of the local atmosphere above the telescope (see Hedelt et al. 2013 for a detailed description of the model).

We define a set of ten filters i) to sample the Chappuis band and the Rayleigh continuum; and ii) to isolate the most interesting narrow features in the profile shown in Fig. 24, that is, the oxygen A and B bands, and the water vapour around 9500 Å. The low spectral resolution 10 Å bin in this figure shows that the oxygen A-band peaks about 20 km above the lowest altitude around 10 200 Å. The Earth radius is thus 0.3% larger and the transit 0.6% deeper through a filter that isolates the oxygen A-band. The Earth transit depth is 8.3 × 10-5, and the transit variation at that spectral resolution for the O2 band represents 7.5 × 10-7 of the stellar flux.

According to Porto de Mello et al. (2006), there are 33 solar-type stars within 10 pc, 13 of which are reasonably similar to the Sun in terms mass, age, metallicity, and evolutionary status. Considering that the transit probability of an Earth twin in front of its star is about 1%, we understand that the probability of detection is low, even if all stars have a planet in their habitable zone. The Darwin All Sky Star Catalogue (Kaltenegger et al. 2010) lists 878 single FGK stars within 30 pc, a number fully consistent with the 33 solar-type stars numbered by Porto de Mello et al. (2006) for a 27 times smaller volume. Nevertheless, although the probability of transit detection significantly increases with 878 targets, the flux for a star at 30 pc will be 9 times fainter than at 10 pc meaning that 9 times more transits will be necessary to reach the same S/N as for a target at 10 pc.

Table 7

Transit of an Earth twin at 10 pc observed for 8 h by the E-ELT, at low spectral resolution and assuming no atmospheric perturbations.

We consider an earth at 10 pc transiting a G2V star observed during a full night of 8 h (an Earth-like 13 h single transit being longer than the night). The simulation assumes that the atmosphere above the telescope is perfectly stable. The simulation results are given in Table 7 and shown in Fig. 27. The result shows that transit photometry through the defined filters allows in principle the detection of the O2A-band with a detection level of 2.3  σ. The oxygen B and water-vapour bands will be more difficult to detect and will require more than one transit to reach a detection level >2σ. This is also true for the ozone Chappuis band, which is blended with the Rayleigh scattering.

It is encouraging to see that the E-ELT will beat the photon noise in only one single transit to identify oxygen in the visible range on a terrestrial planet transiting a solar-type star at 10 pc. But it will be extremely difficult, at least at low spectral resolution, to detect the 7.5 × 10-7 flux variation in the O2 band through the slightly fluctuating atmosphere above the telescope.

We considered broadband and narrow-band photometry, with a spectral resolution of up to 170. But it is worth noting that at a much higher spectral resolution (of the order of 105), the telluric lines become distinct from those of an exoplanet as soon as its radial velocity is sufficiently different from the BERV. Under these conditions, ground-based detection at high spectral resolution of O2 or H2O in an exoplanet atmosphere becomes easier, especially by cross-correlating the transit spectrum and an O2 or H2O spectrum, respectively (Vidal-Madjar et al. 2010; Snellen et al. 2013; Rodler & López-Morales 2014).

thumbnail Fig. 25

Zooms on the UVES altitude profile from method 3 shown in Fig. 24. The UVES oxygen spectrum at arbitrary scale is shown. The three oxygen bands are clearly detected. Note that we were unable to calculate the core of the oxygen line for the A-band because the signal in the line core is too weak. The sodium layer is not detected. Water-vapour absorption is too weak in the Na I region to appear above 19 km.

thumbnail Fig. 26

Zooms on the UVES altitude profile from method 3 shown in Fig. 24. The UVES water-vapour absorption spectrum at arbitrary scale is shown. Water vapour is very well detected, even at an altitude above 20 km at λ > 8000 Å where the water-vapour absorption is the strongest in the UVES wavelength domain.

8. Summary and conclusion

We have described the analysis of the observations of a lunar eclipse used as a proxy to observe Earth in transit. The results fully confirm those obtained with SOPHIE and the observation of the August 2008 eclipse (Vidal-Madjar et al. 2010) and extend the observed h(λ) profile towards the near infrared thanks to the UVES data. They agree well with the models (Ehrenreich et al. 2006; Kaltenegger & Traub 2009; Snellen et al. 2013; Bétrémieux & Kaltenegger 2013). The Rayleigh increase in the blue is visible. Biogenic oxygen and ozone are also well visible in the atmosphere altitude profiles, and are therefore relevant detectable biosignatures to look for during transits of Earth-like exoplanets in the [4000–10 500] Å range. The oxygen A and B bands show up clearly with HARPS or UVES.

Water vapour is well detected at λ > 8000 Å with UVES because water-vapour absorption is stronger at these wavelengths. Therefore even low densities of water vapour at these altitudes are detected. Water vapour is seen at an altitude of up to 20 km around 9500 Å. Moreover, Rayleigh scattering and ozone absorption are essentially absent above 8000 Å, therefore we are also able to probe deeper layers in the Earth atmosphere where low-altitude water vapour is also present at higher densities. In HARPS data at λ < 7000 Å  water vapour appears above 25 km with method 1 but remains below 25 km with method 2 and below 12 km with method 3, because water-vapour absorption is weaker in the visible range.

The results from UVES are nevertheless less reliable than those of HARPS, because the PWV was not stable during the beginning of the night when the calibration spectra were recorded on the full-Moon before the eclipse, which compromised the quality of the eclipse observation. Although we were able to recover the water-vapour content in the calibration spectra a posteriori, the results for the detection of water vapour in the altitude profiles with UVES are not coherent between method 1, 2 and 3. This underlines the difficulties of these observations which were made through the local atmosphere above the telescope, while our goal was to measure atmospheric features at a planetary scale. Lunar eclipse observations done from space, either at low or high spectral resolution, would not suffer from the telluric pollution, with the additional benefit of an access to UV and the prominent Hartley band of ozone.

Another bias in the recovering of the Earth atmosphere profile from a lunar eclipse is the effect of refraction in method 3 that makes use of refracted light from the umbra. We have shown that this effect can be partially corrected with the results from method 1, but the profiles – either from HARPS or UVES – still show a residual negative slope with respect to the profile from method 1 or 2.

thumbnail Fig. 27

Simulated observations with the E-ELT of an Earth twin at 10 pc transiting a G2V star, at low spectral resolution and assuming no atmospheric perturbations (Table 7). The simulated observations are shown in red. The horizontal bars represent the filter bandwidth and the vertical error bars are the ±1σ estimated error from the model. The solid line is the 10 Å bin UVES profile (method 3).

Finally, these results provide input for the preparation of future observations with the next generation of giant telescopes. We showed that the E-ELT observing in the visible range at low spectral resolution and in the absence of atmospheric perturbations is in principle able to detect the O2A-band in the atmosphere of a transiting Earth twin. It will be extremely challenging in practice through the atmosphere. If we overcome this difficulty, especially with high spectral resolution spectroscopy, the O2 signature may then be increasing from year to year (from transit to transit), as will the enthusiasm of the observers.


6

See for example the IMO (International Meteor Organization) data at http://www.imo.net/live/geminids2012/ and http://www.imo.net/live/perseids2012/

Acknowledgments

We thank Gaspare Lo Curto (HARPS), and Claudio Melo and Lorenzo Monaco (UVES) for their help in the preparation of these relatively unusual observations, as well as the teams who operated the telescopes that night. We also thank I. Snellen, and the referee, W. A. Traub, for their constructive comments on the manuscript. R. R. Querel acknowledges funding from Conicyt through Fondecyt grant 3120150.

References

  1. Bouchy, F., Hébrard, G., Udry, S., et al. 2009, A&A, 505, 853 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Brinkmann, R. T. 1968, ApJ, 154, 1087 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bétrémieux, Y., & Kaltenegger, L. 2013, ApJ, 772, L31 [NASA ADS] [CrossRef] [Google Scholar]
  4. Charbonneau, D., Brown, T. M., Noyes, R. W., & Gilliland, R. L. 2002, ApJ, 568, 377 [NASA ADS] [CrossRef] [Google Scholar]
  5. Combaz, J.-P., Dubois, J., & Poumeyrol, F. 1950, L’Astronomie, 64, 457 [NASA ADS] [Google Scholar]
  6. Danjon, A. 1925, L’Astronomie, 39, 272 [NASA ADS] [Google Scholar]
  7. Dekker, H., D’Odorico, S., Kaufer, A., Delabre, B., & Kotzlowski, H. 2000, SPIE, 4008, 534 [Google Scholar]
  8. Dubois, J. 1950, Ciel et Terre, 66, 36 [NASA ADS] [Google Scholar]
  9. Dubois, J. 1953, L’Astronomie, 67, 43 [NASA ADS] [Google Scholar]
  10. Dubois, J. 1962, L’Astronomie, 76, 188 [NASA ADS] [Google Scholar]
  11. Dubois, J., & Link, F. 1970, The Moon, 1, 462 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dubois, J., & Link, F. 1980, Moon and the Planets, 23, 109 [NASA ADS] [Google Scholar]
  13. Dubois, J., & Mianes, P. 1954, L’Astronomie, 68, 381 [NASA ADS] [Google Scholar]
  14. Dubois, J., & Poumeyrol, F. 1952, L’Astronomie, 66, 105 [NASA ADS] [Google Scholar]
  15. Dubois, J., & Poumeyrol, F. 1954, L’Astronomie, 68, 55 [NASA ADS] [Google Scholar]
  16. Dubois, J., & Rousseau, M. 1956, L’Astronomie, 70, 414 [NASA ADS] [Google Scholar]
  17. Dubois, J., & Rousseau, M. 1957, L’Astronomie, 71, 237 [NASA ADS] [Google Scholar]
  18. Eliason, E., Isbell, C., Lee, E., et al. 1999, Mission to the Moon, The Clementine UVVIS Global Mosaic, PDS Volumes USA NASA, PDS CL_4001-4078, http://www.lpi.usra.edu/lunar/tools/clementine/instructions/UVVIS_DIM_Info.html [Google Scholar]
  19. Ehrenreich, D., Tinetti, G., Lecavelier Des Etangs, A., Vidal-Madjar, A., & Selsis, F. 2006, A&A, 448, 379 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Ehrenreich, D., Vidal-Madjar, A., Widemann, T., et al. 2012, A&A, 537, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Grainger, J. F., & Ring, J. 1962, Nature, 193, 762 [NASA ADS] [CrossRef] [Google Scholar]
  22. García-Muñoz, A., & Pallé, E. 2011, J. Quant. Spectr. Rad. Transf., 112, 1609 [Google Scholar]
  23. García-Muñoz, A., Zapatero Osorio, M. R., Barrena, R., et al. 2012, ApJ, 755, 103, 10 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hayes, D. S., & Latham, D. W. 1975, ApJ, 197, 593 [NASA ADS] [CrossRef] [Google Scholar]
  25. Hedelt, P., von Paris, P., Godolt, M., et al. 2013, A&A, 553, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Hestroffer, D., & Magnan, C. 1998, A&A, 333, 338 [NASA ADS] [Google Scholar]
  27. Hunten, D. M. 1967, Space Sci. Rev., 6, 493 [NASA ADS] [CrossRef] [Google Scholar]
  28. Johnston, P. V., & McKenzie, R. L. 1989, J. Geophys. Res., 94, 3473 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kaltenegger, L., & Traub, W.A. 2009, ApJ, 698, 519 [NASA ADS] [CrossRef] [Google Scholar]
  30. Kaltenegger, L., Eiroa, C., & Fridlund, C. V. M. 2010, Astrophys. Space Sci., 326, 233 [NASA ADS] [CrossRef] [Google Scholar]
  31. Kattawar, G. W., Young, A. T., & Humphreys, T. J. 1981, ApJ, 243, 1049 [NASA ADS] [CrossRef] [Google Scholar]
  32. Link, F. 1969, Eclipse phenomena in astronomy (Berlin: Springer) [Google Scholar]
  33. Léger, A., Pirre, M., & Marceau, F. J. 1993, A&A, 277, 309 [NASA ADS] [Google Scholar]
  34. Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 [NASA ADS] [Google Scholar]
  35. McEwen, A. S. 1996, Lunar and Planetary Science, 27, 841 [Google Scholar]
  36. Moussaoui, N., Clemesha, B. R., Holzlöhner, R., et al. 2010, A&A, 511, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Noxon, J. F., Whipple, E. C. Jr.,, & Hyde, R. S. 1979, J. Geophys. Res., 84, 5047 [NASA ADS] [CrossRef] [Google Scholar]
  38. Owen, T. 1980, in Strategies for the Search for Life in the Universe, Proc. Meeting held 15, 16 August 1979 in Montreal, Canada, ed. M. D. Papagiannis, (Dordrecht: D. Reidel Publishing Co.), Astrophys. Space Sci. Libr., 83, 177 [Google Scholar]
  39. Pallamraju, D., Baumgardner, J., & Chakrabarti, S. 2000, Geophys. Res. Lett., 27, 1875 [NASA ADS] [CrossRef] [Google Scholar]
  40. Pallé, E., Osorio, M. R. Z., Barrena, R., Montañés-Rodríguez, P., & Martín, E. L. 2009, Nature, 459, 814 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  41. Perruchot, S., Kohler, D., Bouchy, F., et al. 2008 in Ground-based and Airborne Instrumentation for Astronomy II, eds. McLean, Ian, S., & Casali, Mark, M., Proc. SPIE, 7014, 17 [Google Scholar]
  42. Plane, J. 2003, Chem. Rev., 103, 4963 [Google Scholar]
  43. Porto de Mello, G., Del Peloso, E. F., & Ghezzi, L. 2006, Astrobiology, 6, 308 [NASA ADS] [CrossRef] [Google Scholar]
  44. Querel, R. R., Naylor, D. A., & Kerber, F. 2011, PASP, 123, 222 [NASA ADS] [CrossRef] [Google Scholar]
  45. Randich, S., Gratton, R., Pallavicini, R., Pasquini, L., & Carretta, E. 1999, A&A, 348, 487 [NASA ADS] [Google Scholar]
  46. Rauer, H., Gebauer, S., Paris, P. V., et al. 2011, A&A, 529, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Reaves, G., & Walker, M. F. 1952, PASP, 64, 15 [NASA ADS] [CrossRef] [Google Scholar]
  48. Rodler, F., & López-Morales, M. 2014, ApJ, 781, 54 [NASA ADS] [CrossRef] [Google Scholar]
  49. Schindler, T. L., & Kasting, J. F. 2000, Icarus, 145, 262 [NASA ADS] [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
  50. Schneider, J. 1992, in Frontiers of Life, Colloque interdisciplinaire du CNRS, Blois, France, 14–19 Oct. 1991, eds. J. Trân Thanh Vân, et al. (Éditions Frontières), 471 [Google Scholar]
  51. Schneider, J. 1994, Astrophys. Space Sci., 212, 321 [Google Scholar]
  52. Snellen, I. A. G., de Kok, R. J., le Poole, R., Brogi, M., & Birkby, J. 2013, ApJ, 764, 182 [NASA ADS] [CrossRef] [Google Scholar]
  53. Seager, S., Bains, W., & Hu, R. 2013, ApJ, 775, 104 [NASA ADS] [CrossRef] [Google Scholar]
  54. Sioris, C. E., & Evans, W. F. J. 2000, Geophys. Res. Lett., 64, 15 [Google Scholar]
  55. Vidal-Madjar, A., Lecavelier des Etangs, A., Désert, J.-M., et al. 2003, Nature, 422, 143 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  56. Vidal-Madjar, A., Désert, J.-M., Lecavelier des Etangs, A., et al. 2004, ApJ, 604, L69 [NASA ADS] [CrossRef] [Google Scholar]
  57. Vidal-Madjar, A., Arnold, L., Ehrenreich, D., et al. 2010, A&A, 523, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Vollmer, M., & Gedzelman, S. D. 2008, Appl. Opt., 47, H52 [NASA ADS] [CrossRef] [Google Scholar]
  59. Vountas, M., Rozanov, V. V., & Burrows, J. P. 1998, J. Quant. Spectr. Rad. Transf., 60, 943 [NASA ADS] [CrossRef] [Google Scholar]
  60. Wagner, T. K., Chance, U., Frieß, M., et al. 2001, Proc. of the 1st DOAS Workshop, Heidelberg, 13, 14 Sept., Heidelberg, Germany [Google Scholar]
  61. Wallace, L., Hinkle, K., Livingston, W. C., & Davis, S. P. 2011, ApJS, 195, 6 [NASA ADS] [CrossRef] [Google Scholar]
  62. Walker, M. F., & Reaves, G. 1957, PASP, 69, 153 [NASA ADS] [CrossRef] [Google Scholar]
  63. Winn J. N. 2009, in Transiting Planets, Proceedings of the International Astronomical Union, eds. F. Pont, et al., IAU Symp., 253, 99 [Google Scholar]

All Tables

Table 1

HARPS and UVES spectral ranges.

Table 2

Series of n Moon reference spectra taken with HARPS before the eclipse.

Table 3

HARPS exposures and fiber positions in the penumbra (negative values) and umbra (positive values) with respect to the umbra edge given by Guide-9.

Table 4

Series of three Moon reference spectra taken with UVES before the eclipse.

Table 5

UVES exposures or series of n exposures, and measured slit positions in the penumbra (negative values) and umbra (positive values) with respect to the umbra edge given by Guide-9.

Table 6

Measured PWV before the eclipse and at the end the observation, and reduced PWV for airmass =1.

Table 7

Transit of an Earth twin at 10 pc observed for 8 h by the E-ELT, at low spectral resolution and assuming no atmospheric perturbations.

All Figures

thumbnail Fig. 1

Measured magnitude across the umbra edge with HARPS at 5500 Å. Crosses and diamonds stand for A and B fiber, respectively. The horizontal error bars are 2σ = 20 arcsec long and represent the typical pointing error of the telescope. Vertically, the 2σ error bars are the magnitude errors evaluated from the S/N of the spectra at that wavelength and the evaluated error in the local Moon albedo. The solid line is the model from Eq. (4) fitted through the data. The vertical dotted line marks our definition of the umbra edge. See Sect. 7.3 for details.

In the text
thumbnail Fig. 2

Earth and Sun seen from the Moon deep penumbra where HARPS fiber A was positioned at 07:06:35 UT, at the beginning of the 300 s exposure. Image with Guide-9.

In the text
thumbnail Fig. 3

Example of a path of the HARPS fibers across the reconstructed Moon image, during the 300 s exposure in the penumbra. The darkening from the penumbra and umbra are not superimposed on the image here, since we are only interested in recovering the moon reflectivity. The A and B fibers remained at 232 and 125 arcsec from the umbra edge, while the fibers path across the Moon was 164 arcsec long (length of the yellow lines). The mean reflectivity is measured along the path, taking into account a 10 arcsec pointing error in both RA and Dec directions. The terrain reflectivity fluctuates by σ = 21% along the path for Fiber A and σ = 15% for Fiber B in this highly contrasted lunar region. The error on the mean reflectivity due to pointing uncertainty is estimated to be about σ = 4%. The figure shows that while Fiber A went through the dark region of Mare Nectaris, Fiber B remained on a region measured 39% brighter between Catharina and Cyrillus craters, pointing out the relevance of correcting the flux for lunar local reflectivity.

In the text
thumbnail Fig. 4

La Silla atmosphere transmission function T(λ) for AM = 1 as evaluated from full-Moon spectra taken with HARPS just before the eclipse (Table 2). The spectrum is normalized at 6866 Å. To verify this experimental T(λ) spectrum, a transmission fit based on the model of Hayes & Latham (1975) for the Rayleigh and aerosol components normalized at 6866 Å is displayed (blue dashed line). We also add an ozone transmission spectrum to the model (light-blue solid line) that correctly fits the observed wide Chappuis band. The best least-squares fit for Rayleigh and ozone over La Silla that night is obtained for an airmass of AM = 0.77 for ozone and Rayleigh, and an optical depth of 0.181 at 5320 Å for the aerosols, between 4000 and 6866 Å. The fit for ozone is made only where O2 and H2O absorption bands are absent. The O2γ and B absorption bands are clearly seen, around 6300 and 6900 Å respectively, the other deep bands are due to H2O. The small gap in HARPS data is visible between 5300 and 5340 Å (Table 1).

In the text
thumbnail Fig. 5

Example of a path of the UVES slit across the reconstructed Moon image during the 600 s exposure in the umbra. The slit, 10 arcsec long and perpendicular to the path, remains at 936 arcsec from the umbra edge, while the slit path over the Moon was 330 arcsec long (length of the yellow line). The mean reflectivity is measured along the path, taking into account a 10 arcsec pointing error in both RA and Dec directions. The terrain reflectivity fluctuates by σ = 13% along the path in this bright lunar region, and the error on the mean reflectivity due to pointing uncertainty is estimated to be about σ = 4%.

In the text
thumbnail Fig. 6

Part of the spectrum showing the ratio of a high airmass over a low airmass UVES full-Moon spectra. The oxygen B and A bands at 6900 and 7600 Å  respectively, are in absorption as expected, but the water vapour bands at 7200 and 8200 Å appear in emission, showing that the amount of water vapour along the line of sight was higher when the airmass was lower, indicating that the amount of water vapour above Paranal increased during the first part of the night, before the eclipse. To calculate this ratio, the spectra are aligned by cross-correlation to match the radial velocity of the telluric lines, so the solar lines do not properly cancel out in the ratio, and appear with a wing in emission and the other in absorption.

In the text
thumbnail Fig. 7

Paranal atmosphere transmission function T(λ) for AM = 1 as evaluated from full-Moon spectra taken with UVES before the eclipse, following the reduction process described in Sect. 6.2. The UVES reduction pipeline automatically removes from all individual spectra tabulated values for the Rayleigh, aerosol and ozone signatures at the corresponding airmass of the observation. These signatures therefore do not appear here as for HARPS (Fig. 4). A weak residual curvature of the continuum (5% maximum at 7800 Å, with no resemblance to a residual signature of ozone or Rayleigh) was fitted with a second-order function. The spectrum shown here is normalized to this fit to have a flat continuum at 1 (horizontal line). The O2γ, A, and B bands are clearly seen (in red), the other bands are due to H2O (in blue).

In the text
thumbnail Fig. 8

Altitude profile with HARPS as obtained from Eq. (1) (method 1). A 10 Å bin is superimposed in light-blue over the full-resolution profile. The profile shows the large bump due to ozone Chappuis band around 6000 Å, the oxygen γ- and B-band at 6300 at 6900 Å  respectively, and also the water-vapour band around 6500 Å. The expected increased of h(λ) in the blue due to Rayleigh is visible as a negative slope between 4000 and 5000 Å. Oxygen and water-vapour absorption spectra at arbitrary scale are shown.

In the text
thumbnail Fig. 9

Zoom in the Na I doublet region of the altitude profile from Fig. 8 with HARPS and method 1. The 2 Å-wide photospheric wings of the Na lines are still present, showing that the solar signature is not properly removed from this h(λ) profile. A residual solar line of Ni I is visible at 5892.9 Å (Wallace et al. 2011).

In the text
thumbnail Fig. 10

Altitude profile with HARPS and method 1 without Ring correction (upper panel) and different empirical Ring corrections. The amount of correction is 1% of the continuum at 5880 Å: i) Noxon et al. (1979); ii) Noxon et al. (1979) except in deep telluric bands (T < 0.5) where the Ring effect is negligible because Raman-scattered photons are strongly absorbed (Sioris & Evans 2000); iii) a λ-2 power law to mimic the chromatic dependance observed by Pallamraju et al. (2000) or Brinkmann (1968); iv) a λ-2 power law subtracted except in deep telluric bands (T < 0.5); and v) a subtraction of the reciprocal of the spectrum (Johnston & McKenzie 1989), which mimics the Ring spectrum. Here again, the obtained Ring spectrum is not removed from deep telluric lines (T < 0.5). The Na I doublet peaks at 5890 Å for all types of Ring corrections tested.

In the text
thumbnail Fig. 11

Same as Fig. 10 but zoomed-in on the Na I doublet region of the altitude profile without Ring correction (upper panel) and different empirical Ring corrections.

In the text
thumbnail Fig. 12

Final altitude profile obtained with HARPS, method 1, and with a Noxon correction for the Ring effect, expect in the deep telluric lines (Sioris & Evans 2000). A 10 Å bin is superimposed in light-blue over the full-resolution profile; 30 and 100 Å bins are shown too, shifted upwards for clarity. The Na I peaks are visible near 5900 Å. Oxygen and water-vapour absorption spectra at arbitrary scale are shown. The altitude profile has been shifted to have an altitude of 23.8 km in the [4520–4540] Å domain, as predicted by Ehrenreich et al. (2006) (dots and dashes in green – dashes to improve the visibility of the curve around 6000 Å).

In the text
thumbnail Fig. 13

Zooms on the HARPS final altitude profile from method 1 (Fig. 12) with a Noxon correction for the Ring effect, except in the deep telluric lines (Sioris & Evans 2000). Oxygen and water vapour absorption spectra at arbitrary scale are shown. Oxygen shows up very clearly. Water vapour is barely visible around 6500 Å above 25 km and weak lines are visible in the higher Na I D1 D2 region. Note that our Noxon correction does not completely remove the wings of the Na I solar lines. The wings disappear with a stronger correction, but the two Na I peaks then reach an altitude well above 100 km which is not realistic.

In the text
thumbnail Fig. 14

Final altitude profile obtained with UVES and method 1. No Ring correction is applied. A 10 Å bin is superimposed in light-blue on the full-resolution profile. Reference oxygen and water-vapour absorption spectra obtained with UVES are shown at arbitrary scale. The altitude profile has been shifted to have an altitude of 23.8 km in the [4520–4540] Å domain, as predicted by Ehrenreich et al. (2006) (green dashed line, no water vapour in the model).

In the text
thumbnail Fig. 15

Altitude profile with HARPS and Eq. (3) (method 2). Only the spectrum of the deepest penumbra, taken at 66 arcsec from the umbra, is used here (Table 3). Shallower penumbra spectra give similar but noisier profiles. A 10 Å bin is superimposed in light-blue on the full-resolution profiles. The upper panel shows the profile obtained with the HARPS spectra aligned on the telluric lines. Solar lines are consequently slightly misaligned and generate noise in the profile. The lower panel shows a composite profile showing the profile calculated with the spectra aligned on the telluric lines, expect when a solar line is present, where the profile from spectra aligned on the solar lines is used. Oxygen and water-vapour absorption spectra at arbitrary scale are shown. The green dashed line is the model of Ehrenreich et al. (2006).

In the text
thumbnail Fig. 16

Zooms on the HARPS altitude profile from method 2 shown in Fig. 15. Oxygen and water-vapour absorption spectra at arbitrary scale are shown to facilitate the identification of the spectral lines. Oxygen is correctly detected, while water vapour is not detected. Atmospheric Na I D-lines are not detected either, maybe they are hidden by the strong solar residual lines.

In the text
thumbnail Fig. 17

Altitude profile with UVES and method 2. A 10 Å bin is superimposed in light-blue on the full-resolution profile. The upper panel shows the profile obtained with the UVES spectra aligned on the telluric lines, with solar lines consequently slightly misaligned and thus generating noise in h. The lower panel shows a composite profile showing the profile aligned on the solar lines when no telluric line is present, and the profile aligned on the telluric lines when a telluric line is present. Oxygen and water-vapour absorption spectra at arbitrary scale are shown to facilitate the identification of the spectral lines. The model of Ehrenreich et al. (2006) is shown in green.

In the text
thumbnail Fig. 18

Zooms on the oxygen A and B-band regions of the UVES profil shown in Fig. 17 (method 2). An oxygen absorption spectrum at arbitrary scale is shown.

In the text
thumbnail Fig. 19

Altitude profiles obtained with HARPS and method 3 for different positions across the umbra edge, given by parameter c in Eq. (10) and shown in Fig. 20. The vertical axis is the measured altitude, set to 23.8 km at at [4520–4540] Å for the profile calculated for c = 0. Other profiles are shifted vertically for clarity by steps of 40 km. The profiles in this figure are not built at HARPS native resolution, but at a smaller 0.6 Å/pixel resolution. The light-blue lines are 10 Å bins.

In the text
thumbnail Fig. 20

Positions across the umbra edge used to build the altitude profiles shown in Fig. 19. The positions are defined by parameter c (Eq. (10)) and ranging from c =  −0.6 to 0.9. The definition for the umbra edge is taken for c = 0.6 (vertical solid line), at 1 arcmin in the penumbra side with respect to the inflection point, here fixed at abscissa =0. The umbra edge is only about one magnitude fainter than the penumbra at 4 arcmin from the inflection point.

In the text
thumbnail Fig. 21

Altitude profile obtained with HARPS and method 3. The left vertical axis is the measured altitude, set to 23.8 km at at [4520–4540] Å. The right vertical axis is the corrected altitude, if one considers the lowest region of h close to altitude 12 km as in Fig. 12 and again h = 23.8 km at [4520–4540] Å. The horizontal dotted line represents the altitude 23.8 km.The light-blue line is a 10 Å bin superimposed on the full-resolution profile. Oxygen and water-vapour absorption spectra at arbitrary scale are shown.

In the text
thumbnail Fig. 22

Zooms on the HARPS altitude profile from method 3 shown in Fig. 21. The left vertical axis is the measured altitude, set to 23.8 km at [4520–4540] Å. The right vertical axis is the corrected altitude, if one considers the lowest region of h close to altitude 12 km as in Fig. 12 and again h = 23.8 km at [4520–4540] Å. Oxygen and water-vapour absorption spectra at arbitrary scale are shown. Water vapour is not detected above 16 km, suggesting that the amount of water vapour is too low to absorb enough to be visible at this altitude. The sodium layer Na I is not detected.

In the text
thumbnail Fig. 23

Measured magnitude across the umbra edge with UVES at 10 000 Å. The horizontal error bars are 20 arcsec long and represent the typical pointing error of the telescope. Vertically, the 2σ error bars are the magnitude errors given by the UVES pipeline and the evaluated error in the local Moon albedo. The solid line is the model from Eq. (4) fitted through the data with parameter a3 set to the mean value of a3 found for the HARPS data. The vertical dotted line marks the edge of the umbra following our definition.

In the text
thumbnail Fig. 24

Altitude profile obtained with UVES and method 3. The left vertical axis is the measured altitude, set to 23.8 km at [4520–4540] Å. The right vertical axis is the corrected altitude, if one considers the lowest region of h close to altitude 0 km and again h = 23.8 km at [4520–4540] Å. The horizontal dotted line represents the altitude 23.8 km. A 10 Å bin is superimposed in light-blue on the full-resolution profile, and a 100 Å bin in red is also shown, shifted upwards for clarity. The profile for λ ≤  ≈ 4300 Å is not shown, because the very weak umbra signal induces noise in the altitude profile which becomes noisy and h unexpectedly decreases. Oxygen and water-vapour absorption spectra at arbitrary scale are shown.

In the text
thumbnail Fig. 25

Zooms on the UVES altitude profile from method 3 shown in Fig. 24. The UVES oxygen spectrum at arbitrary scale is shown. The three oxygen bands are clearly detected. Note that we were unable to calculate the core of the oxygen line for the A-band because the signal in the line core is too weak. The sodium layer is not detected. Water-vapour absorption is too weak in the Na I region to appear above 19 km.

In the text
thumbnail Fig. 26

Zooms on the UVES altitude profile from method 3 shown in Fig. 24. The UVES water-vapour absorption spectrum at arbitrary scale is shown. Water vapour is very well detected, even at an altitude above 20 km at λ > 8000 Å where the water-vapour absorption is the strongest in the UVES wavelength domain.

In the text
thumbnail Fig. 27

Simulated observations with the E-ELT of an Earth twin at 10 pc transiting a G2V star, at low spectral resolution and assuming no atmospheric perturbations (Table 7). The simulated observations are shown in red. The horizontal bars represent the filter bandwidth and the vertical error bars are the ±1σ estimated error from the model. The solid line is the 10 Å bin UVES profile (method 3).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.