利用VLBI数据确定"探测一号"卫星的轨道

双星计划的“探测一号”是中国首颗真正严格意义上的科学实验卫星,其运行轨道为中国迄今所发射的卫星中距地球最远,远地点地心距达7．8万公里．采用射电天文的VLBI技术可以对“探测一号”以及更远的深空目标,如探月飞行器实现跟踪．为了验证VLBI技术在我国探月计划中的作用,上海天文台组织了国内目前仅有的上海、乌鲁木齐和昆明3个台站对“探测一号”进行试跟踪,利用对“探测一号”约两天的VLBI观测数据,确定“探测一号”卫星的轨道,对VLBI的定轨能力做初步的探讨．按照测控部门提供的初轨 (其精度仅保证跟踪)推算的轨道与VLBI时延的拟合误差平均约2 km,时延率的拟合误差平均约15 cm／s．而利用VLBI数据定轨后的拟合程度相对于初轨有了很大的改善,结果表明,单独利用VLBI时延定轨,时延的拟合精度约5．5 m,作为外部检核的VLBI时延率的拟合精度在2 cm／s左右．单独利用VLBI时延率定轨,时延率的拟合精度约为1．3 cm／s,作为外部检核的VLBI时延的拟合精度约为29 m．而若将时延和时延率数据联合定轨,采用其内符精度加权,VLBI时延和时延率的残差分别为5．5 m和 2 cm／s．为了合理地评估VLBI定轨的真实精度,利用模拟数据进行误差协方差分析,结果表明VLBI定轨精度受动力学模型误差的影响较大,由于"探测一号”卫星的动力学模型难以精确确定,所以利用两天弧段的VLBI数据确定“探测一号”卫星轨道的位置误差为km量级,而速度误差可达cm／s量级．模拟计算还表明, VLBI和USB数据联合定轨可以大大提高定轨精度．     

The accuracy of geopotential models

Extensive tests of two recent geopotential models (GEM 7 and 8) have been made with observations not used in the solutions. Several other recent models are also evaluated. These tests show the accuracy of the satellite derived model (GEM 7, with 400 coefficients) to be about 4.3 m (r.m.s.) with respect to the global geoid surface. The corresponding accuracy of the combined satellite and surface gravimetry model (GEM 8, with 706 coefficients) is found to be 3.9m (r.m.s.). These results include a calibration for the commission errors of the coefficients in the models and an estimate of the errors from omitted coefficients. For GEM 7, the formal precision (commission errors) of the solution gives 0.7 m for the geoid error which after calibration increases to 2.4 m.

Independent observations used in this assessment include: 159 lumped coefficients from 35 resonant orbits of 1 and 9 through 15 revolutions per day, two sets of (8, 8) fields derived from optical-only and laser-only data, sets of zonal and resonant coefficients derived from largely independent sources and geoid undulations measured by satellite altimetry. In addition, the accuracy of GEM 7 has been judged by the gravimetry in GEM 8. The ratio of estimated commission to formal error in GEM 7 and 8 ranges from 2 to 5 in these tests.     

Upper-atmosphere zonal winds: Variation with height and local time

The average rotation rate of the upper atmosphere can be found by analysis of the changes in the orbital inclinations of satellites, and results previously obtained have indicated that the atmospheric rotation rate appreciably exceeds the Earth's rotation rate at heights between 200 and 400 km.We have examined all such results previously published in the light of current standards of accuracy: some are accepted, some revised, and some rejected as inadequate in accuracy. We also analyse a number of fresh orbits and, adding these to the accepted and revised previous results, we derive the variation of zonal wind speed with height and local time. The rotation rate (rev/day) averaged over all local times increases from near 1.0 at 150 km height to 1.3 near 350 km (corresponding to an average west-to-east wind of 120 m/s), and then decreases to 1.0 at 400 km and probably to about 0.8 at greater heights. The maximum west-to-east winds occur in the evening hours, 18–24 h local time: these evening winds increase to a maximum of about 150 m/s at heights near 350 km and decline to near zero around 600 km. In the morning, 4–12 h local time, the winds are east to west, with speeds of 50–100 m/s above 200 km. We also tentatively conclude that, at heights above 350 km, the average rotation rate is greater in equatorial latitudes (0–25°) than at higher latitudes.     

The determination of gravitational potential differences from satellite-to-satellite tracking

