关于北斗G2卫星的溯往追终

2022年1月,失效的北斗G2卫星被实践21号卫星从地球静止轨道拖入了坟墓轨道.为了这项捕获任务的安全实施,需要预先确定北斗G2的旋转状态.基于过去10 yr的测光观测数据展示了北斗G2卫星自转的演化过程.根据北斗G2的自转速度和轨道的演化,确认了在过去的10 yr里发生的6次异常事件.据推测, 2012年的小碎片碰撞事件,是随后几年燃油泄漏的导火索. 2017年之后剩余燃油完全释放,再也没有出现转速异常.将2014年太阳能帆板损坏和2016年的解体事件后建立的旋转动力学模型外推1 yr,转轴的标准偏差小于3°,转速标准偏差为0.11°·s^-1,能够有效地满足捕获任务时刻旋转状态的精度要求.     

基于ITRF96和ITRF97的全球板块运动模型

利用国际地球自转服务(IERS)发布的国际地球参考架ITRF96和ITRF97的水平运动,分别建立了完全基于现代空间大地测量实测结果的新的现今全球板块运动模型．研究表明,ITRF96和ITRF97与协议地球参考架(CTRF)的无整体旋转(NNRNo-Net-Rota-tion)的定义不符,都存在一个整体旋转．ITRF96参考架的旋转角速度约为0．11°／m．yr,旋转极位于N-4．13°,E89．3°;ITRF97的旋转角速度约为0．10°／m．yr,旋转极位于N-27．2°,E74．7°．这种不一致性对ITRF96和ITRF97的高精度应用和长期维持、地球自转参数的长期变化都将产生一定的影响。     

太阳活动与太阳自转速率变化关系的分析

研究发现,太阳自转速率的变化与太阳活动之间存在一定的联系,但是不同学者的研究结论存在着矛盾:有的认为两者为正相关,而有的却认为是负相关.究竟两者之间是什么关系,需要做进一步深入的分析.利用EEMD (Ensemble Empirical Mode Decomposition)等方法对太阳自转速率和太阳黑子数据序列进行相关关系以及相位关系的计算和分析,以探讨太阳自转速率变化与太阳活动之间的关系.研究发现:两者的长期趋势项分量呈显著负相关;在11 yr左右周期分量上,观测到的太阳自转速率滞后太阳黑子的变化约2 yr时,呈显著负相关关系,超前3 yr时呈现次显著的正相关;对太阳活动第12–23周各周内部太阳黑子与太阳自转速率的相关分析表明,两者的关系比较复杂,但负相关关系更为显著.这为进一步理解太阳活动变化与太阳自转速率变化之间的成因联系提供了新的依据.     

3颗脉冲星自转性质的研究

正常脉冲星和亳秒脉冲星都表现出计时噪声.脉冲星的自转变化是造成计时噪声的重要原因之一.通过联合新疆天文台南山25m射电望远镜和澳大利亚Parkes64m射电望远镜的脉冲到达时间数据，使用脉冲星计时的方法对PSRJ1539--5626、J1832- -0827和J1847- -0402的自转进行分析研究，它们在到达时间残差上均表现出很强的计时噪声，低频噪声的功率谱分别符合谱指数为-6、-6、 -4.5的幂律,功率谱强度分别为1.77 x 10^-17 yr³、4.43x 10^-18 yr³和2.09x 10^-18yr3.这3颗脉冲星在自转频率1阶导数的变化上都表现出较为明显的振荡,振荡幅度分别为0.61(3)x 10^-15s^-2、0.54(5) x 10^-16s^-2和0.11(2) x 10^-15s^-2 (其中括号内数字代表末位数字的有效误差,下同)，自转频率1阶导数振荡变化的相对大小分别是0.75(5)%、0.035(9)%、 0.076(2)%.利用Parkes 64 m射电望远镜的观测数据,分别获得了这3颗脉冲星积分脉冲轮廓及其半高全宽,发现3颗脉冲星的脉冲轮廓的宽度均有一定 的变化，变化幅度分别为0.0028(6)、0.00059(3)和0.00011(4)个相位.没有探测到自转减慢率的变化与辐射变化之间存在明显的相关性.     

一种基于全天相机组网的区域火流星监测网的实现和应用

流星监测网是小尺寸近地小天体撞击监测、判断陨石落点的主要工具. 提出了一种基于多站布局的全天视频相机组网监测系统, 并在江苏及周边构建了一个区域级原型系统, 实现了火流星监测组网控制、视频数据采集、数据处理及流星体定轨的完整流程. 通过1yr的实测运行表明, 该系统可观测流星极限视星等为-1.0等, 可以实现绝对星等-2.5等流星的完备检测; 根据监测数据得到火流星通量为2.68×10^-7km^-2 ·h^-1;群流星和偶发流星占比分别为46%和54%,偶发流星中类小行星轨道和类彗星轨道比例分别为27.1%和72.9,统计结果与国际主要流星监测网相接近,验证了监测网系统在实际组网使用中的监测能力.     

