月球磁异常区的光学特性探讨

月球无全球性磁场,在某些大型撞击坑的对峙区域存在较强的磁场强度。这些磁场异常区域往往与光学异常区域呈地理相关,推断与磁异常诱发的微磁层弱化空间风化效应有关。嫦娥一号卫星携带的干涉成像光谱仪获取了月球表面可见光-近红外波长范围内32个波段的光谱,为研究月球空间风化效应提供了一个可能物理量——月表反射率。通过对几个常见磁异常结构磁场内外区域的平均光谱对比发现,磁场内区域反射率一般高于磁场外区域。这一结果对月球表面空间风化效应的研究产生了新的推动。     

基于空间光学成像的4179号小行星表面反射特性研究

2012年12月13日,嫦娥二号在距地球约700万公里的深空,成功飞掠4179号小行星,获得了最高分辨率优于3 m的系列可见光图像.使用最小二乘拟合方法,对成像相机辐射定标数据进行处理,获取绝对定标系数和相对定标校正矩阵,校正原始小行星图像及定标数据,同时反演成像时刻小行星表面辐射亮度.根据Hapke对Nicodemus反射率定义标准在行星科学应用中的具体描述,获取小行星半球反照率.在R、G、B谱段上,小行星表面平均半球反照率分别为0.208 3、0.126 9和0.134 6.小行星表面半球反照率为0.156 6.首次使用基于空间光学成像的方法获得了4179号小行星表面光谱反照率,并获得了4179号小行星的表面反照率分布图,对行星科学研究有应用价值.     

小倾角LEO多星天基平台光学跟踪GEO精密定轨

针对地基卫星测控系统(Tracking Telemetry and Command, TT&C)系统对地球静止轨道(Geostation-\lk ary Earth Orbit, GEO)卫星在空间和时间覆盖上的局限性, 提出小倾角低地球轨道(Low Earth Orbit, LEO)多星组网天基平台对GEO卫星进行跟踪定轨的方法. 根据空间环境和光学可视条件对仿真数据进行筛选以模拟真实的观测场景, 利用光学测角数据, 使用数值方法对GEO卫星的轨道进行确定. 结果与参考轨道进行重叠对比, 在平台轨道精度5 m、测量精度5rq\rq、 定轨弧长12 h的情况下, 两颗LEO卫星对GEO卫星进行跟踪定轨的精度可达到千米量级, 4颗LEO卫星对GEO目标进行跟踪定轨的精度可达到百米量级. 随着LEO组网卫星数量的增加, 定轨精度得到了较大的提高.     

1m新真空太阳望远镜成像系统观测数据干涉条纹特性分析及消除方法

抚仙湖1 m新真空太阳望远镜(New Vacuum Solar Telescope,简称NVST)是我国新一代地基太阳观测设备之一,其中Hα成像系统是进行太阳色球观测的主要终端.实测Hα观测数据(特别是偏带观测数据)中明显存在由薄膜干涉引起的等厚条纹,而常规平场校正有时无法将其消除,并在后期图像高分辨率重建过程中造成严重影响.为了解决这个问题,利用连续4 h的Hα偏带观测数据对干涉条纹的空间和时间分布特性做了分析并发现:干涉条纹的间距或二维空间分布样式随时间固定不变,但条纹的强度随时间变化明显(4 h左右可见度增加9倍),是导致常规平场校正无法将其消除的主要原因.由此推断造成上述时变特性的主要原因是观测过程中入射光的强度以及入射光与CCD靶面位置的相对角度随时间发生了变化,而且产生条纹的薄膜结构与靶面位置距离较近.尝试了两种降低条纹可见度的方法.首先对探测器(pco.4000型CCD)光敏介质的前封窗进行了改造,通过改造前封窗楔角至2?的方案来调制条纹间距至CCD像元尺度,改造之后小范围视场内可识别出可见度仅为0.6%的少量干涉条纹.其次基于图像滤波技术提取干涉条纹图样并生成所谓条纹平场,分别采用和比较了频域傅里叶变换滤波和空域中值滤波两种方法,结果基本相同.观测数据再次进行条纹平场校正之后,条纹消除效果明显:可见度可降低8倍左右(即由之前的4.7%降低为0.6%以下).同时也指出多幅图像的积分(累加)可以有效降低精细太阳结构在条纹平场提取中的影响,并给出累加时间的经验值约为20 min.     

