Parque Astron´omicode Atacama: An Ideal Site for Millimeter, Submillimeter, and Mid-Infrared Astronomy

R. Bustos1, M. Rubio2,3, A. Ot´arola4,5, and N. Nagar6

ABSTRACT The area of Chajnantor, at more than 5000 m altitude in northern , offers unique at- mospheric and operational conditions which arguably make it the best site in the world for millimeter, submillimeter, and mid-infrared observatories. Long-term monitoring of the precip- itable water vapor (PWV) column on the Chajnantor plateau has shown its extreme dryness with annual median values of 1.1 mm. Simultaneous measurements of PWV on the Chajnantor plateau (5050 m) and on Cerro Chajnantor (5612 m) show that the latter is around 36% lower under normal atmospheric conditions and up to 80% lower than the plateau in the presence of temperature inversion layers. Recently, the Government of Chile has consolidated the creation of the Parque Astron´omicode Atacama (Atacama Astronomical Park), an initiative of the Chilean Commission for Science and Technology (CONICYT). This new park offers an opportunity for national and international projects to settle within its boundaries, gain access to an extremely dry site that is suitable for observations over a broad spectral range, especially in the millimeter to mid-infrared wavelengths, and benefit from operational and logistical support within a secure legal framework. Subject headings: astronomical phenomena and seeing

1. Introduction cal observatory, a high altitude (to reduce the at- mospheric column) and/or extremely dry (low wa- Astronomical objects emit electromagnetic ter vapor column) site is desirable. The rare com- waves over a wide range of frequencies, or wave- bination of both conditions creates an ideal site for lengths. Our atmosphere, mainly due to the ab- millimeter, submillimeter, and mid-infrared ob- sorption of energy by atmospheric water vapor, servations from the ground. The Chajnantor area partially or totally blocks sections of the electro- in Chile (23◦ 1.40 S, 67◦ 45.20 W), at over 5000 m magnetic spectrum that carry information about above sea level on the western side of the the cosmos. To mitigate this natural blocking, mountains and about 50 km east of San Pedro de which is suffered by any ground-based astronomi- Atacama, provides these particular conditions. Since 1997, several state-of-the-art instruments, 1Facultad de Ingenier´ıa, Universidad Cat´olica de la arXiv:1410.2451v2 [astro-ph.IM] 30 Mar 2015 Sant´ısimaConcepci´on,Alonso de Ribera 2850, Concepci´on, covering different wavelength regimes, have settled Chile; [email protected]. at various sites near the Chajnantor area. Some 2Departamento de Astronom´ıa, Universidad de Chile, have settled on the Chajnantor plateau (5050 m): Casilla 36-D, Santiago, Chile the Cosmic Background Imager (CBI; Padin et 3Programa Astronom´ıa,CONICYT, Moneda 1375, Piso al. (2002)), the Atacama Pathfinder Experiment 7, Santiago, Chile (APEX; G¨ustenet al. (2006)), the Q/U Imag- 4TMT Observatory, 1111 South Arroyo Parkway, Pasadena, CA 91105, USA ing Experiment (QUIET; QUIET collaboration 5Atmospheric Sciences, University of Arizona, 1118 E (2013)), and the Atacama Large Millimeter Array 4th St, Tucson, Arizona 85721, USA (ALMA; Wootten & Thompson (2009)). Others, 6Departamento de Astronom´ıa, Universidad de Con- over various time spans, have settled on a plateau cepci´on,Casilla 160-C, Concepci´on,Chile (5200 m) on the nearby Cerro Toco (5604 m): the

1 Fig. 1.— Panoramic view of the Chajnantor area taken from Cerro Honar in 2007. To the left, Volc´an (5920 m). To the center-right, Cerro Chajnantor (5640 m). Cerro Toco (5604 m) is in between. The white dots on Llano de Chajnantor are the CBI and ALMA (under construction). To the right, the Pampa La Bola plateau (4865 m). Photo credit: ALMA (ESO, NAOJ, NRAO).