A new, rigorous model is developed for the difference of gravitational potential between two close earth-orbiting satellites in terms of measured range-rates, velocities and velocity differences, and specific forces. It is particularly suited to regional geopotential determination from a satellite-to-satellite tracking mission. Based on energy considerations, the model specifically accounts for the time variability of the potential in inertial space, principally due to earth’s rotation. Analysis shows the latter to be a significant (±1 m²/s²) effect that overshadows by many orders of magnitude other time dependencies caused by solar and lunar tidal potentials. Also, variations in earth rotation with respect to terrestrial and celestial coordinate frames are inconsequential. Results of simulations contrast the new model to the simplified linear model (relating potential difference to range-rate) and delineate accuracy requirements in velocity vector measurements needed to supplement the range-rate measurements. The numerical analysis is oriented toward the scheduled Gravity Recovery and Climate Experiment (GRACE) mission and shows that an accuracy in the velocity difference vector of 2×10⁻⁵ m/s would be commensurate within the model to the anticipated accuracy of 10⁻⁶ m/s in range-rate. This revised version was published online in July 2006 with corrections to the Cover Date.     

An Orbital Prediction Algorithm for LEO Satellites Based on Optical Observations of Short Arcs at Single Station

The forecasting technique of the target tracking based on the short arcs at single station is an important way to guarantee that high-precision photoelectric theodolites can normally track and capture the targets in unconventional environments. We construct the tracking prediction algorithm based on nonlinear filter, which can provide the guiding data for the closed loop tracking under normal circumstances. At the same time we also construct the target prediction algorithm based on the nonlinear transformation, without valid observational data, which can provide a track guidance for the theodolite and ensure that the targets will not be lost. It is demonstrated that the nonlinear filtering is more effective than the EKF (extended Kalman filter) in the tracking prediction algorithm of the short arcs at single station. The results indicate that the nonlinear filter designed in this paper can be used as the guiding algorithm for the optical tracking equipments. And its guiding accuracy is in the same order of magnitude of the theodolite's random measurement accuracy. When the systematic error of the equipments reaches 50″, the accuracy can achieve 20″ for predictions in 60 s. This still satisfies the requirement of the field of view of the tracking equipments.     

CCD图像数字定心算法的比较

从CCD采集的数字图像中提取天体的位置和光度信息对于天文研究具有基础意义.其中,恒星星像中心像素位置的精确测量对于天体测量至关重要.针对国际上常用的三种位置测量算法:修正矩方法、Gauss拟合法和中值法,利用实测的CCD图像进行了试验研究.具体地,采用云南天文台1m望远镜观测的CCD图像同时使用这些算法进行了实际测量,根据每幅图像中测得的像素位置,对这三类定心算法进行了比较和精度分析.可视化方法被用于不同阈值的选取和设定比较.实验数据表明,无阈值的二维Gauss拟合方法是一种精度相对较高的定心算法.     

Research on the panel adjustment method of an active main reflector for a large radio telescope

The surface accuracy of a large parabolic antenna is an important indicator to evaluate the quality of the antenna.It not only directly affects the antenna’s aperture efficiency,thereby determining the shortest wavelength that the antenna can work,but also affects the main lobe width and side lobe structure of the antenna pattern.Microwave holography is an important method for parabolic antenna profile detection.In this article we adopt a new algorithm to adjust the panels for the large radio telescope with an active main reflector through the TM65 m antenna’s aperture phase profile.The panels of the TM65 m radio telescope is in a radial pattern with 14 rings.Each corner of the panel is fixed on the screw of the actuator to move up and down,and the adjacent corners of the four panels share an actuator.We use the method of plane fitting to calculate the adjustment value of every panel’s corner.But one actuator,which simultaneously controls the common corner of the adjacent panels,will have different adjustment values according to the different plane fitting equation based on adjacent panels.In this paper,the adjustment value of the adjacent panels’crosspoints are constrained to be equal to the constraint condition to calculate each actuator’s adjustment value of the TM65 m radio telescope.Through multiple adjustments and application of the new algorithms,the accuracy of the TM65 m antenna reflector profile has been improved from the original 0.28 mm to the current 0.19 mm.     

Gauss's method for secular dynamics, softened

We show that the algorithm proposed by Gauss to compute the secular evolution of gravitationally interacting Keplerian rings extends naturally to softened gravitational interactions. The resulting tool is ideal for the study of the secular dynamical evolution of nearly Keplerian systems such as stellar clusters surrounding black holes in galactic nuclei, cometary clouds or planetesimal discs. We illustrate its accuracy, efficiency and versatility on a variety of configurations. In particular, we examine a secularly unstable system of counterrotating discs, and follow the unfolding and saturation of the instability into a global, uniformly precessing, lopsided  ( m = 1)  mode.     

Orbit determination of satellite “Tance 1” with VLBI data

