In this contribution, we present the most recent progresses we obtained in the context of a long-term program we undertook since a few years towards the implementation of operational forecast systems (a) on top-class ground-based telescopes assisted by AO systems to support the flexible scheduling of observational scientific programs in night as well in day time and (b) on ground-stations to support free space optical communication. Two topics have been treated and presented in the Conference AO4ELT6: 1. ALTA is an operational forecast system for the OT and all the critical atmospheric parameters affecting the astronomical ground-based observations conceived for the LBT. It operates since 2016 and it is in continuous evolution to match with necessities/requirements of instruments assisted by AO of the LBT (SOUL, SHARK-NIR, SHARK-VIS, LINC-NIRVANA,...). In this contribution, we present a new implemented version of ALTA that, thanks to an auto-regression method making use of numerical forecasts and real-time OT measurements taken in situ, can obtain model performances (for forecasts of atmospherical and astroclimatic parameters) never achieved before on time scales of the order of a few hours. 2. We will go through the main differences between optical turbulence forecast performed with mesoscale and general circulation models (GCM) by clarifying some fundamental concepts and by correcting some erroneous information circulating recently in the literature.

The efficiency of the management of top-class ground-based astronomical facilities supported by Adaptive Optics (AO) relies on our ability to forecast the optical turbulence (OT) and a set of relevant atmospheric parameters. Indeed, in spite of the fact that the AO is able to achieve, at present, excellent levels of wavefront corrections (a Strehl Ratio up to 90% in H band), its performances strongly depend on the atmospheric conditions. Knowing in advance the turbulence conditions allows an optimization of the AO use. It has already been proven that it is possible to provide reliable forecasts of the optical turbulence (CN2 profiles and integrated astroclimatic parameters such as seeing, isoplanantic angle, wavefront coherence time, ...) for the next night. In this paper we prove that it is possible to improve the forecast performances on shorter time scales (order of one or two hours) with consistent gains (order of 2 to 8) using filtering techniques. This has permitted us to achieve forecasts accuracies never obtained before and reach a fundamental milestone for the astronomical applications. The time scale of one or two hours is the most critical one for an efficient management of the ground-based telescopes supported by AO. Results shown here open, therefore, to an important revolution in the field. We implemented this method in the operational forecast system of the Large Binocular Telescope, named ALTA Center that is, at our knowledge, the first operational system providing forecasts of turbulence and atmospheric parameters at short time scales to support science operations.

In this contribution we evaluate the impact of filtering techniques in enhancing the accuracy of forecasts of optical turbulence and atmospheric parameters critical for ground-based telescopes. These techniques make use of the data continuously provided by the telescope sensors and instruments to improve the performances of real-time forecasts which have an impact on the telescope operation. In previous works we have already shown how a mesoscale high-frequency forecast (Meso-NH and Astro-Meso-Nh models) can produce reliable predictions of different atmospheric parameters and the optical turbulence. The mesoscale forecast has an advantage on the global model in having a better implementation of the physical atmospheric processes, including turbulence, and produces an output with greater spatial resolution (up to 100m or beyond). Filtering techniques that make use of the real-time sensor data at the telescope may help in removing potential biases and trends which have an impact on short term mesoscale forecast and, as a consequence, may increase the accuracy of the final output. Given the complexity and cost of present and future top-class telescope installations, each improvement of forecasts of future observing conditions will definitely help in better allocating observing time, especially in queue-mode operation, and will definitely benefit the scientific community in medium-long term.

(35-words maximum) In this talk I present the scientific drivers related to the optical turbulence forecast applied to the ground-based astronomy supported by Adaptive Optics, the state of the art of the achieved results and the most relevant challenges for future progresses.

Forecast of the atmospheric parameters and optical turbulence applied to the ground-based astronomy is very crucial mainly for the queue scheduling. So far, most efforts have been addressed by our group in developing algorithms for the optical turbulence (CN2) and annexed integrated astroclimatic parameters and quantifying the performances of the Astro-Meso-Nh package in reconstructing such parameters. Besides, intensive analyses on the Meso-Nh performances= in reconstructing atmospheric parameters relevant for the ground-based astronomy has been carried out. Our studies referred always to the night time regime. To extend the applications of our studies to the day time regime, we present, in this contribution, preliminary results obtained by comparing model outputs and measurements of classical atmospheric parameter relevant for the ground-based astronomy in night and day time. We chose as a test case, the Roque de los Muchachos Observatory (Canary Islands), that offers a very extended set of measurements provided by different sensors belonging to different telescopes on the same summit/Observatory. The convective regime close to the ground typical of the day time is pretty different from the stable regime characterising the night time. This study aims therefore to enlarge the domain of validity of the Astro-Meso-Nh code to new turbulence regimes and it permits to cover the total 24 hours of a day. Such an approach will permit not only an application to solar telescopes (e.g. EST) but also applications to a much extended set of scientific fields, not only in astronomical context such as satellite communications.

Water vapour in the atmosphere is the main source of the atmospheric opacity in the infrared and sub-millimetric regimes and its value plays a critical role in observations done with instruments working at these wavelengths on ground-based telescopes. The scheduling of scientific observational programs with instruments such as the VLT Imager and Spectrometer for mid Infrared (VISIR) at Cerro Paranal and the Large Binocular Telescope Interferometer (LBTI) at Mount Graham would definitely benefit from the ability to forecast the atmospheric water vapour content. In this contribution we present a study aiming at validating the performance of the non-hydrostatic mesoscale Meso-NH model in reliably predicting precipitable water vapour (PWV) above the two sites. For the VLT case we use, as a reference, measurements done with a Low Humidity and Temperature PROfiling radiometer (LHATPRO) that, since a few years, is operating routinely at the VLT. LHATPRO has been extensively validated on previous studies. We obtain excellent performances on forecasts performed with this model, including for the extremely low values of the PWV (<= 1 mm). For the LBTI case we compare one solar year predictions obtained with the Meso-NH model with satellite estimates again obtaining an excellent agreement. This study represents a further step in validating outputs of atmospheric parameters forecasts from the ALTA Center, an operational and automatic forecast system conceived to support observations at LBT and LBTI.

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).

Wide Field Adaptive Optics (WFAO) systems are among the most sophisticated AO systems available today on large telescopes. The knowledge of the vertical spatio-temporal distribution of the wind speed (WS) and direction (WD) are fundamental to optimize the performance of such systems. Previous studies already proved that the Gemini Multi-Conjugated AO system (GeMS) is able to retrieve measurements of the WS and WD stratification using the SLODAR technique and to store measurements in the telemetry data. In order to assess the reliability of these estimates and of the SLODAR technique applied to such a kind of complex AO systems, in this study we compared WS and WD retrieved from GeMS with those obtained with the atmospherical model Meso-Nh on a rich statistical sample of nights. It has been previously proved that, the latter technique, provided an excellent agreement with a large sample of radiosoundings both, in statistical terms and on individual flights. It can be considered, therefore, as an independent reference. The excellent agreement between GeMS measurements and the model that we find in this study, proves the robustness of the SLODAR approach. To by-pass the complex procedures necessary to achieve automatic measurements of the wind with GeMS, we propose a simple automatic method to monitor nightly WS and WD using the Meso-Nh model estimates. Such a method can be applied to whatever present or new generation facilities supported by WFAO systems. The interest of this study is, therefore, well beyond the optimization of GeMS performance.