Mobile Anisotropy Telescope (MAT; Torbet et Cerro Chajnantor is almost 600 m higher than the al. (1999)), the Q-band Mobile Anisotropy Probe plateau, and thus has a significantly lower atmo- (QMAP; Miller et al. (2002)), the Millimeter-wave spheric moisture column. Both Cerro Toco and Interferometer (MINT; Fowler et al. (2005)), the Cerro Chajnantor, as well as their immediate sur- Atacama Cosmology Telescope (ACT; Swetz et roundings, are inside the boundaries of the Parque al. (2011)), the Radiative Heating in Underex- Astron´omicode Atacama. plored Bands Campaign (RHUBC II; Turner et The extremely dry conditions of the Atacama al. (2012)), Polarbear (Kermish et al. 2012), and Desert result from a combination of several effects, the Atacama B-Mode Search (ABS; Kusaka et al. of which the primary is its subtropical location, (2014)). The plateau (5525 m) on Volc´anSaire- which is right under the descending branch of the cabur hosted the Receiver Lab Telescope (RLT; Hadley cell of dry and warm air, whose subsidence Marrone et al. (2004)). On the nearby Pampa La maintains the southeast Pacific anticyclone. This Bola (4865 m): the Atacama Submillimeter Tele- high pressure keeps the middle and lower tropo- scope (ASTE; Kohno (2005)) and the NANTEN2 sphere warm and provides conditions that explain Observatory (Kawamura et al. 2005). On Cerro both the aridity and dryness of this subtropical Chajnantor (5640 m), a 1 m telescope (mini-TAO) region (Garreaud 2011). The subsiding air that was installed as a pathfinder for a future 6 m mid- creates this high surface pressure is also responsi- infrared telescope (Sako et al. 2008). New projects ble for the strong temperature inversion layer that such as the Cosmology Large Angular Scale Sur- gives origin to the southeast pacific maritime cloud veyor (CLASS), High Altitude Terahertz Solar deck (Mu˜nozet al. 2010). The coastal mountain Telescope (HATS), CCAT, Tokyo Atacama Ob- range on the west side of Chile, at over 2000 m al- servatory (TAO), and Chilean Medidor Aut´onomo titude, combined with the temperature inversion de Radio Interferencia (MARI) are planning to use layer, block the moist air flow to the interior of the Chajnantor area in the near future. this region. To the east side, the Andes mountain Based on transparency measurements (Radford range blocks the humidity flow from the Amazon & Holdaway 1998) and the number of installed and basin, except in Austral summer, the so-called Al- planned telescopes, the Chajnantor area has be- tiplanic Winter (Bustos 2001). These effects, in come a leading global site for millimeter and sub- combination with high elevation plateaus on the millimeter astronomy. Similarly, based on the first western side of the Andes, produce ideal sites to observations by mini-TAO, the area has proven to scientifically exploit the millimeter, submillimeter, be excellent for observations in the mid-infrared and mid-infrared spectral windows. regime (Sako et al. 2008; Motohara et al. 2008). Various methods are used to characterize the Figure 1 shows a panoramic view of the Chaj- atmospheric transmission at a given site. Most are nantor plateau and the nearby mountains, looking based on passive radiometry to monitor the sky north from Cerro Honar in 2007. The summit of brightness within a given spectral band. On Cha-

2 PWV: Llano Chajnantor (5050 m) & Cerro Chajnantor (5612 m) 10 Cerro Chajnantor Llano Chajnantor 8

6

PWV [mm] 4

2

0 Dec 18 Dec 19 Dec 20 Dec 21 Dec 22 Dec 23

PWV decrease between Llano and Cerro Chajnantor in % 100

80

60 [%] 40 36.3 %

20

0 Dec 18 Dec 19 Dec 20 Dec 21 Dec 22 Dec 23 2011