从我国第六颗卫星轨道倾角分析测定高层大气旋转速度

我国第六颗卫星(1976 87A)发射于1976年8月30日,在轨道上存在816天,于1978年11月25日陨落,初始轨道倾角约69°.16,至1978年5月约69°.10,减少0°.06左右;初始周期108~m.7,至1978年5月为98~m.9,约减少10分钟。从轨道倾角分析,测定了卫星平均近地点高度约220公里附近,平均大气旋转速度∧=1.13±0.17(圈/天)。     

北斗亚欧全视时间比对试验

基于卫星导航双频时间传递型接收机的伪码观测量,利用国际全球卫星导航系统服务组织(International Global Navigation Satellite System (GNSS) Service, IGS)提供的高精度卫星轨道和钟差产品,实现了北斗全视法时间比对.以IGS提供的时间尺度为两个待比对站的公共参考时间,首先使用双频组合法消除电离层对伪距观测的影响,然后将对流层和地球自转效应带来的时延利用理论模型在伪码观测量中进行扣除,分别获得两个比对站时间与公共参考时间之差后,将2者再做差,便得到了北斗全视时间比对结果.以中国科学院国家授时中心(NTSC)、德国物理技术研究院(PTB)和西班牙海军天文台(ROA)所保持的国家标准时间作为比对对象,开展了长基线北斗全视时间比对试验,获得北斗全视时间传递结果,最后利用阿伦方差和时间方差两项关键性能指标以及卫星双向时间比对对其进行性能评估.结果表明:北斗全视时间比对的天稳为10-14量级,可以满足国际时间比对需求.     

武汉SLR站地心坐标的精密测定

目前全球激光测距网日臻完善，国内的激光测距网已初具规模，上海、武汉、长春等的SLR站已投入了正常观测。本文利用两批武汉站的1985年激光测LAGEOS卫星的资料，与同期的全球SLR同多的资料一起，主淫地，精密测定了武汉站的地心坐标，在归算中为了减弱地球自转参数的误差和力学模型的不确定性对测定测站坐标的影响，我们设计了多级复弧法，在不同弧段长度内联合解算地球自转参数、卫星轨道和武汉站的坐标，测得的武汉站地心坐标：高度h=38.87m±0.053m；经度λ=114°.3462470±1°.210^-6；纬度φ=30°. 5418007±1°.1×10^ -6。     

基于无量纲的轨道参数开展卫星性质统计研究

目前已发现了285颗围绕太阳系八大行星公转的卫星, 它们的轨道和物理性质呈现了丰富多样性. 目前为止, 几乎所有的卫星研究工作都基于单个卫星系统或者卫星群, 似乎缺少统一的研究. 提出了一个新的与行星性质无关、只与恒星半径有关的轨道参数n, 定义为以太阳半径为单位的轨道半长轴的自然对数. 不同行星的卫星的n值都存在双极分布, 绝大部分卫星在$n\gtrsim2$区间, 其次在$n\lesssim-1$区间, 位于中间区域的行星则很少. 从卫星物理参数和轨道参数与n的关系中发现, 分属六大行星的卫星有明显的共同特征. 首先, 轨道偏心率和轨道倾角偏大的卫星的n值都在3.5左右, 它们都是巨行星的不规则卫星. 其次, n值在-1和1之间的卫星绝大部分体积大、质量大、反照率高、自转速度慢. 从文献中找到11颗系外卫星候选体, 获得了它们轨道n值和卫星质量, 发现后者也是在-1< n< 1区间最大,其他区间偏小.这些统一的 规律暗示,太阳系内不同行星的卫星形成机制以及太阳系外卫星的形成机制可能一样或类似.     

热海王星系统HD 106315的近$2:1$ 平运动共振捕获与轨道演化

通过结合理论分析和数值模拟方法,可以对热海王星系统HD 106315轨道迁移中的近2:1平运动共振捕获机制以及潮汐作用下的演化过程进行研究.在轨道迁移阶段,初始轨道半长径、初始偏心率以及行星c的偏心率衰减系数K会对系统轨道构型产生影响.数值模拟结果显示当初始轨道半长径分别为ab～0.4 au、ac～0.8 au,偏心率eb和ec均小于0.03时, HD 106315b和HD 106315c在中央恒星的引力作用以及原行星盘粘滞作用下向内迁移, 65000 yr左右两颗行星均可迁移至当前观测位置附近并形成近2:1平运动共振捕获.此外,中央恒星的潮汐效应也可能会对行星系统共振构型产生影响,理论分析表明当行星潮汐耗散系数Q=100时,潮汐效应造成的轨道半长径衰减使系统轨道周期比发生的变化可能是系统脱离共振构型的原因.数值模拟结果显示, HD 106315系统内两颗行星Q103时,来自中央恒星的潮汐效应并不会使行星系统产生明显的偏心率和轨道半长径衰减,不足以使HD 106315行星系统在剩余寿命内脱离2:1平运动共振轨道构型.     

