日本SELENE月球探测计划和卫星间多普勒跟踪的数学模型

日本月球探测计划（SELENE）定于2004年夏季利用HIIa火箭发射一组共3颗绕月人造卫星。他们是主卫星、跟踪中断卫星和空间VLBI电波源。其主要科学目标之一是利用对绕月卫星的多普勒跟踪数据精确测定月球重力场,研究月球的起源与演化。SELENE计划中实现这个科学目标的关键技术是引入中继卫星,目的在于当处于低轨道的主卫星飞行到月球背面地面观测站无法观测时,采用卫星间跟踪方法（SST）,建立地面站与主卫星之间的联系,以得到月球背面重力场的直接测量数据。介绍了几种典型的四程卫星间多普勒跟踪模式和相应的数学模型,并针对SELENE计划中采用的特殊四程多普勒跟踪模式建立了卫星相对观测站速度与跟踪信号多普勒频移之间的转换关系。提出了利用GEODYNⅡ定轨分析软件处理SELENE多普勒跟踪数据的流程。     

Orbit Determination Using Satellite—to—Satellite Tracking Data

1 INTRODUCTIONThe tracking arc-length should be increased in order to approve the accuracy in orbitdetermination of LEO (Low Earth Orbit) satellites. The local ground-based tracking networkdoes not provide sufficient orbit coverage for the user satellites. The most promising methodis to use high orbiting satellites, such as GPS and TDRS, as trackers to observe the usersatellites. For examPle, tWO geosynclironous satellites could cover more than 85% of the orbitof any given user sate…     

The lunar gravity mission MAGIA: preliminary design and performances

The importance of an accurate model of the Moon gravity field has been assessed for future navigation missions orbiting and/or landing on the Moon, in order to use our natural satellite as an intermediate base for next solar system observations and exploration as well as for lunar resources mapping and exploitation. One of the main scientific goals of MAGIA mission, whose Phase A study has been recently funded by the Italian Space Agency (ASI), is the mapping of lunar gravitational anomalies, and in particular those on the hidden side of the Moon, with an accuracy of 1 mGal RMS at lunar surface in the global solution of the gravitational field up to degree and order 80. MAGIA gravimetric experiment is performed into two phases: the first one, along which the main satellite shall perform remote sensing of the Moon surface, foresees the use of Precise Orbit Determination (POD) data available from ground tracking of the main satellite for the determination of the long wavelength components of gravitational field. Improvement in the accuracy of POD results are expected by the use of ISA, the Italian accelerometer on board the main satellite. Additional gravitational data from recent missions, like Kaguya/Selene, could be used in order to enhance the accuracy of such results. In the second phase the medium/short wavelength components of gravitational field shall be obtained through a low-to-low (GRACE-like) Satellite-to-Satellite Tracking (SST) experiment. POD data shall be acquired during the whole mission duration, while the SST data shall be available after the remote sensing phase, when the sub-satellite shall be released from the main one and both satellites shall be left in a free-fall dynamics in the gravity field of the Moon. SST range-rate data between the two satellites shall be measured through an inter-satellite link with accuracy compliant with current state of art space qualified technology. SST processing and gravitational anomalies retrieval shall benefit from a second ISA accelerometer on the sub-satellite in order to decouple lunar gravitational signal from other accelerations. Experiment performance analysis shows that the stated scientific requirements can be achieved with a low mass and low cost sub-satellite, with a SST gravimetric mission of just few months.     

卫星跟踪卫星模式中轨道参数需求分析

首次基于半解析法利用GRACE(Gravity Recovery and Climate Experiment)双星K波段星间速度误差、GPS接收机轨道误差和加速度计非保守力误差影响累计大地水准面精度的联合模型开展了卫星跟踪卫星模式中轨道参数的需求分析.建议我国将来首颗重力卫星的平均轨道高度设计为400 km和平均星间距离设计为220 km较优.此研究不仅为我国将来卫星重力测量计划中轨道参数的优化选取以及全球重力场精度的有效和快速估计提供了理论基础和计算保证,同时对将来国际GRACE Follow-On地球重力测量计划和GRAIL(Gravity Recovery and Interior Laboratory)月球重力探测计划的发展方向具有一定的指导意义.     