Fig. 2.— Top: A comparison of the 1-minute averages of the simultaneous PWV measurements at Llano de Chajnantor (red) and Cerro Chajnantor (blue) over a 5 day run. Bottom: The drop in PWV in percent between Llano and Cerro Chajnantor. The dotted line shows the 36.3% decrease in PWV. jnantor, the atmospheric transparency at 225 GHz decay in altitude at a rate given by the scale height was monitored as part of the ALMA site-testing h0 (Ot´arolaet al. 2010), as seen in equation 1. effort using a 225 GHz tipping radiometer (McK- − h innon 1988). The Harvard-Smithsonian Cen- PWV (h) = PWV0 · e h0 (1) ter for Astrophysics (Paine et al. 2000) and the Nobeyama Radio Observatory (Matsushita et al. Here, PWV0 (in millimeters) is the correspond- 1999) monitored the sky transparency using a ing PWV at the surface level, and the scale height Fourier transform spectrometer from 200 GHz up h0 (in meters) corresponds to the e-folding alti- to 3.5 THz and 1.5 THz, respectively. Another tude in the exponential distribution of water va- method to characterize the atmospheric trans- por with altitude. Deviations from the model of mission is to combine measurements of the in- equation 1 are due to horizontal advection of mois- tegrated water vapor of the atmospheric column ture at varying altitudes in the atmosphere and/or with a physical model of the atmospheric absorp- changes in water vapor density mainly due to the tion spectrum as a function of atmospheric mois- presence of temperature inversion layers in the tro- ture. This method was followed by the ALMA posphere (Roosen & Angione 1977; Ot´arolaet al. site-testing team using two passive radiometers to 2010). monitor the intensity of the 183 GHz water vapor For the Chajnantor area, Giovanelli et al. emission line (Delgado et al. 1999). Considering (2001) determined a median h0 of 1135 m from that water vapor pressure is limited by temper- balloon measurements and Bustos et al. (2000) ature through the Clausius-Clayperon equation, estimated h0 ∼ 1500 m from a reanalysis of his- the precipitable water vapor (PWV) content at torical data. Neither of these methods measured any given height is modeled with an exponential h0 continuously. For h0 = 1135 m (1500 m), the expected drop in PWV between Llano de Chaj-

3

00 h nantor (5050 m) and Cerro Chajnantor (5612 m) 100 % 22 h 02 h is 39.1% (31.2%). Here, we present preliminary 80 %

results on simultaneous PWV measurements, di- 60 % 20 h 04 h rectly obtaining the difference in PWV between 40 % both sites. NIGHT−TIME 20 %

2. PWV observations: Llano versus Cerro 18 h 06 h Chajnantor DAY−TIME

