上海天文台GPS掩星技术研究现状

简要地介绍了GPS／LEO卫星无线电掩星技术的基本原理、意义以及国际上的最近动态，列举了几个主要的LEO卫星计划；重点分析了GPS掩星探测地球大气技术中需要注意的若干问题，并推出可能采取的改进方案和相应的技术路线。上海天台在GPS／LEO空基气象学中已经成功地提出：通过通约LEO／GPS轨道方法，实现掩星地面观测点控制的新思路。作为一个应用实例，针对东中国海洋上空大气状况对中国东部经济发达的地区特别是上海地区的天气有着重要影响的背景，提出利用通约LEO／GPS卫星可对该地区大气实施监控。为了充分利用GPS／LEO掩星观测资源，上海天台开展了振幅反演的计算方法研究；提出了在振幅反演中存在着无线电信号几何衰减和物理衰减两种不同的衰减机制，并考虑了它们的数学模型以及对反演结果的影响。简要地叙述了掩星质量因子的定义和计算方法。作为国内天地球动力学研究中心，上海天台将与国内外一些合作单位考虑建立GPS／LEO掩星技术观测处理的软件系统。作为空基GPS气象学的推广，还考虑进行山基和飞机载掩星观测实验，论证用这些方法监测局部大气剖面的可行性以及它们在国民经济和科学研究上的作用。     

上海天文台GPS掩星技术研究现状

简要地介绍了GPS/LEO卫星无线电掩星技术的基本原理、意义以及国际上的最近动态 ,列举了几个主要的LEO卫星计划 ;重点分析了GPS掩星探测地球大气技术中需要注意的若干问题 ,并推出可能采取的改进方案和相应的技术路线。上海天文台在GPS/LEO空基气象学中已经成功地提出 :通过通约LEO/GPS轨道方法 ,实现掩星地面观测点控制的新思路。作为一个应用实例 ,针对东中国海洋上空大气状况对中国东部经济发达的地区特别是上海地区的天气有着重要影响的背景 ,提出利用通约LEO/GPS卫星可对该地区大气实施监控。为了充分利用GPS/LEO掩星观测资源 ,上海天文台开展了振幅反演的计算方法研究 ;提出了在振幅反演中存在着无线电信号几何衰减和物理衰减两种不同的衰减机制 ,并考虑了它们的数学模型以及对反演结果的影响。简要地叙述了掩星质量因子的定义和计算方法。作为国内天文地球动力学研究中心 ,上海天文台将与国内外一些合作单位考虑建立GPS/LEO掩星技术观测处理的软件系统。作为空基GPS气象学的推广 ,还考虑进行山基和飞机载掩星观测实验 ,论证用这些方法监测局部大气剖面的可行性以及它们在国民经济和科学研究上的作用     

中性大气掩星标准反演技术

结合JPL和GFZ的中性大气标准反演算法，发展了上海天文台中性大气掩星标准反演算法流程．对CHAMP观测资料进行中性大气反演，并将反演结果与ECMWF(the European Center for Medium Range Weather Forecasting)的预报大气数据进行比较．作为讨论，指出将来的发展方向．     

GPS/LEO掩星技术中超折射效应的修正

在GPS掩星探测地球大气技术中，Abel积分变换要求大气折射指数n是折射半径a的单值函数，当大气折射率的垂直梯度达到小于一个极限值dN／dr≈-0．16N-unit m^-1上述单值性不成立，称其为超折射，此时Abel积分变换不再适用,如果还是在形式上应用经典的Abel变换，在反演结果中就会产生负大气折射率偏差．描述了低对流层中超折射现象的物理特性和数学表示；在广义Abel积分变换的基础上，讨论了超折射层内和超折射层下的大气剖面反演算法；选择了一个简单的采样间隔内等大气折射率垂直梯度假设，对英国高分辨率无线电探空观测资料进行模拟计算，验证了负超折射与大气折射率偏差的关系，并提出的广义Abel变换的合理性.     

在97.8GHz处测量地球大气的不透明度

用研制的全自动辐射计在97．8GHz上测量南京市内大气的不透明度以及这种不透明度随时间的变化．辐射计记录大气在不同天顶高度的辐射温度,由此测出天顶方向的大气不透明度。在56小时的测量时间内,取得了1400多组数据,拟合得到了在观测时段内天顶方向的大气不透明度在0．2至0．7之间,典型值为0．4,表明测试地点的大气条件允许进行3mm波段的射电天文观测．由于在工作频段内大气不透明度的主要来源是对流层水汽的吸收,测量的不透明度可用来直接反映大气内水汽的含量,并实时校准大气吸收．     