In this paper we quantify the performances of an automated weather forecast system implemented on the Large Binocular Telescope (LBT) site at Mt. Graham (Arizona) in forecasting the main atmospheric parameters close to the ground. The system employs a mesoscale non-hydrostatic numerical model (Meso-Nh). To validate the model we compare the forecasts of wind speed, wind direction, temperature and relative humidity close to the ground with the respective values measured by instrumentation installed on the telescope dome. The study is performed over a large sample of nights uniformly distributed over two years. The quantitative analysis is done using classical statistical operators (bias, RMSE and $\sigma$) and contingency tables, which allows to extract complementary key information, such as the percentage of correct detection (PC) and the probability to obtain a correct detection within a defined interval of values (POD). Results of our study indicate that the model performances in forecasting the atmospheric parameters we have just cited are very good, in some cases excellent: RMSE for temperature is below 1{\deg} C, for relative humidity is 14%, for the wind speed is around 2.5m/s. The relative error of the RMSE for wind direction varies from 9% to 17% depending on the wind speed conditions. This work is performed in the context of ALTA (Advanced LBT Turbulence and Atmosphere) Center project, which final goal is to provide forecasts of all the atmospheric parameters and the optical turbulence to support LBT observations, adaptive optics facilities and interferometric facilities.

One of the main goals of the feasibility study MOSE (MOdellig ESO Sites) is to evaluate the performances of a method conceived to forecast the optical turbulence above the ESO sites of the Very Large Telescope and the European-Extremely Large Telescope in Chile. The method implied the use of a dedicated code conceived for the optical turbulence (OT) called Astro-Meso-Nh. In this paper we present results we obtained at conclusion of this project concerning the performances of this method in forecasting the most relevant parameters related to the optical turbulence (CN2, seeing , isoplanatic angle theta_0 and wavefront coherence time tau_0). Numerical predictions related to a very rich statistical sample of nights uniformly distributed along a solar year and belonging to different years have been compared to observations and different statistical operators have been analyzed such as classical bias, RMSE and and more sophisticated statistical operators derived by the contingency tables that are able to quantify the score of success of a predictive method such as the percentage of correct detection (PC) and the probability to detect a parameter within a specific range of values (POD). The main conclusions of the study tell us that the Astro-Meso-Nh model provides performances that are already very good to definitely guarantee a not negligible positive impact on the Service Mode of top-class telescopes and ELTs. A demonstrator for an automatic and operational version of the Astro-Meso-Nh model will be soon implemented on the sites of VLT and E-ELT.

In this paper we study the abilities of an atmospherical mesoscale model in forecasting the classical atmospherical parameters relevant for astronomical applications at the surface layer (wind speed, wind direction, temperature, relative humidity) on the Large Binocular Telescope (LBT) site - Mount Graham, Arizona. The study is carried out in the framework of the ALTA project aiming at implementing an automated system for the forecasts of atmospherical parameters (Meso-Nh code) and the optical turbulence (Astro-Meso-Nh code) for the service-mode operation of the LBT. The final goal of such an operational tool is to provide predictions with high time frequency of atmospheric and optical parameters for an optimized planning of the telescope operation (dome thermalization, wind-dependent dome orientation, observation planning based on predicted seeing, adaptive optics optimization, etc...). Numerical simulations are carried out with the Meso-Nh and Astro-Meso-Nh codes, which were proven to give excellent results in previous studies focused on the two ESO sites of Cerro Paranal and Cerro Armazones (MOSE Project). In this paper we will focus our attention on the comparison of atmospherical parameters forescasted by the model close to the ground with measurements taken by the observatory instrumentations and stored in the LBT telemetry in order to validate the numerical predictions. As previously done for Cerro Paranal (Lascaux et al., 2015), we will also present an analysis of the model performances based on the method of the contingency tables, that allows us to provide complementary key information with the respect to the bias and RMSE (systematic and statistical errors), such as the percentage of correct detection and the probability to obtain a correct detection inside a defined interval of values.

Wide Field Adaptive Optics (WFAO) systems represent the more sophisticated AO systems available today at large telescopes. A critical aspect for these WFAO systems in order to deliver an optimised performance is the knowledge of the vertical spatiotemporal distribution of the CN2 and the wind speed. Previous studies (Cortes et al., 2012) already proved the ability of GeMS (the Gemini Multi-Conjugated AO system) in retrieving CN2 and wind vertical stratification using the telemetry data. To assess the reliability of the GeMS wind speed estimates a preliminary study (Neichel et al., 2014) compared wind speed retrieved from GeMS with that obtained with the atmospherical model Meso-Nh on a small sample of nights providing promising results. The latter technique is very reliable for the wind speed vertical stratification. The model outputs gave, indeed, an excellent agreement with a large sample of radiosoundings (~ 50) both in statistical terms and on individual flights (Masciadri et al., 2013). Such a tool can therefore be used as a valuable reference in this exercise of cross calibrating GeMS on-sky wind estimates with model predictions. In this contribution we achieved a two-fold results: (1) we extended analysis on a much richer statistical sample (~ 43 nights), we confirmed the preliminary results and we found an even better correlation between GeMS observations and the atmospherical model with basically no cases of not-negligible uncertainties; (2) we evaluate the possibility to use, as an input for GeMS, the Meso-Nh estimates of the wind speed stratification in an operational configuration. Under this configuration these estimates can be provided many hours in advanced with respect to the observations and with a very high temporal frequency (order of 2 minutes or less).

In this contribution we present the most relevant results obtained in the context of a feasibility study (MOSE) undertaken for ESO. The principal aim of the project was to quantify the performances of a mesoscale model (Astro-Meso-NH code) in forecasting all the main atmospherical parameters relevant for the ground-based astronomical observations and the optical turbulence (CN2 and associated integrated astroclimatic parameters) above Cerro Paranal (site of the VLT) and Cerro Armazones (site of the E-ELT). A detailed analysis on the score of success of the predictive capacities of the system have been carried out for all the astroclimatic as well as for the atmospherical parameters. Considering the excellent results that we obtained, this study proved the opportunity to implement on these two sites an automatic system to be run nightly in an operational configuration to support the scheduling of scientific programs as well as of astronomical facilities (particularly those supported by AO systems) of the VLT and the E-ELT. At the end of 2016 a new project for the implementation of a demonstrator of an operational system to be run on the two ESO's sites will start. Our team is also responsible for the implementation of a similar automatic system at Mt.Graham, site of the LBT (ALTA Project). Our system/method will permit therefore to make a step ahead in the framework of the Service Mode for new generation telescopes. Among the most exciting achieved results we cite the fact that we proved to be able to forecast CN2 profiles with a vertical resolution as high as 150 m. Such a feature is particularly crucial for all WFAO systems that require such detailed information on the OT vertical stratification on the whole 20 km above the ground. This important achievement tells us that all the WFAO systems can rely on automatic systems that are able to support their optimized use.

This article aims at proving the feasibility of the forecast of all the most relevant classical atmospherical parameters for astronomical applications (wind speed and direction, temperature) above the ESO ground-base site of Cerro Paranal with a mesoscale atmospherical model called Meso-Nh. In a precedent paper we have preliminarily treated the model performances obtained in reconstructing some key atmospherical parameters in the surface layer 0-30~m studying the bias and the RMSE on a statistical sample of 20 nights. Results were very encouraging and it appeared therefore mandatory to confirm such a good result on a much richer statistical sample. In this paper, the study was extended to a total sample of 129 nights between 2007 and 2011 distributed in different parts of the solar year. This large sample made our analysis more robust and definitive in terms of the model performances and permitted us to confirm the excellent performances of the model. Besides, we present an independent analysis of the model performances using the method of the contingency tables. Such a method permitted us to provide complementary key informations with respect to the bias and the RMSE particularly useful for an operational implementation of a forecast system.