In this experiment, we used two radiometers 16 h 08 h Dec 18 that measured the intensity of the 183 GHz water Dec 19 vapor line. One of these was originally built and Dec 20 14 h 10 h Dec 21 used for the ALMA site-testing campaign (Del- Dec 22 12 h gado et al. 1999). This unit was refurbished by the Laboratory of Radio Astronomy at Universi- Fig. 3.— 24 h polar plot of the drop in PWV dad de Concepci´onin 2009 with the installation between the Llano and Cerro Chajnantor for each of a new backend, replacement of the tipping mir- of the 5 days. Nighttime (23:00-9:30 UTC) and ror motors, and the implementation of new soft- daytime are separated by the solid black line. The ware for observing, control, and communications black circle is the 36.3% level. (Bustos & Nagar 2011). The second 183 GHz ra- diometer used in the experiment is the one that is mounted in the cabin of the APEX radio tele- tor is around 36.3% (dotted line). This ∼36.3% scope located on Llano de Chajnantor, whose data difference is more often seen in the afternoons, is available on the web1. since the atmosphere is generally well-mixed and The refurbished radiometer was first installed no significant inversion layers are expected. On next to APEX for a 10-month, side-by-side cali- the nights of December 18, 19, and 20, the extra bration test in January 2011. This test showed difference in PWV (up to 90% on Dec. 18 and up that both PWV time-series compared well to a to 70% on the other two nights) is likely explained level of ±0.14 mm in 1-minute samples. At the by the formation of an inversion layer between the end of 2011, the refurbished radiometer was relo- Llano and Cerro Chajnantor. From December 21 cated to the CCAT site weather monitoring sta- to 23, no evidence for a significant inversion layer tion, slightly below the summit of Cerro Chaj- is seen. This could be explained by the relatively nantor, starting PWV data acquisition on 2011 high amount of PWV, which is a usual occurrance December 17. From the Llano, the horizontal dis- previous to snowstorms such as the one on Decem- tance of these radiometers is 3 km in the norhteast ber 22 that produced the power failure. direction and the altitude difference is 562 m. A To show that inversion layers between Llano snowstorm that started on 2011 December 22 pro- and Cerro Chajnantor generally occur during duced a power failure on the solar panels of the night-time, Figure 3 separates nighttime and day- CCAT station that was not resolved until Octo- time on a 24 h polar plot. Nighttime is considered ber 2012. Hence, only a week of PWV data was to be between 23:00 and 9:30 UTC (8 p.m. to obtained in this first campaign. 6:30 a.m., local time). For December 18 to 20, the Simultaneous PWV observations on Llano de drop in PWV is significantly larger at nighttime Chajnantor (APEX site) and Cerro Chajnantor compared to daytime. (CCAT site) from 2011 December 18-23 are pre- To determine the contribution and strength of sented in Figure 2. This figure shows that under the inversion layers, the drop in PWV was normal- normal atmospheric conditions, without the pres- ized to the 36.3% level with a standard deviation ence of inversion layers between the two sites, the σ of 13.9%. Figure 4 shows the normalized drop in drop in PWV from the Llano to Cerro Chajnan- PWV between both sites for each hour and each day. As mentioned in Figure 3, the differences are 1http://www.apex-telescope.org stronger and more significant for the first three

4 Normalized difference between APEX and CCAT 4 Dec 18 Dec 19 3 Dec 20 Dec 21 Dec 22 2

1

0

−1 00h 04h 08h 12h 16h 20h 24h UTC

Fig. 4.— Normalized PWV difference between Llano and Cerro Chajnantor, averaged over 1 h for each of the five days. The zero level corresponds to a PWV difference of 36.3%. The standard deviation of all data together is 13.9%.

nights. During all five afternoons and in the last the time, which is twice that of Llano de Chaj- two days, the drop in PWV is within the 1-σ level. nantor and the South Pole. Given the ±0.14 mm error in 1-minute samples Using the ATM model (Pardo et al. 2001) and we obtained between both radiometers in the side- the 350 µm zenith optical depth cumulative distri- by-side calibration test at the APEX site, we ex- butions obtained by Radford (2011), we compared pect a statistical error in 1 h of ±0.018 mm. Thus, the atmospheric transmissions at Llano de Chaj- for a PWV of 1.0 mm, the difference in PWV is ex- nantor and Cerro Chajnantor (Figure 5). Cerro pected to have a measurement error of 1.8%, sig- Chajnantor starts to show significant improve- nificantly lower than the 1-σ level of 13.9% seen ments in atmospheric transmission at frequencies in Figure 4. More details on the calibration cam- above 600 GHz. In the THz frequency regime, the paign and a larger dataset will be presented in a atmospheric transmission on Cerro Chajnantor is following paper. twice that of Llano de Chajnantor. For the mid-infrared to far-infrared (mid-IR 3. Worldwide site comparison to far-IR) spectral regime, Cerro Chajnantor presents an atmospheric transparency of up to Radford (2011) have compared the atmospheric 40% in the 36-38 µm band (Sako et al. 2008) and transparency conditions at several important a median atmospheric seeing of 0.6900 at 500 nm sites worldwide: Mauna Kea, Hawaii (4070 m); (Motohara et al. 2008). In this work, we im- the South Pole (2800 m); Llano de Chajnantor plemented a line-by-line, layer-by-layer radiative (5050 m); and Cerro Chajnantor (5612 m). Cu- transfer model in Matlab to calculate the atmo- mulative distributions through several years were spheric transmission for Cerro Chajnantor (Fig- obtained by measuring the 350 µm and 225 GHz ure 6). The model included a dry atmosphere zenith optical depth at these four sites. They have composition (Table 1) and a wet contribution of shown that Cerro Chajnantor, Llano de Chajnan- 0.70 mm (50% of the time), 0.42 mm (25%), and tor, and the South Pole are usually drier than 0.21 mm (10%) of PWV. The spectral information Mauna Kea. If we consider the superb observ- in the radiative transfer model comes from the HI- ing conditions when 350 µm zenith optical depths TRAN 2008 database at high spectral resolution are below 1.1, the best sites are arguably found (Rothman et al. 2009). in Chile. In this case, Cerro Chajnantor presents 350 µm zenith optical depths below 1.1 for 50% of Figure 6 shows the resulting transmission spec-

