Sky subtraction for LAMOST

Sky subtraction is a key technique in data reduction of multi-fiber spectra. Knowledge of characteristics related to the instrument is necessary to determine the method adopted in sky subtraction.In this study, we describe the sky subtraction method designed for the Large sky Area Multi-Object fiber Spectroscopic Telescope(LAMOST) survey. The method has been integrated into the LAMOST 2D Pipeline v2.6 and applied to data from LAMOST DR3 and later. For LAMOST, calibration using sky emission lines is used to alleviate the position-dependent(and thus time-dependent) ~ 4% fiber throughput uncertainty and small wavelength instability(0.1 A) during observation. Sky subtraction using principal component analysis(PCA) further reduces 25% of the sky line residual from OH lines in the red part of LAMOST spectra after the master sky spectrum, which is derived from a B-spline fit of 20 sky fibers in each spectrograph. Using this approach, values are adjusted by a sky emission line and subtracted from each fiber. Further analysis shows that our wavelength calibration accuracy is about 4.5 km s~(-1), and the averages of residuals after sky subtraction are about 3% for sky emission lines and 3% for the continuum region. The relative sky subtraction residuals vary with moonlight background brightness, and can reach as low as 1.5% for regions that have sky emission lines during a dark night.Tests on F stars with both similar sky emission line strength and similar object continuum intensity show that the sky emission line residual of LAMOST is smaller than that of the SDSS survey.     

Harmonic analysis of cosmic microwave background data – II. From ring-sets to the sky

Despite the fact that the physics of the cosmic microwave background anisotropies is most naturally expressed in Fourier space, pixelized maps are almost always used in the analysis and simulation of microwave data. A complementary approach is investigated here, in which maps are used only in the visualization of the data, and the temperature anisotropies and polarization are only ever expressed in terms of their spherical multipoles. This approach has a number of advantages: there is no information loss (assuming a band-limited observation); deconvolution of asymmetric beam profiles and the temporal response of the instrument are naturally included; correlated noise can easily be taken into account, removing the need for additional 'destriping'; polarization is also analysed in the same framework; and reliable estimates of the spherical multipoles of the sky and their errors are obtained directly for subsequent component separation and power spectrum estimation. The formalism required to analyse experiments which survey the full sky by scanning on circles is derived here, with particular emphasis on the Planck mission. A number of analytical results are obtained in the limit of simple scanning strategies. Although there are non-trivial computational obstacles to be overcome before the techniques described here can be implemented at high resolution, if these can be overcome the method should allow for a more robust return from the next generation of full-sky microwave background experiments.     

Detection of Leonids meteoric dust in the upper atmosphere by polarization measurements of the twilight sky

A new method for the meteoric dust detection in the upper atmosphere based on the polarimetric observations of the twilight sky is proposed. Polarization measurements are effective for detection of the meteoric dust scattering on the background consisting basically of troposphere multiple scattering. The method is based on the observed and explained polarization properties of the sky background during different stages of twilight. It is used to detect the mesosphere dust after the Leonids maximum in 2002, estimate its altitude range and to investigate its evolution—slow decrease of the altitude. The polarization method takes into account the multiple scattering and sufficient contribution of moonlight scattering background and turns out to be more sensitive than existing analogs used in the present time.     

On the reduction and presentation of data in astronomical two-channel photopolarimetry

Many different methods exist for reducing data obtained when an astronomical source is studied with a two-channel polarimeter, such as a Wollaston prism system. This paper presents a rigorous method of reducing the data from raw aperture photometry, and evaluates errors both by a statistical treatment, and by propagating the measured sky noise from each frame. The reduction process performs a hypothesis test for the presence of linear polarization. The probability of there being a non-zero polarization is obtained, and the best method of obtaining the normalized Stokes Parameters is discussed. Point and interval estimates are obtained for the degree of linear polarization, which is subject to positive bias; and the polarization axis is found.     

Directionality in the Wilkinson Microwave Anisotropy Probe polarization data

Polarization is the next frontier of cosmic microwave background analysis, but its signal is dominated over much of the sky by foregrounds which must be carefully removed. To determine the efficacy of this cleaning, it is necessary to have sensitive tests for residual foreground contamination in polarization sky maps. The dominant Galactic foregrounds introduce a large-scale anisotropy on to the sky, so it makes sense to use a statistic sensitive to overall directionality for this purpose. Here, we adapt the rapidly computable     statistic of Bunn and Scott to polarization data, and demonstrate its utility as a foreground monitor by applying it to the low-resolution Wilkinson Microwave Anisotropy Probe 3-yr sky maps. With a thorough simulation of the maps' noise properties, we find no evidence for contamination in the foreground cleaned sky maps.     