The advent of a new generation of Adaptive Optics systems called Wide Field AO (WFAO) mark the beginning of a new era. By using multiple Guide Stars (GSs), either Laser Guide Stars (LGSs) or Natural Guide Stars (NGSs), WFAO significantly increases the field of view of the AO-corrected images, and the fraction of the sky that can benefit from such correction. Different typologies of WFAO have been studied over the past years. They all require multiple GSs to perform a tomographic analysis of the atmospheric turbulence. One of the fundamental aspects of the new WFAO systems is the knowledge of the spatio-temporal distribution of the turbulence above the telescope. One way to get to this information is to use the telemetry data provided by the WFAO system itself. Indeed, it has been demonstrated that WFAO systems allows one to derive the Cn2 and wind profile in the main turbulence layers (see e.g. Cortes et al. 2012). This method has the evident advantage to provide information on the turbulence stratification that effectively affects the AO system, property more difficultly respected by independently vertical profilers. In this paper, we compare the wind speeds profiles of GeMS with those predicted by a non-hydrostatical mesoscale atmospherical model (Meso-NH). It has been proved (Masciadri et al., 2013), indeed, that this model is able to provide reliable wind speed profiles on the whole troposphere and stratosphere (up to 20-25 km) above top-level astronomical sites. Correlation with measurements revealed to be very satisfactory when the model performances are analyzed from a statistical point of view as well on individual nights. Such a system appears therefore as an interesting reference to be used to quantify the GeMS wind speed profiles reliability.

The Multi Aperture Scintillation Sensor (MASS) and the Generalized-Scintillation Detection and Ranging (Generalized SCIDAR) are two instruments conceived to measure the optical turbulence (OT) vertical distribution on the whole troposphere and low stratosphere (~ 20 km) widely used in the astronomical context. In this paper we perform a detailed analysis/comparison of measurements provided by the two instruments and taken during the extended site testing campaign carried out on 2007 at Cerro Paranal and promoted by the European Southern Observatory (ESO). The main and final goal of the study is to provide a detailed estimation of the measurements reliability i.e dispersion of turbulence measurements done by the two instruments at different heights above the ground. This information is directly related to our ability in estimating the absolute value of the turbulence stratification. To better analyse the uncertainties between the MASS and the GS we took advantage of the availability of measurements taken during the same campaign by a third independent instrument (DIMM - Differential Imaging Motion Monitor) measuring the integrated turbulence extended on the whole 20 km. Such a cross-check comparison permitted us to define the reliability of the instruments and their measurements, their limits and the contexts in which their use can present some risk.

We present the overview of the MOSE project (MOdeling ESO Sites) aiming at proving the feasibility of the forecast of the classical atmospherical parameters (wind speed intensity and direction, temperature, relative humidity) and the optical turbulence OT (CN2 profiles and the most relevant integrated astro-climatic parameters derived from the CN2: the seeing, the isoplanatic angle, the wavefront coherence time) above the two ESO ground-based sites of Cerro Paranal and Cerro Armazones. The final outcome of the study is to investigate the opportunity to implement an automatic system for the forecast of these parameters at these sites. In this paper we present results related to the Meso-Nh model ability in reconstructing the vertical stratification of the atmospherical parameters along the 20 km above the ground. The very satisfactory performances shown by the model in reconstructing most of these parameters (and in particular the wind speed) put this tool of investigation as the most suitable to be used in astronomical observatories to support AO facilities and to calculate the temporal evolution of the wind speed and the wavefront coherence time at whatever temporal sampling. The further great advantage of this solution is that such estimates can be available in advance (order of some hours) with respect to the time of interest

This article is the second of a series of articles aiming at proving the feasibility of the forecast of all the most relevant classical atmospherical parameters for astronomical applications (wind speed and direction, temperature, relative humidity) and the optical turbulence (Cn2 and the derived astro-climatic parameters like seeing, isoplanatic angle, wavefront coherence time...). This study is done in the framework of the MOSE project, and focused above the two ESO ground-bases sites of Cerro Paranal and Cerro Armazones. In this paper we present the results related to the Meso-Nh model ability in reconstructing the surface layer atmospherical parameters (wind speed intensity, wind direction and absolute temperature, [0-30] m a.g.l.). The model reconstruction of all the atmospherical parameters in the surface layer is very satisfactory. For the temperature, at all levels, the RMSE (Root Mean Square Error) is inferior to 1{\deg}C. For the wind speed, it is ~2 m/s, and for the wind direction, it is in the range [38-46{\deg}], at all levels, that corresponds to a RMSE_relative in a range [21-26{\deg}]. If a filter is applied for the wind direction (the winds inferior to 3 m/s are discarded from the computations), the wind direction RMSE is in the range [30-41{\deg}], i.e. a RMSE_relative in the range [17-23{\deg}]. The model operational forecast of the surface layer atmospherical parameters is suitable for different applications, among others: thermalization of the dome using the reconstructed temperature, hours in advance, of the beginning the night; knowing in advance the main direction which the strong winds will come from during the night could allow the astronomer to anticipate the occurrence of a good/bad seeing night, and plan the observations accordingly; preventing adaptive secondary mirrors shake generated by the wind speed.

We present very encouraging preliminary results obtained in the context of the MOSE project, an on-going study aiming at investigating the feasibility of the forecast of the optical turbulence and meteorological parameters (in the free atmosphere as well as in the boundary and surface layer) at Cerro Paranal (site of the Very Large Telescope - VLT) and Cerro Armazones (site of the European Extremely Large Telescope - E-ELT), both in Chile. The study employs the Meso-Nh atmospheric mesoscale model and aims at supplying a tool for optical turbulence forecasts to support the scheduling of the scientific programs and the use of AO facilities at the VLT and the E-ELT. In this study we take advantage of the huge amount of measurements performed so far at Paranal and Armazones by ESO and the TMT consortium in the context of the site selection for the E-ELT and the TMT to constraint/validate the model. A detailed analysis of the model performances in reproducing the atmospheric parameters (T, V, p, H, ...) near the ground as well as in the free atmosphere, is critical and fundamental because the optical turbulence depends on most of these parameters. This approach permits us to provide an exhaustive and complete analysis of the model performances and to better define the model operational application. This also helps us to identify the sources of discrepancies with optical turbulence measurements (when they appear) and to discriminate between different origins of the problem: model parameterization, initial conditions, ... Preliminary results indicate a great accuracy of the model in reproducing most of the main meteorological parameters in statistical terms as well as in each individual night in the free atmosphere and in proximity of the surface. The study is co-funded by ESO and INAF-Arcetri (Italy).

In the context of the MOSE project, in this contribution we present a detailed analysis of the Meso-NH mesoscale model performances and their dependency on the model and orography horizontal resolutions in proximity of the ground. The investigated sites are Cerro Paranal (site of the ESO Very Large Telescope - VLT) and Cerro Armazones (site of the ESO European Extremely Large Telescope - E-ELT), in Chile. At both sites, data from a rich statistical sample of different nights are available - from AWS (Automated Weather Stations) and masts - giving access to wind speed, wind direction and temperature at different levels near the ground (from 2 m to 30 m above the ground). In this study we discuss the use of a very high horizontal resolution (dX=0.1 km) numerical configuration that overcomes some specific limitations put in evidence with a standard configuration with dX=0.5 km. In both sites results are very promising. The study is co-funded by ESO and INAF.

