Astrometry.net与SCAMP在天测位置定标中的应用研究

为探究Astrometry.net和SCAMP两款软件在天测位置定标时的应用差异以及以不同方式运行软件所得定标效果的差异,分别介绍了两款软件位置定标的原理,以ZTF(Zwicky Transient Facility)的100张巡天图像为测试图像,以Gaia DR2为参考星表,分别以4种不同的定标流程对测试图像进行了位置定标。从计算总耗时、星表匹配情况和坐标均方根(Root Mean Square,RMS)偏差3个方面,比较了不同流程所得定标的结果,同时也与ZTF给出的定标结果进行了对比。结果表明,Astrometry.net能够快速给出粗略的线性定标结果,SCAMP可以快速在前者的基础上直接进行扭曲修正,两者相结合的定标流程平均耗时仅为1 s,所得定标结果的平均赤经赤纬均方根偏差均小于70毫角秒(小于0.1个像素),优于仅使用Astrometry.net所得定标结果。由此可见,搭配使用Astrometry.net和SCAMP可以实现更快速且更准确的天测位置定标,可推广应用到时域天文巡天项目的数据处理上。

The Possibility of Detecting our Solar System through Astrometry

Searching for exoplanets with different methods has always been the focus of astronomers over the past few years.Among multiple planet detection techniques,astrometry stands out for its capability to accurately determine the orbital parameters of exoplanets.In this study,we examine the likelihood of extraterrestrial intelligent civilizations detecting planets in our solar system using the astrometry method.By conducting injection-recovery simulations,we investigate the detectability of the four giant planets in our solar system under different observing baselines and observational errors.Our findings indicate that extraterrestrial intelligence could detect and characterize all four giant planets,provided they are observed for a minimum of 90 yr with signal-noise ratios exceeding 1.For individual planets such as Jupiter,Saturn,and Neptune,a baseline that surpasses half of their orbital periods is necessary for detection.However,Uranus requires longer observing baselines since its orbital period is roughly half of that of Neptune.If the astrometry precision is equal to or better than 10μas,all 8707 stars located within30 pc of our solar system possess the potential to detect the four giant planets within 100 yr.Additionally,our prediction suggests that over 300 stars positioned within 10 pc from our solar system could detect our Earth if they achieve an astrometry precision of 0.3μas.

Improved astrometry of space debris with image restoration

In order to implement an observing strategy, image degradation that occurs during optical observation of space debris is ineluctable and has distinct characteris- tics. Image restoration is presented as a way to remove the influence of degradation in CCD images of space debris, based on assumed PSF models with the same F-WHM as images of the object. In the process of image restoration, the maximum entropy method is adopted. The results of reduction using observed raw CCD images indi- cate that the precision in estimating positions of objects is improved and the effects of degradation are reduced. Improving the astrometry of space debris using image restoration is effective and feasible.

Astrometry of three near Earth asteroids with the Lijiang 2.4 m telescope

Under the framework of observational campaigns organized by the GAIA Follow Up Network for Solar System Objects, three near Earth asteroids, 367943 Duende, 99942 Apophis and 2013 TV135, were observed with the Lijiang 2.4m telescope administered by Yunnan Observatories. The software package PRISM was used to calibrate the CCD fields and measure the positions of 99942 Apophis and 2013 TV135, and our own software was used for 367943 Duende. A comparison of the results show that the ephemerides of INPOP10a and JPL are consistent for 99942 Apophis and 2013 TV135, however, they are quite inconsistent for 367943 Duende. Moreover, we have found that differences between the mean values in the ephemerides of INPOP10a and JPL are about 72＂ and -199＂ in right ascension and declination respectively for 367943 Duende. Moreover, the ephemeris published by JPL is reliable in terms of the mean observed-minus-calculated （O - C） residuals in right ascension and declination of about 2.72＂ and 1.49＂ respectively.

PyMsOfa:A Python Package for the Standards of Fundamental Astronomy(SOFA)Service

The Standards of Fundamental Astronomy(SOFA)is a service provided by the International Astronomical Union that offers algorithms and software for astronomical calculations,which was released in two versions for FORTRAN77 and ANSI C,respectively.In this work,we implement the Python package PyMsOfa for SOFA service by three ways:(1)a Python wrapper package based on a foreign function library for Python(ctypes),(2)a Python wrapper package with the foreign function interface for Python calling C code(cffi)and(3)a Python package directly written in pure Python codes from SOFA subroutines.The package PyMsOfa has fully implemented 247 functions of the original SOFA routines released on 2023 October 11.In addition,PyMsOfa is also extensively examined,which is exactly consistent with those test examples given by the original SOFA.This Python package can be suitable to not only the astrometric detection of habitable planets from the Closeby Habitable Exoplanet Survey mission,but also for the frontier themes of black holes and dark matter related to astrometric calculations and other fields.The source codes are available via http://pypi.org/project/PyMsOfa/and https://github.com/CHES2023/PyMsOfa.

