Using meteorological forecasts to predict astronomical ‘seeing’
The main limitations affecting visible-light and IR astronomy are related to visibility, atmospheric humidity, and spatial resolution. High-resolution imaging is fundamentally limited by optical turbulence (OT) in the atmosphere. OT is induced by dynamical turbulence developing within a stratified temperature field and depends on weather conditions. The resulting degraded telescope resolution is called ‘seeing.’ At the best ground-based observatories, characterized by a seeing of order 0.5 arcseconds, the resolution of a telescope with an aperture diameter of 10m is no better than 20cm (at ground level) at optical wavelengths.
At the present time, when large and expensive telescope projects are increasingly becoming the norm, observatory managers rely on a proper assessment and understanding of weather conditions for optimal scheduling of scientific operations. Flexible rotas (i.e., planning optimal use of the available instrumentation) are now frequently used at the largest astronomical observatories. Site testing in preparation for new observatory developments is another important challenge. For site characterization for the European 42m (diameter) Extremely Large Telescope (E-ELT)1 or the Antarctic ‘Dome C’ project,2 we designed the DIMM (differential image-motion monitor) and SSS (single-star scidar)3 optical monitors. Since monitoring requires the availability of specific instruments and logistics, and it is limited in place and time, we are investigating other ways to retrieve atmospheric OT conditions from meteorological measurements or forecasts.
Atmospheric conditions are defined by the speed, temperature, and composition of the constituent elements. The corresponding state of the atmosphere is governed by the Navier-Stokes equations prevalent above and around the site. For astronomical seeing prediction, one must obtain the local (microscopic) atmospheric conditions on the metric OT scale. Several models can be used, but most are computationally time consuming, unable to provide a global analysis for a given site, or use the meteorological forecast directly to predict the temporal OT-strength Cn2 profiles, where Cn2 is the OT index.
We analyzed4 the correlation between the Cn2 and macroscopic mean profiles of the horizontal wind-speed component and temperature, with a vertical resolution of 100m to 1km. We developed and validated a parametric model on the basis of our analysis of a large Cn2 and meteorological profile database.5 This model can be implemented easily in a forecasting scheme and used for astronomical site-management purposes (see Figure 1). Any mesoscale model can be used to provide hourly forecasts over 24 hours and deliver the main meteorological parameters of interest. The typical horizontal grid step is 1km, while we simultaneously use more than 30 vertical levels. Higher horizontal resolution may be achieved, depending on the terrain model. Our parametric model is initialized using global data provided by meteorological agencies, complemented by local measurements. The OT model is then applied to the meteorological forecast, thus providing predictions for the seeing, altitude, and intensity of the main OT-producing layers.
Here we present results obtained with the Weather Research and Forecasting Model and applied to La Palma (Canary Islands, Spain), where the Observatorio del Roque de los Muchachos (ORM) is sited on a crater rim in the northern part of the island. We ran the simulations for one month, with hourly sampling. The model was initialized using data from the Global Forecast System. Figure 2 shows a typical forecast, displaying the island's seeing conditions with 1km resolution. Figure 3 represents the temporal evolution of the Cn2 profiles at the ORM over the course of one month.
Meteorological studies and predictions are crucial to assess astronomical observing conditions. Seeing forecasts are of major importance since they enable planning for optimal use of observing time. They also allow tuning of adaptive-optics systems as well as site characterization. Conceptual tools are now available to calculate model results. Our next step is to compare our model results with measurements to validate the models as forecasting products. Once this has been achieved, this approach could be implemented and routinely used at astronomical observatories.
Hervé Trinquet received his PhD from the University of Nice in 2003. His current research relates to E-ELT contracts and site testing at Dome C. He has been involved in OT parametrization and analysis of gravity waves in the atmosphere.
Jean Vernin is a senior scientist involved in astronomical site testing and analysis of turbulent flows in stratified media.