We have performed H and Ks band observations of the planetary system around HR 8799 using the new AO system at the Large Binocular Telescope and the PISCES Camera. The excellent instrument performance (Strehl ratios up to 80% in H band) enabled detection the inner planet HR8799e in the H band for the first time. The H and Ks magnitudes of HR8799e are similar to those of planets c and d, with planet e slightly brighter. Therefore, HR8799e is likely slightly more massive than c and d. We also explored possible orbital configurations and their orbital stability. We confirm that the orbits of planets b, c and e are consistent with being circular and coplanar; planet d should have either an orbital eccentricity of about 0.1 or be non-coplanar with respect to b and c. Planet e can not be in circular and coplanar orbit in a 4:2:1 mean motion resonances with c and d, while coplanar and circular orbits are allowed for a 5:2 resonance. The analysis of dynamical stability shows that the system is highly unstable or chaotic when planetary masses of about 5 MJup for b and 7 MJup for the other planets are adopted. Significant regions of dynamical stability for timescales of tens of Myr are found when adopting planetary masses of about 3.5, 5, 5, and 5 Mjup for HR 8799 b, c, d, and e respectively. These masses are below the current estimates based on the stellar age (30 Myr) and theoretical models of substellar objects.

We present diffraction-limited \ks band and \lprime adaptive optics images of the edge-on debris disk around the nearby F2 star HD 15115, obtained with a single 8.4 m primary mirror at the Large Binocular Telescope. At \ks band the disk is detected at signal-to-noise per resolution element (SNRE) \about 3-8 from \about 1-2\fasec 5 (45-113 AU) on the western side, and from \about 1.2-2\fasec 1 (63-90 AU) on the east. At \lprime the disk is detected at SNRE \about 2.5 from \about 1-1\fasec 45 (45-90 AU) on both sides, implying more symmetric disk structure at 3.8 \microns . At both wavelengths the disk has a bow-like shape and is offset from the star to the north by a few AU. A surface brightness asymmetry exists between the two sides of the disk at \ks band, but not at \lprime . The surface brightness at \ks band declines inside 1\asec (\about 45 AU), which may be indicative of a gap in the disk near 1\asec. The \ks - \lprime disk color, after removal of the stellar color, is mostly grey for both sides of the disk. This suggests that scattered light is coming from large dust grains, with 3-10 \microns -sized grains on the east side and 1-10 \microns dust grains on the west. This may suggest that the west side is composed of smaller dust grains than the east side, which would support the interpretation that the disk is being dynamically affected by interactions with the local interstellar medium.

As the only directly imaged multiple planet system, HR 8799 provides a unique opportunity to study the physical properties of several planets in parallel. In this paper, we image all four of the HR 8799 planets at H-band and 3.3 microns with the new LBT adaptive optics system, PISCES, and LBTI/LMIRCam. Our images offer an unprecedented view of the system, allowing us to obtain H and 3.3$ micron photometry of the innermost planet (for the first time) and put strong upper-limits on the presence of a hypothetical fifth companion. We find that all four planets are unexpectedly bright at 3.3 microns compared to the equilibrium chemistry models used for field brown dwarfs, which predict that planets should be faint at 3.3 microns due to CH4 opacity. We attempt to model the planets with thick-cloudy, non-equilibrium chemistry atmospheres, but find that removing CH4 to fit the 3.3 micron photometry increases the predicted L' (3.8 microns) flux enough that it is inconsistent with observations. In an effort to fit the SED of the HR 8799 planets, we construct mixtures of cloudy atmospheres, which are intended to represent planets covered by clouds of varying opacity. In this scenario, regions with low opacity look hot and bright, while regions with high opacity look faint, similar to the patchy cloud structures on Jupiter and L/T transition brown-dwarfs. Our mixed cloud models reproduce all of the available data, but self-consistent models are still necessary to demonstrate their viability.

Generalized Scidar (GS) measurements taken at the Paranal Observatory in November/December 2007 in the context of a site qualification for the future European Extremely Large Telescope E-ELT are re-calibrated to overcome the bias induced on the CN2 profiles by a not correct normalization of the autocorrelation of the scintillation maps that has been recently identified in the GS technique. A complete analysis of the GS corrected measurements as well as of the corrected errors is performed statistically as well as on individual nights and for each time during all nights. The relative errors of the CN2 profiles can reach up to 60% in some narrow temporal windows and some vertical slabs, the total seeing up to 12% and the total integrated turbulence J up to 21%. However, the statistic analysis tells us that the absolute errors of the median values of the total seeing is 0.06 arcsec (relative error 5.6%), for the boundary seeing 0.05 arcsec (relative error 5.6%) and for the seeing in the free atmosphere 0.04 arcsec (relative error 9%). We find that, in spite of the fact that the relative error increases with the height, the boundary and the free atmosphere seeing contribute in an equivalent way to the error on the total seeing in absolute terms. Besides, we find that there are no correlations between the relative errors and the value of the correspondent seeing. The absolute error of the median value of the isoplanatic angle is 0.13 arcsec (relative error 6.9%).

The mesoscale model Meso-NH is used to simulate the optical turbulence at Mt Graham (Arizona, USA), site of the Large Binocular Telescope. Measurements of the CN2-profiles obtained with a generalized scidar from 41 nights are used to calibrate and quantify the model's ability to reconstruct the optical turbulence. The measurements are distributed over different periods of the year, permitting us to study the model's performance in different seasons. A statistical analysis of the simulations is performed for all the most important astroclimatic parameters: the CN2-profiles, the seeing {\epsilon}, the isoplanatic angle {\theta}0 and the wavefront coherence time {\tau}0. The model shows a general good ability in reconstructing the morphology of the optical turbulence (the shape of the vertical distribution of CN2) as well as the strength of all the integrated astroclimatic parameters. The relative error (with respect to measurements) of the averaged seeing on the whole atmosphere for the whole sample of 41 nights is within 9.0 %. The median value of the relative error night by night is equal to 18.7 %, so that the model still maintains very good performances. Comparable percentages are observed in partial vertical slabs (free atmosphere and boundary layer) and in different seasons (summer and winter). We prove that the most urgent problem, at present, is to increase the ability of the model in reconstructing very weak and very strong turbulence conditions in the high atmosphere. This mainly affects the model's performances for the isoplanatic angle predictions, for which the median value of the relative error night by night is equal to 35.1 %. No major problems are observed for the other astroclimatic parameters. A variant to the standard calibration method is tested but we find that it does not provide better results, confirming the solid base of the standard method.

We present in this study a mapping of the optical turbulence (OT) above different astronomical sites. The mesoscale model Meso-NH was used together with the Astro-Meso-Nh package and a set of diagnostic tools allowing for a full 3D investigation of the Cn2. The diagnostics implemented in the Astro-Meso-Nh, allowing for a full 3D investigation of the OT structure in a volumetric space above different sites, are presented. To illustrate the different diagnostics and their potentialities, we investigated one night and looked at instantaneous fields of meteorologic and astroclimatic parameters. To show the potentialities of this tool for applications in an Observatory we ran the model above sites with very different OT distributions: the antarctic plateau (Dome C, Dome A, South Pole) and a mid-latitude site (Mt. Graham, Arizona). We put particular emphasis on the 2D maps of integrated astroclimatic parameters (seeing, isoplanatic angles) calculated in different slices at different heights in the troposhere. This is an useful tool of prediction and investigation of the turbulence structure. It can support the optimization of the AO, GLAO and MCAO systems running at the focus of the ground-based telescopes.From this studies it emerges that the astronomical sites clearly present different OT behaviors. Besides, our tool allowed us for discriminating these sites.

A characterization of the optical turbulence vertical distribution and all the main integrated astroclimatic parameters derived from the CN2 and the wind speed profiles above Mt. Graham is presented. The statistic includes measurements related to 43 nights done with a Generalized Scidar (GS) used in standard configuration with a vertical resolution of ~1 km on the whole 20-22 km and with the new technique (HVR-GS) in the first kilometer. The latter achieves a resolution of ~ 20-30 m in this region of the atmosphere. Measurements done in different periods of the year permit us to provide a seasonal variation analysis of the CN2. A discretized distribution of the typical CN2 profiles useful for the Ground Layer Adaptive Optics (GLAO) simulations is provided and a specific analysis for the LBT Laser Guide Star system ARGOS case is done including the calculation of the 'gray zones' for J, H and K bands. Mt. Graham confirms to be an excellent site with median values of the seeing without dome contribution equal to 0.72", the isoplanatic angle equal to 2.5" and the wavefront coherence time equal to 4.8 msec. We provide a cumulative distribution of the percentage of turbulence developed below H* where H* is included in the (0,1 km) range. We find that 50% of the whole turbulence develops in the first 80 m from the ground. The turbulence decreasing rate is very similar to what has been observed above Mauna Kea.