A Positioning System based on Communication Satellites and the Chinese Area Positioning System(CAPS)

The Chinese Area Positioning System （CAPS） is a positioning system based on satellite communication that is fundamentally different from the 3＂G＂ （GPS, GLONASS and GALILEO） systems. The latter use special-purpose navigation satellites to broadcast navigation information generated on-board to users, while the CAPS transfers ground-generated navigation information to users via the communication satellite. In order to achieve accurate Positioning, Velocity and Time （PVT）, the CAPS employs the following strategies to over- come the three main obstacles caused by using the communication satellite： （a） by real-time following-up frequency stabilization to achieve stable frequency; （b） by using a single carrier in the transponder with 36 MHz band-width to gain sufficient power; （c） by incorporating Decommissioned Geostationary Orbit communication satellite （DGEO）, barometric pressure and Inclined Geostationary Orbit communication satellite （IGSO） to achieve the 3-D posi- tioning. Furthermore, the abundant transponders available on DGEO can be used to realize the large capacity of communication as well as the integrated navigation and communication. With the communication functions incorporated, five new functions appear in the CAPS： （1） combination of navigation and communication; （2） combination of navigation and high accu- racy orbit measurement; （3） combination of navigation message and wide/local area differen- tial processing; （4） combination of the switching of satellites, frequencies and codes; and （5） combination of the navigation message and the barometric altimetry. The CAPS is thereby labelled a PVT5C system of high accuracy. In order to validate the working principle and the performance of the CAPS, a trial system was established in the course of two years at a cost of about 20 million dollars. The trial constellation consists of two GEO satellites located at E87.5° and E110.5°, two DGEOs located at E130° and E142°, as well as barometric altimetry as a virtual satellite. Static and dynamic performance tests were completed for the Eastern, the Western, the Northern, the Southern and the Middle regions of China. The evaluation results are as follows： （1） land static test, plane accuracy range： C/A code, 15-25 m; P code, 5-10 meters; altitude accuracy range, 1- m; （2） land dynamic test, plane accuracy range, C/A code, 15-25 m; P code, 8-10m; （3） velocity accuracy, C/A code, 0.13-0.3 m s-1, P code, 0.15-0.17 m s- 1; （4） timing accuracy, C/A code, 160 ns, P code, 13 ns; （5） timing compared accuracy of Two Way Satellite Time and Frequency Transfer （TWSTFT）, average accuracy, 0.068 ns; （6） random error of the satellite ranging, 10.7 mm; （7） orbit determination accuracy, better than 2 m. The above stated random error is 1σ error. At present, this system is used as a preliminary operational system and a complete system with 3 GEO, 3 DGEO and 3 IGSO is being established.

The principle of measuring unusual change of underground mass by optical astrometric instrument

Optical observations from a ground-based astrometric instrument provide a unique set of data sensi- tive to variations of local plumb line, which represent the change of local gravitational field on the ground. The local gravitational change could be caused by density variation or movement of underground mass in addition to tidal action. Earthquake was found to be related with gravitational change. A seismic zone could have abnor- mal underground material density, causing the variation of plume line on the ground. In this frame, the astro- metric telescope can measure the variation of local plumb line by observing the local astronomical longitudes and latitudes from stars in the Hipparcos catalog. In this study, we estimate the deflection angle of the plumb line on a ground site, and give a relation be- tween the angle, abnormal mass and site distance （depth and horizontal distance）. Then we derive the abnor- mality of underground material density using the plumb lines measured at different sites, and study the earth- quake gestation, development and occurrence. Using the deflection angles of plumb lines observed at two sites, we give a method to calculate the mass and the center of gravity of underground materials. We also esti- mate the abnormal masses of latent seismic zones with different energy, using thermodynamic relations, and in- troduce a new optical astrometric instrument we had developed.

Experiment on diffuse reflection laser ranging to space debris and data analysis

Space debris poses a serious threat to human space activities and needs to be measured and cataloged. As a new technology for space target surveillance, the measurement accuracy of diffuse reflection laser ranging （DRLR） is much higher than that of microwave radar and optoelectronic measurement. Based on the laser ranging data of space debris from the DRLR system at Shanghai Astronomical Observatory acquired in March-April, 2013, the characteristics and precision of the laser ranging data are analyzed and their applications in orbit determination of space debris are discussed, which is implemented for the first time in China. The experiment indicates that the precision of laser ranging data can reach 39 cm-228 cm. When the data are sufficient enough （four arcs measured over three days）, the orbital accuracy of space debris can be up to 50 m.

