We describe the implementation of a robotic SLODAR instrument at the Cerro Paranal observatory. The instrument measures the vertical profile of the optical atmospheric turbulence strength, in 8 resolution elements, to a maximum altitude ranging between 100 m and 500 m. We present statistical results of measurements of the turbulence profile on a total of 875 nights between 2014 and 2018. The vertical profile of the ground layer of turbulence is very varied, but in the median case most of the turbulence strength in the ground layer is concentrated within the first 50 m altitude, with relatively weak turbulence at higher altitudes up to 500 m. We find good agreement between measurements of the seeing angle from the SLODAR and from the Paranal DIMM seeing monitor, and also for seeing values extracted from the Shack-Hartmann active optics sensor of VLT UT1, adjusting for the height of each instrument above ground level. The SLODAR data suggest that a median improvement in the seeing angle from 0.689 arcsec to 0.481 arcsec at wavelength 500 nm would be obtained by fully correcting the ground-layer turbulence between the height of the UTs (taken as 10 m) and altitude 500 m.

At Paranal Observatory, the least predictable parameter affecting the short-term scheduling of astronomical observations is the optical turbulence, especially the seeing, coherence time and ground layer fraction. These are critical variables driving the performance of the instruments of the Very Large Telescope (VLT), especially those fed with adaptive optics systems. Currently, the night astronomer does not have a predictive tool to support him/her in decision-making at night. As most service-mode observations at the VLT last less than two hours, it is critical to be able to predict what will happen in this time frame, to avoid time losses due to sudden changes in the turbulence conditions, and also to enable more aggressive scheduling. We therefore investigate here the possibility to forecast the turbulence conditions over the next two hours. We call this "turbulence nowcasting", analogously with weather nowcasting, a term already used in meteorology coming from the contraction of "now" and "forecasting". We present here the results of a study based on historical data of the Paranal Astronomical Site Monitoring combined with ancillary data, in a machine learning framework. We show the strengths and shortcomings of such an approach, and present some perspectives in the context of the Extremely Large Telescope.

The Results of the site selection campaign conducted for the proposed 3.4m optical telescope of the Iranian National Observatory are reported. During the first 3 years, among 33 nominated regions throughout the country, the potential regions were confined to four provinces, namely, Southern Khorasan located in the East, Kerman in the South East and Qom and Esfahan in the central part by examining the long-term meteorological and geographical parameters. Over the following 3 years, after astro-climate considerations and short-term atmospheric seeing measurements using differential image motion monitor (DIMM) technique were carried out, four candidate sites were selected among these regions. Simultaneous seeing measurements were carried out on the four candidate sites and finally, the Dinava located at the common border of the Esfahan and Qom provinces and Gargash in Kashan region were selected as the most promising sites. Continuous seeing measurements during 23 months yielded a median value of about 1.0" for both sites. The latitude and longitude of the Dinava are 50, 54 E and 34, 09 N, respectively and is situated at an altitude of about 3000 m and about 57 Km direct from Qom towards the south. The latitude and longitude of the Gargash are 51, 19 E and 33, 40 N, respectively and is situated at an altitude of about 3600 m and about 37km direct from Kashan towards the south. With two sites in hand, further comparative measurements of seeing and weather data continued at the two sites, which resulted in the selection of the Gargash site for the construction of the Iranian National Observatory.

Atmospheric water vapor is the main limiting factor of atmospheric transparency in the mm and submm wavelength spectral windows. Thus, dry sites are needed for the installation and successful operation of radio astronomy observatories exploiting those spectral windows. Temperature and wind are variables of special consideration when planning the installation and operation of large-aperture radio telescopes, as these parameters affect the mechanical response of radio telescopes exposed to the environmental conditions. Temperature, and in particular temperature gradients, induce thermal deformation of mechanical structures, while high wind speeds and gusts induce pointing jitter affecting the tracking accuracy of astronomical sources during the observations. This work summarizes the statistics of precipitable water vapor (PWV), temperature, and wind monitored at sites by the coastal mountain range, as well as on the west slope of the Andes mountain range in the region of Antofagasta, Chile. This information could prove useful for the planning of extended baselines for the Atacama Large Millimeter and Submillimeter Array (ALMA), and/or new radio telescope projects, such as the Atacama Large Aperture Submm/ mm Telescope (AtLAST) initiative.

