Speckle lifetime in XAO coronagraphic images: temporal evolution of SPHERE coronagraphic images
Abstract
The major source of noise in highcontrast imaging is the presence of slowly evolving speckles that do not average with time. The temporal stability of the pointspreadfunction (PSF) is therefore critical to reach a high contrast with extreme adaptive optics (xAO) instruments. Understanding on which timescales the PSF evolves and what are the critical parameters driving the speckle variability allow to design an optimal observing strategy and data reduction technique to calibrate instrumental aberrations and reveal faint astrophysical sources. We have obtained a series of 52 min, AOcorrected, coronagraphically occulted, highcadence (1.6Hz), Hband images of the star HR 3484 with the SPHERE (SpectroPolarimeter Highcontrast Exoplanet REsearch [Beuzit2008]) instrument on the VLT. This is a unique data set from an xAO instrument to study its stability on timescales as short as one second and as long as several tens of minutes. We find different temporal regimes of decorrelation. We show that residuals from the atmospheric turbulence induce a fast, partial decorrelation of the PSF over a few seconds, before a transition to a regime with a linear decorrelation with time, at a rate of several tens parts per million per second (ppm/s). We analyze the spatial dependence of this decorrelation, within the wellcorrected radius of the adaptive optics system and show that the linear decorrelation is faster at short separations. Last, we investigate the influence of the distance to the meridian on the decorrelation.
Speckle lifetime in XAO coronagraphic images: temporal evolution of SPHERE coronagraphic images
Keywords: Extreme Adaptive Optics,Speckles, High Contrast
1 Introduction
Direct imaging of exoplanets and circumstellar environments requires to probe very close separations, within a fraction of an arcsecond from the host star, and to supress the light from the star [Oppenheimer2003]. A new generation of instruments have entered into operations in the last years, known as extreme adaptive optics (xAO) instruments, to push further down the achievable contrast, with the primary goal to better characterize the population of young gas giants orbiting beyond a few astronomical units of their host stars [Beuzit2008, Macintosh2014, Close2012, Guyon2010, Dekany2013, Esposito2011]. These instruments involve a highperformance adaptive optics system with a dense deformable mirror and sensitive wavefront sensors to reach high Strehl ratios. They use in conjunction advanced coronagraphs to reject the starlight [Mawet2012]. However, residual starlight still contaminates the area of interest. It comes from both the diffracted light of the telescope pupil and the wavefront errors due to uncorrected atmospheric perturbations or imperfect optics within the light path. The diffraction pattern can be attenuated by specific diffraction control techniques and is stable in time for pupilstabilized observations. However the second source of residual starlight is more difficult to address. Wavefront errors create speckles in the image plane that mimic point sources and degrade the contrast [Racine1999]. These speckles vary temporally on different timescales, from a few milliseconds for noncorrected atmospheric perturbations up to several minutes, hours or even days for quasistatic speckles slowly variable with the temperature, the telescope flexures or the movement of the optical elements in the light beam such as the derotator or the atmospheric dispersion corrector. In addition they can interfere with the static pattern to create pinned speckles [Bloemhof2001]. Eventually, those quasistatic speckles represent the major source source of noise at short separations because they are correlated in time and do not average out during an observing sequence [Hinkley2007], unlike other sources of noise [Macintosh2005].
To calibrate them, various observation strategies have been designed, based on differential imaging. One of the most successful is known as Angular Differential Imaging [Marois2006] (ADI), and relies on the differential motion between the pupil and the sky to disentangle speckles from an astrophysical signal. A reference image is constructed as a linear combination of frames recorded at different parallactic angles from the science frame. Because this reference frame is built from images taken at different times, this method is intrisically subject to the variability of the speckles and will work best for a highly stable temporal sequence. We therefore aim here at analysing empirically how coronagraphic images evolve with time in the wellcorrected region of an xAO instrument, based on onsky images. This complements the analysis already made base on laboratory measurements[Martinez2012, Martinez2013].
2 Instrumental Setup
2.1 Measurements
The observations were carried out on March 3, 2015 between 03:53 and 04:35UT on the bright star HR 3484 of magnitude and . They used the pupiltracking mode of VLT / SPHERE, to keep the aberrations as stable as possible. To reach a high contrast, we made use of the coronagraphic combination N_ALC_YJH_S corresponding to an apodizer, a Lyot mask of diameter 185 mas and an undersized Lyot stop to block the starlight rejected off the mask and to cover the telescope spiders.
We used the IRDIS subsystem [Dohlen2008] in dualband imaging [Langlois2010]. The IRDIS imager splits the incoming light in two channels where the dualband filters H2 and H3 were inserted (m, nm, m, nm). IRDIS provides a 1111 field of view with a pixel scale of 12.25 mas. The star was observed after meridian passage, during 50 minutes under good and stable atmospheric conditions. The seeing was between 0.3 and 0.5, the coherence time derived from DIMM measurements [Sarazin2002] was between 3 and 4ms. The real time controller of the SPHERE AO system indicated a wind speed of 3 to 4m/s and derived a Strehl ratio of 92% to 93%. The spatial filter upstream the ShackHartmann wavefront sensor was set to its smallest aperture to get the least aliasing. In order to obtain a high frame rate, the detector was windowed to read only 128 rows (out of a total of 1048), which enabled to use a detector integration time (DIT) of 0.1s. 5000 frames were recorded in this setup for a total of 3087s. This corresponds to a parallactic angle variation of . We rejected the first 50 frames for which the AO loop had not reached a stationary regime. Because of the time to read the detector, the actual frame rate is one image per 617ms or 1.62Hz. The median frame of this sequence is shown in Fig. 1 left. This deep coronagraphic sequence was followed by an acquisition of 500 frames on an empty sky region for measurement of the background, and by a pointspread function (hereafter PSF) measurement out of the coronagraphic mask with the neutral density filter ND_2.
To compare these onsky measurements with data obtained without disturbance by the atmosphere, we recorded an additional sequence with the same instrumental setup but with the internal lamp and the instrument shutter closed, on November 24, 2015 between 01:08 and 02:17 UT. The sequence consists of 5740 frames with the AO loop closed. To be closer to onsky conditions, the derotator was turned on to simulate pupiltracking of a star at a declination of above Paranal observed at an hour angle between 20 min and +40 min. The detector was read at the same speed in the same windowed region as for the onsky observations. The median image of the sequence is shown in Fig. 1 right. The only difference with onsky observations in the instrument after the optical fiber injecting the light source is the atmospheric dispersion corrector, which is turned off while using the internal lamp.
2.2 Methodology
 Definition