施密特巡天星表的现状和展望

扼要介绍了施密特巡天底片及基于此编制而成的施密特巡天星表的发展历程,并详细介绍了GSC2.3和USNO-B1.0的情况.分析研究了这类星表存在的问题和原因所在,并提出了一个新的计划:利用现有施密特底片资料并增加新的第3期观测,编制一个具有绝对自行、多色测光、系统均匀的高密度全天星表以满足各方面的需要.     

Digital complex correlator for a C-band polarimetry survey

The international Galactic Emission Mapping project aims to map and characterize the polarization field of the Milky Way. In Portugal it will map the sky polarized emission of the Northern Hemisphere in C-band and provide templates for map calibration and foreground control of microwave space data to be provided by ESA Planck Surveyor mission and later missions. The receiver system is equipped with a novel receiver with a full digital back-end using a low-cost Field Programmable Gate Array without compromising its performance relation. This new digital backend comprises a base-band complex cross-correlator outputting the four Stokes parameters of the incoming polarized radiation. In this document we describe the design and implementation of the complex correlator using the FPGA and the dedicated digitizers at each receiver arm, detailing the method applied at the several algorithm stages. This correlator is suitable for large sky area polarization continuum surveys.     

Harmonic analysis of cosmic microwave background data – I. Ring reductions and point-source catalogue

We present a harmonic model for the data analysis of an all-sky cosmic microwave background survey, such as Planck , where the survey is obtained through ring-scans of the sky. In this model, resampling and pixelization of the data are avoided. The spherical transforms of the sky at each frequency, in total intensity and polarization, as well as the bright-point-source catalogue, are derived directly from the data reduced on to the rings. Formal errors and the most significant correlation coefficients for the spherical transforms of the frequency maps are preserved. A clean and transparent path from the original samplings in the time domain to the final scientific products is thus obtained. The data analysis is largely based on Fourier analysis of rings; the positional stability of the instrument's spin axis during these scans is a requirement for the data model and is investigated here for the Planck satellite. Brighter point sources are recognized and extracted as part of the ring reductions and, on the basis of accumulated data, used to build the bright-point-source catalogue. The analysis of the rings is performed in an iterative loop, involving a range of geometric and detector response calibrations. The geometric calibrations are used to reconstruct the paths of the detectors over the sky during a scan and the phase offsets between scans of different detectors; the response calibrations eliminate short- and long-term variations in detector response. Point-source information may allow the reconstruction of the beam profile. The reconstructed spherical transforms of the sky in each frequency channel form the input to the subsequent analysis stages. Although the methods in this paper were developed with the data processing for the Planck satellite in mind, there are many aspects which have wider implementation possibilities, including the construction of real-space pixelized maps.     

Polarimetry of the twilight sky and stratospheric aerosol

We propose a modification of the method of polarimetric measurements of the twilight sky, traditionally performed in a zenith direction, to study physical properties of the stratospheric aerosol (at altitudes higher than 30 km). The measurements carried out in zenith directions as a rule limit phase angles by values of 80–100°. We suggest setting up the declination of the telescope equal to the declination of the sun and measuring the polarization degree of the twilight sky at different values of the right ascension. It will allow us not only to enhance the range of the phase angles but also to plan observations in a way to obtain data on the phase dependence of the polarization degree of the light scattered by atmospheric layers at different altitudes.     

Predicted interplanetary distribution of Lyman-σ intensity and polarization

Our current understanding of the spatial and temporal distribution of interplan etary neutral hydrogen is currently limited to a comparison of Lyman-σ photometric data with predictions of the solar backscattered radiation using theoretical models. In this paper, how the uncertainties in current model calculations could be reduced through the future use of polarization measurements made from interplanetary spacecraft is investigated. In particular, inquiry into how a mapping of the degree of linear polarization made from a spacecraft at various locations in the Solar System can improve knowledge of the interstellar wind parameters, number density, temperature, and velocity, is made. A polarization measurement can, in principle, be made with very high precision. In this regard, being a relative quantity, a polarization measurement can be made independent of instrumental calibration and long-term sensitivity degredation. Furthermore, the sky distribution of both intensity and polarization has been calculated using a variety of models for the neutral hydrogen. It is found that the polarization distribution over the sky is quite different from that of the intensity distribution. It is also showed that the maximum degree of polarization of the Laymam-σ line increases with heliocentric distance of the spacecraft, varying from 0 up to ～ 18% at 20 AU.     

