非洲磁异常对地磁场结构及其长期变化的控制作用

地球非偶极磁场在主磁场结构及其长期变化中起着重要作用.非偶极磁场主要表现为行星尺度磁异常,它们是南大西洋磁异常、非洲磁异常、欧亚大陆磁异常、澳洲磁异常和北美磁异常.在这5块磁异常中,非洲磁异常对磁赤道的形状和位置以及全球长期变化特征有极大的影响.非洲磁异常的重要性主要表现在3方面：第一,由于异常区位于赤道这一特殊的地理位置,所以它极大地影响磁赤道的形状和位置.相对于偶极场的地磁赤道而言,异常区所在的中北非洲和中大西洋地区的磁赤道向北移动,最大移动量可达约15°.第二,非洲磁异常的快速西漂对全球长期变化的分布起着决定性作用,它在该异常区西边的中美洲形成了全球最主要的长期变化区,在1900~2005年期间,最大年变率Z_max超过200 nT/a.第三,非洲负磁异常区与其南面的南大西洋正磁异常区相结合,〖HJ〗它们的变化使西半球地磁场强度大大减弱,也使全球磁场发生显著畸变.这两块磁异常与深部的反极性斑区有着成因联系.     

地球主磁场的NOC模型

从1900~2000年国际参考地磁场球谐模型系列出发, 建立一种新的地磁场模型──自然正交分量模型, 简称NOC模型. 首先, 由IGRF1900-2000的高斯系数求出地磁场的本征模; 然后, 以此作为基函数系, 将每一年代的磁场展开, 求出各本征模的强度系数, 即可得到表示地球主磁场空间结构和时间变化的NOC模型. 对NOC模型的收敛特征和基函数的稳定性进行了数值检验, 结果表明, 与传统的IGRF球谐模型相比, NOC模型具有级数短、收敛快的特点；在所研究的100年内, 低阶NOC基函数比较稳定, 而高阶基函数有较大的变化. 讨论了NOC模型基函数的物理内涵, 并揭示了地磁场的空间结构及其长期变化之间的内在联系.     

1900-2000年非偶极子磁场的全球变化

根据第７代ＩＧＲＦ模型,计算并绘制了１９００-２０００年（时间间隔为５ａ）非偶极子磁场的全球分布,计算和分析了各个磁异常中心位置和强度的变化,东亚正磁异常、大洋洲负磁异常、南大西洋正磁异常、非洲负磁异常和北美洲正磁异常是分布范围广、异常强度大的５个磁异常。南大西洋正磁异常是强度最大的磁异常。在２０世纪９０年代以前,东亚正磁异常的强度位居第２位,９０年代以后,非洲负磁异常的强度（绝对值）超过东亚正磁异常,成为第２大磁异常。磁异常强度增长最快的是非洲负磁异常、南大西洋正磁异常和大洋洲负磁异常。南大西洋正磁异常和非洲负磁异常是磁异常中心位置变化最快的两个磁异常。     

行星尺度地磁异常的长期变化

为了综合反映地球表面行星尺度磁异常的展布面积、磁场极值以及磁场分布特征等多种因素及其与磁能的关系，本文用穿过各异常区的“无符号磁通量”为特征参数来表征磁异常区强度.用第八代国际参考地磁场模型(IGRF)，分析了1900年到2000年全球最大的5个磁异常区的长期变化，结果表明，在一百年中，南大西洋(SAT)、大洋洲(AUS)和非洲(AF)3个异常区的磁通量均增加了200MWb以上，欧亚异常(EA)磁通量增加幅度稍小(157MWb)，上述4个异常区磁通量增幅为30%-60%，而北美异常(NAM)的磁通量则减小了50MWb.各异常区面积虽有变化，但最大变化仅为%左右.对磁异常区的西向漂移研究表明，地球表面和核幔界面的西漂明显存在差异：地表磁场有持续而稳定的西向漂移，全球平均西漂速度为0.2°/a；但核幔界面磁场的西向漂移速度要小得多，最大不超过0.1°/a.形成这种差异的原因可能是组成地磁场的不同球谐分量有不同的漂移速度；地表磁场的西漂主要决定于占优势的低阶分量，而核幔界面的西漂则受到高阶分量的重大影响.本文指出，在把地表西漂值用作地核磁流体运动速度的典型值时必须十分谨慎.     

地球主磁场的能量密度谱及其长期变化