5 Atmospheric Transmission Llano de Chajnantor (5,050 m) Atmospheric Transmission Cerro Chajnantor (5,612 m) 1 1 50% 1.17 mm 50% 0.70 mm 0.9 25% 0.70 mm 0.9 25% 0.42 mm 0.8 10% 0.40 mm 0.8 10% 0.21 mm

0.7 0.7

0.6 0.6

0.5 0.5

0.4 0.4 Transmission Transmission 0.3 0.3

0.2 0.2

0.1 0.1

0 0 200 400 600 800 1000 1200 1400 1600 200 400 600 800 1000 1200 1400 1600 Frequency [GHz] Frequency [GHz]

Fig. 5.— Atmospheric transmission in the millimeter and submillimeter range for Llano de Chajnantor (left) and Cerro Chajnantor (right). Transmissions are derived from observed 350 µm zenith optical depths and the ATM model (Pardo et al. 2001). The green, red, and blue lines show the transmission for 50%, 25%, and 10% of the time, respectively. The inset lists the corresponding PWV in millimeters. Comparing left and right panels clearly shows the gain in atmospheric transmission by going to a higher site, especially at terahertz frequencies.

Besides the excellent atmospheric conditions, Table 1: Dry atmosphere composition other advantages of the Chajnantor area are its geographic location, accessibility, and logistics. Molecular species Formula Fraction per volume From latitude 23◦ S, the southern sky and a good Nitrogen N 0.7808 2 fraction of the northern sky is available for astro- Oxygen O 0.2095 2 nomical research. It has year-round accessibility Argon Ar 0.0093 via a paved road (the international route connect- Carbon Dioxide CO 0.000387 2 ing Chile and Argentina). It is close to large cities Methane CH 0.000001745 4 (Calama, Antofagasta) with ports, airports, and Carbon Monoxide CO 0.000000750 infrastructure that can provide general services, Nitrous Oxide N O 0.000000314 2 communications, and energy by local companies.

4. Parque Astron´omicode Atacama tra in the mid-IR to far-IR bands. The dryness Given the unique and excellent transparency of Cerro Chajnantor is such that for PWV con- conditions for astronomical observations of the ditions of 0.21 mm (10% of the time), the atmo- Chajnantor area, the Government of Chile has spheric transmission is higher than 95% in the Q secured an area of 363.81 km2 of fiscal land for (17-25 µm) band, higher than 60% in the Z (28- the exclusive use of scientific activities through 40 µm) band, and higher than 25% in several ob- the Parque Astron´omicode Atacama2. The Min- serving windows beyond Z band. Therefore, in istry of National Assets has given a long-term land extreme dry conditions (10% of the time), it is concession (50 years, starting from 2013) to the expected that Cerro Chajnantor offers completely National Commission for Research in Science and unique conditions for mid-IR to far-IR observa- Technology (CONICYT) (Diario Oficial 40.688, tions of cosmic sources from the ground. 2013) and given CONICYT the responsibility for These results are in agreement with previous its administration. simultaneous measurements of 350 µm zenith op- In Figure 7, the left panel shows the loca- tical depths obtained in Radford et al. (2008) and tion of the Parque Astron´omicode Atacama in Radford (2011) and scale height values from bal- loon measurements in Giovanelli et al. (2001). 2http://www.conicyt.cl/astronomia