Quick Detection System for Planck satellite

Q uick D etection S ystem or qds is a software package that has been developed for detecting point sources in the Planck satellite data as soon as the data become available, a few days after transmission to the Earth. Point sources are detected by filtering the data with a filter defined by the Mexican hat wavelet. An alert is generated on those detections that are found to be interesting, such as prominent flaring, according to the criteria specified to the software. The goal is to detect spectral or flux variability in active galactic nuclei so that instant multifrequency follow-up observations with other instruments could be arranged to study the interesting behaviour.     

Calculating the Light Period of BL Lac Object S5 0716+714 by Wavelet Analysis Method

A kind of wavelet analysis method for identifying the light period of the BL Lac object S5 0716+714 is introduced. The rather complete observed data in the four optical wavebands B, V, R, I are collected, and the long-term light curves based on 10-day averaging are obtained. On these light curves, the periodicity analysis is performed by using a wavelet analysis method. The result demonstrates that the wavelet analysis method is preferable for searching and identifying the light periods of BL Lac objects. From the contour map of the real part of the wavelet transform coefficient, the periodical light variations of the BL Lac object S5 0716+714 can be precisely identified. By analyzing the wavelet variance curves at the 4 wavebands, it is found that the BL Lac Object S5 0716+174 has a stable light period of 1160 days. This result is consistent with the 3.3-year period given by Raiteri et al. It is predicted that the next outburst in this object will happen around the August of 2011.     

Solar Energetic Particle Event of 2005 January 20： Release Times and Possible Sources

Based on cosmic ray data obtained by neutron monitors at the Earth's surface, and data on near-relativistic electrons measured by the WIND satellite, as well as on solar X-ray and radio burst data, the solar energetic particle (SEP) event of 2005 January 20 is studied. The results show that this event is a mixed event where the flare is dominant in the acceleration of the SEPs, the interplanetary shock accelerates mainly solar protons with energies below 130 MeV, while the relativistic protons are only accelerated by the solar flare. The interplanetary shock had an obvious acceleration effect on relativistic electrons with energies greater than 2 MeV. It was found that the solar release time for the relativistic protons was about 06:41 UT, while that for the near-relativistic electrons was about 06:39 UT. The latter turned out to be about 2 min later than the onset time of the interplanetary type III burst.     

四种综合脉冲星时算法比较

由单颗脉冲星定义的脉冲星时受多种噪声源的影响,其短期和长期稳定度都不够好.为了削弱这些噪声源对单脉冲星时的影响,可以采取合适的算法对多个单脉冲星时进行综合得到综合脉冲星时,从而提高综合脉冲星时的长期稳定度.文中介绍4种综合脉冲星时算法:经典加权算法、小波分析算法、维纳滤波算法和小波域中的维纳滤波算法,将这4种算法分别应用于Arecibo天文台对两颗毫秒脉冲星PSR B1855+09和PSRB1937+21观测得到的计时残差并作出比较.     

Steps Towards a Fully Automated Classification and Redshiftmeasurement Pipeline for LAMOST Spectra. I. Continuum level and wavelength estimation for galaxies

1 INTRoDUCTIONThere are aiready several wide-field multi-fiber spectrograPhic surveys in various stages ofdevelopment, most notably the stellar popu1ation, galaxy and QSO suxveys using the 2dF faciLity as described by Lewis et al. (1998) and the Sloan Digital Sky Survey (SDSS) aJs summarizedby York et al. (2000), both of which are already underway, and those surveys planning to usethe 6dF facility as described by Wason et al. (2000) and the LAMOST facility as outllnedby Wang et a…     

The ASCA spectrum of the z = 4.72 blazar GB 1428+4217

