arxiv: v1 [astro-ph.im] 13 Jan 2010

Transcription

1 Light pollution at high zenith angles, as measured at Cerro Tololo Inter-American Observatory 1 Kevin Krisciunas, 2 Hector Bogglio, 3 Pedro Sanhueza, 3 and Malcolm G. Smith 4 arxiv: v1 [astro-ph.im] 13 Jan 2010 ABSTRACT On the basis of measurements of the V-band sky brightness obtained at Cerro Tololo Inter-American Observatory in December 2006 and December 2008 we confirm the functional form of the basic model of Garstang (1989, 1991). At high zenith angles we measure an enhancement of a factor of two over Garstang s later model when there is no marine cloud layer over La Serena/Coquimbo. No corresponding enhancement is found in the B-band. Subject headings: Astronomical Phenomena and Seeing 1. Introduction In a previous paper (Krisciunas et al. 2007, hereafter K07) we presented measurements of the BV RI sky brightness at Cerro Tololo Inter-American Observatory, obtained with the CTIO 0.9-m, 1.3-m, and 1.5-m telescopes and CCD detectors. As a sanity check we also obtained B- and V-band data with a 15-cm reflector and single channel photometer previously used by Krisciunas (1997) for an 11-year series of sky brightness measurements obtained at the 2800-m level of Mauna Kea. We found that at a zenith angle of 80 deg we could measure the effect of artificial lights in La Serena some 55 km away. The excess V-band light amounted to an enhancement of 44 to 72 percent above the model of Garstang (1989). The December 2006 single channel B- and V-band measuremens are included in the analysis of the present paper. 1 Based in part on observations taken at the Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. 2 George P. and Cynthia Woods Mitchell Institute for Fundamental Physics & Astronomy, Texas A. & M. University, Department of Physics, 4242 TAMU, College Station, TX 77843; krisciunas@physics.tamu.edu 3 Oficina de Protección de la Calidad del Cielo del Norte de Chile (OPCC), 1606 Cisternas, La Serena, Chile; hbogglio@opcc.cl, psanhueza@opcc.cl 4 Cerro Tololo Inter-American Observatory, Casilla 603, La Serena, Chile; msmith@ctio.noao.edu

2 2 In December of 2008 we obtained more data with the single channel system of Krisciunas (1996). Our purpose was twofold: 1) to confirm the functional form of the basic model of Garstangbytakingdataatthefullrangeofzenithangles; and2)todetermineabaselinelevel of light pollution measurable from Cerro Tololo. As La Serena and Coquimbo s population continues to grow, we can use these measurements to determine the effectiveness of the lighting control measures in place. 2. The Data Our more and more antiquated, but still functional, photometric system consists of a 15-cm f/5.9 Newtonian reflector and a very light photometer whose light sensitive element is an uncooled RCA 931A photomultiplier (pm) tube effectively the same as a 1P21 tube. The V-band filter is a 2 mm thick piece of Schott GG495 glass, which has a short wavelength cut on at about 4800 Å. At 5500 Å it reaches a maximum transmission of 91 percent, which it retains to red end of the visible portion of the electromagnetic spectrum. The effective V-band filter transmission results from the combination of this transmission curve and the quantum efficiency (QE) of the pm tube, which falls to near zero at 7000 Å(Lallemand 1962, Fig. 1). The B-band filter is composed of a 2 mm thick piece of Schott GG385 glass and a 1 mm thick piece of Schott BG-12 glass. Together they give a filter transmission which cuts on at 3600 Å, has peak transmission at 4000 Å, diminishes to zero shortly after 5000 Å, but then has a 15 percent red leak starting at 7000 Å. Given the near zero QE of the tube at the red end, the red leak would have no significant effect on the B-band measurements. From 28 to 30 December 2008 UT our typical observing strategy was as follows: 1) start and end the night with measurements of standard stars chosen from the Bright Star Catalogue (Hoffleit & Jaschek 1982) so that the photometric zeropoints are interpolated and never extrapolated; 2) measure the background sky near the zenith, then very low in the sky, then at an intermediate zenith angle; 3) measure a standard star every 40 minutes or so. Our principal standard stars were BS 1179 and ζ Cae, with ρ For observed less often. The first two stars could be observed on the same gain setting used for the sky brightness measures. This eliminates one source of potential systematic error. The output of the pm tube is amplified by a factor up to 10 6 with an analog amplifier. The output is displayed real-time on a milliammeter and recorded with a Hewlett Packard strip chart recorder. To determine the photometric zeropoints each night we assumed mean extinction coefficients of k v = mag/airmass and k b = mag/airmass, based on several years of CCD photometry at CTIO. Since our standards were observed between 1.0 and 1.1 airmasses, the adopted extinction values hardly mattered, but on any given night the derived

