地球磁场偶极子和非偶极子分量变化特征及其相关的地球深部过程

通过分析１９００年以来每５ａ地球磁场各阶球谐系数的相关系数发现,地球磁场偶极子分量的ρ值在１９０５-１９７０年基本保持不变,１９７０-１９８５年下降,１９８５-１９９５年又上升,呈Ｖ字型变化。而非偶极子分量的　值则具有周期性变化,在１９１５年、１９４５年和１９８０年出现峰值,在１９３５年和１９６５年则处于低谷,变化周期约为３０ａ。通过分析地球磁场非偶极子分量　值周期变化与日长周期变化的相关性,得出两点新的认识：１．非偶极子分量的周期性变化与核幔边界耦合过程相关,而偶极子分量的变化对应外核的深部过程;２．外核顶部流体具有周期性的变化,其变化周期与日长及非偶极子分量的变化周期一致。     

核幔耦合对地球自由核章动的激发影响

地球自由核章动（FCN）是地幔与液核相互作用的重要动力学现象,其激发机制涉及地表流体层、地幔和地核等圈层之间的耦合,此前研究多利用地表流体层角动量数据单独研究其对FCN的激发,对核幔耦合的影响考虑不足.本文基于角动量守恒理论分析了核幔耦合对FCN周期及振幅的影响,并结合多个大气及海洋角动量函数时间序列首次估算了核幔耦合在FCN激发过程中的贡献.结果表明核幔耦合对FCN周期产生的固定和时变影响对FCN激发的作用均不可忽视,尤其时变影响可达几十个微角秒,对于进一步解释FCN时变特征非常重要;核幔耦合对FCN振幅的直接影响是地表流体层的激发与实测FCN不相符的主要原因,黏滞、电磁和地形等耗散耦合的存在对地表流体的激发振幅有67%左右的减弱效果.     

核－幔边界的动力学背景

根据传播矩阵方法，并把由联合反演得到的同时满足长波地形起伏、板块运动速度、重力位异常资料以及地震层析先验知识的全地幔三维异常密度作为载荷，以求取核－幔边界的动力学背景．计算结果显示：１．所求的核－幔边界起伏图像与Ｈａｇｅｒ等根据格林函数方法所求得的核－幔边界的起伏在全球范围内基本相符．２．核－幔边界处的环型场仅在数量上降为地表处的环型场的１／８左右，而极型场较地表处的极型场的流动图像有显著变化，数值也增为地表处极型场的３倍左右．     

Ray parameters of diffracted long period P and S waves and the velocities at the base of the mantle

The derivation of P and S velocities at the core-mantle boundary (CMB) from long-period diffracted waves by the use of the simple ray-theoretical formulav _CMB=r _c /p (v _CMB=velocity at the CMB;r _c =core radius;p=ray parameter) yields apparent velocity values which differ from the true velocities. Using a dominant period of about 20 sec for calculating theoretical seismograms, we found a linear relation between the apparent velocity and the average velocity in a transition zone at the base of the mantle with fixed velocity on top.The ray parameters determined from long-period earthquake data are found to be 4.540±0.035 and 8.427±0.072 sec/deg for P_diff and S_diff, respectively. These values yield apparent velocities of 13.378±0.103 for P and 7.207±0.062 km/sec for S waves. By means of the theoretical relation between apparent and average velocity and under the assumption of linear variation of velocity with depth, one can invert the apparent velocities into true CMB velocities of 13.736±0.170 and 7.320±0.124 km/sec. These results imply positive velocity gradients at the base of the mantle and hence no significant departures from adiabaticity and homogeneity.Contribution No. 211 of the Geophysical Institute, University of Karlsruhe.     

A comparison of the reflectivity method and full wave theory as applied to SH diffraction around the core

Summary To derive the structure of the Lower Mantle, observations of P and S waves at epicentral distances beyond 90 degrees have been used. The arrival times and amplitude decay with distance proved to be a source of information to derive the structure near the core-mantle interface. These observations were compared with results of theoretical calculations based on full wave theory [1, 2, 5, 6]. Another important theoretical method is the reflectivity method [4]. Differences between the two methods have a direct bearing on the interpretation of the observations and therefore need to be investigated.In the two methods the wave propagation involved in diffraction is simulated in different ways. In this paper the importance of approximations involved in each of the methods will be discussed. For this purpose the calculation of the frequency-dependent decay with distance of diffracted SH waves has been chosen. The approximations are considered to be important both in P and SH diffraction. The advantage of studying SH waves is the simpler and more straightforward formulation.     

Core mantle boundary topography from short period PcP, PKP, and PKKP data