Atmospherical mesoscale models can offer unique potentialities to characterize and discriminate potential astronomical sites. Our team has recently completely validated the Meso-Nh model above Dome C (Lascaux et al. 2009, 2010). Using all the measurements of CN2 profiles (15 nights) performed so far at Dome C during the winter time (Trinquet et al. 2008) we proved that the model can reconstruct, on rich statistical samples, reliable values of all the three most important parameters characterizing the turbulence features of an antarctic site: the surface layer thickness, the seeing in the free atmosphere and in the surface layer. Using the same Meso-Nh model configuration validated above Dome C, an extended study is now on-going for other sites above the antarctic plateau, more precisely South Pole and Dome A. In this contribution we present the most important results obtained in the model validation process and the results obtained in the comparison between different astronomical sites above the internal plateau. The Meso-Nh model confirms its ability in discriminating between different optical turbulence behaviors, and there is evidence that the three sites have different characteristics regarding the seeing and the surface layer thickness. We highlight that this study provides the first homogeneous estimate, done with comparable statistics, of the optical turbulence developed in the whole 20-22 km above the ground at Dome C, South Pole and Dome A.

In two recent papers the mesoscale model Meso-NH, joint with the Astro-Meso-NH package, has been validated at Dome C, Antarctica, for the characterization of the optical turbulence. It has been shown that the meteorological parameters (temperature and wind speed, from which the optical turbulence depends on) as well as the Cn2 profiles above Dome C were correctly statistically reproduced. The three most important derived parameters that characterize the optical turbulence above the internal antarctic plateau: the surface layer thickness, the seeing in the free-atmosphere and in the total atmosphere showed to be in a very good agreement with observations. Validation of Cn2 has been performed using all the measurements of the optical turbulence vertical distribution obtained in winter so far. In this paper, in order to investigate the ability of the model to discriminate between different turbulence conditions for site testing, we extend the study to two other potential astronomical sites in Antarctica: Dome A and South Pole, which we expect to be characterized by different turbulence conditions. The optical turbulence has been calculated above these two sites for the same 15 nights studied for Dome C and a comparison between the three sites has been performed.

The characterization of the wind speed vertical distribution V(h) is fundamental for an astronomical site for many different reasons: (1) the wind speed shear contributes to trigger optical turbulence in the whole troposphere, (2) a few of the astroclimatic parameters such as the wavefront coherence time (tau_0) depends directly on V(h), (3) the equivalent velocity V_0, controlling the frequency at which the adaptive optics systems have to run to work properly, depends on the vertical distribution of the wind speed and optical turbulence. Also, a too strong wind speed near the ground can introduce vibrations in the telescope structures. The wind speed at a precise pressure (200 hPa) has frequently been used to retrieve indications concerning the tau_0 and the frequency limits imposed to all instrumentation based on adaptive optics systems, but more recently it has been proved that V_200 (wind speed at 200 hPa) alone is not sufficient to provide exhaustive elements concerning this topic and that the vertical distribution of the wind speed is necessary. In this paper a complete characterization of the vertical distribution of wind speed strength is done above Mt.Graham (Arizona, US), site of the Large Binocular Telescope. We provide a climatological study extended over 10 years using the operational analyses from the European Centre for Medium-Range Weather Forecasts (ECMWF), we prove that this is representative of the wind speed vertical distribution at Mt. Graham with exception of the boundary layer and we prove that a mesoscale model can provide reliable nightly estimates of V(h) above this astronomical site from the ground up to the top of the atmosphere (~ 20 km).

The atmospheric properties above three sites (Dome A, Dome C and the South Pole) are investigated for astronomical applications using the monthly median of the analyses from the ECMWF (European Centre for Medium-Range Weather Forecasts). Radiosoundings extended on a yearly time-scale at the South Pole and Dome C are used to quantify the reliability of the ECMWF analyses in the free atmosphere as well as in the boundary and surface layers, and to characterize the median wind speed in the first 100 m above the two sites. Thermodynamic instabilities in the free atmosphere above the three sites are quantified with monthly median values of the Richardson number. We will present a ranking of the sites with respect to the thermodynamic stability, using the Richardson number, and with respect to the wind speed, in the free atmosphere (using ECMWF analyses) as well as in the surface layer (using radiosoundings).

During winter and springtime, the flow above Antarctica at high altitude (upper troposphere and stratosphere) is dominated by the presence of a vortex centered above the continent. It lasts typically from August to November. This vortex is characterized by a strong cyclonic jet centered above the polar high. In a recent study of our group (Hagelin et al., 2008) of four different sites in the Antarctic internal plateau (South Pole, Dome C, Dome A and Dome F), it was made the hypothesis that the wind speed strength in the upper atmosphere should be related to the distance of the site to the center of the Antarctic polar vortex. This high altitude wind is very important from an astronomical point of view since it might trigger the onset of the optical turbulence and strongly affect other optical turbulence parameters. What we are interested in here is to localize the position of the minimum value of the wind speed at high altitude in order to confirm the hypothesis of Hagelin et al. (2008).

It has already been demonstrated that a mesoscale meteorological model such as Meso-NH is highly reliable in reproducing 3D maps of optical turbulence. Preliminary measurements above the Antarctic Plateau have so far indicated a pretty good value for the seeing: around 0.3" at Dome C. However some uncertainties remain. That's why our group is focusing on a detailed study of the atmospheric flow and turbulence in the internal Antarctic Plateau. Our intention is to use the Meso-NH model to do predictions of the atmospheric flow and the corresponding optical turbulence in the internal plateau. The use of this model has another huge advantage: we have access to informations inside an entire 3D volume which is not the case with observations only. Two different configurations have been used: a low horizontal resolution (with a mesh-size of 100 km) and a high horizontal resolution with the grid-nesting interactive technique (with a mesh-size of 1 km in the innermost domain centered above the area of interest). We present here the turbulence distribution reconstructed by Meso-NH for 16 nights monitored in winter time 2005, looking at the the seeing and the surface layer thickness.

The atmospheric properties above three sites on the Internal Antarctic Plateau are investigated for astronomical applications calculating the monthly median of the analysis-data from ECMWF (European Centre for Medium-Range Weather Forecasts) for an entire year (2005) thus covering all seasons. Radiosoundings extended on a yearly time scale from Dome C and the South Pole are used to verify the reliability of the analyses in the free atmosphere and to study the wind speed in the first 100 m as the analysis-data are not optimized for this altitude-range. The wind speed in the free atmosphere is obtained from the ECMWF analyses from all three sites. It appears that the strength of the wind speed in the upper atmosphere in winter is correlated to the distance of the site from the centre of the polar high. The Richardson number is employed to investigate the stability of the free atmosphere and, consequently, the probability to trigger thermodynamic instabilities above the three sites. We find that, in a large majority of the cases, the free atmosphere over the Internal Antarctic Plateau is more stable than at mid-latitude sites. Given these data we can obtain a ranking of the three sites with respect to wind speed, in the free atmosphere as well as in the surface layer, and with respect to the stability of the atmosphere, using the Richardson number.