The spiral structure of the Milky Way

The morphology and kinematics of the spiral structure of the Milky Way are long-standing problems in astrophysics.In this review we firstly summarize various methods with different tracers used to solve this puzzle.The astrometry of Galactic sources is gradually alleviating this difficult situation caused mainly by large distance uncertainties, as we can currently obtain accurate parallaxes(a few μas) and proper motions(≈1 km s-1) by using Very Long Baseline Interferometry(VLBI).On the other hand, the Gaia mission is providing the largest, uniform sample of parallaxes for O-type stars in the entire Milky Way.Based upon the VLBI maser and Gaia O-star parallax measurements, nearby spiral structures of the Perseus, Local, Sagittarius and Scutum Arms are determined in unprecedented detail.Meanwhile, we estimate fundamental Galactic parameters of the distance to the Galactic center,R0, to be 8.35 ± 0.18 kpc, and circular rotation speed at the Sun, Θ0, to be 240±10 km s-1.We found kinematic differences between O stars and interstellar masers: the O stars, on average, rotate faster,>8 km s-1than maser-traced high-mass star forming regions.

The use of laser ranging to measure space debris

Space debris is a major problem for all the nations that are currently active in space. Adopting high-precision measuring techniques will help produce a reliable and accurate catalog for space debris and collision avoidance. Laser ranging is a kind of real-time measuring technology with high precision for space debris observation. The first space-debris laser-ranging experiment in China was performed at the Shanghai Observatory in July 2008 with a ranging precision of about 60-80 cm. The experi- mental results showed that the return signals from the targets with a range of 900 km were quite strong, with a power of 40W （2J at 20 Hz） using a 10ns pulse width laser at 532 nm wavelength. The performance of the preliminary laser ranging system and the observed results in 2008 and 2010 are also introduced.

CHES: A Space-borne Astrometric Mission for the Detection of Habitable Planets of the Nearby Solar-type Stars

The Closeby Habitable Exoplanet Survey(CHES) mission is proposed to discover habitable-zone Earth-like planets of nearby solar-type stars(~10 pc away from our solar system) via microarcsecond relative astrometry.The major scientific objectives of CHES are:to search for Earth Twins or terrestrial planets in habitable zones orbiting100 FGK nearby stars;further to conduct a comprehensive survey and extensively characterize nearby planetary systems.The primary payload is a high-quality,low-distortion,high-stability telescope.The optical subsystem is a coaxial three-mirror anastigmat(TMA) with a 1.2 m-aperture,0°.44 × 0°.44 field of view and 500 nm-900 nm working wave band.The camera focal plane is composed of a mosaic of 81 scientific CMOS detectors each with4 k × 4 k pixels.The heterodyne laser interferometric calibration technology is employed to ensure microarcsecond level(1 μas) relative astrometry precision to meet the requirements for detection of Earth-like planets.The CHES satellite operates at the Sun-Earth L2 point and observes all the target stars for 5 yr.CHES will offer the first direct measurements of true masses and inclinations of Earth Twins and super-Earths orbiting our neighbor stars based on microarcsecond astrometry from space.This will definitely enhance our understanding of the formation of diverse nearby planetary systems and the emergence of other worlds for solar-type stars,and finally provide insights to the evolution of our own solar system.

Searching for exoplanets with HEPS:I.detection probability of Earth-like planets in multiple systems

The astrometry method has great advantages in searching for exoplanets in the habitable zone around solar-like stars. However, the presence of multiple planets may cause a problem with degeneracy when trying to compute accurate planet parameters from observation data and reduce detectability. The degeneracy problem is extremely critical, especially in a space mission which has limited observation time and cadence. In this series of papers, we study the detectability of habitable Earth-mass planets in different types of multi-planet systems, aiming to find the most favorable targets for the potential space mission–Habitable ExoPlanet Survey(HEPS). In the first paper, we present an algorithm to find planets in the habitable zone around solar-like stars using astrometry. We find the detectability can be well described by planets' signal-to-noise ratio(SNR) and a defined parameter S = M2/(T1-T2)2, where M2 and T2are the mass and period of the second planet, respectively. T1 is the period of the planet in the habitable zone. The parameter S represents the influence of planetary architectures. We fit the detectability as a function of both the SNR of the planet in the habitable zone and the parameter S. An Earth-like planet in a habitable zone is harder to detect(with detectability PHP< 80%) in a system with a hot Jupiter or warm Jupiter(within2 AU), in which the parameter S is large. These results can be used in target selections and to determine the priority of target stars for HEPS, especially when we select and rank nearby planet hosts with a single planet.

