基于ITRF2000的全球板块运动模型

利用最新国际地球自转服务 (IERS)发表的国际地球参考架ITRF2 0 0 0速度场 ,建立一个完全基于现代空间大地测量实测结果的现今全球板块运动模型ITRF2 0 0 0VEL。板块的总角动量 |L| =0 .12 7sr M·a-1,即ITRF2 0 0 0不满足无整体旋转的要求 ,与协议地球参考架CTRF定义不符 ,存在整体性旋转     

分析昆明GPS站相对欧亚板块的运动

利用国际联测和综合处理 ,得到昆明GPS站地壳形变速率每年以( -4 .3± 0 .5)mm下沉趋势 ,水平分别以 3 7.5mm ,方位 1 47°± 1 .°5运动。利用最新国际地球自转服务 (IERS)发布的国际地球参考架ITRF2 0 0 0速度场 ,建立最新反映全球板块模型 ,基于最新的ITRF2 0 0 0地球参考架和欧亚板块的欧拉矢量 ,估计得到昆明站的地壳垂直形变速率为每年 -0 .0 1mm呈下降趋势 ,水平形变速率为每年 ( 8.1± 0 .2 0 )mm ,方位 1 3 7°± 1 .°5;并进一步分别基于几百万年地质地磁模型NNR -NUVEL1A和ITRF96、ITRF97模型的欧亚板块的欧拉矢量 ,得到较一致的结果 ,说明全球板块运动稳定性与模型的可靠性。本文基于不同板块模型分析了昆明GPS站相对于欧亚板块运动     

基于ITRF97的上海和鹿岛VLBI站地壳垂直形变的VLBI测定

利用美国宇航局(NASA)在1999年的VLBI全球解(glb1123)中关于上海VLBI站与国际上13个VLBI站的基线长变化率高精度测定结果,以及日本鹿岛Kashima和Kashima 34两VLBI站分别与国际VLBI站间的相应结果,基于最近的国际地球参考架ITRF97,估计得到了上海VLBI站和日本鹿岛两VLBI站相对于板块稳定区的地壳垂直形变率分别为-1.91±0.56 mm/yr、-3.72±0.74 mm/yr和-8.81±0.84mm/yr,表明上海VLBI站存在每年近2毫米的下降,讨论并证实日本鹿岛两VLBI站的地壳垂直形变速率是不一致的;基于不同板块运动模型的欧拉旋转矢量,估计得到了上述3站的水平形变速率;还估计给出了太平洋板块中部的形变率结果.     

ITRF2005简介和评析

该文对最新的国际地球参考框架 ITRF2005进行了简单介绍,指出了IFRF2000与ITRF2005之间在解的生成,基准的定义和实现等方面的差别,特别评析了ITRF2005所作改进的理由和合理性.同时也指出了ITRF2005与ITRF2000一样,没有满足协议地球参考系(CTRS)的定向随时间的演化需遵循相对于地壳无整体旋转(NNR)的约束条件,并就此提出了解决办法.     

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

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

MORVEL构造板块的转动张量（英文）

NNR-MORVEL56板块运动模型描述了全球56个构造板块在无整体旋转参考架下的角速度运动参数。这些板块可以近似描述为单位球上的无重叠球面多边形区域。用ITRF速度场计算这56个板块相对于无整体旋转参考架下的绝对运动时,板块的几何参数起着至关重要的作用。详细给出了计算板块几何参数的方法并且编写了FORTRAN90程序以供参考,使得计算单位球上板块的面积和转动惯性张量得以实现。文中的计算方法和程序主要采用球面三角算法和自适应辛普森双积分算法,并对全球56个板块的几何参数进行了计算,得到了较为可靠的计算结果。     

并置VLBI站的相对形变及其与全球解结果的比较

基于ITRF2000地球参考架的三维VLBI站速度矢量和NNR-NUVEL-1A地球板块运动模型，采用观测的VLBI基线长变化率作为约束，估计和讨论了全球6个并置VLBI站的局部或区域性地壳的相对形变，并与国际地球参考架ITRFs解和VLBI全球解glb2003以及VTRF2003的结果进行了比较．结果表明日本Kashima与Kashim34VLBI站之间，美国GREENBANK的NRAO20与NRAO85-3站之间可能存在每年约近4毫米的相对垂直形变；WESTFORD与HAYSTACK二站间每年近2毫米的垂直形变率差基本可得到肯定，KOKEE和KAUAI二站的垂直形变率约有每年1至3毫米的差．利用如下方法能较可靠地检测出并置VLBI站间的局部或区域性的地壳垂直形变，改进ITRF解的可能不足．     