根据BJ地磁场模型和第8代国际地磁参考场(IGRF)模型，计算并分折了1690～2000年期间地球表面主磁场和分量的能量密度谱及其长期变化。结果表明1690～2000年期间主磁场的能量密度谱一直在减小，其主要原因是由于偶板子磁场的衰减产生的，而非偶板子磁场的能量密度谱在1690～1780年期间减小，1780～1890年增大，1890年以后快速增大。在地磁场总能谱的变化中，Z分量的贡献起主要作用。能量密度随谐波阶数的变化在半对数坐标中近似线性，非偶板子场的等效磁源深度位于核幔边界附近，n＝1，2谐波项的等效磁源位于地球内核边缘，其位置随时间变化。1780和1890年前后是地磁场变化的转折时期。能量密度谱的长期变化存在大约60a的周期规律。     

地磁场的漂移运动和强度变化

在修正了Briggs提出的移动变形图案相关分析法基础上,对全球非偶极子磁场以及6个行星尺度磁异常区的漂移特性和强度变化进行了研究.结果表明,在1900-2000年期间,全球非偶极子磁场以0.15°/a的平均速度向西漂移,强度累计增长了29%;6个行星尺度异常区的西漂运动存在明显差异,其中漂移最快的是赤道附近的非洲异常,平均西漂速度为0.26°/a,其次是南半球的澳大利亚异常(0.23°/a),最慢的是欧亚异常(0.09°/a).除西漂外,大多数异常区还有较小的北向漂移.在1940-1955年期间北半球的欧亚、北美和北大西洋3个异常区以及赤道地区的非洲异常几乎同时由西漂或西南漂转为西北向漂移;紧接着,南半球的南大西洋和澳大利亚两个异常区的漂移特征也发生明显变化,主要是漂移明显减慢,而不是漂移方向的转折.在6个异常区中,澳大利亚、南大西洋、非洲和欧亚4个异常区的强度有明显增加,而北美和北大西洋两个异常区的强度则显示了减小的趋势.     

1900—2010年地磁场水平分量梯度的全球变化

本文根据最新的国际地磁参考场模型(IGRF11),结合梯度法分析能够很好的表现细微变化部分特征的优点,计算和绘制了1900—2010年期间间隔5年的全球地磁场水平分量H的水平梯度和垂直梯度,分析其时空分布,各磁异常中心位置和强度的变化规律。结果表明,地磁场的各个异常中心位置在水平方向和垂直方向上都有变化,主要以西漂为主,南大西洋正磁异常区水平分量梯度异常中心点强度值变化最大。     

1900—2010年地磁场水平分量梯度的全球变化

本文根据最新的国际地磁参考场模型（IGRF11）,结合梯度法分析能够很好的表现细微变化部分特征的优点,计算和绘制了1900—2010年期间间隔5年的全球地磁场水平分量H的水平梯度和垂直梯度,分析其时空分布,各磁异常中心位置和强度的变化规律。结果表明,地磁场的各个异常中心位置在水平方向和垂直方向上都有变化,主要以西漂为主,南大西洋正磁异常区水平分量梯度异常中心点强度值变化最大。     

地磁场长期变化和日长十年尺度变化的周期特征

根据历史地磁场模型GUFM1、第10代国际参考地磁场（IGRF10）模型和日长资料,采用小波变换方法,分析了地磁场磁矩、能量、西向漂移等参数的长期变化和日长十年尺度变化的周期分量及其时变特征.结果表明,1800～2005年期间,偶极子磁场长期变化有82年和48年准周期分量,它们与日长变化的周期没有直接关系.非偶极子磁场参数的长期变化与日长变化有66年和32年准周期分量,66年准周期比32年准周期强.在66年准周期分量,西向漂移比日长变化超前8.8年,非偶极子磁场能量比日长变化滞后15.6年.日长十年尺度波动和地磁场长期变化的起源不存在因果关系.     

冀鲁豫交界区流动地磁场总强度异常特征及地震活动性分析

运用2014年7月至2016年4月冀鲁豫交界区开展的流动地磁9期次测量资料,分析相邻期、半年期、一年期的地磁场变化特征。结果表明,除2015年11月、12月地磁场总强度相邻期异常变化较大外,其余期次的相邻期地磁场变化较小,在10nT左右;半年、一年尺度磁异常整体表现为上升。将范县附近8个测点磁异常变化与地震序列进行对比分析发现,磁异常变化梯度较大时,范县易发生小震群活动。     

NOC model of the earth's main magnetic field

The international geomagnetic reference field (IGRF) is a standard model for describing the spatial structure and temporal variation of the earth抯 main magnetic field[1—3]. The first IGRF model, designated IGRF 1965, was adopted by IAGA in 1968[4]. In l…     

Secular variations in the main geomagnetic field within the scope of the dynamic model of field sources

The effect of different-level sources on the spatial structure of the secular variations has been considered based on the dynamic model of sources of the main geomagnetic field developed by us. It has been obtained that the development of 13 most powerful dipoles only roughly characterizes global anomalies of the secular variations, and each anomaly results from the superposition of the dynamics of several sources. The model secular variations have been compared with the data from the observatories. It has been obtained that it is impossible to describe local anomalies of the observed secular variations ignoring sources of the third order of smallness as compared to the main dipole. It has been assumed that topographic vortices, originating around inhomogeneities of the core-mantle boundary, can be physical sources responsible for dipoles of the third order.     