基于精密光程差调制的时域干涉信号闭合相位检测方法研究

闭合相位法是实现长基线恒星光干涉高分辨成像的重要技术手段之一,获得精确的闭合相位信息是进行光干涉图像重构的先决条件.提出一种基于精密光程差调制的时域干涉信号闭合相位检测方法,在3路干涉臂上进行非冗余精密光程调制,并通过多次干涉测量结合数据拟合的方法消除光程差调制中存在的正弦误差,使得光程调制的精度达到20 nm以内.引入高速探测器件提升时域干涉信号的采样频率,对探测器上获得的时域干涉信号进行傅立叶变换处理,获得3路干涉臂精确的闭合相位信息.室内实验结果表明,基于精密光程调制的时域信号闭合相位计算精度可以达到1/50波长以内.     

GPS星历对LEO星载GPS精密定轨精度的影响

目前,越来越多的低轨卫星上都搭载了用于精密定轨的星载GPs接收机,星载GPS已成为低轨卫星精密定轨的主要手段之一.星载GPS精密定轨精度依赖于GPS星历及钟差精度.基于SHORDE-Ⅲ非差动力学定轨功能,以2005年8月1日至8月7日一周的GRACE卫星实测数据为例,采用事后精密轨道(igs)、快速轨道(igr)和超快速轨道(igu)三种GPS星历在同等条件下定轨,估计GPS星历精度对低轨卫星定轨精度的影响,实际计算结果表明igs和igr两类GPS星历定轨精度相当,约为9.5 cm,igu星历定轨精度略低于igs和igr星历,约为10.5cm:高频GPS卫星钟差数据对定轨精度会产生1-6cm影响.     

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

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

环火卫星运动的坐标系附加摄动及相应坐标系的选择

与处理地球卫星相关问题类似,在研究和处理环火卫星(尤其是低轨卫星)的轨道问题时,宜采用火心历元平赤道坐标系,即火心天球坐标系,其xy坐标面和x轴方向就是相应的平赤道面和平春分点方向.与地球的岁差章动现象类似,在该坐标系中,火星赤道面在空间的摆动同样会引起坐标系附加摄动.采用类似对地球岁差章动的处理方法,在一定精度前提下,基于IAU2000火星定向模型,处理了火星赤道面摆动中的岁差效应,并在此基础上,研究岁差对环火卫星轨道的影响,给出了相应的火星非球形引力位的变化及其导致的卫星轨道的坐标系附加摄动解,其表达形式简单,引用方便.与高精度数值解的比对表明,该分析解能够满足通常的精度要求.因此,在处理环火卫星(即使是低轨卫星)轨道及其相关问题时,可以采用统一坐标系:火心天球坐标系.而不必像当初处理地球卫星那样,为了避免计算坐标系附加摄动而引进一种混合型赤道坐标系,即采用瞬时真赤道面和历元平春分点方向作为其xy坐标面和x轴方向.在统一坐标系的选择下,实际工作中就不会存在坐标系转换的麻烦.     

连线干涉系统测量共位双星宽带信号处理与分析

随着人类空间活动的急剧增加,同步轨道资源已变得十分稀缺。为此人们采用了多星共位技术,来充分挖掘同步轨道资源。多星共位中,由于卫星存在位置漂移的情况,安全问题成为维护卫星正常工作必须考虑的问题之一。考虑到同步卫星的轨道特性,常采用三程测距的方式对卫星进行测轨,而且需要较长时间的测轨弧段。连线干涉测量可以获得高精度的差分相位时延,这对卫星在垂直视线方向上构成了很强的约束。给出了连线干涉测量系统设计及测量数据处理方法,并对测量中的关键技术与设备特性、测量数据精度进行了初步分析。     

嫦娥二号伽玛射线谱仪

嫦娥二号卫星是我国继嫦娥一号卫星成功发射后的第2颗探月卫星,其上搭载的伽玛射线谱仪(Gamma-Ray Spectrometer,GRS)的主要科学目标是探测月表O、Si、Fe、Ti、U、Th、K、Mg、Al、Ca等主量元素.相比嫦娥一号伽玛射线谱仪,嫦娥二号的能量分辨率和探测效率都大大提高.描述了嫦娥二号伽玛射线谱仪的设计和性能测试结果以及在轨飞行的初步探测成果.     

Photometric correction and in-flight calibration of Chang’ E-1 Interference Imaging Spectrometer (IIM) data