A survey of polarization in the JVAS/CLASS flat-spectrum radio source surveys – II. A search for aligned radio polarizations

We have used the very large Jodrell Bank VLA Astrometric Survey/Cosmic Lens All-Sky Survey 8.4-GHz surveys of flat-spectrum radio sources to test the hypothesis that there is a systematic alignment of polarization position angle vectors on cosmological scales of the type claimed by Hutsemékers et al. The polarization position angles of 4290 sources with polarized flux density ≥1 mJy have been examined. They do not reveal large-scale alignments either as a whole or when split in half into high-redshift  ( z ≥ 1.24)  and low-redshift subsamples. Nor do the radio sources which lie in the specific areas covered by Hutsemékers et al. show any significant effect. We have also looked at the position angles of parsec-scale jets derived from very long baseline interferometry observations and again find no evidence for systematic alignments. Finally, we have investigated the correlation between the polarization position angle and those of the parsec-scale jets. As expected, we find that there is a tendency for the polarization angles to be perpendicular to the jet angles. However, the difference in jet and polarization position angles does not show any systematic trend in different parts of the sky.     

埃文斯目视日晕光度计测量原理研究

埃文斯目视日晕光度计（Evans Visual Sky Photometer,EVSP）是应用于日冕仪选址的重要仪器,从20世纪40年代一直使用至今,它为现代日晕光度计的定标提供了参考标准。通过使用云南天文台现存的一架EVSP研究了它的工作原理,并重点介绍了所利用的简易定标方法。给出了国际现有的多台EVSP日晕亮度定标曲线。由于EVSP内部光学元件反射率和透过率,以及中性渐变光楔的光学密度等存在未知的时间缓变特点,因此利用这种新定标手段可以高效经济地获得各自的定标曲线。     

The Rhodes/HartRAO 2326-MHz radio continuum survey

We have made scanning observations with the HartRAO 26-m radio telescope to obtain a pencil-beam map of 67 per cent of the sky at 2326 MHz. This is the highest resolution and highest frequency radio continuum map of this type made of such a large area of sky.   In this paper we describe the observations and data reduction procedures used to produce the survey. The resulting map has an angular resolution (HPBW) of 20 arcmin, and the rms pointing accuracy is 1.2 arcmin. The rms noise fluctuations are less than 30-mK T _FB over the whole map. We estimate that the uncertainty in the temperature scale is less than 5 per cent, and that the error in the absolute zero level is better than 80-mK T _FB in any direction.   High-contrast half-tone images of the data with a model of the diffuse galactic background subtracted are presented. These images show many complex emission structures up to and beyond 50° latitude, and illustrate the quality of the data. Extracts from the survey data are available via FTP by arrangement with the authors.     

Comment on: “Electrical coupling of the E- and F-regions”: by R. A. Heelis,P. C. Kendall,R. J. Moffett,D. W. Windle and H. Rishbeth

New measurements have been made of the reduced polarization in a Fraunhofer line (Hβ) of blue sky light. The value of the depolarization agrees well with a previously reported intensity filling-in effect. There is no evidence of a change of the azimuth of vibration of the polarization across the line profile.     

Simulation of a zenith-field sky survey on RATAN-600

In this paper we simulate a deep multi-frequency zenith-field sky survey on RATAN-600 (the RZF survey). In our simulations we use the 1.4-GHz sky images obtained in the NVSS survey. We convolved NVSS images with the two-dimensional power beam pattern of RATAN-600 and obtain simulated 24-hour scans of sky transits at all wavelengths of the RZF survey. For the 7.6-cm wavelength we analyze the effect of the image area size on the results of the simulation. We estimate the accuracy of the determination of source fluxes on simulated scans and derive the distributions of the spectral indices of the sources. We use the simulated scans to clean real records of the RZF survey at 7.6 cm. The standard error of the residual noise at this wavelength is about 1 mJy.     

Eigenvector Sky Subtraction

We develop a new method for estimating and removing the spectrum of the sky from deep spectroscopic observations. Our method does not rely on simultaneous measurement of the sky spectrum with the object spectrum. The technique is based on the iterative subtraction of continuum estimates and eigenvector sky models derived from singular value decompositions of sky spectra and sky spectra residuals. Using simulated data derived from small-telescope observations, we demonstrate that the method is effective for faint objects on large telescopes. We discuss simple methods to combine our new technique with the simultaneous measurement of sky to obtain sky subtraction very near the Poisson limit.     

Dust coma of comet C/2009 P1 (Garradd) by imaging polarimetry