In addition to astro-meteorological parameters, such as seeing, coherence time and isoplanatic angle, the vertical profile of the Earth's atmospheric turbulence strength and velocity is important for instrument design, performance validation and monitoring, and observation scheduling and management. Here we compare these astro-meteorological parameters as well as the vertical profile itself from a forecast model based on a General Circulation Model from the European Centre for Median range Weather Forecasts and the stereo-SCIDAR, a high-sensitivity turbulence profiling instrument in regular operation at Paranal, Chile. The model is fast to process as no spatial nesting or data manipulation is performed. This speed enables the model to be reactive based on the most up to date forecasts. We find that the model is statistically consistent with measurements from stereo-SCIDAR. The correlation of the median turbulence profile from the model and the measurement is 0.98. We also find that the distributions of astro-meteorological parameters are consistent. We compare contemporaneous measurements and show that the free atmosphere seeing, isoplanatic angle and coherence time have correlation values of 0.64, 0.40 and 0.63 respectively. We show and compare the profile sequences from a large number of trial nights. We see that the model is able to forecast the evolution of dominating features. In addition to smart scheduling, ensuring that the most sensitive astronomical observations are scheduled for the optimum time, this model could enable remote site characterisation using a large archive of weather forecasts and could be used to optimise the performance of wide-field AO system.

Knowledge of the optical turbulence profile is important in adaptive optics (AO) systems, particularly tomographic AO systems such as those to be employed by the next generation of 40 m class extremely large telescopes (ELTs). Site characterisation and monitoring campaigns have produced large quantities of turbulence profiling data for sites around the world. However AO system design and performance characterisation is dependent on Monte-Carlo simulations that cannot make use of these large datasets due to long computation times. Here we address the question of how to reduce these large datasets into small sets of profiles that can feasibly be used in such Monte-Carlo simulations, whilst minimising the loss of information inherent in this effective compression of the data. We propose hierarchical clustering to partition the dataset according to the structure of the turbulence profiles and extract a single profile from each cluster. This method is applied to the Stereo-SCIDAR dataset from ESO Paranal containing over 10000 measurements of the turbulence profile from 83 nights. We present two methods of extracting turbulence profiles from the clusters, resulting in two sets of 18 profiles providing subtly different descriptions of the variability across the entire dataset. For generality we choose integrated parameters of the turbulence to measure the representativeness of our profiles and compare to others. Using these criterion we also show that such variability is difficult to capture with small sets of profiles associated with integrated turbulence parameters such as seeing.

SPHERE is the high-contrast exoplanet imager and spectrograph installed at the Unit Telescope 3 of the Very Large Telescope. After more than two years of regular operations, we analyse statistically the performance of the adaptive optics system and its dependence on the atmospheric conditions above the Paranal observatory, as measured by the suite of dedicated instruments which are part of the Astronomical Site Monitor and as estimated by the SPHERE real-time calculator. We also explain how this information can be used to schedule the observations in order to yield the best data quality and to guide the astronomer when processing his/her data.

Knowledge of the Earth's atmospheric optical turbulence is critical for astronomical instrumentation. Not only does it enable performance verification and optimisation of existing systems but it is required for the design of future instruments. As a minimum this includes integrated astro-atmospheric parameters such as seeing, coherence time and isoplanatic angle, but for more sophisticated systems such as wide field adaptive optics enabled instrumentation the vertical structure of the turbulence is also required. Stereo-SCIDAR is a technique specifically designed to characterise the Earth's atmospheric turbulence with high altitude resolution and high sensitivity. Together with ESO, Durham University has commissioned a Stereo-SCIDAR instrument at Cerro Paranal, Chile, the site of the Very Large Telescope (VLT), and only 20~km from the site of the future Extremely Large Telescope (ELT). Here we provide results from the first 18 months of operation at ESO Paranal including instrument comparisons and atmospheric statistics. Based on a sample of 83 nights spread over 22 months covering all seasons, we find the median seeing to be 0.64" with 50% of the turbulence confined to an altitude below 2 km and 40% below 600 m. The median coherence time and isoplanatic angle are found as 4.18 ms and 1.75" respectively. A substantial campaign of inter-instrument comparison was also undertaken to assure the validity of the data. The Stereo-SCIDAR profiles (optical turbulence strength and velocity as a function of altitude) have been compared with the Surface-Layer SLODAR, MASS-DIMM and the ECMWF weather forecast model. The correlation coefficients are between 0.61 (isoplanatic angle) and 0.84 (seeing).