Test area of the SAGE survey

Sky surveys represent one of the most important efforts to improve developments in astrophysics,especially when using new photometric bands. We are performing the Stellar Abundance and Galactic Evolution(SAGE) survey with a self-designed SAGE photometric system, which is composed of eight photometric bands. The project mainly aims to study the stellar atmospheric parameters of ～0.5 billion stars in ～12 000 deg2 of the northern sky, which mainly focuses on Galactic astronomy, as well as some aspects of extragalactic astronomy. This work introduces the detailed data reduction process of the test field NGC 6791, including the data reduction of single-exposure images and stacked multi-exposure images, and properties of the final catalog.

Digitizer of astronomical plates at Shanghai Astronomical Observatory and its performance test

Before CCD detectors were widely employed in observational astronomy, the main method of detection was the use of glass astrophotographic plates. Astronomical plates have been used to record information on the position and activity of celestial bodies for more than 100 years. There are about 30 000 astronomical plates in China, and the digitization of astronomical plates is of great significance for permanent preservation and to make full use of these valuable observation data. A digitizer with high precision and high measuring speed is a key piece of equipment for carrying out the task of dig- itizing these astronomical plates. A digitizer for glass astrophotographic plates was developed jointly by Shanghai Astronomical Observatory and Nishimura Co., Ltd of Japan. The digitizer＇s hardware was manufactured by Nishimura Co., Ltd, and the performance test, error corrections as well as image pro- cessing of the digitizer were carried out by Shanghai Astronomical Observatory. The main structure and working mode of the digitizer are introduced in this paper. A performance test shows that brightness uniformity of illumination within the measuring area is better than 0.15%, the repeatability of digitized positions is better than 0.2 μm and the repeatability of digitized brightness is better than 0.01 instrumen- tal magnitude. The systematic factors affecting digitized positions, such as lens distortion, the actual optical resolution, non-linearity of guide rails, non-uniformity of linear motors in the mobile platform, deviation of the image mosaic, and non-orthogonality between the direction of scanning and camera lin- ear array, are calibrated and evaluated. Based on an astronomical plate with a size of 300 mm × 300 mm, which was digitized at different angles, the conversion residuals of positions of common stars on dif- ferent images were investigated. The results show that the standard deviations of the residuals are better than 0.9μm and the residual distribution is almost random, which demonstrates the digitizer has a higher precision for digitization.

An Approach of Tropospheric Correction for VLBI Phase-Referencing using GPS Data

The dominant source of error in VLBI phase-referencing is the troposphere at observing frequencies above 5 GHz. We compare the tropospheric zenith delays derived from VLBI and GPS data at VLBA stations collocated with GPS antennas. The systematic biases and standard deviations both are at the level of sub-centimeter. Based on this agreement, we suggest a new method of tropospheric correction in phase-referencing using combined VLBI and GPS data.

A new source extraction algorithm for optical space debris observation

Specific challenges arise in the task of real-time automatic data reduction of optical space debris observations. Here we present an automatic technique that optimally detects and measures the sources from images of optical space debris ob- servations. We show that highly reliable and accurate results can be obtained on most images produced by our specific sensors, and due to optimizations, the whole pipeline works fast and efficiently. Tests demonstrate that the technique performs better than SExtractor from the point of view of fast and accurate detection, therefore it is well suited for data reduction of optical space debris observations.

Preliminary results of solving the problem of geometric distortion for the 2.4 m telescope at Yunnan Observatory

We observed the open clusters NGC 1664 （43 exposures） and M35 （42 ex- posures） by the Yunnan Faint Object Spectrograph and Camera in the 2.4 m telescope at Yunnan Astronomical Observatory on 2011 January 3, and processed them by a method recently proposed by us. The result shows that there is a geometric distortion effect in the field of view and the maximum distortion is - 0.25 （i.e. 1 pixel）. After correcting the geometric distortion, the precision of stellar positional measurement is significantly improved. The best precision in each direction is 6 mas for well-exposed stars.

Eclipses and Occultations of Galilean Satellites Observed at Yunnan Observatory in 2003

We describe and analyze observations of mutual events of Galilean satellites made at the Yunnan Observatory in February 2003 from CCD imaging for the first time in China. Astrometric positions were deduced from these photometric observations by modelling the relative motion and the photometry of the involved satellites during each event.