大气对地球自转变化的影响-方法和数据的深入分析

大气相对固体地球的运动产生大气相对角动量，它的变化可以激发地球自转多时间尺度的变化．计算大气相对角动量现在采用两种不同的垂直积分高度，一种为从地形表面积分到顶层大气，称之为SP方法；另一种为从1000hPa积分到顶层大气，称之为BP方法，对采用这两种方法所得到的大气相对角动量进行了详细的比较．应用欧洲中距气象预报中心（ECMWF）和美国国家环境预报中心（NCEP）大气再分析数据，重新研究了大气相对角动量变化的时空特征．通过对大气相对角动量季节平均，季节振幅和时空特征的分析，得出ECMWF和NCEP的大气相对角动量变化对地球自转周年极移的影响，在亚洲季风区域和南极洲区域差别最为明显．     

视宁仪与长焦拖影的对比观测

本文给出了大气视宁仪和长焦拖影两种方法、两种仪器在同一圆顶内的同时观测结果及相关分析。结果表明，两种方法所观测到的结果有明显的相关性，相关系数γ＝０．９１。本文还论证了安徽光机所宋正芳等人所研制的大气视宁仪系统是可靠的。     

大气折射对大视场多目标光纤光谱观测的影响

大气折射从两个方面影响大视场多目标光纤光谱的观测：一是大气较差折射，二是由大气折射引起的大气色散．本文定量地分析了这些影响，给出其对一具有５°视场、位于北纬４０．４°的望远镜在不同纬度时影响的大小．     

大气湍流理论的进展及其天文学意义

传统的大气湍流理论始于大约３０年前。湍流理论的发展使得天文学家对大气湍流对天文观测的影响有了很好的认识，尤其是它为天文台选址和天文高分辨率技术提供了基础。近年来新的天文观测对已传统的大气湍流理论提出了挑战，并可能给地面天文观测带来一场革命。在回顾了大气湍流理论的发展历史后对传统大气湍流理论的基本特性及其应用作了系统的综述，并介绍了新的天文观测事实以及为此而提出的新的大气湍流理论模型。     

GPS掩星折射率剖面一维变分同化

近年来,GPS／LEO(全球定位系统／低地球轨道)卫星无线电掩星技术给出了地球大气探测的新途径．从LEO卫星观测到的掩星数据可以反演的地球大气的气压、水汽、温度等剖面;它们对气象和大气科学研究,是具有潜在价值的数据资源．掩星数据资料的同化技术可以有效地改进这些气象参数的剖面,从而改进目前的数值天气预报模式．在当前采用的一维变分同化反演技术中,可以用掩星观测资料的大气折射率或弯曲角剖面进行同化,来反演大气水汽和温度剖面以及海平面压强．作为独立自主开发的GPD／LEO掩星技术系统的一部分,以欧洲中尺度天气预报分析(ECMWF)资料为背景场,CHAMP 掩星观测得到的折射率剖面为观测值,采用Levenberg—Marquardt方法实行GPS掩星资料一维变分同化．在讨论中,用掩星观测点附近相应的探空气球资料来检验CHAMP掩星资料变分同化的结果．     

Optical multiband observations of BL Lacertae during the outburst of 2005

The aim of our observations is to investigate the intranight variability properties and the spectral variability of BL Lacertae. 799 optical multiband observations were intensively made with the Beijing–Arizona–Taiwan–Connecticut (BATC) 60/90 cm Schmidt telescope during the outburst composed of two subsequent flares in 2005. The second flare, whose rising phase lasted at least 44 d, was observed with amplitudes of more than 1.1 mag in three BATC optical bands. In general, on intranight time-scale the amplitude of variability and the variation rate are larger at the shorter wavelength, and the variation rate is comparable in the rising and decaying phases within each band. A possible time-lag between the light curves in bands e and m , around 11.6 min, was obtained. Based on the analysis of the colour index variation with the source brightness, the variability of BL Lacertae can be interpreted as having two components: a 'strongly chromatic' intranight component and a 'mildly chromatic' internight component, which may be the results of both intrinsic physical mechanism and geometric effects.     

GNSS空间信号质量评估系统接收通道性能测试