6 Atmospheric Transmission Cerro Chajnantor (5,612 m) 1 50% 0.70 mm 0.9 25% 0.42 mm 0.8 10% 0.21 mm

0.7

0.6

0.5

0.4 Transmission 0.3

0.2

0.1

0 5 10 15 20 25 30 35 40 45 50 Wavelength [µm]

Fig. 6.— Atmospheric transmission on Cerro Chajnantor in the mid-IR to far-IR bands. Green, red, and blue lines are 50%, 25%, and 10% of the time, respectively, with the corresponding PWV in millimeters. the II Region of Antofagasta in Chile, and the tions, security procedures, coordination between right panel shows its extension in the marked different projects, and the provision of common area. This area has been protected by the Chilean services, among others. It will also maintain a re- government for scientific research and astronom- lationship with local native communities and the ical observations since 1998 by Presidential De- local administration. cree No 185. Two short-term land concessions (5 The area has also been declared a radio years each) to CONICYT allowed the establish- zone (no radio emissions above 31.3 GHz) over ment of the CBI telescope in 1999, followed by a 30 km radius and has a coordination zone of ASTE, APEX, NANTEN-2, ACT, mini-TAO, and 120 km radius to limit the granting of radio emis- Polarbear. CONICYT transferred part of the pro- sion licenses by the Ministry of Telecommunica- tected area to the ALMA project in 2003. tions around the park. This radio protection adds The main purpose of this park is to promote the to the existing light pollution protection by the installation of new astronomical facilities, to help Chilean government for the II, III and IV regions provide necessary infrastructure, and to coordi- of the country issued by the Ministry of Environ- nate the development of the area with the current ment (Decree No. 043, 2013) to preserve the dark and future observatories or telescopes. The park skies over this large territory of the country. is open to host big international projects as well as national or international university telescopes that 5. Conclusions could be operated in person or remotely. Other types of projects which do not interfere with as- Simultaneous measurements of PWV at the tronomical observations, i.e. meteorological mon- Llano de Chajnantor (5050 m) and Cerro Chaj- itoring, seismic stations, solar research, etc. are nantor (5612 m) were analyzed using 183 GHz ra- also welcome. diometers between 2011 December 18-23. We find that PWV decreases by ∼36% under normal atmo- To achieve this purpose, in October 2013, the spheric conditions on Cerro Chajnantor as com- National Commission for Research in Science and pared to Llano de Chajnantor. This decrease is Technology (CONICYT) created a foundation as consistent with a well developed boundary layer the entity responsible for managing and develop- of at least 600 m depth above the plateau. In the ing the Parque Astron´omicoAtacama. The foun- presence of inversion layers, the drop in PWV can dation will be responsible for preserving and pro- be as high as 90%, improving the observing condi- tecting the area from mining, geothermic, or other tions for the millimeter to mid-infrared astronomy. industrial activities that could interfere with as- tronomical activities, as well as generating gen- These measurements can constrain the scale eral norms and regulations for telescope installa- height continuously, and therefore they provide a