球状星团M3的绝对自行测定及其空间运动

使用上海天文台佘山40厘米折射望远镜拍摄的两期共14张球状星团M3底片的PDS扫描结果,以ACT星表中的24颗恒星的位置和绝对自行为参考架,用中心重叠法进行天体测量归算,得到了这个星团中心附近1°5．5x1°．5天区内534颗恒星的位置和绝对自行,并利用这些自行对成员概率进行估算,得到这个星团的赤经方向绝对自行为-0．3±0．3mas/yr,赤纬方向绝对自行为-3．1±0．3mas/yr．使用该自行值,结合已知的M3的距离和视向速度,计算了其在给定的银河系引力势中的运动轨道．     

利用地球定向参数确定依巴谷输入星表与射电天球参考架之间的旋转

通过比较地球定向参数（ＥＯＰ）序列和对应的地固参考架，本文得到了依巴谷输入星表（ＨＩＣ）和国际地球自转服务（ＩＥＲＳ）基于甚长基线干涉测量（ＶＬＢＩ）建立的河外天球参考架（ＩＣＲＦ）之间的旋转参数，在Ｊ２０００．０历元，前者至后者的旋转参数为［８．４±１．１ｍａｓ，４６．７±１．１ｍａｓ，４５．５±１．５ｍａｓ］，其时间变率为［－０．１４±０．０５ｍａｓ／ａ，３．２２±０．０５ｍａｓ／ａ，５．７０±０．０７ｍａｓ／ａ］。依巴谷输入星表是完全基于ＦＫ５系统的，本文得到的旋转参数与ＩＣＲＦ和ＦＫ５之间的系统差异情况相符     

H37地极坐标序列的低频变化及其与POLE95相应结果的比较

利用依巴谷暂定星表（H37）为统一参考架解算出的地极坐标序列的POLE95序列，分析了地极价值的长期变化、Markowitz变化和其它低频变化。从两个序列得到地极移动长期漂移速度的平均值为3.86mas/yr，方向为西经82°.7；对Markowitz波动的分析表明，其周期约为25－30年，周期和振幅的变化程度分别为10％和40％。     

Determination of the regional deformation rates of Shanghai and Kashima VLBI stations based on ITRF97

The vertical deformation rates (VDRs) and horizontal deformation rates (HDRs) of Shanghai VLBI station in China and Kashima and Kashima34 VLBI stations in Japan were re-analysed using the baseline length change rates from Shanghai to 13 global VLBI stations, and from Kashima to 27 stations and from Kashima34 to 12 stations, based on the NASA VLBI global solution glb1123 (Ma, 1999). The velocity vectors of the global VLBI stations were referred to the ITRF97 reference frame, and the Eulerian vectors of different models of plate motion were used for comparative solutions. The VDR of Shanghai station is estimated to be −1.91±0.56 mm/yr, and those of Kashima and Kashima34 stations, −3.72±0.74mm/yr and −8.81±0.84mm/yr, respectively. The difference between the last two was verified by further analysis. Similar estimates were also made for the Kokee, Kauai and MK_VLBA VLBI stations in mid-Pacific.     

A thermodynamic model of high temperature lava vaporization on Io