GNSS空间信号质量评估系统的接收通道是评估导航信号质量的主要误差源之一。通过分析接收通道特性引起的信号幅度误差、相位失真、频率偏移等现象,研究了接收通道增益平坦度、幅度误差、相位误差、频率偏移、矢量误差幅度、群时延等机理,提出了相应的测试方法,测试结果表明接收通道性能能够满足空间信号质量评估的要求。     

Models of young powerful radio sources

Observations of compact symmetric double sources suggest that these objects with physical scales of order tens of parsecs to kiloparsecs are young radio active galactic nuclei. There is, in general, a striking similarity between the structures of these compact objects and the structures of large evolved radio galaxies although the latter are two to three orders of magnitude larger. This has led to the use of self-similar models of the evolution of radio sources as a framework for interpreting observational data. However, the assumptions on which the self-similar models are based become increasingly less valid on the physical scales which are probed by the observations of the smallest sources. In this paper, the dynamics of sources on these small scales is examined and a model developed which extends earlier work in a self-consistent way to small physical scales. The limit of applicability of the models is identified as is the transition from an early evolutionary phase to the self-similar phase of expansion.     

Discovery of the 1.80 h spin period of the white dwarf of the symbiotic system BF Cyg

We report on the discovery of a coherent periodicity in the B light curve of the symbiotic star BF Cyg. The signal was detected in some sections of the light curve of the star recorded in the year 2003 as double-hump periodic variations with an amplitude of ≃7 mmag. In the year 2004, the signal was also present in only a subsection of the light curve. In that year, the system was about twice as bright and the amplitude of the oscillations was about half of what it was in 2003. In 2004, the cycle structure was of a single hump, the phase of which coincided with the phase of one of the humps in the 2003 cycle. No periodic signal was detected in a third, short series of observations performed in the year 2007, when the star was 3 mag brighter than in 2003. We interpret the periodicity as the spin period of the white dwarf component of this interacting binary system. We suggest that the signal in 2003 originated in two hotspots on or near the surface of the white dwarf most likely around the two antipodes of an oblique dipole magnetic field of this star. Magnetic field lines funnelled accreted matter from the wind of the cool component to the pole areas, where the falling material created the hotspots. This process is apparently intermittent in its nature. In 2004, the activity near only one pole was enhanced enough to raise the signal above the threshold of our detection ability.     

Multiwavelength modelling the SED of supersoft X-ray sources I. The method and examples

Radiation of supersoft X-ray sources (SSS) dominates both the supersof X-ray and the far-UV domain. A fraction of their radiation can be reprocessed into the thermal nebular emission, seen in the spectrum from the near-UV to longer wavelengths. In the case of symbiotic X-ray binaries (SyXBs) a strong contribution from their cool giants is indicated in the optical/near-IR. In this paper I introduce a method of multiwavelength modelling the spectral energy distribution (SED) of SSSs from the supersoft X-rays to the near-IR with the aim to determine the physical parameters of their composite spectra. The method is demonstrated on two extragalactic SSSs, the SyXB RX J0059.1-7505 (LIN 358) in the Small Magellanic Cloud (SMC), RX J0439.8-6809 in the Large Magellanic Cloud (LMC) and two Galactic SSSs, the classical nova RX J2030.5+5237 (V1974 Cyg) during its supersoft phase and the classical symbiotic star RX J1601.6+6648 (AG Dra) during its quiescent phase. The multiwavelength approach overcomes the problem of the mutual dependence between the temperature, luminosity and amount of absorption, which appears when only the X-ray data are fitted. Thus, the method provides an unambiguous solution. It was found that selection of the model (a blackbody or an atmospheric model) is not of crucial importance in fitting the global X-ray/IR SED. The multiwavelength modelling of the SED of SSSs is essential in determining their physical parameters.     

X-ray variability in a deep, flux-limited sample of QSOs

We present an analysis of X-ray variability in a flux-limited sample of quasi-stellar objects (QSOs). Selected from our deep ROSAT survey, these QSOs span a wide range in redshift (0.1< z <3.2) and are typically very faint, so we have developed a method to constrain the amplitude of variability in ensembles of low signal-to-noise ratio light curves. We find evidence for trends in this variability amplitude with both redshift and luminosity. The mean variability amplitude declines sharply with luminosity, as seen in local active galactic nuclei (AGN), but with some suggestion of an upturn for the most powerful sources. We find tentative evidence that this is caused by redshift evolution, since the high-redshift QSOs ( z >0.5) do not show the anticorrelation with luminosity seen in local AGN. We speculate on the implications of these results for physical models of AGN and their evolution. Finally, we find evidence for X-ray variability in an object classified as a narrow-emission-line galaxy, suggesting the presence of an AGN.     