A characterization of the optical turbulence vertical distribution (Cn2 profiles) and all the main integrated astroclimatic parameters derived from the Cn2 and the wind speed profiles above the site of the Large Binocular Telescope (Mt. Graham, Arizona, US) is presented. The statistic includes measurements related to 43 nights done with a Generalized Scidar (GS) used in standard configuration with a vertical resolution Delta(H)~1 km on the whole 20 km and with the new technique (HVR-GS) in the first kilometer. The latter achieves a resolution Delta(H)~20-30 m in this region of the atmosphere. Measurements done in different periods of the year permit us to provide a seasonal variation analysis of the Cn2. A discretized distribution of Cn2 useful for the Ground Layer Adaptive Optics (GLAO) simulations is provided and a specific analysis for the LBT Laser Guide Star system ARGOS (running in GLAO configuration) case is done including the calculation of the 'gray zones' for J, H and K bands. Mt. Graham confirms to be an excellent site with median values of the seeing without dome contribution epsilon = 0.72", the isoplanatic angle theta0 = 2.5" and the wavefront coherence time tau0= 4.8 msec. We find that the optical turbulence vertical distribution decreases in a much sharper way than what has been believed so far in proximity of the ground above astronomical sites. We find that 50% of the whole turbulence develops in the first 80+/-15 m from the ground. We finally prove that the error in the normalization of the scintillation that has been recently put in evidence in the principle of the GS technique, affects these measurements with an absolutely negligible quantity (0.04").

The internal antarctic plateau revealed in the last years to be a site with interesting potentialities for the astronomical applications due to the extreme dryness and low temperatures, the typical high altitude of the plateau, the weak level of turbulence in the free atmosphere down to a just few tens of meters from the ground and the thin optical turbulence layer developed at the ground. The main goal of a site testing assessment above the internal antarctic plateau is to characterize the site (optical turbulence and classical meteorological parameters) and to quantify which is the gain we might obtain with respect to equivalent astronomical observations done above mid-latitude sites to support plans for future astronomical facilities. Our group is involved, since a few years, in studies related to the assessment of this site for astronomical applications that include the characterization of the meteorological parameters and optical turbulence provided by general circulation models as well as mesoscale atmospherical models and the quantification of the performances of Adaptive Optics (AO) systems. In this talk I will draw the status of art of this site assessment putting our studies in the context of the wide international site testing activity that has been done in Antarctica. I will focus on the site assessment relevant for astronomical applications to be done in the visible up to the near infrared ranges, i.e. those ranges for which the optical turbulence represents a perturbing element for the quality of the images and the AO techniques an efficient tool to correct these wavefront perturbations.

In a recent paper the authors presented an extended study aiming at simulating the classical meteorological parameters and the optical turbulence at Dome C during the winter with the atmospherical mesoscale model Meso-NH. A statistical analysis has been presented and the conclusions of that paper have been very promising. Wind speed and temperature fields revealed to be very well reconstructed by the Meso-NH model with better performances than what has been achieved with the European Centre for Medium-Range Weather Forecast (ECMWF) global model, especially near the surface. All results revealed to be resolution-dependent and it has been proved that a grid-nesting configuration (3 domains) with a high horizontal resolution (1km) for the innermost domain is necessary to reconstruct all the optical turbulence features with a good correlation to measurements. High resolution simulations provided an averaged surface layer thickness just ~14 m higher than what is estimated by measurements, the seeing in the free atmosphere showed a dispersion from the observed one of just a few hundredths of an arcsecond (~0.05"). The unique limitation of the previous study was that the optical turbulence in the surface layer appeared overestimated by the model in both low and high resolution modes. In this study we present the results obtained with an improved numerical configuration. The same 15 nights have been simulated, and we show that the model results now match almost perfectly the observations in all their features: the surface thickness, the seeing in the free atmosphere as well as in the surface layer. This result permits us to investigate now other antarctic sites using a robust numerical model well adapted to the extreme polar conditions (Meso-NH).

Mesoscale model such as Meso-Nh have proven to be highly reliable in reproducing 3D maps of optical turbulence (see Refs. 1, 2, 3, 4) above mid-latitude astronomical sites. These last years ground-based astronomy has been looking towards Antarctica. Especially its summits and the Internal Continental Plateau where the optical turbulence appears to be confined in a shallow layer close to the icy surface. Preliminary measurements have so far indicated pretty good value for the seeing above 30-35 m: 0.36" (see Ref. 5) and 0.27" (see Refs. 6, 7) at Dome C. Site testing campaigns are however extremely expensive, instruments provide only local measurements and atmospheric modelling might represent a step ahead towards the search and selection of astronomical sites thanks to the possibility to reconstruct 3D Cn2 maps over a surface of several kilometers. The Antarctic Plateau represents therefore an important benchmark test to evaluate the possibility to discriminate sites on the same plateau. Our group8 has proven that the analyses from the ECMWF global model do not describe with the required accuracy the antarctic boundary and surface layer in the plateau. A better description could be obtained with a mesoscale meteorological model. In this contribution we present the progress status report of numerical simulations (including the optical turbulence - Cn2) obtained with Meso-Nh above the internal Antarctic Plateau. Among the topic attacked: the influence of different configurations of the model (low and high horizontal resolution), use of the grid-nesting interactive technique, forecasting of the optical turbulence during some winter nights.

Due to the recent interest shown by astronomers towards the Antarctic Plateau as a potential site for large astronomical facilities, we assisted in the last years to a strengthening of site testing activities in this region, particularly at Dome C. Most of the results collected so far concern meteorologic parameters and optical turbulence measurements based on different principles using different instruments. At present we have several elements indicating that, above the first 20-30 meters, the quality of the optical turbulence above Dome C is better than above whatever other site in the world. The challenging question, crucial to know which kind of facilities to build on, is to establish how much better the Dome C is than a mid-latitude site. In this contribution we will provide some complementary elements and strategies of analysis aiming to answer to this question. We will try to concentrate the attention on critical points, i.e. open questions that still require explanation/attention.

The image quality from Ground-Layer Adaptive Optics (GLAO) can be gradually increased with decreased contiguous field of view. This trade-off is dependent on the vertical profile of the optical turbulence (Cn2 profiles). It is known that the accuracy of the vertical distribution measured by existing Cn2 profiling techniques is currently quite uncertain for wide field performance predictions 4 to 20 arcminutes. With assumed uncertainties in measurements from Generalized-SCIDAR (GS), SODAR plus MASS we quantify the impact of this uncertainty on the trade-off between field of view and image quality for photometry of science targets at the resolution limit. We use a point spread function (PSF) model defined analytically in the spatial frequency domain to compute the relevant photometry figure of merit at infrared wavelengths. Statistics of this PSF analysis on a database of Cn2 measurements are presented for Mt. Graham, Arizona and Dome C, Antarctica. This research is part of the activities of ForOT (3D Forecasting of Optical Turbulence above astronomical sites).

The simulation of the optical turbulence (OT) for astronomical applications obtained with non-hydrostatic atmospherical models at meso-scale presents, with respect to measurements, some advantages. The future of the ground-based astronomy relies upon the potentialities and feasibility of the ELTs. Our ability in knowing, controlling and 'managing' the effects of the turbulence on such a new generation telescopes and facilities are determinant to assure their competitiveness with respect to the space astronomy. In the past several studies have been carried out proving the feasibility of the simulation of realistic Cn2 profiles above astronomical sites. The European Community (FP6 Program) decided recently to fund a Project aiming, from one side, to prove the feasibility of the OT forecasts and the ability of meso-scale models in discriminating astronomical sites from optical turbulence point of view and, from the other side, to boost the development of this discipline at the borderline between the astrophysics and the meteorology. In this contribution I will present the scientific and technological goals of this project, the challenges for the ground-based astronomy that are related to the success of such a project and the international synergies that have been joint to optimize the results.

These last years ground-based astronomy has been looking towards Antarctica, especially its summits and the internal continental plateau where the optical turbulence (OT) appears to be confined in a shallow layer close to the surface. Preliminary measurements have so far indicated pretty good value for the seeing above 30-35 m: 0.36" (Agabi et al. 2006), 0.27" (Lawrence et al. 2004) and 0.3" (Trinquet et al. 2008) at Dome C. Site testing campaigns are however extremely expensive, instruments provide only local measurements and atmospheric modeling might represent a step ahead towards the search and selection of astronomical sites thanks to the possibility to reconstruct 3D Cn2 maps over a surface of several kilometers. The Antarctic Plateau represents therefore an important benchmark test to evaluate the possibility to discriminate sites on the same plateau. Our group (Hagelin et al. 2008) has proven that the analyses from the ECMWF global model do not describe with the required accuracy the antarctic boundary and surface layers in the plateau. A better description could be obtained with a mesoscale model. The Meso-Nh model has proven to be reliable in reproducing 3D maps of OT above mid-latitude astronomical sites (Masciadri et al. 1999ab, 2004, Masciadri and Jabouille 2005). In this paper we study the ability of the Meso-Nh model in reconstructing the meteorological parameters as well as the OT above Dome C with different model configurations. We concentrate our attention on the model abilities in reproducing the OT surface layer thickness Hsl and the integral of the Cn2 in the free atmosphere and in the surface layer. It is worth to highlight that these are the first estimates ever done so far with a mesoscale model of the optical turbulence above the internal Antarctic Plateau.

The atmospheric properties above three sites (Dome C, Dome A and the South Pole) on the Internal Antarctic Plateau are investigated for astronomical applications using the monthly median of the analyses from ECMWF (the European Centre for Medium-Range Weather Forecasts). Radiosoundings extended on a yearly time scale at the South Pole and Dome C are used to quantify the reliability of the ECMWF analyses in the free atmosphere as well as in the boundary and surface layers, and to characterize the median wind speed in the first 100 m above the two sites. Thermodynamic instability properties in the free atmosphere above the three sites are quantified with monthly median values of the Richardson number. We find that the probability to trigger thermodynamic instabilities above 100 m is smaller on the Internal Antarctic Plateau than on mid-latitude sites. In spite of the generally more stable atmospheric conditions of the Antarctic sites compared to mid-latitude sites, Dome C shows worse thermodynamic instability conditions than those predicted above the South Pole and Dome A above 100 m. A rank of the Antarctic sites done with respect to the strength of the wind speed in the free atmosphere (ECMWF analyses) as well as the wind shear in the surface layer (radiosoundings) is presented.

We examine the implications for the distribution of extrasolar planets based on the null results from two of the largest direct imaging surveys published to date. Combining the measured contrast curves from 22 of the stars observed with the VLT NACO adaptive optics system by Masciadri et al. (2005), and 48 of the stars observed with the VLT NACO SDI and MMT SDI devices by Biller et al. (2007) (for a total of 60 unique stars; the median star for our survey is a 30 Myr K2 star at 25 pc), we consider what distributions of planet masses and semi-major axes can be ruled out by these data, based on Monte Carlo simulations of planet populations. We can set this upper limit with 95% confidence: the fraction of stars with planets with semi-major axis from 20 to 100 AU, and mass >4 M_Jup, is 20% or less. Also, with a distribution of planet mass of dN/dM ~ M^-1.16 between 0.5-13 M_Jup, we can rule out a power-law distribution for semi-major axis (dN/da ~ a^alpha) with index 0 and upper cut-off of 18 AU, and index -0.5 with an upper cut-off of 48 AU. For the distribution suggested by Cumming et al. (2007), a power-law of index -0.61, we can place an upper limit of 75 AU on the semi-major axis distribution. At the 68% confidence level, these upper limits state that fewer than 8% of stars have a planet of mass >4 M_Jup between 20 and 100 AU, and a power-law distribution for semi-major axis with index 0, -0.5, and -0.61 cannot have giant planets beyond 12, 23, and 29 AU, respectively. In general, we find that even null results from direct imaging surveys are very powerful in constraining the distributions of giant planets (0.5-13 M_Jup) at large separations, but more work needs to be done to close the gap between planets that can be detected by direct imaging, and those to which the radial velocity method is sensitive.

The commonality of collisionally replenished debris around main sequence stars suggests that minor bodies are frequent around Sun-like stars. Whether or not debris disks in general are accompanied by planets is yet unknown, but debris disks with large inner cavities - perhaps dynamically cleared - are considered to be prime candidates for hosting large-separation massive giant planets. We present here a high-contrast VLT/NACO angular differential imaging survey for eight such cold debris disks. We investigated the presence of massive giant planets in the range of orbital radii where the inner edge of the dust debris is expected. Our observations are sensitive to planets and brown dwarfs with masses >3 to 7 Jupiter mass, depending on the age and distance of the target star. Our observations did not identify any planet candidates. We compare the derived planet mass upper limits to the minimum planet mass required to dynamically clear the inner disks. While we cannot exclude that single giant planets are responsible for clearing out the inner debris disks, our observations constrain the parameter space available for such planets. The non-detection of massive planets in these evacuated debris disks further reinforces the notion that the giant planet population is confined to the inner disk (<15 AU).

We present the results of a survey of 45 young (<250 Myr), close (<50 pc) stars with the Simultaneous Differential Imager (SDI) implemented at the VLT and the MMT for the direct detection of extrasolar planets. As part of the survey, we observed 54 objects total, of which 45 were close, young stars, 2 were somewhat more distant (<150 pc), very young (<10 Myr) stars, 3 were stars with known radial velocity planets, and 4 were older, very nearby (<20 pc) solar analogues. Our SDI devices use a double Wollaston prism and a quad filter to take images simultaneously at 3 wavelengths surrounding the 1.62 um methane absorption bandhead found in the spectrum of cool brown dwarfs and gas giant planets. In our VLT data, we achieved H band contrasts > 10 mag (5 sigma) at a separation of 0.5" from the primary star on 45% of our targets and H band contrasts of > 9 mag at a separation of 0.5'' on 80% of our targets. With this degree of attenuation, we should be able to image (5sigma detection) a 5 M_{Jup} planet 15 AU from a 70 Myr K1 star at 15 pc or a 5 M_{Jup} planet at 2 AU from a 12 Myr M star at 10 pc. We believe that our SDI images are the highest contrast astronomical images ever made from ground or space for methane rich companions <1'' from their star. For the best 20 of our survey stars, we attained 50% 5 sigma completeness for 6-10 M_Jup planets at semi-major axes of 20-40 AU. Thus, our completeness levels are sufficient to significantly test theoretical planet distributions. From our survey null result, we can rule out (at the 98% confidence/2.0sigma level) a model planet population using a planet distribution where N(a) $\propto$ constant out to a distance of 45 AU (further model assumptions discussed within).

The active K2V star $\epsilon$ Eri hosts the most nearby known extrasolar planet. With an angular separation of about 1" on average, and an age of a few to several hundred Myrs, $\epsilon$ Eri b is one of the prime candidates for becoming the first definitive extrasolar planet imaged directly. We present a multi-epoch deep differential imaging survey performed with NACO-SDI at the VLT with the aim of finding the planet. The results are combined with recent astrometry in an attempt to further constrain the detection limits. No convincing candidate is found among the many coherent structures that constitute the residual speckle noise, which is the dominant noise at small angular scales. We present our detection limits, compare them with the estimated brightness of $\epsilon$ Eri b, and analyze how the limits can be improved further. It is found that integration time remains a very important parameter for achieving good results, even in the speckle-dominated regimes. The results yield new, improved upper 3$\sigma$ limits on the absolute H-band (1.6 $\mu$m) brightness of the 1.55 $M_{\rm jup}$ companion of 19.1 to 19.5 mag, depending on the specific age of the system.

We have carried out a direct and differential imaging search for sub-stellar companions to eps Indi A using the adaptive optics system NACO at the ESO VLT. The observations were carried out in September 2004 with NACO/SDI as well as with NACO's S27 camera in the H and Ks filters. The SDI data cover an area of \~2.8" around eps Indi A. No detection was achieved in the inner neighbourhood down to 53 Mj (5 sigma confidence level) at a separation > 0.4" (1.45 AU) and down to 21 Mj for separations > 1.3" (4.7 AU). To cover a wider field of view, observations with the S27 camera and a coronagraphic mask were obtained. We detected a faint source at a separation of (7.3 +/- 0.1)" and a position angle of (302.9 +/- 0.8) degree. The photometry for the candidate companion yields m(H)=(16.45 +/- 0.04)mag and m(Ks) = (15.41 +/- 0.06)mag, respectively. Those magnitudes and the resulting color (H-Ks) = (1.04 +/- 0.07)mag fit best to a spectral type of L5 - L9.5 if it is bound. Observations done with HST/NICMOS by M. Endl have shown the source to be a background object.

In this paper we present the characterization of all the principal meteorological parameters (wind speed and direction, pressure, absolute and potential temperature) extended over 25 km from the ground and over two years (2003 and 2004) above the Antarctic site of Dome C. The data set is composed by 'analyses' provided by the General Circulation Model (GCM) of the European Center for Medium Weather Forecasts (ECMWF) and they are part of the catalog MARS. A monthly and seasonal (summer and winter time) statistical analysis of the results is presented. The Richardson number is calculated for each month of the year over 25 km to study the stability/instability of the atmosphere. This permits us to trace a map indicating where and when the optical turbulence has the highest probability to be triggered on the whole troposphere, tropopause and stratosphere. We finally try to predict the best expected isoplanatic angle and wavefront coherence time employing the Richardson number maps, the wind speed profiles and simple analytical models of CN2 vertical profiles.

We discuss the preliminary results of a survey of young (<300 Myr), close (<50 pc) stars with the Simultaneous Differential Extrasolar Planet Imager (SDI) implemented at the VLT and the MMT. SDI uses a double Wollaston prism and a quad filter to take 4 identical images simultaneously at 3 wavelengths surrounding the 1.62 um methane bandhead found in the spectrum of cool brown dwarfs and gas giants. By performing a difference of images in these filters, speckle noise from the primary can be significantly attenuated, resulting in photon noise limited data. In our survey data, we achieved H band contrasts >25000 (5 sigma Delta F1(1.575 um) > 10 mag, Delta H > 11.5 mag for a T6 spectral type) at a separation of 0.5" from the primary star. With this degree of attenuation, we should be able to image (5 sigma detection) a 2-4 Jupiter mass planet at 5 AU around a 30 Myr star at 10 pc. We are currently completing our survey of young, nearby stars. We have obtained complete datasets for 35 stars in the southern sky (VLT) and 7 stars in the northern sky (MMT). We believe that our SDI images are the highest contrast astronomical images ever made from ground or space for methane rich companions.

We report on a survey devoted for the search of exo-planets around young and nearby stars carried out with NACO at the VLT. The detection limit for 28 among the best available targets vs. the angular separation from the star is presented. The non-detection of any planetary mass companion in our survey is used to derive, for the first time, the frequency of the upper limit of the projected separation planet-stars. In particular, we find that in 50 % of cases, no 5 MJ planet (or more massive) has been detected at projected separations larger than 14 AU and no 10 MJ planet (or more massive) has been detected at projected separations larger than 8.5 AU. In 100 % of cases, these values increase to 36 AU and 65 AU respectively. The excellent sensitivity reached by our study leads to a much lower upper limit of the projected planet-star separation compared with previous studies. For example, for the Beta Pictoris group, (~ 12 Myr), we did not detect any 10 MJ planet at distances larger than 15 AU. A previous study carried out with 4 m class telescopes put an upper limit for 10 MJ planets at ~ 60 AU. For our closest target (V 2306 Oph - d = 4.3 pc) it is shown that it would be possible to detect a 10 MJ planet at a minimum projected separation from the star of 1 AU and a 5 MJ planet at a minimum projected separation of 3.7 AU. Our results are discussed with respect to mechanisms explaining planet formation and migration and to forthcoming observational strategies and future planet finder observations from the ground.

First evidences of IR excess and disk mass accretion (strong H$\alpha$ emission) around brown dwarfs seem to indicate the existence of circumstellar disks around these sub-stellar objects. Nothing is known at the present time about outflows which potentially might be launched from brown dwarfs, although jets are typically associated with the accretion in standard T Tauri star disks. In this paper we calculate the H$\alpha$ emission of internal working surfaces produced by a radiative jet in a neutral and in a photoionized environment as a function of the jet parameters (the ejection velocity $v_j$, shock velocity $v_s$, mass loss rate ${\dot M}$ and radius $r_j$ of the jet) and we provide estimates of the H$\alpha$ luminosity for the parameters of ``standard'' Herbig-Haro (HH) jets from T Tauri stars and for the parameters expected for jets from BDs. Interestingly, we find that while the mass loss rates associated with jets from BDs are found to be two orders of magnitude lower than the mass loss rates associated with ``standard'' HH jets (from T Tauri stars), their velocities are likely to be similar. Based on our calculations, we discuss the conditions in which jets from BDs can be detected and we conclude that the H$\alpha$ luminosities of internal working surfaces of jets from BDs in a photoionized environment should have only one order of magnitude lower than the H$\alpha$ luminosities of T Tauri jets in a neutral environment.

Observations of the standard late type M8 star LHS 2397aA were obtained at the ESO-VLT 8m telescope ``Yepun'' using the NAOS/CONICA Adaptive Optics facility. The observations were taken during the NACO commissioning, and the infrared standard star LHS 2397aA was observed in the H, and Ks broad band filters. In both bands the brown dwarf companion LHS2397aB was detected. Using a program recently developed (Bouy et al., 2003) for the detection of stellar binaries we calculated the principal astrometric parameters (angular binary separation and position angle P.A.) and the photometry of LHS 2397aA and LHS 2397aB. Our study largely confirms previous results obtained with the AO-Hokupa'a facility at Gemini-North (Freed et al., 2003); however a few discrepancies are observed.

The giant jets represent a fundamental trace of the historical evolution of the outflow activity over timescales which are comparable to the accretion time of the outflow sources in their main protostellar phase. The study of such huge jets provides the possibility of retrieving important elements related to the life of the outflow sources. In this paper, we study the role of precession (combined with jet velocity-variability and the resulting enhanced interaction with the surrounding environment) as a deceleration mechanism for giant jets using a numerical approach. We obtain predictions of H alpha intensity maps and position-velocity diagrams from 3D simulations of the giant HH 34 jet (including an appropriate ejection velocity time-variability and a precession of the outflow axis), and we compare them with previously published observations of this object. Our simulations represent a step forward from previous numerical studies of HH objects, in that the use of a 7-level, binary adaptive grid has allowed us to compute models which appropiately cover all relevant scales of a giant jet, from the ~ 100 AU jet radius close to the source to the ~ 1 pc length of the outflow. A good qualitative and quantitative agreement is found between the model predictions and the observations. Moreover, we show that a critical parameter for obtaining a better or worse agreement with the observations is the ratio rho_j/rho_a between the jet and the environmental densities. The implications of this result in the context of the current star formation models are discussed (ABRIDGED).