Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

ViSiON: Visibility Service for Observing Nights

ViSiON: Visibility Service for Observing Nights Preparation of detailed night schedule prior to an observing run can be tedious, especially for solar system objects which coordinates are epoch-dependent. We aim at providing the community with a Web service compliant with Virtual Observatory (VO) standards, to create tables of observing conditions, together with airmass and sky charts, for an arbitrary list of targets. We take advantage of available VO services such as the SIMBAD astronomical database, the Aladin sky atlas, and the Miriade ephemerides generator to build a new service dedicated to the planning of observations. The requests for ephemerides charts and tables are handled by a VO-compliant Web service. For each date, and each target, coordinates in local and equatorial frames are computed, and used to select targets accordingly to user’s criteria for visibility. This new service, dubbed ViSiON for Visibility Service for Observing Night, is a new method of Miriade Web service hosted at IMCCE. It allows anyone to create graphics of observing conditions and tables summarizing them, provided as PDF, VOTable, and xHTML documents. Keywords: ephemerides, virtual observatory tools, methods: numerical 1. Introduction standards, to produce tables and charts of visibility for both fixed coordinates targets and moving solar system objects. To prepare and conduct observations from a ground-based This service is dubbed Visibility Service for Observing Night telescope, many astronomers create airmass charts, also called (ViSiON). Targeted users are observers willing to prepare an visibility plots. These graphics represent the elevation above upcoming night, or to check the observability of some targets the horizon, or altitude, of each target during the night, provid- on-the-fly, and developers of night scheduler and telescope con- ing a graphical and easy way to schedule the sequence of ob- trol software willing to include airmass charts in graphical in- servations, in which each target should be observed as closely terfaces. as possible to its highest altitude. Another type of visual guide- The present article is structured as following: we quickly de- lines are fish-eye sky charts, in which the successive positions scribe the core of ephemeris computation in Sect. 2. The al- of the targets are displayed by tracks on a stereographic projec- gorithm of ViSiON is detailed in Sect. 3, and the criteria which tion of the local frame. It also helps optimizing observations, define observability are presented in Sect. 4. The different ways by avoiding pointing the telescope back and forth from different to access to the service are provided in Sect. 5 and its links with regions of the sky. other VO services in Sect. 6 If many astronomers have developed their own piece of code for that purpose, several efforts have been made to make 2. Ephemeris computation such tools freely available to the community. The staralt pages , by J. Me´ndez, allows productions of such graphics. The The ephemerides of planets and small Solar system ob- skycalc software , maintained by J. Thorstensen, represents jects are computed into a quasi-inertial reference frame, tak- an easily-portable and interactive solution to these points. Yet, ing into account the post-Newtonian approximations. The ge- these tools have been developed for objects with fixed coordi- ometric positions of the major planets and the Moon are pro- nates, and the eight planets. Planning observations of moving vided by INPOP planetary theory (Fienga et al., 2014). Those Solar System Objects (SSOs) may therefore result tedious: one of small SSOs (asteroids, comets, Centaurs, trans-neptunian should compute first the coordinates of each object, and then objects) are calculated by numerical integration of the N- enter them into the software interface. This often results im- body perturbed problem (Gragg-Bulirsch-Stoer algorithm, see practicable if the list of targets is long, or many dates are envi- Bulirsch and Stoer, 1966; Stoer and Bulirsch, 1980), using sioned. their latest published osculating elements, from the astorb We aim here at providing the community with a simple (Bowell et al., 1993) and cometpro (Rocher and Cavelier, Web service, fully compliant with Virtual Observatory (VO) 1996) databases. The overall accuracy of asteroid and comet http://catserver.ing.iac.es/staralt/ 2 3 http://www.dartmouth.edu/ physics/labs/skycalc/flyer.html http://vo.imcce.fr/webservices/miriade/?vision Preprint submitted to PSS October 9, 2018 arXiv:1806.09425v1 [astro-ph.IM] 25 Jun 2018 ephemerides provided by ViSiON is at the level of tens of milli- Table 1: List of parameters reported for each visible target, for each requested arcseconds, mainly depending on the accuracy of their osculat- night. The values reported in ViSiON correspond to the geometry at the time of ing elements. The positions of natural satellites are obtained the highest altitude on sky. thanks to dedicated solutions of their motion, e.g., Lainey et al. Quantity Description Units (2004a,b, 2007) for Mars and Jupiter, Vienne and Duriez Target Target designation – (1995) for Saturn, Laskar and Jacobson (1987) for Uranus, and m Apparent magnitude mag Le Guyader (1993) for Neptune satellites. φ Apparent diameter arcsec The ephemerides of stars, and more generally of all D Duration of visibility window h:m SIMBAD objects (Wenger et al., 2000), are also computed into a Alt Altitude deg quasi-inertial reference frame, taking into account their proper Az. Azimuth deg motion, parallax and radial velocity, as long as the corre- RA Right Ascension h:m:s sponding input parameters are provided by the catalogue. The DEC Declination d:m:s ephemeris accuracy of SIMBAD objects is thus at the level of the Rate Apparent non-sidereal motion arcsec/h accuracy of their catalogued position. λ Galactic longitude deg All the ephemerides in ViSiON are computed relatively to a β Galactic latitude deg given observer, in the ICRF reference frame (mean J2000 equa- r Range to observer au tor and equinox), and expressed in the UTC time scale. The Δ Heliocentric distance au coordinates of the observer (the topocenter) can be either pro- α Solar phase angle deg vided directly by users (longitude, latitude, altitude), or by us- SEO Solar elongation deg ing the observatory code provided by IAU Minor Planet Center MEO Moon elongation deg (MPC) for listed observatories. following parameters can be specified (see Table 2 for details): 3. Description of ViSiON’s algorithm altitude, apparent magnitude, solar phase angle, solar elonga- For each requested date, a series of computations are per- tion, lunar elongation, angular diameter, duration of visibil- formed prior to selection of targets suitable for observation. ity window, and events limiting the “night” (sunrise/sunset, or First, the ephemerides of the Sun are computed, to determine civil/nautical/astronomical twilights). the instants of sunset and sunrise. The limits of the civil, nau- ViSiON proceeds as following: for each target, it computes tical, and astronomical twilights, corresponding to the times its ephemerides for each observing date, as described above ◦ ◦ ◦ when the altitude of the Sun is 6 , 12 , and 18 below the hori- (Sect. 3). Each parameter (apparent magnitude, elevation, etc) zon, are also computed at this stage. is checked against the selection criteria at every 5 min time step. Second, the phase and position of the Moon throughout the A time step is considered valid only if the tests against all the night are computed and stored. Finally, the topocentric astro- criteria have passed. Finally, the target will be considered valid metric coordinates of each target are computed, with a time step for the date if the number of time steps exceeds the specified of 5 minutes, and expressed as equatorial (RA, Dec), horizontal minimum duration of visibility (30 min by default) between the (Az, Alt), and galactic coordinates (λ , β ). G G two events defining the “night”. Additionally, both solar and lunar elongations are computed For each date, all the targets that pass the selection are drawn at each time step, together with (for Solar system objects only) on the graphics, and included in the summary table. Their or- the solar phase angle, the heliocentric distance and the range to der can also be chosen by users, among all the parameters listed the observer. These detailed coordinates are used to determine in Table 1. A double sorting is even possible, with objects be- for how long each target is visible during the night. The entire ing grouped by type first (asteroid, comet, planet, satellite, and list of criteria defining observability is provided in Sec. 4. If fixed-coordinate) then by the requested parameter. there is at least one target observable, the night is included in the result, and these detailed coordinates are the inputs to the 5. Access to the service airmass and sky charts (see Figs. 1 and 2 for examples). For each date, a table lists all the targets that are visible (see Table 4 In 2005, the IMCCE started to implement Virtual Obser- for an example), and summarizes their observing conditions at vatory (VO) compliant interfaces to its ephemerides services the time of highest altitude (see Table 1 for a description of each (Thuillot et al., 2005). A Web portal describe the various parameter). services, such as solar system objects identification (Skybot, Berthier et al., 2006), or general ephemerides computation 4. Selection and sorting criteria (Miriade, Berthier et al., 2009). All our services are acces- sible via Web services (based on SOAP and HTTP POST verb) A suite of parameters are tested to determine if a target shall which allow anyone to interface her/his own software with the be considered as visible. The thresholds for selection are pro- services, via HTTP request and Web forms, and are integrated vided as inputs, accordingly to each observer’s preference. The 4 5 http://www.minorplanetcenter.net/iau/lists/ObsCodes.html http://vo.imcce.fr/ 2 Sunset Twilight Twilight Sunrise h m h m h m h m 21 59 23 21 09 58 11 20 h m h m h m h m h m h m h m h m LST 9 01 11 01 13 02 15 02 17 02 19 03 21 03 23 03 Mars M70 Camilla Triton TYC5597−234−1 67P Moon SMC 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 21:00 23:00 01:00 03:00 05:00 07:00 09:00 11:00 Time (UT) st Figure 1: Example of airmass chart produced by ViSiON for a request a La Silla observatory (ESO) on June the 1 , 2016. The different levels of gray shading on each side correspond to the day light, civil, nautical, and astronomical twilights. Both universal time (UT) and local sidereal time (LST) are reported on x-axes, and altitude and airmass are provided on y-axes. Each curve represents a target, color-coded as a function of type: planet (Mars), asteroid (12 116 and Camilla), comet (67P), satellite (Triton), and fixed-coordinates such as galaxy (Small Magellanic Cloud: SMC), star (TYC5597-234-1), or clusters (M70). in several VO-compliant software such as the well-spread Al- or “s:” suffixes, e.g., p:Phobos, s:Enceladus for the corre- adin Sky Atlas (Bonnarel et al., 2000). sponding satellites of Mars and Saturn. There are two ways to use the ViSiON web service: by writ- For objects with fixed coordinates, ViSiON queries the ting a client to send requests to the Miriade server and to re- SIMBAD service. The targets are therefore inherited from ceive and analyze the response (as for instance the ephemerides SIMBAD, which currently contains above 8 millions entries. tool for observations in Service Mode hosted by ESO, see These targets are submitted to ViSiON using the “u:” prefix, Carry and Berthier, 2016), or by using a command line inter- e.g., u:Sirius, u:HIP 87937, u:M 31 for Sirius, the Hippar- face and a data transfer program such as curl or wget. We also cos star 87937, and Messier 31. For objects without a catalog propose a simple query form on our SSO VO portal . In any entry, it is always possible to directly submit equatorial coor- case, the service accepts nine parameters, described in Table 3. dinates (in decimal hours and degrees in astrometric J2000) Only the list of targets is mandatory, all the other parameters to ViSiON, with the “e:” prefix. For instance, the equatorial being optional. We describe them in detail here. North pole can be queried via e:0.0+90.0. In all cases, ViSiON accepts aliases for each target, that can Targets (name). Planets can be requested by setting the name be used to simplify the labels on the graphics, with the syntax parameter to their name or number, preceeded by the “p:” pre- prefix:name=alias. For instance in Figs. 1 and 2, the Small fix, for instance p:Mars and p:5 to request visibility of Mars Magellanic Cloud was requested as e:0.9-72=SMC, while and Jupiter. Any entry in astorb and cometpro can be re- Mars was requested using p:4 and M70 with u:M 70. This trieved, using the “a:” and “c:” prefixes, e.g., a:Ceres, a:21, option may be useful for minor planets with high numbers and or c:67P for the asteroids (1) Ceres and (21) Lutetia, and the provisional designation: to avoid the display of 2000 WO137 comet 67P/Churyumov-Gerasimenko. All the satellites present on the graphics when querying for a:123456, one can force the in the references listed in Sect. 2 are available, though the “p:” usage of the IAU number by specifying a:123456=123456. http://vo.imcce.fr/webservices/miriade/?forms#Vision Altitude ( ) airmass 5.60 2.90 2.00 1.55 1.30 1.15 1.06 1.02 1.0 N 67P Moon E W TYC5597−234−1 Triton Camilla Mars M70 SMC 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 st Figure 2: Example of sky chart produced by ViSiON for a request a La Silla observatory (ESO) on June the 1 , 2016. As in Fig. 1, each curve represents a target, color-coded as a function of its type. Topocenter (observer). The location on Earth of the observer and tag :{max:value } to set only one extrema, or max can be specified by providing the longitude (0,360 West), lat- tag :{min:value , max:value } to set both. For in- max i min itude (-90,+90 North), and altitude (meters above sea level) stance, to limit targets to sources brighter than V = 20, one or the IAU MPC observatory code, in which case ViSiON re- should specify mag:{max:20}. To select only objects matching trieves the coordinates of the observer from the MPC list. For all criteria during 2 hours between the sunrise and sunset, the instance, the Paris observatory can be specified as -2.336537 two parameters duration:120, event:sun must be speci- +48.8364 132.0 or 007. fied. Sorting options (sort). The entries in the summary tables (see Dates of computation (ep, nbd, step). The ep parameter Table 4 for an example) are sorted in increasing right ascen- sets the first (or only) date at which the visibility is sought sion by default. Changing this behavior is achieved by speci- (YYYY-MM-DD). If more than one night is desired, their total fying sort=key, in which key can take any value among the number can be specified with nbd (integer ≥1) and the num- reported parameters (Table 1). For instance, to sort object by ber of skipped nights between each by step (integer ≥0). For increasing apparent magnitude, one should specify sort=mv. It instance, four dates, one each week, starting from June, the st is also possible to first group targets by type before applying the 1 , 2016, these parameters should be set to ep=2016-06-01, sorting, using sort=type>key. nbd=4, and step=7. Output formats (mime). The results can be generated in three Selection criteria (cuts). Any number of selection crite- formats: VOTable, xHTML pages, and PDF file by specifying ria can be specified, using the syntax tag :value , ..., 1 1 the mime parameter. Most users will likely use the PDF and tag :value , where the tag are listed in Table 2. Tags N N i xHTML outputs, for convenient display and interaction. The accepting both a minimum and a maximum value shall VOTable output has been designed for a more advanced us- be set using the following syntax tag :{min:value } i min age, into a workflow for instance. It allows to directly query ViSiON from, e.g., a software to control observing queues, re- http://www.minorplanetcenter.net/iau/lists/ObsCodes.html trieving the values in the table, and the graphics for display 4 Table 2: Selection criteria in ViSiON. For each, we indicate if a minimum and/or a maximum value can be specified, and which are the default limits. The phase angle and apparent diameter selection only applies to solar system objects, and is disregarded for the other targets. Quantity Description Minimum Maximum Defaults Alt. Highest altitude reached!! [30 to 90 ] m Apparent magnitude!! [No default] α Solar phase angle!! [0 to 180 ] SEO Solar elongation! [No default] MEO Moon elongation! [No default] φ Apparent diameter!! [No default] D Duration of visibility window! [30 min] Event Start/End of the “night”! [civil] Table 3: List of input parameters to ViSiON (see Sect. 5 for the list of access points). For each, the default (between brackets) and range are specified. Parameter Description Range, or example values name Designation of the requested targets p:Mars, a:Pluto, u:M 31 observer Observer’s location (IAU code or geographic coordinates) [Paris, France] ep Starting epoch [now] — JD — YYYY-MM-DD nbd Number of date of computation [1] ≤ nbd < 366 step Time interval (days) between two dates [1] cuts Criteria for visibility selection (see Sect. 4) sort Parameter used to sort entries [RA] mime Mime type of the results votable — html — [pdf] from Name of the caller application within a graphical interface. 6.2. Target information at CDS The orbital elements of the asteroids and comets from the User identity (from). We encourage users to provide their astorb and cometpro catalogs hosted at CDS, used in the identity using the from parameter when they query ViSiON via computation of ephemerides, are displayed in a Web browser. its Web service (option not available in the query form). We For targets with fixed-coordinates, the SIMBAD entry corre- only keep this information for statistical purposes. It also al- sponding to their identifier is opened if resolved as a celestial lows to identify frequent users and their requests, in the case object. support is required. 6. Bridges to other services 6.3. Miriade at IMCCE In the case of moving targets, the coordinates reported in the To help further the observers in planning or conducting their summary table at the time of highest altitude may not be suffi- observations, ViSiON dynamically builds links to three VO ser- cient to point the targets. This is the case for particularly fast- vices for each target, provided in the last three columns of the moving targets, such as near-Earth asteroids (NEAs), or for in- summary table (see the example in Table 4). They provide, with struments with small fields of view. Ephemerides with a time a simple click, access to a view of the sky around the target in step of 10 min can therefore be displayed in a Web browser, Aladin (Sec. 6.1), information on the target at CDS (Sec. 6.2), computed on the fly at the IMCCE by the Miriade service and detailed ephemerides at IMCCE (Sec. 6.3). (Berthier et al., 2009). These ephemerides are dedicated to ob- servation, and local, horizontal, and equatorial coordinates are 6.1. Aladin at CDS listed for the entire night. An interactive view of the sky around the target position is loaded in an instance of the Aladin Sky Atlas (Bonnarel et al., 2000), host at the Centre de Donne´es astronomiques de Stras- 7. Conclusion bourg (CDS). By default, ViSiON requests the Digital Sky Sur- vey (DSS) to be displayed (Fig. 3), and overplots the known We have created a freely accessible service, named ViSiON, moving targets (at the time of passage at the maximum alti- compliant with Virtual Observatory standards, which allows tude) from the Skybot service (Berthier et al., 2006). Users anyone to create graphics and data tables of observing con- can then take advantage of all the tools offered by Aladin, such ditions of solar system and (extra)galactic objects. The ac- as overlaying additional catalogs and fields of view for instance, cess is standardized, and both command-line and Web form are around their target position. available. Results can be returned in human-readable format 5 Lainey, V., Dehant, V., Pa¨tzold, M., Apr. 2007. First numerical ephemerides of the Martian moons. A&A465, 1075–1084. Lainey, V., Duriez, L., Vienne, A., Jun. 2004b. New accurate ephemerides for the Galilean satellites of Jupiter. I. Numerical integration of elaborated equa- tions of motion. A&A420, 1171–1183. Laskar, J., Jacobson, R. A., Dec. 1987. GUST86 - an analytical ephemeris of the Uranian satellites. A&A188, 212–224. Le Guyader, C., May 1993. Solution of the N-Body Problem Expanded Into Taylor Series of High Orders. Applications to the Solar System Over Large Time Range. A&A272, 687–694. Rocher, P., Cavelier, C., 1996. Production d’elements orbitaux de cometes sur PC. In: Ferraz-Mello, S., Morando, B., Arlot, J.-E. (Eds.), Dynamics, Ephemerides, and Astrometry of the Solar System. Vol. 172 of IAU Sympo- sium. p. 357. Stoer, J., Bulirsch, R., 1980. Introduction to numerical analysis. Springer- Verlag New York Inc. Thuillot, W., Berthier, J., Vachier, F., Lainey, V., Arlot, J.-E., Aug. 2005. Vir- tual Observatory tools for Solar System bodies. In: AAS/Division for Plan- etary Sciences Meeting Abstracts #37. Vol. 37 of Bulletin of the American Astronomical Society. p. 638. Vienne, A., Duriez, L., May 1995. TASS1.6: Ephemerides of the major Satur- nian satellites. A&A297, 588. Wenger, M., Ochsenbein, F., Egret, D., Dubois, P., Bonnarel, F., Borde, S., Genova, F., Jasniewicz, G., Laloe¨, S., Lesteven, S., Monier, R., Apr. 2000. The SIMBAD astronomical database. The CDS reference database for as- tronomical objects. A&AS143, 9–22. Figure 3: Example of the Aladin instance launched by ViSiON. Here the target was the comet 67P/ChuryumovGerasimenko on June 4, 2016. (HTML, PDF) and VOTable for insertion in other programs, like telescope scheduler for instance. The results are tight with three other VO services, hosted at CDS and IMCCE to provide further information on targets and help in preparing the obser- vations. Acknowledgements This development of this service was partly funded by Eu- roplanet 2020 Research Infrastructure project of the European Union’s Horizon 2020 research and innovation programme un- der grant agreement No 654208. References Berthier, J., Hestroffer, D., Carry, B., Vachier, F., Lainey, V., Emelyanov, N. V., Thuillot, W., Arlot, J.-E., Ephemerides Service, I., Sep. 2009. Miriade: A Service for Solar Sytem Ojects Ephemerides in the VO Framework. In: Eu- ropean Planetary Science Congress 2009. p. 676. Berthier, J., Vachier, F., Thuillot, W., Fernique, P., Ochsenbein, F., Genova, F., Lainey, V., Arlot, J.-E., Jul. 2006. SkyBoT, a new VO service to identify Solar System objects. In: Gabriel, C., Arviset, C., Ponz, D., Enrique, S. (Eds.), Astronomical Data Analysis Software and Systems XV. Vol. 351 of Astronomical Society of the Pacific Conference Series. p. 367. Bonnarel, F., Fernique, P., Bienayme´, O., Egret, D., Genova, F., Louys, M., Ochsenbein, F., Wenger, M., Bartlett, J. G., Apr. 2000. The ALADIN inter- active sky atlas. A reference tool for identification of astronomical sources. Astronomy & Astrophysics Supplement 143, 33–40. Bowell, E., Muinonen, K. O., Wasserman, L. H., 1993. In: Asteroids, Comets, Meteors 1993. Vol. 810 of LPI Contributions. p. 44. Bulirsch, R., Stoer, J., 1966. Numerische Mathematik 8, 1. Carry, B., Berthier, J., Mar. 2016. A Simpler Procedure for Specifying Solar System Objects in Phase 2. The Messenger 163, 12–13. Fienga, A., Manche, H., Laskar, J., Gastineau, M., Verma, A., May 2014. IN- POP new release: INPOP13b. Scientific notes. Lainey, V., Arlot, J. E., Vienne, A., Nov. 2004a. New accurate ephemerides for the Galilean satellites of Jupiter. II. Fitting the observations. A&A427, 371–376. 6 7 d [ Target m φ D Alt Az. RA DEC Rate λ β r Δ α SEO MEO Links V G G h m h m s ◦ ′ ′′ Moon -7.03 1997 12 30 29 61 2 17 50 10 08 50 1772 335 -46 0.002 1.012 143.3 36.6 0.0 – h m h m s ◦ ′ ′′ 67P Churyumov-. . . 18.54 0.11 5 49 55 359 11 20 43 10 04 6 23 248 62 2.86 3.12 18.9 95.3 139.8 h m h m s ◦ ′ ′′ 12116 18.09 0.00 7 50 78 356 12 40 44 -13 37 24 8 299 48 2.46 3.12 15.8 122.6 162.3 h m h m s ◦ ′ ′′ TYC5597-234-1 11.85 0.00 10 50 72 2 15 46 30 -7 32 31 0 360 34 1.00 1.00 0.0 160.2 153.4 – h m h m s ◦ ′ ′′ Mars -1.97 18.56 11 00 86 348 15 40 47 -21 19 39 53 347 25 0.50 1.51 9.4 165.9 152.5 – h m h m s ◦ ′ ′′ M70 9.06 0.00 10 30 82 181 18 43 12 -32 17 31 0 363 -13 1.00 1.00 0.0 151.5 112.9 – h m h m s ◦ ′ ′′ (107) Camilla 13.07 0.10 8 30 74 2 19 57 9 -9 05 38 10 392 -19 2.96 3.71 11.8 131.4 95.0 h m h m s ◦ ′ ′′ Triton 13.52 0.12 5 40 73 359 22 54 2 -7 54 59 0 243 -56 29.93 29.96 1.9 90.3 53.8 – h m h m s ◦ ′ ′′ SMC 0.00 0.00 7 30 39 168 0 53 59 -72 00 0 0 302 -45 1.00 1.00 0.0 101.9 83.3 – – Table 4: Example of the summary table associated with Figs. 1 and 2. Values are reported at the time of the highest altitude. See Table 1 for a description of each column. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Astrophysics arXiv (Cornell University)