10-GHz Tenerife cosmic microwave background observations at 8° resolution and their analysis using a new maximum entropy method

The complete set of data from the Tenerife 10-GHz (8° FWHM) twin-horn, drift scan experiment is described. These data are affected by both long-term atmospheric baseline drifts and short-term noise. A new maximum entropy procedure, utilizing the time invariance and spatial continuity of the astronomical signal, is used to achieve a clean separation of these effects from the astronomical signal, and to deconvolve the effects of the beam-switching. We use a fully positive/negative algorithm to produce two-dimensional maps of the intrinsic sky fluctuations. Known discrete sources and Galactic features are identified in the deconvolved map. The data from the 10-GHz experiment, after baseline subtraction with MEM, are then analysed using conventional techniques, and new constraints on Galactic emission are made.     

Photometric and spectroscopic investigation of the oscillating Algol type binary: EW Boo

We obtained the physical and geometrical parameters of the EW Boo system, which exhibits short period and small amplitude pulsations as well as brightness variations due to orbital motion of components. Towards this end we carried out photometric observations at Ankara University Kreiken Observatory (AUKO) as well as spectroscopic observations at TUBITAK National Observatory (TNO). The light and radial velocity curves obtained from these observations have been simultaneously analyzed with PHOEBE and the absolute parameters of the system along with the geometric parameters of the components have been determined. Using model light curves of EW Boo, light curve regions in which the pulsations are active have been determined and as a result of analyses performed in the frequency region, characteristic parameters of pulsations have been obtained. We find that the results are compatible with current parameters of similar systems in the literature. The evolutionary status of the components is propounded and discussed.     

Flare-induced signals in polarization measurements during the X2.6 flare on 2005 January 15

Flare-induced signals in polarization measurements which were manifested as apparent polarity reversal in magnetograms have been reported since 1981. We are motivated to further quantify the phenomenon by asking two questions: can we distinguish the flare-induced signals from real magnetic changes during flares, and what we can learn about flare energy release from the flare-induced signals? We select the X2.6 flare that occurred on 2005 January 15, for further study. The flare took place in NOAA active re-gion (AR) 10720 at approximately the central meridian, which makes the interpretation of the vector magnetograms less ambiguous. We have identified that flare-induced signals during this flare appeared in six zones. The zones are located within an average distance of 5 Mm from their weight center to the main magnetic neutral line, have an average size of (0.6±0.4)×1017 cm2, duration of 13±4 min, and flux density change of 181±125 G in the area of reversed polarity. The following new facts have been revealed by this study: (1) the flare-induced signal is also seen in the transverse magnetograms but with smaller magnitude, e.g., about 50 G; (2) the flare-induced signal mainly manifests itself as apparent polarity reversal, but the signal starts and ends as a weakening of flux density; (3) The flare-induced signals appear in phase with the peaks of hard X-ray emission as observed by the Ramaty High Energy Solar Spectroscopic lmager (RHESSI), and mostly trace the position of RHESSI hard X-ray footpoint sources. (4) in four zones, it takes place cotemporally with real magnetic changes which persist after the flare. Only for the other two zones does the flux density recover to the pre-flare level immediately after the flare.The physical implications of the flare-induced signal are discussed in view of its relevance to the non-thermal electron precipitation and primary energy release in the flare.     

The diurnal variation of polarization characteristics of the Earth treated as an exoplanet

The diurnal polarization variation of the Earth, treated as an exoplanet with an unresolved disk but resolvable from its host star, is presented in three wavelength bands centered at 490 nm, 670 nm and865 nm respectively according to French satellite-borne PARASOL data. We aim to estimate disk-integrated polarization of the Earth with a phase angle of 55?. It is shown that:(1) the linear polarization signal and its variation are ascribed to the combination of surface feature distribution and atmospheric conditions acting as a variable polarimetric modulator; and(2) the polarimetric wavelength dependence is strong due to the atmospheric wavelength sensitivity. During the period when the PARASOL data were acquired, the cloud coverage ranged from 44.0% to 57.7%, and the polarimetric diurnal variation amplitude was within 1.8%in the 490 nm band, 1.3% in the 670 nm band and 1.5% in the 865 nm band.