基于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站的水平形变速率;还估计给出了太平洋板块中部的形变率结果.     

以基线长度变化率为约束估算VLBI站的垂直形变率

基于NNR-NUVEL-1A地球板块运动模型和ITRF2000地球参考架的三维VLBI站速度矢量,采用实测的VLBI基线长度变化率作为约束,重新估计了部分国际VLBI站的局部或区域性地壳的垂直形变,并与国际地球参考架ITRFs解和VLBI全球解GLB2003,VTRF2003和VTRF2005的结果进行了比较。结果表明,欧亚板块的URUMQI站和太平洋板块的KWAJAL26站,南极OHIGGINS站的垂直形变率、ITRFs解和VLBI全球解存在6-15mm/a的差异,北美YUMA站可能有15-31mm/a 的垂直形变率,而美国西部太平洋板块的San Francisco(PRESIDIO)站的垂直形变率还有待进一步的研究。此外,SC-VLBA,CRIMEA和EFLSBERG站的垂直形变率、ITRFs解和VLBI全球解的差约为1-6mm/a。用不同方法得到的VLBI站的水平形变率解有较好的一致性。     

并置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解的可能不足．     

基于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的高精度应用和长期维持、地球自转参数的长期变化都将产生一定的影响。     

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

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

亚洲区域构造块体运动模型的建立及其运动学分析

利用中国地壳运动监测网、亚太地球动力学计划(APSG)联测以及IGS速度场建立亚洲构造块体运动模型,分析其板内形变,并求出各构造块体相对欧亚板块运动的欧拉矢量,进而求出各构造块体绝对运动速度,分析各构造块体绝对运动,揭示亚洲区域构造块体的运动规律和特征．     

对2004年苏门答腊-安达曼岛地震的测地学机制讨论

采用GPS测地学观测资料,基于NUVEL-1A板块运动模型,分别估计和讨论了INDI板块与Burma微板块之间、AUST板块与Burma/SE EURA板块之间的聚敛率,以及INDI与AUST板块间的相对速率。结果显示,沿苏门答腊-尼科巴群岛和沿安达曼群岛两个阶段的地壳破裂分别与AUST板块和Burma板块间、INDI板块和Burma板块间的聚敛率密切相关,从板块运动学角度提供了解释整个苏门答腊-安达曼群岛地壳破裂带上板块间耦合由南向北递减观测结论的佐证;在NUVEL-1A板块运动模型的背景下,2004年苏门答腊-安达曼岛主震应当为AUST板块相对于Burma/SE EURA板块以近60 mm/a,方向NE 30°～35°的运动所致,而非INDI板块或印-澳(INDI-AUST)板块运动所致;根据2004年苏门答腊-安达曼岛地壳破裂特征并结合其余震的分布,进一步肯定了AUST板块西北边界带的地质位置分布。     

ITRF2005简介和评析

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

印、澳板块边界分布与GPS资料用于对2004年Sumatra-Andaman岛地壳破裂测地学机制的再讨论

基于作者此前对印度尼西亚苏门答腊-安达曼岛地震测地学机制的讨论，进一步根据McCloskey等人应用的INDI板块与AUST板块之间大致呈SW（≈S5）-NE（≈N2）走向的板块边界划分模型，利用GPS站速度资料，估计了沿整个苏门答腊-尼科巴岛地壳破裂带INDI板块与Burma／SEEURA板块之间的板块聚敛率，并与基于NUVEL-1A模型的INDI、AUST板块边界的分布所估计的AUST板块与Burma／SEEURA板块问聚敛率结果作了比较；根据2004年印尼苏门答腊-安达曼岛地壳破裂特征，对该地壳破裂的测地学机制作了进一步讨论，认为基于NUVEL-1A模型的INDI、AUST板块边界分布条件下地壳破裂的测地学机制讨论和2004年主震由AUST板块运动触发的结论可能更符合2004年地壳破裂事件的观测结论。目前的Burma微板块的欧拉旋转参数可能不适合用于对该次地壳破裂事件的测地学估计和机制的讨论；INDI、AUST板块间的相对运动是这两者板块间造山带的主要动力源。     

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

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

Impact of the Quadrupole Moment of the Sun on the Dynamics of the Earth-Moon System

Range of values of the Sun's mass quadrupole moment of coefficient J₂ arising both from experimental and theoretical determinations enlarge across literature on two orders of magnitude, from around 10^-7 until to 10^-5. The accurate knowledge of the Moon's physical librations, for which the Lunar Laser Ranging data reach an outstanding precision level, prove to be appropriate to reduce the interval of J₂ values by giving an upper bound of J₂. A solar quadrupole moment as high as 1.1 10^-5 given either from the upper bounds of the error bars of the observations, or from the Roche's theory, is not compatible with the knowledge of the lunar librations accurately modeled and observed with the LLR experiment. The suitable values of J₂ have to be smaller than 3.0 10^-6. As a consequence, this upper bound of 3.0 10^-6 is accepted to study the impact of the Sun's quadrupole moment of mass on the dynamics of the Earth-Moon system. Such as effect (with J₂ = 5.5±1.3 × 10^-6) has been already tested in 1983 by Campbell & Moffat using analytical approximate equations, and thus for the orbits of Mercury, Venus, the Earth and Icarus. The approximate equations are no longer sufficient compared with present observational data and exact equations are required. As if to compute the effect on the lunar librations, we have used our BJV relativistic model of solar system integration including the spin-orbit coupled motion of the Moon. The model is solved by numerical integration. The BJV model stems from general relativity by using the DSX formalism for purposes of celestial mechanics when it is about to deal with a system of n extended, weakly self-gravitating, rotating and deformable bodies in mutual interactions. The resulting effects on the orbital elements of the Earth have been computed and plotted over 160 and 1600 years. The impact of the quadrupole moment of the Sun on the Earth's orbital motion is mainly characterized by variations of , , and . As a consequence, the Sun's quadrupole moment of mass could play a sensible role over long time periods of integration of solar system models. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the role of the motion of the photospheric footpoints of the magnetic field lines