3 3 sky brightness values are systematically in error by the difference of the true extinctions and the assumed values. Of greater consequence is the drift of the zeropoints with time. We found that the drift was linear and amounted to a maximum of 0.03 mag/hour. The sky brightness measures presented here take into account these linear drifts. It should be noted that extinction corrections are not made to sky brightness measures, only to measures of standard stars. Once we determine the equivalent B- and V-band magnitudes of the sky in our square arc minute beam, we obtain the number of magnitudes per square arc second by adding to the single-beam magnitudes. This additive factor is just 2.5 log 10 A, where A is the beam size in square arc seconds. Almost all of our high zenith angle measurements made in December 2006 and December 2008 were taken over La Serena. On nights when there is no marine cloud layer, the city is clearly visible down below. In Fig. 1 we show a photo of this situation, obtained by one of us (MGS). The sky-brightness data are found in Table 1. For the V-band we also give the flux in nanolamberts, using Eq. 3 of K07 or its equivalent (Garstang 1989, Eq. 27). Fluxes in nanolamberts are technically only relevant for the V-band, as it approximates the response of the human eye. We can convert the B-band sky brightness in mag/sec 2 to the flux b B in photons/cm 2 /sec/steradian by using the inverse of Eq. 39 of Garstang (1989), namely: [ ] B b B = exp The observed flux of the sky near the zenith can effectively be corrected to the zenith by dividing it by the secant of the zenith angle. This is basically the same as Eq. 19 of Schaefer (1990): where Z rad is the zenith angle in radians. B zen = B obs /(1 + Z 2 rad/2), 3. Discussion The mean sky brightness at the zenith in December 2008 was V = mag/sec 2. An individual V-band measurement at that flux level was accurate to ± mag/sec 2.

4 4 This is based on 17 observations within 27 degrees of the zenith and corrected to the zenith using the formula of Schaefer (1990) given above. The mean zenith V-band flux was nl, with a 1-σ scatter of ± 3.11 nl (5.5 percent). The mean B-band zenith sky brightness was mag/sec 2. An individual B-band measurement at that flux level was accurate to ± 0.10 mag/sec 2. These mean values are amongst the faintest levels we have ever measured, indicating that the solar minimum was still in effect in December The individual V-band data points from December 2006 and December 2008 are shown in Fig. 2. We have included the locus of Garstang (1989) derived for Mt. Graham, offset by +3.0 nl, and the locus of Garstang (1991) for Junipero Serra (California), similarly offset. Since we did not correct our measurements for the presence of faint stars in the beam, some offset is understandable. Garstang s improved 1991 model includes the addition of an ozone layer, a more accurate representation of the atmospheric molecular density variation as a function of height, a better mathematical representation of the scattering angular function of aerosols, and a thin layer of dust of arbitrary optical thickness and height above sea level. These modifications lead to a reduction of the predicted brightness of the night sky. The ratio of our V-band data to the slightly offset locus of Garstang (1991) is shown in Fig. 3. This shows that except at high zenith angle we do not measure any significant enhancement of the V-band sky brightness at CTIO, even when we can clearly see La Serena down below. For comparison the reader is directed to the interesting article by Luginbuhl, Walker, & Wainscoat (2009), which includes pictorial and graphical data along with modeling based on the Garstang(1991) model and an approximate correction for double scattering. OurB-banddata, converted tophotons/cm 2 /sec/steradian, areshown infig. 4. Unlike our V-band data, there is no strong enhancement of the B-band light at high zenith angle when the lights of La Serena/Coquimbo are visible to the naked eye. This may be because the artificial lights presently in operation do not put out much light in the B-band. Spectra of the scattered light should be obtained to check this. There is one curious anomaly in Fig. 2. Our final high zenith angle measurements on 30 December 2008 were obtained over the town of Andacollo a much more southerly azimuth than La Serena. Though Andacollo was clearly visible (i.e. not covered by cloud), we measured no enhancement of the V-band flux in this direction compared to the model of Garstang (1991). We surmise that this is due to the smaller population in Andacollo and a smaller amount of airborne dust there. A combination of dust and light emitted above the horizontal gives us enhanced light levels over La Serena and Coquimbo when they are visible at Cerro Tololo.