Combination of terrestrial reference frames based on space geodetic techniques in SHAO:methodology and main issues

Based on years of input from the four geodetic techniques （SLR, GPS, VLBI and DORIS）, the strategies of the combination were studied in SHAO to generate a new global terrestrial reference frame as the material realization of the ITRS defined in IERS Conventions. The main input includes the time series of weekly solutions （or fortnightly for SLR 1983-1993） of observational data for satellite techniques and session-wise normal equations for VLBI. The set of estimated unknowns includes 3- dimensional Cartesian coordinates at the reference epoch 2005.0 of the stations distributed globally and their rates as well as the time series of consistent Earth Orientation Parameters （EOPs） at the same epochs as the input. Besides the final solution, namely SOL-2, generated by using all the inputs before 2015.0 obtained from short-term observation processing, another reference solution, namely SOL- 1, was also computed by using the input before 2009.0 based on the same combination of procedures for the purpose of comparison with ITRF2008 and DTRF2008 and for evaluating the effect of the latest six more years of data on the combined results. The estimated accuracy of the x-component and y-component of the SOL- 1 TRF-origin was better than 0.1 mm at epoch 2005.0 and better than 0.3 mm yr- 1 in time evolution, either compared with ITRF2008 or DTRF2008. However, the z-component of the translation parameters from SOL-1 to ITRF2008 and DTRF2008 were 3.4 mm and -1.0 ram, respectively. It seems that the z-component of the SOL-1 TRF-origin was much closer to the one in DTRF2008 than the one in ITRF2008. The translation parameters from SOL-2 to ITRF2014 were 2.2, -1.8 and 0.9 mm in the x-, y- and z-components respectively with rates smaller than 0.4 mmyr-1. Similarly, the scale factor transformed from SOL-1 to DTRF2008 was much smaller than that to ITRF2008. The scale parameter from SOL-2 to ITRF2014 was -0.31 ppb with a rate lower than 0.01 ppb yr-1. The external precision （WRMS） compared with IERS EOP 08 C04 of the combined EOP series was smaller than 0.06 mas for the polar motions, smaller than 0.01 ms for the UT1-UTC and smaller than 0.02 ms for the LODs. The precision of the EOPs in SOL-2 was slightly higher than that of SOL-1.

Micro-arcsecond Celestial Reference Frames:definition and realization——Impact of the recent IAU Resolutions

The adoption of the International Celestial Reference System （ICRS）, based on Very Long Baseline Interferometry （VLBI） observations of extragalactic radiosources by the International Astronomical Union （IAU） since 1998 January 1, opened a new era for astronomy. The ICRS and the corresponding frame, the International Celestial Reference Frame （ICRF）, replaced the Fundamental Catalog （FK5） based on positions and proper motions of bright stars, with the Hipparcos cat- alog being adopted as the primary realization of the ICRS in optical wavelengths. According to its definition, the ICRS is such that the barycentric directions of distant extragalactic objects show no global rotation with respect to these objects; this pro- vides a quasi-inertial reference for measuring the positions and angular motions of the celestial objects. Other resolutions on reference systems were passed by the IAU in 2000 and 2006 and endorsed by the International Union of Geodesy and Geophysics （IUGG） in 2003 and 2007, respectively. These especially concern the definition and re- alization of the astronomical reference systems in the framework of general relativity and transformations between them. First, the IAU 2000 resolutions refined the con- cepts and definition of the astronomical reference systems and parameters for Earth＇s rotation, and adopted the IAU 2000 precession-nutation. Then, the IAU 2006 resolutions adopted a new precession model that is consistent with dynamical theories; they also addressed definition, terminology or orientation issues relative to reference systems and time scales that needed to be specified after the adoption of the IAU 2000 resolutions. An additional IUGG 2007 resolution defined the International Terrestrial Reference System （ITRS） so that it strictly complies with the IAU recommendations. Finally, the IAU 2009 resolutions adopted a new system of astronomical constants and an improved realization of the ICRF. These fundamental changes have led to significant improvements in the fields of astrometry, celestial mechanics, geodynam- ics, geodesy, etc. Of special interest are the improvements in the model for variations in Earth＇s rotation, which, in turn, can provide better knowledge of the dynamics of the Earth＇s interior. These have also contributed to a significant improvement in the accuracy of the ephemerides of the solar system bodies as determined from modern measurements, with a large number of scientific applications. This paper recalls the main aspects of the recent IAU resolutions on reference systems as well as their con- sequences on the concepts, definitions, nomenclature and models that are suitable for the definition, realization and transformation of reference frames at a microarcsecond level.