Scintillation noise is a major limitation of ground base photometric precision. An extensive dataset of stellar scintillation collected at 11 astronomical sites world-wide with MASS instruments was used to estimate the scintillation noise of large telescopes in the case of fast photometry and traditional long-exposure regime. Statistical distributions of the corresponding parameters are given. The scintillation noise is mostly determined by turbulence and wind in the upper atmosphere and comparable at all sites, with slightly smaller values at Mauna Kea and largest noise at Tolonchar in Chile. We show that the classical Young's formula under-estimates the scintillation noise.The temporal variations of the scintillation noise are also similar at all sites, showing short-term variability at time scales of 1 -- 2 hours and slower variations, including marked seasonal trends (stronger scintillation and less clear sky during local winter). Some correlation was found between nearby observatories.

This is the second article of a series devoted to European Extremely Large Telescope (E-ELT) site characterization. In this article we present the main properties of the parameters involved in high angular resolution observations from the data collected in the site testing campaign of the E-ELT during the Design Study (DS) phase. Observations were made in 2008 and 2009, in the four sites selected to shelter the future E-ELT (characterized under the ELT-DS contract): Aklim mountain in Morocco, Observatorio del Roque de los Muchachos (ORM) in Spain, Mac\'on range in Argentina, and Cerro Ventarrones in Chile. The same techniques, instruments and acquisition procedures were taken on each site. A Multiple Aperture Scintillation Sensor (MASS) and a Differential Image Motion Monitor (DIMM) were installed at each site. Global statistics of the integrated seeing, the free atmosphere seeing, the boundary layer seeing and the isoplanatic angle were studied for each site, and the results are presented here. In order to estimate other important parameters such as the coherence time of the wavefront and the overall parameter "coherence \'etendue" additional information of vertical profiles of the wind speed was needed. Data were retrieved from the National Oceanic and Atmospheric Administration (NOAA) archive. Ground wind speed was measured by Automatic Weather Stations (AWS). More aspects of the turbulence parameters such as their seasonal trend, their nightly evolution and their temporal stability were also obtained and analyzed.

Real-time seeing estimation at the focus of a telescope is nowadays strongly emphasized as this knowledge virtually drives the dimensioning of adaptive optics systems and instrument operational aspects. In this context we study the interest of using active optics Shack-Hartmann (AOSH) sensor images to provide accurate estimate of the seeing. The AOSH practically delivers long exposure spot PSFs -- at the critical location of the telescope focus -- being directly related to the atmospheric seeing in the line of sight. Although AOSH sensors are not specified to measure spot sizes but slopes, we show that accurate seeing estimation from AOSH images can be obtained with a dedicated algorithm. The sensitivity and comparison of two algorithms to various parameters is analyzed in a systematic way, demonstrating that efficient estimation of the seeing can be obtained by adequate means.

The precision of radial velocity (RV) measurements depends on the precision attained on the wavelength calibration. One of the available options is using atmospheric lines as a natural, freely available wavelength reference. Figueira et al. (2010) measured the RV of O2 lines using HARPS and showed that the scatter was only of ~10 m/s over a timescale of 6 yr. Using a simple but physically motivated empirical model, they demonstrated a precision of 2 m/s, roughly twice the average photon noise contribution. In this paper we take advantage of a unique opportunity to confirm the sensitivity of the telluric absorption lines RV to different atmospheric and observing conditions: by means of contemporaneous in-situ wind measurements by radiosondes. The RV model fitting yielded similar results to that of Figueira et al. (2010), with lower wind magnitude values and varied wind direction. The probes confirmed the average low wind magnitude and suggested that the average wind direction is a function of time as well. The two approaches deliver the same results in what concerns wind magnitude and agree on wind direction when fitting is done in segments of a couple of hours. Statistical tests show that the model provides a good description of the data on all timescales, being always preferable to not fitting any atmospheric variation. The smaller the timescale on which the fitting can be performed (down to a couple of hours), the better the description of the real physical parameters. We conclude then that the two methods deliver compatible results, down to better than 5 m/s and less than twice the estimated photon noise contribution on O2 lines RV measurement. However, we cannot rule out that parameters alpha and gamma (dependence on airmass and zero-point, respectively) have a dependence on time or exhibit some cross-talk with other parameters (abridged).