5 5 As the population of La Serena and Coquimbo continues to grow, the observations presented here may serve as a reference. If the lighting ordinances in place are effective, we may see only minimally increased light pollution at CTIO. We thank Oscar Saa for technical support at the mountain. Garstang, R. H. 1989, PASP, 101, 306 Garstang, R. H. 1991, PASP, 103, 1109 REFERENCES Hoffleit, D., & Jaschek, C. 1982, The Bright Star Catalogue, 4th ed. (New Haven, Yale Univ. Obs.) Krisciunas, K. 1996, S&T, 91, No. 2, 91 Krisciunas, K. 1997, PASP, 109, 1181 Krisciunas, K., Semler, D. R., Richards, J., Schwarz, H. E., Suntzeff, N. B., Vera, S., & Sanhueza, P. 2007, PASP, 119, 687 (K07) Lallemand, A. 1962, in Astronomical Techniques, W. A. Hiltner, ed., Chicago and London: Univ. of Chicago Press, 126 Luginbuhl, C. B., Walker, C. E., & Wainscoat, R. J. 2009, Physics Today, 62, No. 12, 32 Schaefer, B. E. 1990, PASP, 102, 212 This preprint was prepared with the AAS L A TEX macros v5.2.

6 6 Table 1. CTIO Sky Brightness Values a UT Date UT Azimuth Zenith Angle Filter Sky Brightness Flux (nl) Dec 28 1: V Dec 28 2: V Dec 28 2: V Dec 28 3: V Dec 28 3: V Dec 28 3: V Dec 28 4: V Dec 28 4: V Dec 28 4: V Dec 28 4: V Dec 28 4: V Dec 28 1: B Dec 28 2: B Dec 28 2: B Dec 28 3: B Dec 28 3: B Dec 28 3: B Dec 28 4: B Dec 28 4: B Dec 28 4: B Dec 28 4: B Dec 28 4: B Dec 29 3: V Dec 29 3: V Dec 29 3: V Dec 29 3: V Dec 29 3: V Dec 29 4: V Dec 29 4: V

7 7 Table 1 Continued UT Date UT Azimuth Zenith Angle Filter Sky Brightness Flux (nl) Dec 29 4: V Dec 29 4: V Dec 29 4: V Dec 29 3: B Dec 29 3: B Dec 29 3: B Dec 29 3: B Dec 29 3: B Dec 29 4: B Dec 29 4: B Dec 29 4: B Dec 29 4: B Dec 29 4: B Dec 30 2: V Dec 30 2: V Dec 30 2: V Dec 30 2: V Dec 30 2: V Dec 30 3: V Dec 30 3: V Dec 30 3: V Dec 30 3: V Dec 30 4: V Dec 30 4: V Dec 30 4: V Dec 30 4: V Dec 30 4: V Dec 30 4: V Dec 30 5: V

8 8 Table 1 Continued UT Date UT Azimuth Zenith Angle Filter Sky Brightness Flux (nl) Dec 30 2: B Dec 30 2: B Dec 30 2: B Dec 30 2: B Dec 30 3: B Dec 30 3: B Dec 30 3: B Dec 30 3: B Dec 30 3: B Dec 30 3: B Dec 30 4: B Dec 30 4: B Dec 30 4: B Dec 30 4: B Dec 30 4: B Dec 30 4: B a Year is UT is in hours and minutes. Azimuth and zenith angle are in degrees, and are probably accurate to ± 0.2 deg. Sky brightness is measured in magnitudes per square arc second. For V-band values the sky brightness is also converted to flux in nanolamberts.

9 9 Fig minute exposure over La Serena and Coquimbo, 28 December 2008 UT, taken just south of the #1 40-cm dome at CTIO. Fig. 2. Observed V-band sky brightness at Cerro Tololo during solar minimum and at a variety of zenith angles. The upper (dashed) locus is based on the model of Garstang (1989) for Mt. Graham. The lower (solid) locus is based on the model of Garstang (1991) for Junipero Serra, California. Fig. 3. Ratio of observed V-band sky brightness and Garstang (1991) model shown in Fig. 2. Essentially, only at high zenith angles do we find observed sky brightness in excess of Garstang s model. Fig. 4. Observed B-band sky brightness at Cerro Tololo, converted to flux in cgs units. A third order polynomial fit is also shown.

10 10 1 See: arxiv did not like the giant size of the eps conversion Krisciunas et al. Fig. 1

11 Dec 2006 Observed V-band sky brightness (nl) Dec Dec Dec Dec 2008 Garstang (1989) model Garstang (1991) model V-band sky brightness at Cerro Tololo over La Serena over Andacollo Zenith angle (degrees) Krisciunas et al. Fig. 2

12 V-band flux ratio (observed/model) Dec Dec Dec Dec Dec Zenith angle (degrees) Krisciunas et al. Fig. 3

13 13 Observed B-band sky brightness (photons/cm 2 /sec/steradian) 1e+08 9e+07 8e+07 7e+07 6e+07 5e+07 4e+07 3e+07 2e+07 1e Dec Dec Dec Dec Dec 2008 B-band sky brightness at Cerro Tololo Zenith angle (degrees) Krisciunas et al. Fig. 4