The satellite “Tance 1” of the “Double-Star Program” is the first truly scientific experimentation satellite of China. Its orbit is the farthest so far launched in China, with a geocentric apogee reaching 78 thousand kilometers. The tracking of “Tance 1” and of more distant space targets, such as the lunar exploration craft, can be realized with the VLBI technique of radio astronomy. In order to test and verify the role which the VLBI technique plays in the lunar exploration program of China, Shanghai Astronomical Observatory organized the only 3 tracking stations in China (located at Shanghai, Urumqi and Kunming), to carry out test tracking of “Tance 1,” and used the time delay data obtained to determine the orbit of “Tance 1” over a two-day period, so providing a preliminary assessment of the possibility of VLBI orbit determination. The fitting error of the orbit so obtained is about 5.5 m in the time delay and about 2 cm/s in the delay rate (this for checking only), much better than is provided by the preliminary orbit (used merely for ensuring tracking) in which the corresponding figures are around 2 km and 15 cm/s. Further, if the orbit is determined by using both the time delay and time delay rate data (with weights according to their internal accuracies), then the residuals are 5.5 m in the time delay and 2 cm/s in the delay rate. For an appreciation of the true accuracy of the VLBI orbit determination, we used simulation data (of the observed two-day VLBI data) and found the results depended greatly on the error in the dynamic model of the satellite which, however, is difficult to assess, while the formal residuals are of the order of 1 kin in the delay and of cm/s in the delay rate. The simulation computation also indicates that a joint determination using both VLBI and USB data will have an improved accuracy.     

NAHUAL: a near‐infrared high‐resolution spectrograph for the GTC optimized for studies of ultracool dwarfs

We present the status of an ongoing study for a high‐resolution near‐infrared echelle spectrograph for the 10.4‐m GTC (Gran Telescopio de Canarias) which will soon start operating at the Observatorio del Roque de los Muchachos on the island of La Palma. The main science driver of this instrument, which we have baptized NAHUAL, is to carry out a high precision radial velocity survey of exoplanets around ultracool dwarfs. NAHUAL is being especially designed to achieve the highest possible accuracy for radial velocity measurements. The goal is to reach an accuracy of a few m/s. It is thus required that the instrument is cross‐dispersed and that it covers simultaneously a wide wavelength range. Absorption cells will be placed in front of the slit which will allow a simultaneous self‐reference similar to an iodine‐cell in the optical regime. It is planned to place the instrument at one of the Nasmyth platform of the GTC behind the Adaptive Optics system. Our current design reaches a maximum spectral resolution of λ/Δλ = 50000 with a slit width of 0.175 arcsec, and gives nearly complete spectral coverage from 900 to 2400 nm. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Phase retrieval holography and surface adjustment of SEST using the LES-8 satellite

Phase retrieval holography has been used to determine the surface profile of the Swedish-ESO Sub-millimetre Telescope (SEST). The measurement was performed using a signal transmitted from the geosynchronous LES-8 satellite over a range of elevation angles. The results of the measurements have been used to reset the reflector thus reducing the surface errors from the initial 80 m r.m.s. to 51 m r.m.s. The improvement in the surface accuracy has been confirmed by telescope efficiency measurements at 232–250 GHz.     

Comments about the direct perturbations of Venus and Mars on the Moon's motion

In a previous paper (Standaert, 1980) we have described an algorithm to compute the direct perturbation of the planets on the Moon's motion. A short summary of this algorithm is presented in Section 2 of this paper. Our first results permit us to present some complements and comments about these computations.The algorithm is based upon the Lie transform method and is implemented using Chapront's ELP as solution of the main problem with the partial derivatives of Henrard's Semi-Analytical Lunar Ephemeris (SALE), and Bretagnon's mean Keplerian orbit.An analysis of truncation errors in intermediate results is presented including the resonance effects. The final accuracy of the solution is intended to be about 0.0005 for terms of period up to 2000 yr in the case of Venus and up to 5000 yr in the case of Mars.The effects of second-order terms in the masses are investigated. Only those depending upon the second derivatives of the mean motions are found to be significant to the given accuracy and are included.Proceedings of the Conference on Analytical Methods and Ephemerides: Theory and Observations of the Moon and Planets. Facultés universitaires Notre Dame de la Paix, Namur, Belgium, 28–31 July, 1980.     

SARG: The High Resolution Spectrograph of TNG