We modified the MAGMA chemical equilibrium code developed by Fegley and Cameron (1987, Earth Planet. Sci. Lett. 82, 207-222) and used it to model vaporization of high temperature silicate lavas on Io. The MAGMA code computes chemical equilibria in a melt, between melt and its equilibrium vapor, and in the gas phase. The good agreement of MAGMA code results with experimental data and with other computer codes is demonstrated. The temperature-dependent pressure and composition of vapor in equilibrium with lava is calculated from 1700 to 2400 K for 109 different silicate lavas in the ONaKFeSiMgCaAlTi system. Results for five lavas (tholeiitic basalt, alkali basalt, Barberton komatiite, dunite, and a molten type B1 Ca, Al-rich inclusion) are discussed in detail. The effects of continuous fractional vaporization on chemistry of these lavas and their equilibrium vapor are presented. The predicted abundances (relative to Na) of K, Fe, Si, Al, Ca, and Ti in the vapor equilibrated with lavas at 1900 K are lower than published upper limits for Io's atmosphere (which do not include Mg). We predict evaporative loss of alkalis, Fe, and Si during volcanic eruptions. Sodium is more volatile than K, and the Na/K ratio in the gas is decreased by fractional vaporization. This process can match Io's atmospheric Na/K ratio of 10±3 reported by Brown (2001, Icarus 151, 190-195). Silicon monoxide is an abundant species in the vapor above lavas. Spectroscopic searches are recommended for SiO at IR and mm wavelengths. Reactions of metallic vapors with S- and Cl-bearing volcanic gases may form other unusual gases including MgCl₂, MgS, MgCl, FeCl₂, FeS, FeCl, and SiS.     

Born–Infeld type phantom model in the ω–ω′ plane

In this paper, we investigate the dynamics of Born–Infeld (B–I) phantom model in the ω–ω′ plane, which is defined by the equation of state parameter for the dark energy and its derivative with respect to N (the logarithm of the scale factor a). We find the scalar field equation of motion in ω–ω′ plane, and show mathematically the property of attractor solutions which correspond to ω _φ ∼−1, Ω _φ =1, which avoid the “Big rip” problem and meets the current observations well.      

Shear and vorticity in a combined Einstein-Cartan-Brans-Dicke inflationary lambda-Universe

A combined BCDE (Brans-Dicke and Einstein-Cartan) theory with lambda-term is developed through Raychaudhuri’s equation, for inflationary scenario. It involves a variable cosmological constant, which decreases with time, jointly with energy density, cosmic pressure, shear, vorticity, and Hubble’s parameter, while the scale factor, total spin and scalar field increase exponentially. The post-inflationary fluid resembles a perfect one, though total spin grows, but the angular speed does not (Astrophys. Space Sci. 312: 275, 2007d).      

A Comparison of Several Image Locating Algorithms in Astronomical Instruments

In astronomical observations at optical wavelengths, a fast image tracking system can be adopted to reduce the effects of the atmospheric seeing and telescopic tracking error, and therefore improve the observing efficiency. Aiming at the need of astronomical observations, totally 5 kinds of algorithms in two categories were selected to make a comparative study on their accuracies and stabilities under different noise conditions by both numerical experiment and laboratory test. The results indicate that the normalized cross-correlation method and barycenter method have not only a higher accuracy but also a better reliability against interferences, they will be applied to the high-resolution spectrograph of the Xinglong 2.16 m telescope and the scienti?c instruments of the SONG (Stellar Observations Network Group) project, respectively.     

A Statistical Study of the Mass Distribution of Neutron Stars

By reviewing the methods of mass measurements of neutron stars in four different kinds of systems, i.e., the high-mass X-ray binaries (HMXBs), low-mass X-ray binaries (LMXBs), double neutron star systems (DNSs) and neutron star-white dwarf (NS-WD) binary systems, we have collected the orbital parameters of 40 systems. By using the boot-strap method and the Monte-Carlo method, we have rebuilt the likelihood probability curves of the measured masses of 46 neutron stars. The statistical analysis of the simulation results shows that the masses of neutron stars in the X-ray neutron star systems and those in the radio pulsar systems exhibit different distributions. Besides, the Bayes statistics of these four different kind systems yields the most-probable probability density distributions of these four kind systems to be (1.340 ± 0.230)M8, (1, 505 ± 0.125)M8,(1.335 ± 0.055)M8 and (1.495 ± 0.225)M8, respectively. It is noteworthy that the masses of neutron stars in the HMXB and DNS systems are smaller than those in the other two kind systems by approximately 0.16M8. This result is consistent with the theoretical model of the pulsar to be accelerated to the millisecond order of magnitude via accretion of approximately 0.2M8. If the HMXBs and LMXBs are respectively taken to be the precursors of the BNS and NS-WD systems, then the influence of the accretion effect on the masses of neutron stars in the HMXB systems should be exceedingly small. Their mass distributions should be very close to the initial one during the formation of neutron stars. As for the LMXB and NS-WD systems, they should have already under- gone the process of suffcient accretion, hence there arises rather large deviation from the initial mass distribution.