The main objective of this study is to develop a new photometric correction that is suitable for global Chang’ E-1 (CE-1) Interference Imaging Spectrometer (IIM) data. We considered two improvements in the accuracy of the photometric correction: (1) classifying the whole Moon’s surface into four classes (very bright rays, mature highlands, low FeO basalts, and high FeO basalts) based on the FeO contents; (2) decoupling the combined effects of the solar photometry and the detector responsivity along the cross-track. The results showed that our decoupling method could correct the non-uniformity of detector response but retain the cross-track photometry. Based on these in-flight calibrated data, spectrally continuous photometric functions of the four classes were obtained. By eliminating the contamination from highlands ejecta and low FeO basalts, the opposition surge of the high FeO maria was evident for the first time. The amplitude and width of the opposition surge presents wavelength dependence, which suggests that shadow hiding is the mechanism responsible for the opposition surge. A global lunar reflectance map whose boundaries between adjacent orbits are invisible was generated using the new four-type photometric functions. The quality of the global mosaic and the consistency of the reflectance spectra of the same area obtained in different orbits indicate the effectiveness of our photometric correction method. Although our method was specifically applied to the IIM data, the photometric correction method and the parameter values derived in this study is expected to be suitable for other optical instruments.     

A preliminary experience in the use of Chang’E-1 IIM data

The Interference Imaging Spectrometer (IIM) onboard the first lunar satellite of China, Chang’E-1, has acquired 84% of the area between south and north latitude 70°. To contribute to its usability, this paper presents our preliminary experience in the use of IIM data. Firstly, we provide one practicable method for the on-orbit correction of the inhomogeneity of sensor response. Secondly, aiming at the problem that the spectral range of IIM does not cover the absorption peak of the mafic mineral completely, we explore a method to approximate the absorption band center for IIM data. A strong correlation between the absorption band center and the wavelength at which the first derivative equals to 0 (i.e., stagnation point) was revealed. Based on the corrected data and the correlation, the absorption band center of several large craters was mapped. The distribution of rocks and minerals shown in the map of absorption band center for Aristarchus and Copernicus is in agreement with previous studies but with much finer structure. Horizontal and vertical lithologic diversity was detected in Zucchius crater. This paper demonstrates the potential of IIM data for the identification of lunar rocks due to its high spatial and spectral resolution. In a future study we will produce a global map of the absorption band center with greater accuracy and it is expected that this global map will provide complementary information for other hyperspectral data such as SP on KAGUYA or M³ on Chandrayaan-1.     

倾角函数展开及其在分析法轨道预报中的应用

田谐项摄动是分析法轨道预报中的重要部分,其中包含大量倾角函数及其偏导数的计算.由于具有精度更高、速度更快的优点,倾角函数一般通过递推方法计算.以文献中提出的改进Gooding方法为基础,将其给出的程序稍加改进,在计算2–50阶倾角函数时缩短了约24%的计算时间.考虑到分析法预报过程中轨道平倾角变化很小,以泰勒展开式计算倾角函数,可极大提高计算速度,较大程度地减小分析法预报耗时,且引力场阶次越高,减小幅度越大,取50阶时预报耗时缩短了48%.另一方面,以2阶展开式计算倾角函数时,与改进Gooding法相比,分析法预报星历偏差很小.对于500 km高度的低轨卫星,分别以改进Gooding法和2阶泰勒展开式计算倾角函数,预报3天,当地球引力场阶次不高于50时,二者预报星历偏差RMS (Root Mean Square)低于1 mm,且随着轨道高度的增加,预报星历偏差RMS逐渐减小.     

Saturn's icy satellites investigated by Cassini-VIMS: I. Full-disk properties: 350-5100 nm reflectance spectra and phase curves