SARG is a cross dispersed echelle spectrograph in operation since late spring 2000 at the Italian Telescopio Nazionale Galileo (TNG) 3.5 m telescope, La Palma. SARG offers both single object and long slit (up to 26 arcsec) observing modes covering a spectral range from λ = 0.37 up to1 μm, with resolution ranging from R = 29,000 up to R = 164,000. Cross dispersion is provided by means of a selection of four grisms; interference filters may be used for the long slit mode (up to 26 arcsec). A dioptric camera images the cross dispersed spectra onto a mosaic of two 2048 × 4096 EEV CCDs (pixel size: 13.5 μm) allowing complete spectral coverage at all resolving power for λ < 0.8 μm. In order to reach a high wavelength calibration precision an iodine-absorbing cell is provided. A Distributed Active Temperature Control System (DATCS) maintains constant the temperature of all spectrograph components at a preset value. Early results show that SARG works according to original specifications in terms of wavelength coverage, efficiency (measured peak efficiency is about 13%),resolution (maximum resolution R = 164,000 using a 0.3 arcsec slit, R ∼144,000 using an image slicer), and stability (preliminary estimates of radial velocity accuracy is ∼3 m/s using the iodine cell and ±150 m/s without the iodine cell). This revised version was published online in July 2006 with corrections to the Cover Date.     

基于SLR精密轨道的天文定位精度分析

利用全球卫星激光测距服务系统(ILRS,International Laser Ranging Service)标准点资料对Ajisai卫星进行精密定轨,残差均方根(RMS)优于3 cm,得到该星的精密轨道.进而对长春站40 cm空间碎片光电望远镜获得的Ajisai卫星的天文定位资料进行精度分析,外符合精度约3″左右.单独利用天文定位数据进行轨道改进,内符合精度优于3″.改进轨道的x、y、z坐标3分量在观测数据覆盖范围内的精度在100 m之内.同样地对Jason-1卫星作数据分析,结果和Ajisai卫星精度相当.分析各个弧段的精度变化,发现定标星个数减少,会导致天文定位精度下降.据此提出可以把最少定标星比例作为评定数据质量的参考指标之一.     

Concatenation of HRSC colour and OMEGA data for the determination and 3D-parameterization of high-altitude CO2 clouds in the Martian atmosphere

We used Mars Express HRSC and OMEGA data to investigate mesospheric cloud features observed in the equatorial belt of Mars from December 2007 until early March 2008. This period corresponds to early northern spring of Martian year 29. The reflection peak at 4.26 μm in OMEGA data identifies the clouds as CO₂ ice clouds. HRSC observed the clouds together with OMEGA in five orbits. Cloud features are most prominent in the shortwave HRSC colour channels with wavelength centers at 440 and 530 nm, but rarely visible in all other channels. In the period of Ls 0-36°, OMEGA and HRSC together detected mesospheric CO₂ ice clouds in 40 orbits. They occur in a latitude belt of ±20° around the equator and at longitudes between 240°E (Tharsis) in the West and 30°E (Sinus Meridiani) in the East. The clouds were observed between 3 and 5 p.m. local time with mainly ripple-like to filamentary cloud forms. The viewing angles of the HRSC blue and green colour channels differ by 6.6° and the resulting parallax can be used to directly measure cloud heights by means of ray intersection. 17 HRSC data takes were found to exhibit clouds with heights from 66 to 83 km with an accuracy of 1-2 km. The pushbroom imaging technique also yields a time delay for the two observations in the order of 5-15 s close to periapsis, and therefore time-related cloud movements can be detected. A method was developed to determine the across-track cloud displacements, which can directly be translated to wind velocities. Zonal cloud movements could be measured in 13 cases and were oriented from East to West. Related wind speeds range between 60 and 93 m/s with an accuracy of 10-13 m/s.     

Potential for stratospheric Doppler windspeed measurements of Jupiter by sub-millimetre spectroscopy

The sub-millimetre/microwave range of the spectrum has been exploited in the field of Earth observation by many instruments over the years and has provided a plethora of information on atmospheric chemistry and dynamics—however, this spectral range has not been fully explored in planetary science, having been exclusively employed to carry out ground-based measurements. To this end, a sub-millimetre instrument, the Orbiter Terahertz Infrared Spectrometer (ORTIS), is studied by the University of Oxford and the Rutherford Appleton Laboratory, to meet the requirements of the European Space Agency's Cosmic Visions 2015-2025 programme—in particular, the Europa Jupiter System Mission (EJSM), which has the European Space Agency and the National Aeronautics and Space Administration as partners. ORTIS is designed to measure atmospheric temperature, the abundance of stratospheric water vapour and other jovian gases, and is intended to be capable of retrieving vertical profiles of horizontal windspeed in the stratosphere for the first time, from Doppler-shifted emission lines measured at high spectral resolution. In this work, a preliminary study and implementation of the estimation of windspeed profiles on simulated spectra representative of Jupiter is presented, detailing the development of the retrieval algorithm, showing that a sub-millimetre instrument such as ORTIS should be able to retrieve windspeed profiles to an accuracy of about 15 m/s between 70 and 200 km/0.1-10 mb using a single near-limb measurement, for expected noise amplitudes.     