Global upper ocean heat content and climate variability

Observational data from the Global Temperature and Salinity Profile Program were used to calculate the upper ocean heat content (OHC) anomaly. The thickness of the upper layer is taken as 300 m for the Pacific/Atlantic Ocean and 150 m for the Indian Ocean since the Indian Ocean has shallower thermoclines. First, the optimal spectral decomposition scheme was used to build up monthly synoptic temperature and salinity dataset for January 1990 to December 2009 on 1° × 1° grids and the same 33 vertical levels as the World Ocean Atlas. Then, the monthly varying upper layer OHC field (H) was obtained. Second, a composite analysis was conducted to obtain the total-time mean OHC field ([`([`(H)])] \bar{\bar{H}} ) and the monthly mean OHC variability ( [(\textH)\tilde] \widetilde{\text{H}} ), which is found an order of magnitude smaller than [^(\textH)] \widehat{\text{H}} . Third, an empirical orthogonal function (EOF) method is conducted on the residue data ( [^(\textH)] \widehat{\text{H}} ), deviating from [(\textH)\tilde] \widetilde{\text{H}}  +  [(\textH)\tilde] \widetilde{\text{H}} , in order to obtain interannual variations of the OHC fields for the three oceans. In the Pacific Ocean, the first two EOF modes account for 51.46% and 13.71% of the variance, representing canonical El Nino/La Nina (EOF-1) and pseudo-El Nino/La Nina (i.e., El Nino Modoki; EOF-2) events. In the Indian Ocean, the first two EOF modes account for 24.27% and 20.94% of the variance, representing basin-scale cooling/warming (EOF-1) and Indian Ocean Dipole (EOF-2) events. In the Atlantic Ocean, the first EOF mode accounts for 49.26% of the variance, representing a basin-scale cooling/warming (EOF-1) event. The second EOF mode accounts for 8.83% of the variance. Different from the Pacific and Indian Oceans, there is no zonal dipole mode in the tropical Atlantic Ocean. Fourth, evident lag correlation coefficients are found between the first principal component of the Pacific Ocean and the Southern Oscillation Index with a maximum correlation coefficient (0.68) at 1-month lead of the EOF-1 and between the second principal component of the Indian Ocean and the Dipole Mode Index with maximum values (around 0.53) at 1–2-month advance of the EOF-2. It implies that OHC anomaly contains climate variability signals.     

Variability of chlorophyll-a from ocean color images in the La Plata continental shelf region

Satellite-derived chlorophyll-a fields have been used to investigate temporal and spatial variability of chlorophyll-a concentration over the continental shelf zone (25–40°S and 60–45°W) close to the La Plata River estuary. Ocean color data used in this study were obtained by the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) and consisted of 368 weekly averaged Standard Mapped Images (SMI), from October 1997 to September 2005. Fourier harmonic and EOF analyses were used to study the variability of log-transformed chlorophyll-a concentration in the region. The harmonic analysis has shown that the annual cycle was the most dominant signal followed by the semi- and quadri-annual cycles, in certain areas. The strong annual cycle is mainly present in latitudes lower than 34°S where relatively high amplitudes (∼1.9 mg m⁻³) in pigment variation are seen over the southern Brazilian continental shelf. The semi-annual cycle is mainly associated with the Brazil–Malvinas frontal zone oscillation while the 4-year signal is related to positive La Plata discharge anomalies influenced by El Niño events. After removing the annual signal from the log-transformed chlorophyll anomalies, the EOF results showed that the first three modes captured 85.1% of the variability associated with the regional mean phytoplankton chlorophyll pattern in our smoothed data set. The first three modes explained, respectively, 63.4%, 14.1% and 7.6%. The EOF results showed that the long-term chlorophyll time/space patterns are associated with both La Plata discharge anomaly (mode 1) and alongshore wind stress (mode 2). A reconstruction of the chlorophyll anomaly fields has been made using the two leading EOF modes over two periods of high La Plata River discharge, during ENSO events. In the first event, the spatial patterns of high chlorophyll anomaly were confined to the southern portion of the region, associated with NE winds, which push the plume near the estuary mouth. The second period revealed an elongated tongue of positive chlorophyll anomalies over the Uruguayan and Brazilian middle continental shelves, associated with favorable SW winds. The analyses performed in this study allowed identification of the main modes of variability in SeaWiFS-derived chlorophyll in the region, which were consistent with modulations of important regional environmental forcing mechanisms.     

Structure of the Indian Ocean Meridional Overturning Circulation and its relationship with the zonal wind stress