Saturn's icy satellites are among the main scientific objectives of the Cassini-VIMS (Visual and Infrared Mapping Spectrometer) experiment. This paper contains a first systematic and comparative analysis of the full-disk spectral properties of Dione, Enceladus, Epimetheus, Hyperion, Iapetus, Mimas, Phoebe, Rhea and Tethys as observed by VIMS from July 2004 to June 2005. The disk integrated properties (350-5100 nm reflectance spectra and phase curves at 550-2232 nm) and images of satellites are reported and discussed in detail together with the observed geometry. In general, the spectra in the visible spectral range are almost featureless and can be classified according to the spectral slopes: from the bluish Enceladus and Phoebe to the redder Iapetus, Hyperion and Epimetheus. In the 1000-1300 nm range the spectra of Enceladus, Tethys, Mimas and Rhea are characterized by a negative slope, consistent with a surface largely dominated by water ice, while the spectra of Iapetus, Hyperion and Phoebe show a considerable reddening pointing out the relevant role played by darkening materials present on the surface. In between these two classes are Dione and Epimetheus, which have a flat spectrum in this range. The main absorption bands identified in the infrared are the 1520, 2020, 3000 nm H₂O/OH bands (for all satellites), although Iapetus dark terrains show mostly a deep 3000 nm band while the 1520 and 2020 nm bands are very faint. In this spectral range, the Iapetus spectrum is characterized by a strong reddening. The CO₂ band at 4260 nm and the Fresnel ice peak around 3100 nm are evident only on Hyperion, Phoebe and Iapetus. The phase curves at 550 and at 2232 nm are reported for all the available observations in the 0°-144° range; Rhea shows an opposition surge at visible wavelengths in the 0.5°-1.17° interval. The improvement on the retrieval of the full-disk reflectance spectra can be appreciated by a direct comparison with ground-based telescopic data available from literature. Finally, data processing strategies and recent upgrades introduced in the VIMS-V calibration pipeline (flat-field and destriping-despiking algorithm) are discussed in appendices.     

Accurate Mars Express orbits to improve the determination of the mass and ephemeris of the Martian moons

The determination of the ephemeris of the Martian moons has benefited from observations of their plane-of-sky positions derived from images taken by cameras onboard spacecraft orbiting Mars. Images obtained by the Super Resolution Camera (SRC) onboard Mars Express (MEX) have been used to derive moon positions relative to Mars on the basis of a fit of a complete dynamical model of their motion around Mars. Since, these positions are computed from the relative position of the spacecraft when the images are taken, those positions need to be known as accurately as possible. An accurate MEX orbit is obtained by fitting two years of tracking data of the Mars Express Radio Science (MaRS) experiment onboard MEX. The average accuracy of the orbits has been estimated to be around 20–25 m. From these orbits, we have re-derived the positions of Phobos and Deimos at the epoch of the SRC observations and compared them with the positions derived by using the MEX orbits provided by the ESOC navigation team. After fit of the orbital model of Phobos and Deimos, the gain in precision in the Phobos position is roughly 30 m, corresponding to the estimated gain of accuracy of the MEX orbits. A new solution of the GM of the Martian moons has also been obtained from the accurate MEX orbits, which is consistent with previous solutions and, for Phobos, is more precise than the solution from the Mars Global Surveyor (MGS) and Mars Odyssey (ODY) tracking data. It will be further improved with data from MEX-Phobos closer encounters (at a distance less than 300 km). This study also demonstrates the advantage of combining observations of the moon positions from a spacecraft and from the Earth to assess the real accuracy of the spacecraft orbit. In turn, the natural satellite ephemerides can be improved and participate to a better knowledge of the origin and evolution of the Martian moons.     

Deconvolution of ASCA X-ray data – I. Spectral-imaging method

In this paper we describe a self-contained method for performing the spectral-imaging deconvolution of X-ray data on clusters of galaxies observed by the ASCA satellite. Spatially resolved spectral studies of data from this satellite require such a correction, because its optics redistribute photons over regions that are of comparable size to the angular scales of interest in clusters. This scattering is a function not only of spatial position but also of energy. To perform a correction for these effects we employ maximum-likelihood deconvolution of the image (within energy bands of 1 keV) to determine the spatial redistribution, followed by a Monte Carlo energy reassignment of photon energies with position to determine the spectral redistribution. We present tests on simulated cluster data, convolved with the various instrumental characteristics and the X-ray background, which show that our methodology can successfully recover a variety of intrinsic temperature profiles in typical observational circumstances. In Paper II (this issue) we shall apply our spectral-imaging deconvolution procedure to a large sample of galaxy clusters to determine temperature profiles.     

Speckle observations with PISCO in Merate: IV. Astrometric measurements of visual binaries in 2005