We use a total of 839,369 PcP, PKPab, PKPbc, PKPdf, PKKPab, and PKKPbc residual travel times from [Bull. Seism. Soc. Am. 88 (1998) 722] grouped in 29,837 summary rays to constrain lateral variation in the depth to the core-mantle boundary (CMB). We assumed a homogeneous outer core, and the data were corrected for mantle structure and inner-core anisotropy. Inversions of separate data sets yield amplitude variations of up to 5 km for PcP, PKPab, PKPbc, and PKKP and 13 km for PKPdf. This is larger than the CMB undulations inferred in geodetic studies and, moreover, the PcP results are not readily consistent with the inferences from PKP and PKKP. Although the source-receiver ambiguity for the core-refracted phases can explain some of it, this discrepancy suggest that the travel-time residuals cannot be explained by topography alone. The wavespeed perturbations in the tomographic model used for the mantle corrections might be too small to fully account for the trade off between volumetric heterogeneity and CMB topography. In a second experiment we therefore re-applied corrections for mantle structure outside a basal 290 km-thick layer and inverted all data jointly for both CMB topography and volumetric heterogeneity within this layer. The resultant CMB model can explain PcP, PKP, and PKKP residuals and has approximately 0.2 km excess core ellipticity, which is in good agreement with inferences from free core nutation observations. Joint inversion yields a peak-to-peak amplitude of CMB topography of about 3 km, and the inversion yields velocity variations of ±5% in the basal layer. The latter suggests a strong trade-off between topography and volumetric heterogeneity, but uncertainty analyses suggest that the variation in core radius can be resolved. The spherical averages of all inverted topographic models suggest that the data are best fit if the actual CMB radius is 1.5 km less than in the Earth reference model used (i.e. the average outer core radius would be 3478 km).     

Secular variation in numerical geodynamo models with lateral variations of boundary heat flow

We study magnetic field variations in numerical models of the geodynamo, with convection driven by nonuniform heat flow imposed at the outer boundary. We concentrate on cases with a boundary heat flow pattern derived from seismic anomalies in the lower mantle. At a Rayleigh number of about 100 times critical with respect to the onset of convection, the magnetic field is dominated by the axial dipole component and has a similar spectral distribution as Earth’s historical magnetic field on the core-mantle boundary (CMB). The time scales of variation of the low-order Gauss coefficients in the model agree within a factor of two with observed values. We have determined the averaging time interval needed to delineate deviations from the axial dipole field caused by the boundary heterogeneity. An average over 2000 years (the archeomagnetic time scale) is barely sufficient to reveal the long-term nondipole field. The model shows reduced scatter in virtual geomagnetic pole positions (VGPs) in the central Pacific, consistent with the weak secular variation observed in the historical field. Longitudinal drift of magnetic field structures is episodic and differs between regions. Westward magnetic drift is most pronounced beneath the Atlantic in our model. Although frozen flux advection by the large-scale flow is generally insufficient to explain the magnetic drift rates, there are some exceptions. In particular, equatorial flux spot pairs produced by expulsion of toroidal magnetic field are rapidly advected westward in localized equatorial jets which we interpret as thermal winds.     

地球早期演化的Hf-W同位素制约

~(182)Hf-~(182)W作为短周期放射性衰变体系,可有效约束地球早期演化吸积增生和内部分异过程。本文通过系统总结、归纳太阳系形成初期地球核幔分异过程中Hf-W同位素变化规律、月球与地球硅酸盐的W同位素组成,提出利用~(182)Hf-~(182)W体系测定地球核幔分异时间不确定性的主要原因是地球核幔分异的持续性及开放性,大碰撞时间的~(182)Hf-~(182)W同位素限定主要受控于硅酸盐地球和硅酸盐月球的Hf/W比值,讨论了地幔W同位素不均一性的形成机制,与现代地幔不同的~(182)W/~(184)W组成可能代表了后增生作用之前整体硅酸盐地球的W同位素组成,也可能是~(182)Hf未完全灭绝时形成的区域性Hf/W比值差异经~(182)Hf衰变形成的结果。这些结论为探索类地行星形成与演化提供了重要制约。     

利用VLBI和超导重力资料研究自由核章动周期时变特征

本文根据自由核章动在章动观测和时变重力观测中引起的与其频率相近的受迫章动项或固体潮潮波的共振特性,选取6个不同机构解算的VLBI天极偏差序列和全球超导重力仪网络提供的7个台站高精度时变重力观测资料,根据加权最小二乘方法拟合地球自由核章动参数,研究其本征周期的时间变化特征.VLBI资料获得的本征周期变化幅度在1天之内,存在十年尺度的时间变化特征,迭积多台站同时段重力资料获得的自由核章动本征周期时变特征与VLBI获得的结果相比变化幅度较大,但趋势大致符合.在此基础上,通过自由核章动的理论模型探讨了影响其本征周期的几个主要因素,结果表明FCN周期变化与电磁耦合存在相关性,核幔边界磁感应强度径向分量变化导致的电磁耦合的变化可能是造成FCN周期时变性的原因.