We studied the structure of the Indian Ocean(IO)Meridional Overturning Circulation(MOC)by applying a nonlinear inertia theory and analyzed the coupled relationship between zonal wind stress and MOC anomalies.Our results show that the inertia theory can represent the main characteristics of the IO MOC:the subtropical cell(STC)and cross-equator cell(CEC).The stream function in equatorial and northern IO changes a sign from winter to summer.The anomalies of the zonal wind stress and stream function can be decomposed into summer monsoon mode,winter monsoon mode,and abnormal mode by using the singular vector decomposition(SVD)analysis.The first two modes correlate with the transport through 20°S and equator simultaneously whereas the relationship obscures between the third mode and transports across 20°S and equator,showing the complex air-sea interaction process.The transport experiences multi-time scale variability according to the continuous power spectrum analysis,with major periods in inter-annual and decadal scale.     

Anisotropic tomography of the Atlantic Ocean

We present the first regional three-dimensional model of the Atlantic Ocean with anisotropy. The model, derived from Rayleigh and Love wave phase velocity measurements, is defined from the Moho down to 300 km depth with a lateral resolution of about 500 km and is presented in terms of average isotropic S-wave velocity, azimuthal anisotropy and transverse isotropy.The cratons beneath North America, Brazil and Africa are clearly associated with fast S-wave velocity anomalies. The mid-Atlantic ridge (MAR) is a shallow structure in the north Atlantic corresponding to a negative velocity anomaly down to about 150 km depth. In contrast, the ridge negative signature is visible in the south Atlantic down to the deepest depth inverted, that is 300 km depth. This difference is probably related to the presence of hot-spots along or close to the ridge axis in the south Atlantic and may indicate a different mechanism for the ridge between the north and south Atlantic. Negative velocity anomalies are clearly associated with hot-spots from the surface down to at least 300 km depth, they are much broader than the supposed size of the hot-spots and seem to be connected along a north-south direction.Down to 100 km depth, a fast S-wave velocity anomaly is extenting from Africa into the Atlantic Ocean within the zone defined as the Africa superswell area. This result indicates that the hot material rising from below does not reach the surface in this area but may be pushing the lithosphere upward.In most parts of the Atlantic, the azimuthal anisotropy directions remain stable with increasing depth. Close to the ridge, the fast S-wave velocity direction is roughly parallel to the sea floor spreading direction. The hot-spot anisotropy signature is striking beneath Bermuda, Cape Verde and Fernando Noronha islands where the fast S-wave velocity direction seems to diverge radially from the hot-spots.The Atlantic average radial anisotropy is similar to that of the PREM model, that is positive down to about 220 km, but with slightly smaller amplitude and null deeper. Cratons have a lower than average radial anisotropy. As for the velocities, there is a difference between north and south Atlantic. Most hot-spots and the south-Atlantic ridge are associated with positive radial anisotropy perturbation whereas the north-Atlantic ridge corresponds to negative radial anisotropy perturbation.     

The relationship between significant wave height and Indian Ocean Dipole in the equatorial North Indian Ocean

Based on reanalysis data, we find that the Indian Ocean Dipole (IOD) plays an important role in the variability of wave climate in the equatorial Northern Indian Ocean (NIO). Significant wave height (SWH) in the equatorial NIO, especially over the waters southeast to Sri Lanka, exhibits strong interannual variations. SWH anomalies in the waters southeast to Sri Lanka correlate well with dipole mode index (DMI) during both summer and autumn. Negative SWH anomalies occur over the oceanic area southeast to Sri Lanka during positive IOD events and vary with different types of IOD. During positive prolonged (unseasonable) IOD, the SWH anomalies are the strongest in autumn (summer); while during positive normal IOD, the SWH anomalies are weak in both summer and autumn. Strong easterly wind anomalies over the southeast oceanic area of Sri Lanka during positive IOD events weaken the original equatorial westerly wind stress, which leads to the decrease in wind-sea waves. The longer wave period during positive IOD events further confirms less wind-sea waves. The SWH anomaly pattern during negative IOD events is nearly opposite to that during positive IOD events.     

20世纪地磁长期变化场分析

利用IAGA（国际地磁与高空物理学协会）编制的IGRF(国际地磁参考场)研究了20世纪地磁场变化规律. 20世纪地磁长期变化场的四极子（n=2的高斯系数所表示）变化最为显著,与主磁场相比长期变化场的球谐级数收敛较慢,利用追踪异常焦点位置随时间的变化的方法,发现地磁非偶极子长期变化场的垂直分量Z的等值图上有五大异常,其漂移情况不太统一,但是基本上是西向漂移.这种西向漂移的不一致性,在某种程度上证明了地磁场模型的正确性. 20世纪地磁场长期变化的能量谱与主磁场的不同,偶极子、 四极子和八极子的变化较明显.