The X-ray-luminous quasar GB 1428+4217 at redshift 4.72 has been observed with ASCA . The observed 0.5–10 keV flux is 3.2 Å– 10⁻¹² erg cm⁻² s⁻¹. We report here on the intrinsic 4 − 57 keV X-ray spectrum, which is very flat (photon index 1.29). We find no evidence for flux variability within the ASCA data set or between it and ROSAT data. We show that the overall spectral energy distribution of GB 1428+4217 is similar to that of lower redshift MeV blazars, and present models that fit the available data. The Doppler beaming factor is likely to be at least 8. We speculate on the number density of such high-redshift blazars, which must contain rapidly formed massive black holes.     

Time Scales and Tidal Effects in Minor Mergers

We use controlled N-body simulation to investigate the dynamical processes (dynamical friction, tidal truncation, etc.) involved in the merging of small satellites into bigger halos. We confirm the validity of some analytic formulae proposed earlier based on simple arguments. For rigid satellites represented by softened point masses, the merging time scale depends on both the orbital shape and concentration of the satellite. The dependence on orbital ellipticity is roughly a power law, as suggested by Lacey & Cole, and the dependence on satellite concentration is similar to that proposed by White. When merging satellites are represented by non-rigid objects, Tidal effects must be considered. We found that material beyond the tidal radius are stripped off. The decrease in the satellite mass might mean an increase in the merging time scale, but in fact, the merging time is decreased, because the stripped-off material carries away a proportionately larger amount of of orbital energy and angular momentum.     

TESS空间卫星单扇区变星的周期性研究

TESS (Transiting Exoplanet Survey Satellite)空间卫星提供的短曝光、高精度光度测量为寻找并区分变星与搜寻行星提供了良好的数据.利用变星源的光变曲线,使用周期频谱分析与光变折叠等一系列方法分析了TESS空间卫星21扇区19995颗拥有高质量光变数据的目标源,并对这些源进行了分类,共获得4624颗变星,其中食双星322颗、脉动变星470颗、行星凌星37颗.所得变星结果与VSX (The International Variable Star Index)变星表进行了交叉比较,共交叉匹配了625颗变星源,这些交叉源中共有131颗为食双星系统、31颗为脉动变星,并通过周期频谱分析获取了双星绕转以及脉动周期.另外在59颗变星中发现明显耀发现象,交叉源中有8颗变星为行星凌星并同样通过周期频谱分析获取了行星绕转周期,从而验证了TESS空间卫星数据对变星分析的可行性.通过利用TESS空间卫星21扇区获得的变星周期结果与VSX变星表中提供的变星周期对比,发现与VSX变星表中绝大部分变星的周期一致,有一部分结果与VSX变星表中的结果差别较大,对这些变星周期结果做了进一步修正,并给出了变星表未列出的变星周期结果.     

A high-significance detection of non-Gaussianity in the Wilkinson Microwave Anisotropy Probe 1-yr data using directional spherical wavelets

A directional spherical wavelet analysis is performed to examine the Gaussianity of the Wilkinson Microwave Anisotropy Probe ( WMAP ) 1-yr data. Such an analysis is facilitated by the introduction of a fast directional continuous spherical wavelet transform. The directional nature of the analysis allows us to probe orientated structure in the data. Significant deviations from Gaussianity are detected in the skewness and kurtosis of spherical elliptical Mexican hat and real Morlet wavelet coefficients for both the WMAP and Tegmark, de Oliveira-Costa & Hamilton foreground-removed maps. The previous non-Gaussianity detection made by Vielva et al. using the spherical symmetric Mexican hat wavelet is confirmed, although their detection at the 99.9 per cent significance level is only made at the 95.3 per cent significance level using our most conservative statistical test. Furthermore, deviations from Gaussianity in the skewness of spherical real Morlet wavelet coefficients on a wavelet scale of 550 arcmin (corresponding to an effective global size on the sky of ∼26° and an internal size of ∼3°) at an azimuthal orientation of 72°, are made at the 98.3 per cent significance level, using the same conservative method. The wavelet analysis inherently allows us to localize on the sky those regions that introduce skewness and those that introduce kurtosis. Preliminary noise analysis indicates that these detected deviation regions are not atypical and have average noise dispersion. Further analysis is required to ascertain whether these detected regions correspond to secondary or instrumental effects, or whether in fact the non-Gaussianity detected is due to intrinsic primordial fluctuations in the cosmic microwave background.     