Comet C/2009 P1 (Garradd) was observed by imaging polarimetry for nearly 5 months from October 2011 to March 2012, over an intermediate phase angle range (28°–35°). Two months before perihelion and one month after, dust particles seem to be ejected all around the optocenter and jets extend to distances greater than 40,000 km. An increase of activity is noticed in intensity and polarization after perihelion. Two months before perihelion and one month after, the dust emission seems to be all around the optocenter. Two and three months after perihelion the jets are mainly toward the solar direction with an extension of more than 20,000 km projected on the sky. The values of the aperture polarization are comparable to those of other comets. On the polarization maps in October 2011 and January 2012 the higher polarization zones extend in large regions perpendicularly to the solar direction where jets are also observed. In February and March 2012, the polarization in the jets is larger in the solar direction than in the surrounding coma. By its activity visible on intensity images and polarization maps at large distances from the nucleus, comet Garradd probably belongs to the high‐P_max class of comets.     

A preferred-direction statistic for sky maps

Large patterns could exist on the microwave sky as a result of various non-standard possibilities for the large-scale Universe – rotation or shear, non-trivial topology, and single topological defects are specific examples. All-sky (or nearly all-sky) CMB data sets allow us, uniquely, to constrain such exotica, and it is therefore worthwhile to explore a wide range of statistical tests. We describe one such statistic here, which is based on determining gradients and is useful for assessing the level of 'preferred directionality' or 'stripiness' in the map. This method is more general than other techniques for picking out specific patterns on the sky, and it also has the advantage of being easily calculable for the mega-pixel maps which will soon be available. For the purposes of illustration, we apply this statistic to the four-year COBE DMR data. For future CMB maps, we expect this to be a useful statistical test of the large-scale structure of the Universe. In principle, the same statistic could also be applied to sky maps at other wavelengths, to CMB polarization maps, and to catalogues of discrete objects. It may also be useful as a means of checking for residual directionality (e.g. from Galactic or ecliptic signals) in maps.     

The First Solar Seeing Profile Measurement with Two Apertures and Multiple Guide Regions

The sky brightness is a critical parameter for estimating the coronal observation conditions for a solar observatory. As part of a site-survey project in Western China, we measured the sky brightness continuously at the Lijiang Observatory in Yunnan province in 2011. A sky brightness monitor (SBM) was adopted to measure the sky brightness in a region extending from 4.5 to 7.0 apparent solar radii based on the experience of the Daniel K. Inouye Solar Telescope (DKIST) site survey. Every month, the data were collected manually for at least one week. We collected statistics of the sky brightness at four bandpasses located at 450, 530, 890, and 940 nm. The results indicate that aerosol scattering is of great importance for the diurnal variation of the sky brightness. For most of the year, the sky brightness remains under 20 millionths per airmass before local Noon. On average, the sky brightness is less than 20 millionths, which accounts for 40.41% of the total observing time on a clear day. The best observation time is from 9:00 to 13:00 (Beijing time). The Lijiang Observatory is therefore suitable for coronagraphs investigating the structures and dynamics of the corona.     

Cosmic bandits: Exploration versus exploitation in CMB B-mode experiments

A preferred method to detect the curl-component, or B-mode, signature of inflationary gravitational waves (IGWs) in the cosmic microwave background (CMB) polarization, in the absence of foregrounds and lensing, is a prolonged integration over a single patch of sky of a few square degrees. In practice, however, foregrounds abound and the sensitivity to B modes can be improved considerably by finding the region of sky cleanest of foregrounds. The best strategy to detect B modes thus involves a tradeoff between exploration (to find lower-foreground patches) and exploitation (through prolonged integration). This problem is akin to the multi-armed bandit (MAB) problem in probability theory, wherein a gambler faces a series of slot machines with unknown winning odds and must develop a strategy to maximize his/her winnings with some finite number of pulls. While the optimal MAB strategy remains to be determined, a number of algorithms have been developed in an effort to maximize the winnings. Here, based on this resemblance, we tackle the search for IGW B modes with single frequency experiments in the presence of spatially varying foregrounds by developing adaptive survey strategies to optimize the sensitivity to IGW B modes. We demonstrate, using realistic foreground models and taking lensing-induced B modes into account, that adaptive experiments can substantially improve the upper bound on the tensor-to-scalar ratio (by factors of 2 and 3 in single frequency experiments, and possibly even more). Similar techniques can be applied to other surveys, including 21-cm measurements of signatures of the epoch of reionization, searches for a stochastic primordial gravitational wave background, deep-field imaging by the James Webb Space Telescope or various radio interferometers, and transient follow-up searches.