Seeing-limited resolution in large telescopes working over wide wavelength range depends substantially on the turbulence outer scale and cannot be adequately described by one "seeing" value. We attempt to clarify frequent confusions on this matter. We study the effects of finite turbulence outer scale and partial adaptive corrections by means of analytical calculations and numerical simulations. If a von Karman turbulence model is adopted, a simple approximate formula captures the dependence of atmospheric long-exposure resolution on the outer scale over the entire practically interesting range of telescope diameters and wavelengths. In the infrared (IR), the difference with the standard Kolmogorov seeing formula can exceed a factor of two. We find that low-order adaptive turbulence correction produces residual wave-fronts with effectively small outer scale, so even very low compensation order leads to a substantial improvement in resolution over seeing, compared to the standard theory. Seeing-limited resolution of large telescopes, especially in the IR, is currently under-estimated by not accounting for the outer scale. On the other hand, adaptive-optics systems designed for diffraction-limited imaging in the IR can improve the resolution in the visible by as much as two times.

We probe turbulence structure from centimetric to metric scales by simultaneous imagery at mid-infrared and visible wavelengths at the VLT telescope and show that it departs significantly from the commonly used Kolmogorov model. The data can be fitted by the von Karman model with an outer scale of the order of 30 m and we see clear signs of the phase structure function saturation across the 8-m VLT aperture. The image quality improves in the infrared faster than the standard lambda^{-1/5} scaling and may be diffraction-limited at 30-m apertures even without adaptive optics at wavelengths longer than 8 micron.

Temporal fluctuations of the atmospheric piston are critical for interferometers as they determine their sensitivity. We characterize an instrumental set-up, termed the piston scope, that aims at measuring the atmospheric time constant, tau0, through the image motion in the focal plane of a Fizeau interferometer. High-resolution piston scope measurements have been obtained at the observatory of Paranal, Chile, in April 2006. The derived atmospheric parameters are shown to be consistent with data from the astronomical site monitor, provided that the atmospheric turbulence is displaced along a single direction. Piston scope measurements, of lower temporal and spatial resolution, were for the first time recorded in February 2005 at the Antarctic site of DomeC. Their re-analysis in terms of the new data calibration sharpens the conclusions of a first qualitative examination.

We analyze the monthly average wind velocity at about 12 km above sea level, between 1980 to 1995, for San Pedro Martir, Mauna Kea, another existing observatorie and some sites of interest. We compare the results obtained from two different data sets, the GGUAS and NCEP. Our results show that San Pedro Martir and Mauna Kea are comparable and are amongst the most suitable sites to apply slow wavefront corrugation correction techniques

The region surrounding the Llano de Chajnantor, a high altitude plateau in the Atacama Desert in northern Chile, has caught the attention of the astronomical community for its potential as an observatory site. Combining high elevation and extremely low atmospheric water content, the Llano has been chosen as the future site of the Atacama Large Millimeter Array. We have initiated a campaign to investigate the astronomical potential of the region in the optical/infrared. Here, we report on an aspect of our campaign aimed at establishing a seeing benchmark to be used as a reference for future activities in the region. After a brief description of the region and its climate, we describe the results of an astronomical seeing campaign, carried out with a Differential Image Motion Monitor that operates at 0.5 micron wavelength. The seeing at the Llano level of 5000 m, measured over 7 nights in May 1998, yielded a median FWHM of 1.1". However, the seeing decreased to 0.7" at a modest 100 m gain above the plateau (Cerro Chico), as measured over 38 nights spread between July 1998 and October 2000. Neither of these represents the best seeing expected in the region; the set of measurements provides a reference base for simultaneous dual runs at Cerro Chico and at other sites of interest in the region, currently underway. A comparison between simultaneous measurements at Cerro Chico and Cerro Paranal indicates that seeing at Cerro Chico is about 12% better than at Paranal. The percentage of optically photometric nights in the Chajnantor region is about 60%, while that of nights useful for astronomical work is near 80%.

The performance of an adaptive optics (AO) system on a 100m diameter ground based telescope working in the visible range of the spectrum is computed using an analytical approach. The target Strehl ratio of 60% is achieved at 0.5um with a limiting magnitude of the AO guide source near R~10, at the cost of an extremely low sky coverage. To alleviate this problem, the concept of tomographic wavefront sensing in a wider field of view using either natural guide stars (NGS) or laser guide stars (LGS) is investigated. These methods use 3 or 4 reference sources and up to 3 deformable mirrors, which increase up to 8-fold the corrected field size (up to 60\arcsec at 0.5 um). Operation with multiple NGS is limited to the infrared (in the J band this approach yields a sky coverage of 50% with a Strehl ratio of 0.2). The option of open-loop wavefront correction in the visible using several bright NGS is discussed. The LGS approach involves the use of a faint (R ~22) NGS for low-order correction, which results in a sky coverage of 40% at the Galactic poles in the visible.