By means of a simple relation between the velocity v of the fluid particle and the velocity v_f of the photospheric footpoint of the magnetic field line v_z and B_z being respectively the components of v and the magnetic field B normal to the photospheric surface, it is shown formally that through the phtospheric surface the transport of all the quantities attributed to the magnetic field, such as the magnetic flux, the magnetic energy and the helicity, is independent of v_z, and v_f is the only kinematical quantity on which the transport depends. In addition, in the neighborhood of the neutral line the velocity v_l of the moving curve of constant B_z is found to be equal approximately to the component of v or v_f in the direction of v_l. Since v_l can be measured or extimated, so can the components of v and v_f near the neutral line.     

在后牛顿引力理论(PPN形式)中引力常数变化对地球自转产生的效应

介绍和论述了在后牛顿引力理论(PPN形式)中在优越参考系和非优越参考系中经过参数化后引力常数变化对地球自转产生的效应,其中特别重点介绍了年周期变化的效应。此外也将理论结果同观测结果相对比。     

On the theory of the oblateness of the Sun

In this paper we first emphasize why it is important to know the successive zonal harmonics of the Sun's figure with high accuracy: mainly fundamental astrometry, helioseismology, planetary motions and relativistic effects. Then we briefly comment why the Sun appears oblate, going back to primitive definitions in order to underline some discrepancies in theories and to emphasize again the relevant hypotheses. We propose a new theoretical approach entirely based on an expansion in terms of Legendre's functions, including the differential rotation of the Sun at the surface. This permits linking the two first spherical harmonic coefficients (J ₂ and J ₄) with the geometric parameters that can be measured on the Sun (equatorial and polar radii). We emphasize the difficulties in inferring gravitational oblateness from visual measurements of the geometric oblateness, and more generally a dynamical flattening. Results are given for different observed rotational laws. It is shown that the surface oblateness is surely upper bounded by 11 milliarcsecond. As a consequence of the observed surface and sub-surface differential rotation laws, we deduce a measure of the two first gravitational harmonics, the quadrupole and the octopole moment of the Sun: J ₂=−(6.13±2.52)×10⁻⁷ if all observed data are taken into account, and respectively, J ₂=−(6.84±3.75)×10⁻⁷ if only sunspot data are considered, and J ₂=−(3.49±1.86)×10⁻⁷ in the case of helioseismic data alone. The value deduced from all available data for the octopole is: J ₄=(2.8±2.1)×10⁻¹². These values are compared to some others found in the literature. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1005238718479     

Third Integral and the Dynamics of the Galaxy in the Neighborhood of the Sun

Observational data on the dynamics of stars in the neighborhood of the sun indicate the existence of a third integral besides the integrals of the angular momentum and energy. The Poincaré integral is proposed as a third integral. The consequences of this assumption are derived and compared with available astrophysical data.     

射电日像仪灵敏度测量

灵敏度是射电望远镜的一个重要性能指标,它反映了望远镜监测弱信号的能力。基于明安图射电频谱日像仪(Mingantu Ultrawide Spectral Radioheliograph,MUSER)的调试观测,给出了日像仪灵敏度的测量方法,对天线系统以及整个阵列的灵敏度进行测量分析,得到了日像仪系统整体的灵敏度性能参数。测量同时给出了天线系统的效率以及接收机系统的增益,这将为下一步日像仪展开常规的科学观测提供参考。     

Accounting for the viscosity of the fluid core in the problem of the physical libration of the moon

A two-component theoretical model of the physical libration of the Moon in longitude is constructed with account taken of the viscosity of the core. In the new version, a hydrodynamic problem of motion of a fluid filling a solid rotating shell is solved. It is found that surfaces of equal angular velocity are spherical, and a velocity field of the fluid core of the Moon is described by elementary functions. A distribution of the internal pressure in the core is found. An angular momentum exchange between the fluid core and solid mantle is described by a third-order differential equation with a right-hand side. The roots of a characteristic equation are studied and the stability of rotation is proved. A libration angle as a function of time is found using the derived solution of the differential equation. Limiting cases of infinitely large and infinitely small viscosity are considered and an effect of lag of a libration phase from a phase of action of an external moment of forces is ascertained. This makes it possible to estimate the viscosity and sizes of the lunar fluid core from data of observations.     

Calculation of the global system of the field-aligned currents on the dayside of the magnetosphere

Using the well-known equation for the normal component of the current which exist near the tangential discontinuity in the plasma in the case of the frozen-in magnetic field, and supposing that the current closes in the ionosphere in the auroral oval in the region 1, one calculates and compares with the data of observations the dependence of the density of the field-aligned current at the level of the ionosphere on the local time.     

Rotation of the elastic Earth: the role of the angular-velocity-dependence of the elasticity-caused perturbation

We calculate the so-called convective term, which shows up in the expression for the angular velocity of the elastic Earth, within the Andoyer formalism. The term emerges due to the fact that the elasticity-caused perturbation depends not only on the instantaneous orientation of the Earth but also on its instantaneous angular velocity. We demonstrate that this term makes a considerable contribution into the overall angular velocity. At the same time the convective term turns out to be automatically included into the correction to the nutation series due to the elasticity, if the series is defined by the perturbation of the figure axis (and not of the rotational axis) in accordance with the current IAU resolution. Hence it is not necessary to take the effect of the convective term into consideration in the perturbation of the elastic Earth as far as the nutation is related to the motion of the figure axis.