For each pair of coronagraphic images , we computed the correlation coefficient as follow:
(1) where and are the mean and standard deviation of the region S for the frame taken at time . This is also known as the Pearson’s correlation coefficient. The normalisation by ensures a unit value for . One of the challenge of this analysis is get a Signal to Noise ratio (S/N) large enough to allow to draw meaningfull conclusions. This was done by chosing regions S encompassing a sufficient number of pixels and by resorting to temporal binning in order to average the coeffficient correlations, as detailed below. In Appendix C) we show that such a definition of the correlation is directly related to the contrast obtained after subtracting one image of the pair from the other.
 Global approach over

In the first place, we computed the correlation coefficients between each pair of images in a single region , encompassing all the pixels shown in Fig. 1, exluding those within 8px ( mas) of the central source, where the light is blocked by the coronagraphic mask. These coefficients can be displayed in the form of a symmetric matrix in Fig. 3 (left panel).
Figure 3: Left: Visualisation of the correlation coefficients computed between each pair of images in the form of the squared matrix of dimension . The matrix is symmetric, therefore only the elements below the diagonal are shown here (). When one moves towards the lower right corner, the time ellapsed between the 2 images increases. The dark lines (horizontal or vertical) correspond to frames poorly correlated with the others. The colour scale is nonlinear (squareroot) from 0 (black) to 1 (white). These coefficients are computed in the whole region . Right: Histogram of the correlation coefficients, normalized to an integral of one, the red vertical line is the median (0.7). All correlation coefficients span between 0.2 and 1, with a median of 0.70, as shown in the histogram of Fig. 3 (right). Still, to retrieve a significant trend from this histogram, additional processing and signal extraction is needed and this is the goal of the temporal binning explained below.
 Temporal binning