No Higher Criticism of the Bianchi-corrected Wilkinson Microwave Anisotropy Probe data

Higher Criticism (HC) has been proposed by Donoho & Jin as an effective statistic to detect deviations from Gaussianity. Motivated by the success of the Bianchi VII_h model in addressing many of the anomalies observed in the Wilkinson Microwave Anisotropy Probe ( WMAP ) data (Jaffe et al.), we present calculations in real and in wavelet space of the HC statistic of the Bianchi-corrected WMAP first-year data. At the wavelet scale of 5°, the HC of the WMAP map drops from a value above the 99 per cent confidence level (c.l.) to a value below the 68 per cent CL when corrected by the Bianchi template. An important property of the HC statistic is its ability to locate the pixels that account for the deviation from Gaussianity. The analysis of the uncorrected WMAP data pointed to a cold spot in the Southern hemisphere, centred at   l ∼ 209°, b ∼−57°  . The HC of the Bianchi-corrected map indicates that this spot remains prominent, albeit at a level completely consistent with Gaussian statistics. Consequently, it is debatable how much emphasis should be placed on this residual feature, but we consider the effect of modestly increasing the scaling of the template. A factor of only 1.2 renders the spot indistinguishable from the background level, with no noticeable impact on the results published in Jaffe et al. for the low- l anomalies, large-scale power asymmetry or wavelet kurtosis. A trivial interpretation would be that the Bianchi template may require a small enhancement of power on scales corresponding to the wavelet scale of 5°.     

Ha Line Polarization in the 2B/X4.8 Flare of 2002 July 23

On 2002 July 23, a 2B/X4.8 flare was observed in the Ha line spec-tropolarimetrically by the Large Solar Vacuum Telescope of Baikal Astrophysical Observatory. Linear polarization of 3%-10% was detected in the Ha line, particularly where the line showed central reversal. The linear polarization is mainly radial on the solar disk and appears at the impulsive phase of the hard X-ray and 7-ray bursts. It is limited to some relatively small regions of the flare. The polarization in a limited small region (~ 4" - 5") changed its direction within a short period of time (~ 10s).     

A method for using a fully coupled climate system model to generate detailed surface boundary conditions for paleoclimate modeling investigations: an early Paleogene example

Although fully coupled models of the earth system are now common, simpler model architectures maintain significant utility, and scientific investigations aimed at understanding paleoclimates are frequently conducted with fixed sea surface temperature (SST) or slab ocean modeling experiments. One of the challenges facing the paleoclimate community is that the proxy data used to generate SST boundary conditions exist at a finer resolution and with very limited spatial coverage when compared to a climate model. In addition, SST proxy estimates often represent a single season or annual average conditions. This mismatch in coverage and resolution frequently results in paleoclimate modelers using SST distributions that have very limited spatial and temporal variability. In many regions, a spatially and temporally detailed SST distribution may be necessary for the accurate reproduction of paleoclimatic conditions. Here we borrow from the concept of flux correction and, using available proxy estimates of SST as our guide, force a fully coupled earth system model to produce a spatially and temporally detailed SST distribution for the paleoclimate of the early Paleogene (45–65 Ma). The SST values we produce represent a conservative estimate of early Paleogene high latitude SSTs and match tropical temperatures for this time period well. In addition to matching proxy estimates, our model-derived SST distribution has spatial and temporal variability that meshes well with global climate model resolution. This detailed SST distribution is now available to us as we investigate the causes and sensitivities of early Paleogene climate in fixed SST and slab ocean modeling experiments. The method we used to generate this spatially and temporally detailed SST distribution may prove useful for those investigating other time periods in the past, or the future, for which detailed model boundary conditions are unavailable.     

Time-spectra of chromospheric activity of old solar-type stars: detection of rotational signals from double wavelet analysis