Rotation speed and stellar axis inclination from p modes: how CoRoT would see other suns

In the context of future space-based asteroseismic missions, we have studied the problem of extracting the rotation speed and the rotation-axis inclination of solar-like stars from the expected data. We have focused on slow rotators (at most twice solar rotation speed), first, because they constitute the most difficult case and, secondly, because some of the Convection Rotation and planetary Transits ( CoRoT ) main targets are expected to have slow rotation rates. Our study of the likelihood function has shown a correlation between the estimates of inclination of the rotation axis i and the rotational splitting δν of the star. By using the parameters, i and  δν^⋆=δν sin  i   , we propose and discuss new fitting strategies. Monte Carlo simulations have shown that we can extract a mean splitting and the rotation-axis inclination down to solar rotation rates. However, at the solar rotation rate we are not able to correctly recover the angle i , although we are still able to measure a correct  δν^⋆  with a dispersion less than 40 nHz.     

日月轨道计算对TLE深空编目目标预报精度的影响

双行根数(Two Line Element, TLE)作为一类广泛使用的空间物体编目数据, 其预报精度和误差特性是TLE编目 在空间碎片研究中所要关注的问题之一. TLE编目需要配合SGP4/SDP4 (Simplified General Perturbations 4/Simplified Deep Space 4)模型进行轨道预报, 对深空物体来说, 主要考虑带谐项$J_2$、$J_3$、$J_4$摄动、 第三体日月摄动和特殊轨道共振问题修正等. 其中, SGP4/SDP4模型第三体摄动计算时, 对日月轨道近似采用了长期进动根数和 简单平运动的方式, 在外推10d时存在约2$^\circ$--3${^\circ     

Axial rotation, orbital revolution and solar spin–orbit coupling

The orbital motion of the Sun has been linked with solar variability, but the underlying physics remains unknown. A coupling of the solar axial rotation and the barycentric orbital revolution might account for the relationships found. Some recent published studies addressing the physics of this problem have made use of equations from rotational physics in order to model particle motions. However, our standard equations for rotational velocity do not accurately describe particle motions due to orbital revolution. The Sun's orbital motion is a state of free fall; in consequence, aside from very small tidal motions, the associated particle velocities do not vary as a function of position on or within the body of the Sun. In this note, I describe and illustrate the fundamental difference between particle motions in rotation and revolution, in order to dispel some part of the confusion that has arisen in the past and that which may yet arise in the future. This discussion highlights the principal physical difficulty that must be addressed and overcome by future dynamical spin–orbit coupling hypotheses.     

Upper-atmosphere features revealed by the orbit of 1980-43A

The satellite NOAA-B (1980-43A) was launched in May 1980 into an orbit with perigee height near 260 km and apogee height 1440 km, at an inclination of 92.2°.The lifetime was 11 months. The orbit has been determined at 40 epochs between October 1980 and May 1981 from about 3000 radar and optical observations. The average orbital accuracy, radial and cross-track, was about 100 m, with rather better accuracy in the final 14 days.The variation of orbital inclination has been analysed to determine two good values of atmospheric rotation rate, namely 1.10 ± 0.10 rev day^?1 at 300 km (average local time) and 1.15 ± 0.06 rev day^?1 at 225 km (evening).The effect of atmospheric rotation on the precession of the orbital plane of an actual satellite has never previously been detected; it is clearly apparent for 1980-43A in its last days and conforms to the expected theoretical change.The variation of perigee height has been analysed to determine ten values of atmospheric density scale height, for heights of 280–370 km. These values, accurate to about 3%, exceed by 15% the values indicated by the COSPAR International Reference Atmosphere. Solar activity was higher in the years 1980–1981 than at any time since early 1958 and it appears that the CIRA model underestimates the density and density scale height at high levels of solar activity.     

The old calcium plages differential rotation latitudinal profile: Its gross and fine structure evolution with the solar cycle

This work is a study of the rotational properties of the solar calcium plages, during the time interval 1967–1977; only plages older than 4 days have been the object of this research. We have looked systematically for any significant change occurring during the course of the solar cycle, and any kind of ‘anomaly’ or fine structure in the differential rotation latitudinal profile. We find that such a profile undergoes a cyclic transformation, making it assume the highest steepness at the solar maximum; a sudden flattening then occurs in the first years of declining activity; the last years of the cycle, as the first years of the next one, are characterized by intermediate steepness values. Moreover, we find that, in spite of the general belief that the angular rotation rate is continuously decreasing with increasing heliographic latitude, at least two inversions do exist of such an overall tendency:

1. A narrow, minimal angular-rotation-rate strip lies very close to the equatorward margin of the plage production band; this feature shifts continuously, in a wave-like manner, throughout the solar cycle, from 15/18° to 3/6° latitude.
2. A narrow, maximal angular-rotation-rate strip has been observed lying in the neighbourhood of the poleward margin of the activity band; a process of continuous transformation of the rotation rate profile is always active, in a narrow latitude strip on the equatorward side of such a feature, generating new features of the same kind, which replace the older ones, that disappeared due to the equatorward shift of the plage zone. All that simulates an equatorward shift of the observed ‘anomalies’; we observed them until the minimum activity epoch (1976), at 15/18° latitude. Some relations of these features with both torsionai waves (Howard and LaBonte, 1980) and magnetic activity are briefly discussed.

     

Samos 2 (1961 α 1): Orbit determination and analysis at 31 : 2 resonance

Samos 2, 1961 α 1, launched on 31 January 1961, was the first satellite to enter a sun-synchronous orbit at an inclination of 97.4°. The initial perigee and apogee heights were 474 km and 557 km respectively, the initial period was 94.97 min and the satellite decayed on 21 October 1973 after more than 12 years in orbit.Samos 2 passed through the condition of 31 : 2 resonance in June 1971 and orbital parameters have been determined at 22 epochs from 1674 observations using the RAE orbit refinement program, PROP, between mid-April and Mid-September 1971. The variations of inclination and eccentricity during this time have been analysed and values for six lumped 31st-order harmonic coefficients in the geopotential have been obtained. These have been compared with those derived from the individual coefficients, of order 31 and appropriate degrees, from the most recent Goddard Earth Model, GEM 10C.The decrease in inclination between launch and 1971 has been examined: it is found to be caused mainly by a near-resonant solar gravitational perturbation.     

《大气一号》气球卫星轨道倾角变化分析

引起《大气一号》两颗气球卫星（ＤＱ－１Ａ和ＤＱ－１Ｂ）轨道倾角变化的摄动因素主要是太阳光压摄动、大气旋转和日月引力摄动。太阳光压摄动引起气球卫星轨道倾角增大，平均每天变化约０．００１７，大气旋转引起轨道倾角减小，平均每天变化不到０．０００１，但随着高度下降，变化量亦增大，陨落前达０．００２。本文根据卫星轨道摄动理论，给出气球卫星轨道倾角变化的一种定量分析方法，得到的分析结果为：（１）由太阳光压摄动     

Signature of differential rotation in solar disk‐integrated chromospheric line emission

UARS SOLSTICE data have been subjected to Fourier and wavelet analyses in order to search for the signature of the solar rotation law in the disk‐integrated irradiance of UV lines. Lyman‐α, Mg II, and Ca II data show a different behaviour. In the SOLSTICE data there are significant temporal variations of the rotation rate of the UV tracers over 5—6 years. Often several distinct rotation periods appear almost simultaneously. Beside the basic period around 27 days there are signals at 32—35 days corresponding to the rotation rate at very high latitudes. For more than 5 years during another period of the solar cycle the rotational behaviour is quite different; there is an indication of differential rotation of active regions in these Ca II ground‐based data. The data contain a wealth of information about the solar differential rotation, but it proves difficult to disentangle the effects of the different emitting sources.     

On the solar rotation and activity

The interaction between differential rotation and magnetic fields in the solar convection zone was recently modelled by Brun (2004). One consequence of that model is that the Maxwell stresses can oppose the Reynolds stresses, and thus contribute to the transport of the angular momentum towards the solar poles, leading to a reduced differential rotation. So, when magnetic fields are weaker, a more pronounced differential rotation can be expected, yielding a higher rotation velocity at low latitudes taken on the average. This hypothesis is consistent with the behaviour of the solar rotation during the Maunder minimum. In this work we search for similar signatures of the relationship between the solar activity and rotation determined tracing sunspot groups and coronal bright points. We use the extended Greenwich data set (1878–1981) and a series of full-disc solar images taken at 28.4 nm with the EIT instrument on the SOHO spacecraft (1998–2000). We investigate the dependence of the solar rotation on the solar activity (described by the relative sunspot number) and the interplanetary magnetic field (calculated from the interdiurnal variability index). Possible rotational signatures of two weak solar activity cycles at the beginning of the 20th century (Gleissberg minimum) are discussed. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)