To get a onedimensional evolution of the correlation coefficients with time, we eventually averaged all coefficients calculated beween pairs separated by an equal time step. This enables to ease the interpretation by having a onedimensional curve instead of a matrix and more importantly to get a higher S/N on the correlation by averaging a large number of values for 4950 frames. There are 4948 pairs separated by the elementary time step of 0.1 s and this number decreases linearily as time goes by. Visually, this means that the twodimensional matrix of Fig. 3 (left and middle panels) was averaged along the diagonal direction to be turned into a onedimensional array of correlation coefficients that can directly be plotted against the time ellapsed between the frames. This is shown in Fig. 4.
Figure 4: Evolution of the correlation coefficients calculated over the whole region , with time. A linear fit was performed between 5 and 15 min.  Spatial analysis in subregions

In a second time, to get some spatial resolution for a zonal analysis, the image was divided in 12 by 12 squared regions of 10px wide (122mas or ), as shown in Fig. 5, written with ,. For each subimage, we computed the correlation matrix and averaged it to get a onedimensional plot similar to that shown in Fig. 4 for the whole image. These curves are shown in App. A.
3 Results
3.1 Different regimes of decorrelation
The evolution of the correlation over time (Fig. 4) shows a linear regime of decorrelation over a timescale up to 20 min. The speed of decorrelation, computed as a linear fit between 5 and 15 min reaches 25 parts per million per second (ppm/s) when the analysis is done in the whole region . After 25 min, the correlation increases again. To understand this behaviour, we carried out the same analysis after subtracting the median frame from the sequence in order to remove the static pattern of the coronagraphic image that creates the overall offset of the correlation by about 0.7. This is shown in Fig. 6 (red data points). In such a case, the correlation coefficients are much lower and the decrease in the correlation is much steeper (73 ppm/s) because this is a direct measure of the decorrelation of the residual, quasistatic, speckles. This confirms that the decorrelation observed in the nonmediansubtracted sequence is real. In addition, we checked that we are not dominated by the detector or background noise (see App/B.1 for a verification). The regime where the correlation increases again after 25min corresponds to a regime where the correlation computed on the mediansubtracted sequence is negative. This means that the mediansubtracted frames are anticorrelated. This is also visible in the individual frames of the sequence, such as that shown in Fig. 2 : because the speckle intensity varies over time, negative regions appear after subtraction of the temporal median of the sequence, creating an overall negative correlation coefficient. We refrained from interpretating this regime because the speckle pattern have evolved signifcantly on the globabl scale of the region .
Interestingly, there is a second decorrelation regime, corresponding to a much steeper decrease in the correlation coefficient over a smaller timescale of a few seconds. This is shown in Fig. 7. An exponential decay was adjusted to the measurements in the form of the equation
(2) 
with s, and .
3.2 Zonal analysis: a faster linear decorrelation regime at small separations
The same analysis of the evolution of the correlation coefficient with time was made for the subregions for . We also see a linear decorrelation regime for ellapsed times between a few minutes and 30 minutes. An overall picture of the decorrelation plots is shown in Appendix A. In a second step, we performed a linear fit, as was done previously. We did not subtract the temporal median of each subregion. The parameters of the fit are plotted in Fig. 8 as a function of the separation between the star and the center of the subregion . This plot shows that in the linear regime, the decorrelation is faster close to the optical axis, with a speed exceeding 30 ppm/s for some subregions. We illustrate this trend by the black line in Fig. 8 (top), with a rate of decrease of ppm/s/. Beyond 500 mas no significant change is observerd and the speckles seem to evolve at a rate of ppm/s. The intercept of the linear fit (Fig. 8 bottom) is overall following the profile of the coronagraphic image, which traces the fact that in the nonmedian subtracted sequence, the static part of the sequence has a major influence on the absolute value of the correlation coefficient (in case we subtract the temporal median of the sequence, the top curve shows a very similar trend but the bottom curve is flatter).
3.3 Zonal analysis: a steep decorrelation independant of the separation
In a third step, we analyzed the evolution of the correlation in the first minute. We can also notice in each subregion the fast decorrelation regime visible in the whole region on Fig. 7. Therefore, we fitted for each subregion an exponential decay law, parametrized by Eq. 2. The parameters of the fit are plotted in Fig. 9 as a function of the separation between the star and the center of the subregion . In this case, we conclude that the characterstic time does not depend significantly on the separation, and is about 3 ms (with slightly more scatter at larger separations because the S/N on speckles is smaller). A linear fit (black line in Fig. 9 top) is indeed compatible with a zero slope. A similar characteristic time was also found for the whole region , meaning that there is no zonal effect for the time it takes for the correlation to drop and reach its value after one minute. On the other hand, the amplitude and offset of the exponential law do significantly vary with the separation, which is easily explained because the correlations start with higher values close to the optical axis, as already stated. The same effect is also visible in the intercept of the linear regime (bottom plot of Fig. 8).
3.4 Comparison with the data obtained with the internal lamp
To disentangle the effects of atmospheric perturbations from the aberrations internal to the instrument, we performed the same analysis with an internal light source, as explained in section 2.1. In this case, the only moving part in the instrument was the derotator which simulated pupil tracking of a star observed at an hour angle between 20 and +40 min. We performed the same data analysis as done for the onsky data, namely a calculation of the correlation coefficient in the whole region and in subregions .
The comparison of the decorrelation curves with the onsky data for the whole region is shown in Fig. 10. The more meaningful comparison of the two data sets can be done for sequence that has not been subtracted by the temporal median (right panels). Otherwise, the changes reflect the evolution of the residuals after subtraction of the median, which in the case of the internal data, are relatively low since the sequence is very stable. We however provide also the correlation coefficients in the median subtracted sequence for completeness.
As expected, because the sequence is very stable in time without any perturbation from the atmosphere, the correlation is higher, and stays above 95% during one hour. We can define the same two regimes of decorrelation: a linear regime extending all the way from a few minutes to one hour, and, surprisingly, a steep decorrelation during the first seconds, shown in Fig. 11 . As expected, the linear speed of decorrelation (computed for both internal and onsky data between 5 and 15min) is smaller for the internal data by almost a factor three, indicating that the PSF decorrelates more slowly in the absence of atmospheric perturbations. Surprisingly, the same steep decorrelation occurs on a time scale of a few seconds. The characteristic time is slightly larger for the internal data than for the onsky data (s vs s, right panel of Fig. 11) but this is is less obvious once the temporal median of the sequence has been removed to enhance the relative variation of the quasistatic speckles disregarding the static part of the PSF (s vs s, left panel of Fig. 11).
A zonal analysis of the internal data shows that the characteristic time of the steep exponential decay in the first seconds is only marginally dependent on the separation to the central source, as visible in Fig. 12 left (the black line shows a linear fit only marginally decreasing with the separation, within error bars).
4 Interpretations and additional tests
We will here review some possible explanations for the behaviours presented in section 3.
4.1 A linear decorrelation with time
Between 5 and 20 min, the correlation decreases in a linear way. The atmospheric parameters on such a timescale were stable, we can therefore reasonably rule out external perturbations as the origin of this trend. This regime is also seen with the internal source, although the decorrelation happens at a slower rate. This confirms that the source of this regime is internal to the instrument, and also unrelated to the mechanical flexures of the telescope mirrors. Moreover, we carried out a control experiment by computing the decorrelation in an empty sky region of the detector, and we confirmed the absence of such a linear decrease in this case, detailed in Appendix B.1. We showed that this linear regime evolves faster at short separations, as visible in Fig. 8. This means that the low spatial frequency aberrations (which correspond to regions at small separations in the focal plane) evolve faster than the small scale aberrations in the pupil (corresponding to regions at larger spearations in a focal plane). An additional test carried out with the internal lamp with the AO loop closed but all components fixed in the instrument (including the derotator, that was moving in the measurements presented in section 3.4), shows in this case that the decorrelation speeds are scattered about a zero mean, for the different subregions . This tends to show that this slow linear regime is related to the moving components inside the instrument. A contribution from the thermal expansion of these components might also be involved and could explain why we observe an evolution of small spatial frequencies: the thermal evolution inside the instrument enclosure is of the order of several tens of minutes and induce largescale aberrations on the optics.
If this explanation holds true, this linear regime should evolve slower for an observation carried out about meridian and faster far away after or before meridian. Indeed, for a Nasmith instrument such as VLT / SPHERE, in pupilstabilized mode, the derotator speed cancels at meridian passage and is large far way from meridian. We used therefore our onsky data set to test this hypothesis, but could not come to a clearcut conclusion on this question. The analysis is shown in Appendix D. Our observations were performed after meridian passage, and to study this effect specifically, longer observations about meridian in stable conditions should be repeated. A similar analysis was done with VLT / NaCo around meridian for a longer temporal sequence and could confirm this trend[Milli2016_HIRES].
4.2 A steep decorrelation over a few seconds
A steep regime of decorrelation occurs in the first 30 s, well explained by an exponential decay with a characteristic time of about four seconds. For such a short time scale, ruling out the atmospheric turbulence requires more analysis. Indeed the coherence time was about 3.5 ms in the optical (as estimated from a combination of DIMM measurements and wind speed predictions at 200mbar[Sarazin2002]), which can be converted to a value of 15 ms in the H band. This is too small to explain the steep decorrelation by the presence of atmospherical perturbations uncorrected by the AO system. Macintosh et al. [Macintosh2005] showed that residual speckles caused by noncorrected atmospheric perturbations are refreshed on the pupil every , where is the telescope diameter and the wind speed, inducing a typical decorrelation time of . The wind speed during the observations was 3 to 4 m/s on the ground and higher in altitude, meaning that the resulting time scale is s. This is about the same order of magnitude of what is measured onsky. However, this fast decorrelation regime is also seen with the internal lamp. It is therefore not an effect of the atmospheric perturbations, or the non corrected atmospheric residual phase screen refreshed every second. Moreover, all spatial scales evolve at the same speed because the characteristic times inferred for each individual region do not depend on the separation to the central source. We can also safely rule out any detector persistance effect, as shown in Appendix B.2. In an additional attempt to characterize this effect, we investigated if this can be related to the temporal sampling of the sequence. Frames were recorded every 0.6 s, which is close to the time scale revealed here. We recorded at a later epoch a sequence of frames every 0.18 s and describe the analysis in Appendix B.3. This 4 s timescale is still visible (Fig. 17), despite much different atmospheric conditions, in particular high wind. We can therefore conclude from these analysis that this fast decorrelation over a few seconds is real, related to an internal effect in the intrument independent of the atmospheric conditions or the telescope.
5 Conclusion
We studied the decorrelation timescale of a deep sequence of SPHERE coronagraphic images on a bright star. This sequence is unique because we recorded frames at a high rate of 1.6 Hz by windowing the SPHERE / IRDIS camera in stable atmospheric conditions with a very good AO correction. We used the Pearson correlation coefficients between pairs of frames to evaluate the relative change in the frame speckle pattern during the time ellapsed between the two frames acquisitions. Our analysis highlights two distinct regimes of decorrelation. On the one hand, we observe a linear decorrelation happening for ellapsed times between a few minutes and 20 min. This regime is also observed with an internal lamp, the instrument shutter closed and the derotator simulating pupiltracking, showing that this effect is likely related to the mechanical movements and thermal expansion within the instrument. Additional tests are needed to confirm the origin of this regime, especially a sensitivity analysis with the hour angle, to show that the decorrelation speed is smaller during meridian passage. On the other hand, we also observe a steeper regime of decorrelation taking place during the first few seconds and well fit by an exponential decay. Again this effect is also seen with the internal lamp, shutter closed, showing that this is also an effect internal to the instrument, although the exact cause of this decorrelation is not clear yet.
A Evolution of the correlation coefficient with time in the subregions
Fig. 13 shows a global image of 8 by 12 correlation coefficient curves for the zonal analysis in the subregions . The correlation is higher in the center, close to the coronagraph, as also shown by Fig. 12.
B Verifications
b.1 Analysis of the correlation in an empty sky region
To verify the validity of our measurement procedure, we proceeded to various control experiments. On of them consists in applying the same procedure to a sequence of frames containing only background and thermal noise, without any stellar signal. This was done by using a region of the 128 px wide stripe read on the detector at several arcseconds from the central star, with the same size as the region . The result is shown in Fig. 14. The mean correlation is 0.157 between 5 and 15min and a linear fit to the data yields a zero slope, confirming that the decorrelation seen in the region wellcorrected by the AO system is related to the starlight, its propagation in the atmosphere and in the instrument and its correction by the AO system. This analysis also enables us to get the uncertainty associated to the estimation of the correlation coefficient, measured as the dispersion of the values. It increases towards larger ellapsed time because the number of pairs of frames separated by decreases as increases.
b.2 Remanence of the detector
The SPHERE / IRDIS detector is subject to remanence when a strong illumination is applied. Although the detector was not saturated during the tests presented in this analysis, we verified that the fast decorrelation seen over a few seconds could be due to remanence effects in a region strongly illuminated just outside the coronagraph. To do so, we recorded a sequence of images with a uniform illumination of the detector, e.g. a flat field sequence, and closed the IRDIS shutter during the sequence to study the decay of the detector counts. We windowed the detector to read only 20 rows and be able to have a detector integration time (DIT) of only 25ms. Due to the readout time, we could record a frame every 169ms. For such a short DIT, the illumination was about 1150 ADU, about ten times higher than the level recorded in the baseline test scenario presented here. We computed the mean number of counts in the detector stripe, and show it in Fig. 15 in a window centered on the shutter closure (redshaded area). The count level drops to zero within to readouts, meaning that we can rule out any persistence effect in our analysis for a timescale greater than 340ms.
b.3 Effect of the time sampling
In section 3.1 and 3.3, we show that a fast decorrelation of the images occur in a timescale of about 4s, both for the whole region and in each subregion . The sampling frequency of 1.6 Hz is relatively close to the time scale, and the first data point in Fig. 7 is higher than the best fit value (red curve), therefore we investigated whether a faster image frame rate would lead to a different estimate of . To do so, we windowed the detector to read only 40 raws (equivalent to 0.49on sky), as visible in Fig. 16. This way, we could get a detector integration time (DIT) of 52 ms, and record a frame every 177 ms, equivalent to a frame rate of 5.6 Hz. This sequence was recorded under poor conditions (coherence time of 2.4ms, DIMM seeing 0.7but very high altitude wind and high wind speed on the ground of 12m/s) on the star Crux of magnitude R=1.4, on the night of June, 4 2016. Despite those adverse wind conditions, the Strehl was 80%, but a significant elongation of the PSF is visible along the diagonal of the detector, and induced significant light leak outside the coronagraph radius. We computed the correlation in subregions with the same size as in our baseline sequence. The parameters derived from a fit of an exponential decay (Eq. 2) are shown in Fig. 17. The values of the characteristic time are spread about 4 ms, a value almost identical as that found for the baseline sequence sampled at 1.6 Hz, with some scatter at larger separation where the S/N on the speckles is lower. Therefore we can conclude that the timescale of about 4 s derived for the fast decorrelation regime is robust and does not depend on the frame rate at which the sequence is acquired. Shorter DITs on a brighter star lead to the same estimate of . In addition, this additional observation also shows that does not seem to be related to the quality of the AO correction or the atmospheric conditions, because the Strehl was in this case more than 10% lower than the baseline scenario and the strong wind induced a low coherence time.
C Link with the contrast
The decorrelation speeds might seem small, nevertheless, the contrast which is the relevent quantity that one aims at minimizing does not evolve linearily with the correlation. For a pair of images reduced and centered (e.g subtracted by the spatial mean in the region § and divided by the standard deviation over S), the instantaneous contrast as defined by the norm of the simple difference can be expressed as
(3) 
where is the correlation coefficient as defined in Eq. 1. In such a case, the rate of change of the contrast with time is related to the correlation by
(4) 
This means that for an instantaneous contrast and a decorrelation speed , the contrast decreases at a very high speed of 73 000 ppm/s.
D Influence of the Hour Angle on the Decorrelation
For a telescope with an alt/az mount observing in pupil tracking mode, the quasistatic speckles are expected to be the most stable at meridian passage. At this specific moment, the tracking speed of the altitude axis is indeed going to zero, meaning that the aberrations of the optical parts induced by mechanical flexures of the telescope are minimzed. In addition, the derotator speed, which for a Nasmyth instrument follows the altitude axis, also is moving at its smallest derotation speed, and the ADC movements are minimzed because the target is the highest in the sky. For all these reasons, we tried to detect any influence of the hour angle in our analysis, as it was done for VLT/NaCo data[Milli2016_HIRES]. We therefore binned the onsky data of 3087s into four batches of 13min each and calculated the evolution of the correlation coefficient over time for each batch. The result is shown in Fig. 18. It does not show any significant influence of the hour angle. However, the star was observed after meridian passage, therefore the analysis should be repeated on a sequence centered around meridian to get a higher sensitivity to influence of the hour angle. In addition it should be highlighted that this kind of sensitivity analysis requires a very stable sequence of observations to be able to make a meaningful comparison of the different batches and rule out any effect of the evolution of the atmospheric conditions between the batches.
Acknowledgments
J.M. acknowledges support from the ESO fellowship program and we thank ESO staff and technical operators at the Paranal Observatory.