7 Fig. 7.— Left: The Chajnantor area in northern Chile, close to the border of Bolivia and Argentina. Right: Close-up showing the area of the Parque Astron´omicode Atacama and the ALMA concession (dark brown). The colored dots show potential sites where new telescopes can be installed in the park. very useful tool for site characterization and ra- QUIMAL grant No 130005. M.R. wishes to dio telescope operations, and are an alternative acknowledge support from CONICYT (Chile) method for determining the strength of inversion through FONDECYT grant No 1140839. layers. They also help to illustrate the advantage M.R. would like to thank the presidents of of going to higher sites for astronomical observa- CONICYT, Dr. Eric Goles (2000-2006), Mrs. Vi- tions at submillimeter and infrared wavelengths. vian Heyl (2006-2010), and Dr. Jos´eMiguel Aguil- Technical aspects such as extremely low water era (2010-2013) for the permanent support of the vapor content, high transparency, and sky cover- Parque Astron´omicoAtacama project, and Mrs. age as well as operational aspects such as acces- Alicia Norambuena, Coordinator of the Parque sibility, communications, and services makes the Astron´omicoAtacama, for her strong commitment Chajnantor area unique for astronomical observa- to the project and for ensuring the protection of tories. the area from activities that could interfere with The government of Chile, considering these as- the astronomical projects. We would also like to pects, has established the Parque Astron´omico de thank the Ministry of National Assets for its sup- Atacama administered by CONICYT to protect port during the administrative processes of the and develop this area, identified as the best site in land concessions. the world for millimeter, submillimeter, and mid- We thank Dr. Martin Shepherd and an anony- infrared astronomy. mous reviewer for useful comments that con- tributed to improve the reading of the manuscript. We acknowledge support from Dr. Simon Rad- ford and the CCAT project to install and operate REFERENCES the 183 GHz radiometer. We thank support from Bustos, R. et al., 2000, ALMA Memo Series, 333 the APEX staff and for the APEX PWV data. Bustos, R. 2001, ALMA Memo Series, 379 We also thank Jos´eCort´esand Hugo Pacheco for their help at the site and the CONICYT Astron- Bustos, R. & Nagar, N., 2011, in Astronomical omy Program. Site Testing Data in Chile, RevMexAA, 41, 104 R.B. and N.N. work was supported by ALMA- Delgado, G. et al., 1999, ALMA Memo Series, 271 CONICYT grants 31080022 and 31070015. R.B. Fowler, J. et al., 2005, ApJS, 156, 1 acknowledges support from CONICYT through

8 Garreaud, R. 2011, in Astronomical Site Testing Data in Chile, RevMexAA, 41, 5 Giovanelli, R. et al. 2001, PASP, 113, 803 G¨usten,R. et al., 2006, A&A, 454, 13 Kawamura, A. et al., 2005, IAUS, 235, P275 Kermish, Z. et al., 2012, Proc. SPIE, 8452, 84521C Kohno, K. 2005, in The Cool Universe: Observing Cosmic Dawn, ASP Conference Series Proceed- ings, 344, 242 Kusaka, A. et al., 2014, Rev. Sci. Instrum., 85, 024501 Marrone, D. et al. 2004, ApJ, 612, 940 Matsushita, S. et al., 1999, PASJ, 51, 603 McKinnon, M. 1988, ALMA Memo Series, 49 Miller, A. et al., 2002, ApJS, 140, 115 Motohara, K. et al., 2008, Proc. SPIE 7012, 701244 Mu˜noz,R. et al. 2010, J. Climate, 24, 1013 Ot´arola,A. et al., 2010, PASP, 122, 470 Padin, S. et al., 2002, PASP, 114, 83 Paine, S. et al., 2000, PASP, 112, 108 Pardo, J. et al. 2001, IEEE Trans. Antennas Propag., 49, 1683 QUIET collaboration, 2013, ApJ, 768, 9 Radford, S. & Holdaway, M., 1998, Advanced Technology MMW, Radio, and Terahertz Tele- scopes SPIE, 3357, 486 Radford, S. et al. 2008, Proc. SPIE, 7012, 70121Z Radford, S. 2011, in Astronomical Site Testing Data in Chile, RevMexAA, 41, 87 Roosen, R. & Angione, R., 1977, PASP, 89, 814 Rothman, L. et al., 2009, J. Quant. Spec. Ra- diat. Transf., 110, 533 Sako, S. et al. 2008, Proc. SPIE, 7012, 70122T Swetz, D. et al., 2011, ApJS, 194, 41 Torbet, E. et al., 1999, ApJ, 521, L79 Turner, D. et al., 2012, GeoRL, 39, 10801 Wootten, A. & Thompson, A. R., 2009, Proc. of the IEEE, 97, 8, 1463

This 2-column preprint was prepared with the AAS LATEX macros v5.2.

9