We introduce a novel technique, called the double wavelet analysis (DWA), for the determination of stellar rotation periods from time serial data. This first paper aims narrowly at the discussion, introduction and application of the DWA technique to records of surface magnetism in solar-type (relatively old) lower main sequence stars that are obtained by the Mount Wilson Observatory (MWO) HK Project. The technique takes a series of careful steps that seek to optimize wavelet parameters and normalization schemes, ultimately allowing fine-tuned, arguably more accurate, estimates of rotation-modulated signals (with, e.g., periods of days to months) in records that contain longer periodicities such as stellar magnetic activity cycles (with, e.g., period of years). The apparent rotation periods estimated from the DWA technique are generally consistent with results from both “first-pass” (i.e., ordinary) global wavelet spectrum and earlier classical periodogram analyses. But there are surprises as well. For example, the rotation period of the ancient subdwarf Goombridge 1830 (HD 103095), previously identified as ≈31 days, suggests under the DWA technique a significantly slower period of 60 days. DWA spectra also generally reveal a shift in the cycle period toward high frequencies (hence shorter periods) compared to the first-pass wavelet spectrum. For solar-type stars analyzed here, the character of the DWA spectrum and slope of the first-pass global wavelet spectrum produce a classification scheme that allows a star's record to be placed into one of three categories.     

ROSAT PSPC detection of soft X-ray absorption in GB 1428+4217: the most distant matter yet probed with X-ray spectroscopy

We report on a ROSAT PSPC observation of the highly luminous z =4.72 radio-loud quasar GB 1428+4217 obtained between 1998 December 11 and 17, the final days of the ROSAT satellite. The low-energy sensitivity of the PSPC detector was employed to constrain the intrinsic X-ray absorption of the currently most distant X-ray detected object. Here we present the detection of significant soft X-ray absorption towards GB 1428+4217, making the absorbing material the most distant matter yet probed with X-ray spectroscopy. X-ray variability by 25±8 per cent is detected on a time-scale of 6500 s in the rest frame. The X-ray variation requires an unusually high radiative efficiency of at least 4.2, further supporting the blazar nature of the source.     

Research on Periodicity of Single Sector Variable Star of TESS Space Satellite

The short exposure and high precision photometry of TESS (Transiting Exoplanet Survey Satellite) space satellite provide good data for distinguishing variable stars and searching planets. In this paper, a series of methods, such as periodic spectrum and phase folding, are used to analyze 19995 target sources with high-quality variable data in TESS sector 21, and these sources are classified. A total of 4624 variable stars are obtained, including 322 binary stars, 470 pulsating variable stars, and 37 planetary transients. A total of 625 variable sources were cross matched with VSX (The International Variable Star Index). Among them, 131 were eclipsing binary systems and 31 were pulsating variables, whose periods are obtained by periodic spectrum. The other 59 variable stars have flare phenomena, and 8 stars have transiting planets, whose rotation periods are also obtained through the periodic spectrum. The feasibility of the analysis of variable stars by TESS space satellite data is verified. By comparing the period results obtained by TESS space satellite sector 21 with the variable star period provided by VSX variable star table, we found that the periods of most variable stars are consistent with those in this paper, and a few of them need to be corrected. Finally, the verified and corrected variable star tables are provided in this paper.     

使用小波技术自动搜寻天体谱线

使用小波技术对包括恒星、近邻星系和ＡＧＮ等不同的天体光谱进行了自动处理．用小波滤波的方法将光谱中的连续谱与诸线分离;然后使用小波域隐含马尔可夫模型（ＨＭＭ）,对已去掉连续谱的光谱进行降噪,同时得到了噪声在每个光谱中的分布;在训练ＨＭＭ的过程中,使用改进的Ｔｙｉｎｇ方法增加训练数据以保证训练的可靠性;利用噪声分布确定出谱线信号的局部阈值,在已经降噪的光谱中找到吸收线和发射线;用高斯函数拟会出谱线的形状,标出线心的波长值,作为自动证认的基础．