Experimental Research of First Feed-Support System of FAST

First feed support system, which is one of the key technologies of Five hundred meter Aperture Spherical Telescope (FAST) was discussed in this paper. A 1:30 scaled model (20 m × 20 m × 6 m) has been built to make systemic research on the feasibility of the system. The root mean square of pointing error of moving cable car was measured and the displacements of cable car within the wind speed of 8 m/s, which is the maximal wind speed of the place where the FAST will be located, were evaluated. Then the total maximal pointing error of cable car was computed. Experimental results demonstrate that the total pointing error of the cable car is less than 50 cm, the required pointing accuracy error of the system. This revised version was published online in July 2006 with corrections to the Cover Date.     

An Optical,Dual-Beam,Automated Medium Resolution Spectropolarimeter for the Vainu Bappu Telescope

Fabrication of an optical dual-beam spectropolarimeter as an add-on facility to an existing astronomical spectrograph at the Vainu Bappu Telescope is described. The polarimetric optics consists of a superachromatic Pancharatnam design half-wave plate and a modified Glan-Taylor polarizing beam splitter. Instrumental response, calibration and characterization of the system are presented. Performance of the spectropolarimeter has been assessed based on the results of observations of polarized and unpolarized standard stars. The attainable accuracy of the instrument is found to be dependent on the S/N of the data. The present data set yields an accuracy of ±0.5% at 4000 Å and ±0.3% at 7500 Å, at a spectral resolution of 7.2 Å. Analysis of spectropolarimetric data is systematically covered for removal of instrumental errors. The spectropolarimetric reduction software (SPRS), extremely versatile, user friendly and compatible with the IRAF image processing package, was developed for reducing spectropolarimetric data. The empirical relation by Serkowski (IAU Symposium 52, Interstellar Dust and Related Topics, J. M. Greenberg and H. C. van de Hulst (eds.), Dordrecht, Reidel, 1975, p. 145) for wavelength dependence of polarization due to interstellar medium has been fitted to the data for the star HD 43384 (9 Gem). Our result shows a polarization larger by 0.49% from Hsu and Berger (ApJ 262, 1982, 732). We attribute this difference to the long term variation in P for this star.     

平均风作用对天马望远镜面形精度和指指向精度的影响

天马望远镜的最高工作频段为43 GHz。为保证高质量的观测结果,需要研究风载荷对天线精度的影响。首先对观测站实测的风速风向数据做了统计分析,结果显示:10 m高度处10min时距平均风速小于4 m·s^-1的占比超过80%,主导风向为北-西北方向。然后,通过将倾斜仪实测结果与有限元模拟结果进行对比,验证了模拟的有效性,并进一步分析了在不同迎风姿态、不同风速下天线结构的平均风荷载响应,以及天线面形精度和指向精度的变化。结果表明,平均风荷载对天线指向精度,尤其是俯仰角指向精度的影响较大,对面形精度的影响较小;在弹性范围内,天线面形精度和指向精度与风速间均为二次关系。研究结果可为天线面形精度和指向精度的评估提供参考。     

基于FPGA的空间太阳磁像仪自校正稳像系统研制

太阳磁像仪是开展太阳磁场观测研究的核心仪器,其中的稳像系统是空间太阳磁像仪的关键技术之一,针对深空探测卫星系统对载荷重量、尺寸限制严苛的要求,设计了基于图像自校正方法的稳像观测系统.介绍了一套基于现场可编程门阵列(Field-Programmable Gate Array, FPGA)和数字信号处理器(Digital Signal Processor, DSP),通过基于自相关算法的高精度稳像方法设计,并结合精确偏振调制、准确交替采样控制等系统软硬件设计,克服由于卫星平台抖动、指向误差等因素造成的图像模糊,实现实时相关、校正、深积分的稳像观测系统.针对像素尺寸为1 K×1 K、帧频为20 fps的CMOS (Complementary Metal Oxide Semiconductor)探测器,实现了1像元以内的实时稳像观测精度.在完成实验室测试后, 2021年6月18日在国家天文台怀柔太阳观测基地35 cm太阳磁场望远镜上开展了实测验证,结果表明该系统能够有效地完成太阳磁像仪自校正稳像观测,获得了更高分辨率的太阳磁场数据.稳像系统的成功研制不仅可以为深空太阳磁像仪的研制提供轻量化、高...