We present relative astrometric measurements of visual binaries made during the second semester of 2005, with the speckle camera PISCO at the 102 cm Zeiss telescope of Brera Astronomical Observatory, in Merate. Our sample contains orbital couples as well as binaries whose motion is still uncertain. The purpose of this long term program is to improve the accuracy of the orbits and determine the masses of the components. We performed 130 new observations of 120 objects, with most of the angular separations in the range 0″.1–4″, and with an average accuracy of 0″.01. Most of the position angles could be determined without the usual 180. ambiguity with the application of triple‐correlation techniques, and their mean error is 0°8. We have found a possible new triple system: ADS 11077. The measurements of the closest binaries were made with a new data reduction procedure, based on model fitting of the background of the auto‐correlations. As this procedure proved to be very efficient, we have re‐processed the old observations of close binaries made with PISCO in Merate since 2004. We thus improved 20 measurements already published and obtained 7 new measurements for observations that were previously reported as “unresolved”. We finally present revised orbits for ADS 684, MCA 55Aac (in the Beta 1 Cyg–Albireo multiple system) and ADS 14783 for which the previously published orbits led to large residuals with our measurements and for which the new observations made since their computation allowed a significant improvement of those old orbits. The sum of the masses that we derived for those systems are consistent with the spectral type of the stars and the dynamic parallaxes are in good agreement with the parallaxes measured by Hipparcos. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Astrometric studies of the results of a new reduction of old photographic observations of the Saturnian System based on the comparison with the modern theories of satellite motion

The paper shows the possibility of increasing the accuracy of the results of photographic observations of Saturn and its moons made in the 1970s and reduced using the old reference star catalogues and semiautomatic measurements. New celestial coordinates of the moons (from the third to the eighth), “satellite minus satellite” relative moon coordinates, and Saturn coordinates by positions of satellites are obtained without measuring its images. The results are stored in the Pulkovo Observatory database on the Solar System bodies and are available online at www.puldb.ru. The efficiency of the reduction method based on digitizing of astronegatives using 21 Mpx Canon digital camera and IZMCCD software is shown. The comparison of new results of old observations with the latest theories of moon motion has revealed a significant increase in satellite positioning accuracy. The investigation of the differences (O–C) of celestial coordinates from satellite positions in their apparent Saturn-centric orbits has revealed a noticeable motion of the differences (O–C) in right ascension depending on their distances from Saturn for all moons.     

China's Lunar Exploration Program: Present and future

China launched its first lunar probe, Chang’E-1, at 6:05 p.m. (10:05 GMT), October 24, 2007. Chang’E-1 blasted off on a Long March 3A carrier rocket from the No. 3 launch tower in the Xichang Satellite Launch Center of southwestern Sichuan Province. China National Space Administration performed the lunar orbit injection maneuver for Chang’E-1 at 11:25 a.m. on November 5, 2007 (China Standard Time). Chang’E-1 was injected into the lunar orbit after the maneuver, and will begin to explore the moon in the following 1 year. It is the first step into its ambitious three-phase moon mission, marking a new milestone in the Chinese space exploration history.     

Mapping global lunar abundance of plagioclase,clinopyroxene and olivine with Interference Imaging Spectrometer hyperspectral data considering space weathering effect

Mineral mapping of the lunar surface is critical to understanding the Moon’s geological diversity and history, yet the global lunar abundance of minerals has not been mapped using hyperspectral data. The Interference Imaging Spectrometer (IIM) of Chang’E-1 mission obtained hyperspectral data of the global lunar surface within the wavelength of 480–960 nm in which major minerals can be discriminated by faint differences in 32 contiguous hyperspectral bands. The effect of space weathering produces multiple endmembers of lunar minerals by obscuring the pure spectra of minerals in different levels. In this study, the distributions of plagioclase, clinopyroxene and olivine on the global lunar surface were mapped with IIM hyperspectral data based on the modified Multiple Endmember Spectral Mixture Analysis (MESMA) method considering the space weathering effect. The distribution of lunar space weathering levels was retrieved as a byproduct of mineral mapping. The mineral mapping results were compared with recent mapping results. Although the wavelength of IIM is limited, it shows that our results are basically consistent with the recent research at both global and local scales. The distribution of space weathering levels is also consistent with the map of optical maturity parameter (OMAT) in most parts of the global lunar surface, especially in the highlands. This study demonstrates that the modified MESMA method is an effective approach to quantitative mapping of the lunar minerals and space weathering levels using hyperspectral data. In the future, more minerals can be mapped with higher accuracy if hyperspectral data with a wider spectral range are used based on the method proposed in this study.