ViSiON: Visibility Service for Observing Nights

Carry, Benoit; Berthier, Jerome
Astrophysics , Volume 2018 (1806) – Jun 25, 2018
Download PDF
7 pages

Loading next page...
 
/lp/arxiv-cornell-university/vision-visibility-service-for-observing-nights-Qe8q8QCHgU

References (19)

ISSN
0032-0633
eISSN
ARCH-3330
DOI
10.1016/j.pss.2018.06.012
Publisher site
See Article on Publisher Site

Abstract

Preparation of detailed night schedule prior to an observing run can be tedious, especially for solar system objects which coordinates are epoch-dependent. We aim at providing the community with a Web service compliant with Virtual Observatory (VO) standards, to create tables of observing conditions, together with airmass and sky charts, for an arbitrary list of targets. We take advantage of available VO services such as the SIMBAD astronomical database, the Aladin sky atlas, and the Miriade ephemerides generator to build a new service dedicated to the planning of observations. The requests for ephemerides charts and tables are handled by a VO-compliant Web service. For each date, and each target, coordinates in local and equatorial frames are computed, and used to select targets accordingly to user’s criteria for visibility. This new service, dubbed ViSiON for Visibility Service for Observing Night, is a new method of Miriade Web service hosted at IMCCE. It allows anyone to create graphics of observing conditions and tables summarizing them, provided as PDF, VOTable, and xHTML documents. Keywords: ephemerides, virtual observatory tools, methods: numerical 1. Introduction standards, to produce tables and charts of visibility for both fixed coordinates targets and moving solar system objects. To prepare and conduct observations from a ground-based This service is dubbed Visibility Service for Observing Night telescope, many astronomers create airmass charts, also called (ViSiON). Targeted users are observers willing to prepare an visibility plots. These graphics represent the elevation above upcoming night, or to check the observability of some targets the horizon, or altitude, of each target during the night, provid- on-the-fly, and developers of night scheduler and telescope con- ing a graphical and easy way to schedule the sequence of ob- trol software willing to include airmass charts in graphical in- servations, in which each target should be observed as closely terfaces. as possible to its highest altitude. Another type of visual guide- The present article is structured as following: we quickly de- lines are fish-eye sky charts, in which the successive positions scribe the core of ephemeris computation in Sect. 2. The al- of the targets are displayed by tracks on a stereographic projec- gorithm of ViSiON is detailed in Sect. 3, and the criteria which tion of the local frame. It also helps optimizing observations, define observability are presented in Sect. 4. The different ways by avoiding pointing the telescope back and forth from different to access to the service are provided in Sect. 5 and its links with regions of the sky. other VO services in Sect. 6 If many astronomers have developed their own piece of code for that purpose, several efforts have been made to make 2. Ephemeris computation such tools freely available to the community. The staralt pages , by J. Me´ndez, allows productions of such graphics. The The ephemerides of planets and small Solar system ob- skycalc software , maintained by J. Thorstensen, represents jects are computed into a quasi-inertial reference frame, tak- an easily-portable and interactive solution to these points. Yet, ing into account the post-Newtonian approximations. The ge- these tools have been developed for objects with fixed coordi- ometric positions of the major planets and the Moon are pro- nates, and the eight planets. Planning observations of moving vided by INPOP planetary theory (Fienga et al., 2014). Those Solar System Objects (SSOs) may therefore result tedious: one of small SSOs (asteroids, comets, Centaurs, trans-neptunian should compute first the coordinates of each object, and then objects) are calculated by numerical integration of the N- enter them into the software interface. This often results im- body perturbed problem (Gragg-Bulirsch-Stoer algorithm, see practicable if the list of targets is long, or many dates are envi- Bulirsch and Stoer, 1966; Stoer and Bulirsch, 1980), using sioned. their latest published osculating elements, from the astorb We aim here at providing the community with a simple (Bowell et al., 1993) and cometpro (Rocher and Cavelier, Web service, fully compliant with Virtual Observatory (VO) 1996) databases. The overall accuracy of asteroid and comet http://catserver.ing.iac.es/staralt/ 2 3 http://www.dartmouth.edu/ physics/labs/skycalc/flyer.html http://vo.imcce.fr/webservices/miriade/?vision Preprint submitted to PSS October 9, 2018 arXiv:1806.09425v1 [astro-ph.IM] 25 Jun 2018 ephemerides provided by ViSiON is at the level of tens of milli- Table 1: List of parameters reported for each visible target, for each requested arcseconds, mainly depending on the accuracy of their osculat- night. The values reported in ViSiON correspond to the geometry at the time of ing elements. The positions of natural satellites are obtained the highest altitude on sky. thanks to dedicated solutions of their motion, e.g., Lainey et al. Quantity Description Units (2004a,b, 2007) for Mars and Jupiter, Vienne and Duriez Target Target designation – (1995) for Saturn, Laskar and Jacobson (1987) for Uranus, and m Apparent magnitude mag Le Guyader (1993) for Neptune satellites. φ Apparent diameter arcsec The ephemerides of stars, and more generally of all D Duration of visibility window h:m SIMBAD objects (Wenger et al., 2000), are also computed into a Alt Altitude deg quasi-inertial reference frame, taking into account their proper Az. Azimuth deg motion, parallax and radial velocity, as long as the corre- RA Right Ascension h:m:s sponding input parameters are provided by the catalogue. The DEC Declination d:m:s ephemeris accuracy of SIMBAD objects is thus at the level of the Rate Apparent non-sidereal motion arcsec/h accuracy of their catalogued position. λ Galactic longitude deg All the ephemerides in ViSiON are computed relatively to a β Galactic latitude deg given observer, in the ICRF reference frame (mean J2000 equa- r Range to observer au tor and equinox), and expressed in the UTC time scale. The Δ Heliocentric distance au coordinates of the observer (the topocenter) can be either pro- α Solar phase angle deg vided directly by users (longitude, latitude, altitude), or by us- SEO Solar elongation deg ing the observatory code provided by IAU Minor Planet Center MEO Moon elongation deg (MPC) for listed observatories. following parameters can be specified (see Table 2 for details): 3. Description of ViSiON’s algorithm altitude, apparent magnitude, solar phase angle, solar elonga- For each requested date, a series of computations are per- tion, lunar elongation, angular diameter, duration of visibil- formed prior to selection of targets suitable for observation. ity window, and events limiting the “night” (sunrise/sunset, or First, the ephemerides of the Sun are computed, to determine civil/nautical/astronomical twilights). the instants of sunset and sunrise. The limits of the civil, nau- ViSiON proceeds as following: for each target, it computes tical, and astronomical twilights, corresponding to the times its ephemerides for each observing date, as described above ◦ ◦ ◦ when the altitude of the Sun is 6 , 12 , and 18 below the hori- (Sect. 3). Each parameter (apparent magnitude, elevation, etc) zon, are also computed at this stage. is checked against the selection criteria at every 5 min time step. Second, the phase and position of the Moon throughout the A time step is considered valid only if the tests against all the night are computed and stored. Finally, the topocentric astro- criteria have passed. Finally, the target will be considered valid metric coordinates of each target are computed, with a time step for the date if the number of time steps exceeds the specified of 5 minutes, and expressed as equatorial (RA, Dec), horizontal minimum duration of visibility (30 min by default) between the (Az, Alt), and galactic coordinates (λ , β ). G G two events defining the “night”. Additionally, both solar and lunar elongations are computed For each date, all the targets that pass the selection are drawn at each time step, together with (for Solar system objects only) on the graphics, and included in the summary table. Their or- the solar phase angle, the heliocentric distance and the range to der can also be chosen by users, among all the parameters listed the observer. These detailed coordinates are used to determine in Table 1. A double sorting is even possible, with objects be- for how long each target is visible during the night. The entire ing grouped by type first (asteroid, comet, planet, satellite, and list of criteria defining observability is provided in Sec. 4. If fixed-coordinate) then by the requested parameter. there is at least one target observable, the night is included in the result, and these detailed coordinates are the inputs to the 5. Access to the service airmass and sky charts (see Figs. 1 and 2 for examples). For each date, a table lists all the targets that are visible (see Table 4 In 2005, the IMCCE started to implement Virtual Obser- for an example), and summarizes their observing conditions at vatory (VO) compliant interfaces to its ephemerides services the time of highest altitude (see Table 1 for a description of each (Thuillot et al., 2005). A Web portal describe the various parameter). services, such as solar system objects identification (Skybot, Berthier et al., 2006), or general ephemerides computation 4. Selection and sorting criteria (Miriade, Berthier et al., 2009). All our services are acces- sible via Web services (based on SOAP and HTTP POST verb) A suite of parameters are tested to determine if a target shall which allow anyone to interface her/his own software with the be considered as visible. The thresholds for selection are pro- services, via HTTP request and Web forms, and are integrated vided as inputs, accordingly to each observer’s preference. The 4 5 http://www.minorplanetcenter.net/iau/lists/ObsCodes.html http://vo.imcce.fr/ 2 Sunset Twilight Twilight Sunrise h m h m h m h m 21 59 23 21 09 58 11 20 h m h m h m h m h m h m h m h m LST 9 01 11 01 13 02 15 02 17 02 19 03 21 03 23 03 Mars M70 Camilla Triton TYC5597−234−1 67P Moon SMC 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 21:00 23:00 01:00 03:00 05:00 07:00 09:00 11:00 Time (UT) st Figure 1: Example of airmass chart produced by ViSiON for a request a La Silla observatory (ESO) on June the 1 , 2016. The different levels of gray shading on each side correspond to the day light, civil, nautical, and astronomical twilights. Both universal time (UT) and local sidereal time (LST) are reported on x-axes, and altitude and airmass are provided on y-axes. Each curve represents a target, color-coded as a function of type: planet (Mars), asteroid (12 116 and Camilla), comet (67P), satellite (Triton), and fixed-coordinates such as galaxy (Small Magellanic Cloud: SMC), star (TYC5597-234-1), or clusters (M70). in several VO-compliant software such as the well-spread Al- or “s:” suffixes, e.g., p:Phobos, s:Enceladus for the corre- adin Sky Atlas (Bonnarel et al., 2000). sponding satellites of Mars and Saturn. There are two ways to use the ViSiON web service: by writ- For objects with fixed coordinates, ViSiON queries the ting a client to send requests to the Miriade server and to re- SIMBAD service. The targets are therefore inherited from ceive and analyze the response (as for instance the ephemerides SIMBAD, which currently contains above 8 millions entries. tool for observations in Service Mode hosted by ESO, see These targets are submitted to ViSiON using the “u:” prefix, Carry and Berthier, 2016), or by using a command line inter- e.g., u:Sirius, u:HIP 87937, u:M 31 for Sirius, the Hippar- face and a data transfer program such as curl or wget. We also cos star 87937, and Messier 31. For objects without a catalog propose a simple query form on our SSO VO portal . In any entry, it is always possible to directly submit equatorial coor- case, the service accepts nine parameters, described in Table 3. dinates (in decimal hours and degrees in astrometric J2000) Only the list of targets is mandatory, all the other parameters to ViSiON, with the “e:” prefix. For instance, the equatorial being optional. We describe them in detail here. North pole can be queried via e:0.0+90.0. In all cases, ViSiON accepts aliases for each target, that can Targets (name). Planets can be requested by setting the name be used to simplify the labels on the graphics, with the syntax parameter to their name or number, preceeded by the “p:” pre- prefix:name=alias. For instance in Figs. 1 and 2, the Small fix, for instance p:Mars and p:5 to request visibility of Mars Magellanic Cloud was requested as e:0.9-72=SMC, while and Jupiter. Any entry in astorb and cometpro can be re- Mars was requested using p:4 and M70 with u:M 70. This trieved, using the “a:” and “c:” prefixes, e.g., a:Ceres, a:21, option may be useful for minor planets with high numbers and or c:67P for the asteroids (1) Ceres and (21) Lutetia, and the provisional designation: to avoid the display of 2000 WO137 comet 67P/Churyumov-Gerasimenko. All the satellites present on the graphics when querying for a:123456, one can force the in the references listed in Sect. 2 are available, though the “p:” usage of the IAU number by specifying a:123456=123456. http://vo.imcce.fr/webservices/miriade/?forms#Vision Altitude ( ) airmass 5.60 2.90 2.00 1.55 1.30 1.15 1.06 1.02 1.0 N 67P Moon E W TYC5597−234−1 Triton Camilla Mars M70 SMC 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 2016−06−01 st Figure 2: Example of sky chart produced by ViSiON for a request a La Silla observatory (ESO) on June the 1 , 2016. As in Fig. 1, each curve represents a target, color-coded as a function of its type. Topocenter (observer). The location on Earth of the observer and tag :{max:value } to set only one extrema, or max can be specified by providing the longitude (0,360 West), lat- tag :{min:value , max:value } to set both. For in- max i min itude (-90,+90 North), and altitude (meters above sea level) stance, to limit targets to sources brighter than V = 20, one or the IAU MPC observatory code, in which case ViSiON re- should specify mag:{max:20}. To select only objects matching trieves the coordinates of the observer from the MPC list. For all criteria during 2 hours between the sunrise and sunset, the instance, the Paris observatory can be specified as -2.336537 two parameters duration:120, event:sun must be speci- +48.8364 132.0 or 007. fied. Sorting options (sort). The entries in the summary tables (see Dates of computation (ep, nbd, step). The ep parameter Table 4 for an example) are sorted in increasing right ascen- sets the first (or only) date at which the visibility is sought sion by default. Changing this behavior is achieved by speci- (YYYY-MM-DD). If more than one night is desired, their total fying sort=key, in which key can take any value among the number can be specified with nbd (integer ≥1) and the num- reported parameters (Table 1). For instance, to sort object by ber of skipped nights between each by step (integer ≥0). For increasing apparent magnitude, one should specify sort=mv. It instance, four dates, one each week, starting from June, the st is also possible to first group targets by type before applying the 1 , 2016, these parameters should be set to ep=2016-06-01, sorting, using sort=type>key. nbd=4, and step=7. Output formats (mime). The results can be generated in three Selection criteria (cuts). Any number of selection crite- formats: VOTable, xHTML pages, and PDF file by specifying ria can be specified, using the syntax tag :value , ..., 1 1 the mime parameter. Most users will likely use the PDF and tag :value , where the tag are listed in Table 2. Tags N N i xHTML outputs, for convenient display and interaction. The accepting both a minimum and a maximum value shall VOTable output has been designed for a more advanced us- be set using the following syntax tag :{min:value } i min age, into a workflow for instance. It allows to directly query ViSiON from, e.g., a software to control observing queues, re- http://www.minorplanetcenter.net/iau/lists/ObsCodes.html trieving the values in the table, and the graphics for display 4 Table 2: Selection criteria in ViSiON. For each, we indicate if a minimum and/or a maximum value can be specified, and which are the default limits. The phase angle and apparent diameter selection only applies to solar system objects, and is disregarded for the other targets. Quantity Description Minimum Maximum Defaults Alt. Highest altitude reached!! [30 to 90 ] m Apparent magnitude!! [No default] α Solar phase angle!! [0 to 180 ] SEO Solar elongation! [No default] MEO Moon elongation! [No default] φ Apparent diameter!! [No default] D Duration of visibility window! [30 min] Event Start/End of the “night”! [civil] Table 3: List of input parameters to ViSiON (see Sect. 5 for the list of access points). For each, the default (between brackets) and range are specified. Parameter Description Range, or example values name Designation of the requested targets p:Mars, a:Pluto, u:M 31 observer Observer’s location (IAU code or geographic coordinates) [Paris, France] ep Starting epoch [now] — JD — YYYY-MM-DD nbd Number of date of computation [1] ≤ nbd < 366 step Time interval (days) between two dates [1] cuts Criteria for visibility selection (see Sect. 4) sort Parameter used to sort entries [RA] mime Mime type of the results votable — html — [pdf] from Name of the caller application within a graphical interface. 6.2. Target information at CDS The orbital elements of the asteroids and comets from the User identity (from). We encourage users to provide their astorb and cometpro catalogs hosted at CDS, used in the identity using the from parameter when they query ViSiON via computation of ephemerides, are displayed in a Web browser. its Web service (option not available in the query form). We For targets with fixed-coordinates, the SIMBAD entry corre- only keep this information for statistical purposes. It also al- sponding to their identifier is opened if resolved as a celestial lows to identify frequent users and their requests, in the case object. support is required. 6. Bridges to other services 6.3. Miriade at IMCCE In the case of moving targets, the coordinates reported in the To help further the observers in planning or conducting their summary table at the time of highest altitude may not be suffi- observations, ViSiON dynamically builds links to three VO ser- cient to point the targets. This is the case for particularly fast- vices for each target, provided in the last three columns of the moving targets, such as near-Earth asteroids (NEAs), or for in- summary table (see the example in Table 4). They provide, with struments with small fields of view. Ephemerides with a time a simple click, access to a view of the sky around the target in step of 10 min can therefore be displayed in a Web browser, Aladin (Sec. 6.1), information on the target at CDS (Sec. 6.2), computed on the fly at the IMCCE by the Miriade service and detailed ephemerides at IMCCE (Sec. 6.3). (Berthier et al., 2009). These ephemerides are dedicated to ob- servation, and local, horizontal, and equatorial coordinates are 6.1. Aladin at CDS listed for the entire night. An interactive view of the sky around the target position is loaded in an instance of the Aladin Sky Atlas (Bonnarel et al., 2000), host at the Centre de Donne´es astronomiques de Stras- 7. Conclusion bourg (CDS). By default, ViSiON requests the Digital Sky Sur- vey (DSS) to be displayed (Fig. 3), and overplots the known We have created a freely accessible service, named ViSiON, moving targets (at the time of passage at the maximum alti- compliant with Virtual Observatory standards, which allows tude) from the Skybot service (Berthier et al., 2006). Users anyone to create graphics and data tables of observing con- can then take advantage of all the tools offered by Aladin, such ditions of solar system and (extra)galactic objects. The ac- as overlaying additional catalogs and fields of view for instance, cess is standardized, and both command-line and Web form are around their target position. available. Results can be returned in human-readable format 5 Lainey, V., Dehant, V., Pa¨tzold, M., Apr. 2007. First numerical ephemerides of the Martian moons. A&A465, 1075–1084. Lainey, V., Duriez, L., Vienne, A., Jun. 2004b. New accurate ephemerides for the Galilean satellites of Jupiter. I. Numerical integration of elaborated equa- tions of motion. A&A420, 1171–1183. Laskar, J., Jacobson, R. A., Dec. 1987. GUST86 - an analytical ephemeris of the Uranian satellites. A&A188, 212–224. Le Guyader, C., May 1993. Solution of the N-Body Problem Expanded Into Taylor Series of High Orders. Applications to the Solar System Over Large Time Range. A&A272, 687–694. Rocher, P., Cavelier, C., 1996. Production d’elements orbitaux de cometes sur PC. In: Ferraz-Mello, S., Morando, B., Arlot, J.-E. (Eds.), Dynamics, Ephemerides, and Astrometry of the Solar System. Vol. 172 of IAU Sympo- sium. p. 357. Stoer, J., Bulirsch, R., 1980. Introduction to numerical analysis. Springer- Verlag New York Inc. Thuillot, W., Berthier, J., Vachier, F., Lainey, V., Arlot, J.-E., Aug. 2005. Vir- tual Observatory tools for Solar System bodies. In: AAS/Division for Plan- etary Sciences Meeting Abstracts #37. Vol. 37 of Bulletin of the American Astronomical Society. p. 638. Vienne, A., Duriez, L., May 1995. TASS1.6: Ephemerides of the major Satur- nian satellites. A&A297, 588. Wenger, M., Ochsenbein, F., Egret, D., Dubois, P., Bonnarel, F., Borde, S., Genova, F., Jasniewicz, G., Laloe¨, S., Lesteven, S., Monier, R., Apr. 2000. The SIMBAD astronomical database. The CDS reference database for as- tronomical objects. A&AS143, 9–22. Figure 3: Example of the Aladin instance launched by ViSiON. Here the target was the comet 67P/ChuryumovGerasimenko on June 4, 2016. (HTML, PDF) and VOTable for insertion in other programs, like telescope scheduler for instance. The results are tight with three other VO services, hosted at CDS and IMCCE to provide further information on targets and help in preparing the obser- vations. Acknowledgements This development of this service was partly funded by Eu- roplanet 2020 Research Infrastructure project of the European Union’s Horizon 2020 research and innovation programme un- der grant agreement No 654208. References Berthier, J., Hestroffer, D., Carry, B., Vachier, F., Lainey, V., Emelyanov, N. V., Thuillot, W., Arlot, J.-E., Ephemerides Service, I., Sep. 2009. Miriade: A Service for Solar Sytem Ojects Ephemerides in the VO Framework. In: Eu- ropean Planetary Science Congress 2009. p. 676. Berthier, J., Vachier, F., Thuillot, W., Fernique, P., Ochsenbein, F., Genova, F., Lainey, V., Arlot, J.-E., Jul. 2006. SkyBoT, a new VO service to identify Solar System objects. In: Gabriel, C., Arviset, C., Ponz, D., Enrique, S. (Eds.), Astronomical Data Analysis Software and Systems XV. Vol. 351 of Astronomical Society of the Pacific Conference Series. p. 367. Bonnarel, F., Fernique, P., Bienayme´, O., Egret, D., Genova, F., Louys, M., Ochsenbein, F., Wenger, M., Bartlett, J. G., Apr. 2000. The ALADIN inter- active sky atlas. A reference tool for identification of astronomical sources. Astronomy & Astrophysics Supplement 143, 33–40. Bowell, E., Muinonen, K. O., Wasserman, L. H., 1993. In: Asteroids, Comets, Meteors 1993. Vol. 810 of LPI Contributions. p. 44. Bulirsch, R., Stoer, J., 1966. Numerische Mathematik 8, 1. Carry, B., Berthier, J., Mar. 2016. A Simpler Procedure for Specifying Solar System Objects in Phase 2. The Messenger 163, 12–13. Fienga, A., Manche, H., Laskar, J., Gastineau, M., Verma, A., May 2014. IN- POP new release: INPOP13b. Scientific notes. Lainey, V., Arlot, J. E., Vienne, A., Nov. 2004a. New accurate ephemerides for the Galilean satellites of Jupiter. II. Fitting the observations. A&A427, 371–376. 6 7 d [ Target m φ D Alt Az. RA DEC Rate λ β r Δ α SEO MEO Links V G G h m h m s ◦ ′ ′′ Moon -7.03 1997 12 30 29 61 2 17 50 10 08 50 1772 335 -46 0.002 1.012 143.3 36.6 0.0 – h m h m s ◦ ′ ′′ 67P Churyumov-. . . 18.54 0.11 5 49 55 359 11 20 43 10 04 6 23 248 62 2.86 3.12 18.9 95.3 139.8 h m h m s ◦ ′ ′′ 12116 18.09 0.00 7 50 78 356 12 40 44 -13 37 24 8 299 48 2.46 3.12 15.8 122.6 162.3 h m h m s ◦ ′ ′′ TYC5597-234-1 11.85 0.00 10 50 72 2 15 46 30 -7 32 31 0 360 34 1.00 1.00 0.0 160.2 153.4 – h m h m s ◦ ′ ′′ Mars -1.97 18.56 11 00 86 348 15 40 47 -21 19 39 53 347 25 0.50 1.51 9.4 165.9 152.5 – h m h m s ◦ ′ ′′ M70 9.06 0.00 10 30 82 181 18 43 12 -32 17 31 0 363 -13 1.00 1.00 0.0 151.5 112.9 – h m h m s ◦ ′ ′′ (107) Camilla 13.07 0.10 8 30 74 2 19 57 9 -9 05 38 10 392 -19 2.96 3.71 11.8 131.4 95.0 h m h m s ◦ ′ ′′ Triton 13.52 0.12 5 40 73 359 22 54 2 -7 54 59 0 243 -56 29.93 29.96 1.9 90.3 53.8 – h m h m s ◦ ′ ′′ SMC 0.00 0.00 7 30 39 168 0 53 59 -72 00 0 0 302 -45 1.00 1.00 0.0 101.9 83.3 – – Table 4: Example of the summary table associated with Figs. 1 and 2. Values are reported at the time of the highest altitude. See Table 1 for a description of each column.

Journal

Astrophysics – arXiv (Cornell University)

Published: Jun 25, 2018

There are no references for this article.