地面VLF波穿透电离层的能量衰减变化

地面VLF人工源可以导致辐射带高能粒子沉降.为研究辐射带粒子加速和沉降机制,乃至实现对辐射带电子"人工控制"的设想,需要在精确计算得到VLF人工源在电离层激发电磁场能量分布的基础上,计算哨声波对辐射带电子的调制作用.以往计算VLF人工源在电离层中激发的电磁场能量分布多基于Cary给出的波导中VLF波衰减率和Helliwell提出的经典电离层吸收曲线,但近期研究表明,这些模型结果存在较为明显的误差.本研究建立了地面VLF信号穿透电离层传播的全波计算模型,将计算结果与DEMETER卫星记录的NWC通讯台激发的电磁响应进行对比.虽然模型没有考虑电离层参数水平不均匀性,但模型计算结果与卫星观测值也有较好的一致性.利用经过验证的全波模型计算了不同地磁参数、电离层参数情况下,不同辐射特性的地面VLF辐射在波导中的衰减和穿透电离层时D/E区的吸收值,探讨了上述参数对电磁波能量在波导和D/E区中衰减的影响规律.     

On the nonlinear triggering of VLF emissions by power line harmonic radiation

VLF ground data from Porojarvi in N. Finland has been presented. Spectrograms reveal frequent occurrence of power line harmonic radiation (PLHR), originating from the Finnish power system and from heavy industrial plant. The radiation is seen to penetrate the magnetosphere since numerous occurrences of PLHR triggered emissions are seen. Risers predominate but fallers and hooks are also observed. A well-established 1D Vlasov simulation code has been used to simulate these emissions, using plausible magnetospheric data for a range of L values from L = 4 to L = 5.5. The code is able to reproduce risers fallers and hooks in close agreement with observations. The results shed considerable insight into the generation structure of both risers and fallers.     

基于人工电离层扰动对地-电离层波导传播的VLF波的影响

基于低电离层加热理论和甚低频电波在地-电离层波导中传播理论,建立低电离层扰动对甚低频电波传播影响的分析模型,并利用实验数据验证了该模型的正确性.据此模型,研究了加热功率、加热波极化以及背景参数所导致的低电离层扰动对不同频率甚低频电波传播的影响.结果表明,低电离层扰动越强,则通过该区域内甚低频波幅度和相位的相对变化越强,通过研究地-电离层波导甚低频信号通过人工扰动区域后幅度和相位的变化,可望用于诊断人工电离层扰动强度.     

全球电离层VLF电场功率谱特征

分析2006年3月至2009年2月DEMETER卫星VLF电场功率谱数据发现,卫星高度上全球电离层电场功率谱有以下特征:高纬度区域电场辐射强度平均高于低纬度区域,几个地磁场异常区相应的电场辐射增强;全球电场功率谱,向阳侧强于夜侧,大陆强于海洋,夏季强于分点季,冬季最弱;不同频段的电场功率谱特征有显著差异,某些频段的功率...     

基于NWC甚低频信号的日出效应的观测与分析

频率在3~30 kHz的甚低频(VLF,Very Low Frequency)波具有较小的传播损耗和较高的趋肤深度,可以在地球-低电离层波导中实现长距离传输,广泛应用于航海导航、对潜通信等领域,且在电离层遥测方面具有十分重要的意义.基于武汉大学自主研发的VLF接收机在武汉接收的NWC(North West Cape)台站信号,本文通过分析2018年4月23日—2020年7月22日的观测数据研究了日出期间NWC信号的幅度响应及其特点和规律.结果表明NWC信号日出期间的幅度响应主要包括两种极小值结构：2个幅度极小值(SR1、SR2)的Type Ⅰ结构和3个幅度极小值(SR1、SR2、SR3)的Type Ⅱ结构.在以SR1出现时间为时间零点进行时序叠加分析后发现,Type Ⅰ结构比Type Ⅱ具有更强的规律性和稳定性.在Type Ⅰ结构下,SR2出现时间的波动范围、平均值、标准差分别为43~65 min、54.2 min、4.4 min,而在Type Ⅱ结构下,SR2和SR3出现时间的波动范围分别为48~93 min、80~120 min,平均值分别为64.7 min、96.4 min,标准差分别为10.2 min、11.7 min.在27个月的观测期内,3—7月份Type Ⅰ结构的出现概率100%,未出现Type Ⅱ结构,而在1—2月和8—12月Type Ⅰ结构出现的概率明显下降,最低降至1月份的20.7%,而Type Ⅱ在1月、2月、11月的出现概率均高于70%.按春秋分交替变化(周期1和周期2)的统计结果,在周期1内Type Ⅰ和Type Ⅱ结构出现的概率分别为91.5%、8.5%,而在周期2内Type Ⅰ结构出现的概率降至41.9%,Type Ⅱ结构出现概率则升至58.1%,这表示观测期间内Type Ⅱ结构主要出现在秋冬季,春夏季发生概率较低.

     

基于人工源的甚低频电磁波空间传播特征统计分析

利用2005年1月至2010年11月DEMETER卫星记录的NWC发射站的VLF电场功率谱数据,采用指数拟合的方法,分析了VLF电磁波在卫星高度激发的电场空间分布和衰减特征.研究结果表明: (1) VLF电场在发射站上空及其磁共轭区有着很强的对应关系,存在南、北2个强电场中心涡; (2) 相对于发射站的位置,VLF电场中心点具有经度和纬度偏移,日侧地磁经度偏移均值大于夜侧,而地磁纬度偏移均值则小于夜侧; (3) 日侧VLF电场强度呈现出周期性的年变化; (4) 在VLF电场中心10°范围内,电场强度随距离快速衰减,衰减常数b在长达6年的时间内保持稳定.在以上研究结果基础上初步构建的卫星高度人工源电磁波空间分布特征,将为研究地表-电离层电磁波传播机理提供基础技术支撑.     

基于NWC甚低频信号的日出效应的观测与分析

频率在3~30 kHz的甚低频（VLF,Very Low Frequency）波具有较小的传播损耗和较高的趋肤深度,可以在地球-低电离层波导中实现长距离传输,广泛应用于航海导航、对潜通信等领域,且在电离层遥测方面具有十分重要的意义.基于武汉大学自主研发的VLF接收机在武汉接收的NWC（North West Cape）台站信号,本文通过分析2018年4月23日—2020年7月22日的观测数据研究了日出期间NWC信号的幅度响应及其特点和规律.结果表明NWC信号日出期间的幅度响应主要包括两种极小值结构：2个幅度极小值（SR1、SR2）的Type I结构和3个幅度极小值（SR1、SR2、SR3）的Type II结构.在以SR1出现时间为时间零点进行时序叠加分析后发现,Type I结构比Type II具有更强的规律性和稳定性.在Type I结构下,SR2出现时间的波动范围、平均值、标准差分别为43~65 min、54.2 min、4.4 min,而在Type II结构下,SR2和SR3出现时间的波动范围分别为48~93 min、80~120 min,平均值分别为64.7 min、96.4 min,标准差分别为10.2 min、11.7 min.在27个月的观测期内,3—7月份Type I结构的出现概率100%,未出现Type II结构,而在1—2月和8—12月Type I结构出现的概率明显下降,最低降至1月份的20.7%,而Type II在1月、2月、11月的出现概率均高于70%.按春秋分交替变化（周期1和周期2）的统计结果,在周期1内Type I和Type II结构出现的概率分别为91.5%、8.5%,而在周期2内Type I结构出现的概率降至41.9%,Type II结构出现概率则升至58.1%,这表示观测期间内Type II结构主要出现在秋冬季,春夏季发生概率较低.     

强震前ELF/VLF磁场的扰动特征统计研究

本文利用DEMETER卫星记录的变化磁场数据统计研究了2005—2009年北半球7级以上强震前后空间磁场的扰动特征.在震中±10°范围内利用震前2个月至震后1个月的5年同期观测的ELF/VLF[370~897 Hz]频段磁场功率谱密度数据构建了背景场,定义了表征空间磁场在地震时段相对于背景场的扰动幅度指标.统计研究发现26个强震中42%的地震主要表现为震前磁场扰度幅度逐渐上升,超过3倍标准差,随后在扰度幅度下降过程中发震.35%强震在地震前10天内最大扰动幅度超过3倍标准差,在扰度幅度处于最高值期间发震,震后磁场扰动幅度逐渐回落.23%的强震在震前震后无明显规律.随机事件测试发现在空间天气平静且无强震条件下随机地点上空的磁场无明显的时空演化规律,且随机地点上空的空间磁场相对于背景场的平均扰度幅度很小,最大扰度幅度从未超过2倍标准差,并且没有明显的随时间变化特征.     

玉树地震前电离层VLF信号信噪比变化

2010年4月13日23点49分38秒(世界时),青海省玉树县发生了Ms7.1地震.根据法国DEMETER卫星电场探测仪(ICE)实验数据,研究震区上空半径500 km范围内电离层10-20 kHz甚低频(VLF)电场频谱信噪比发现,玉树地震前3个不同频率的地面VLF发射站信号对应的信噪比均出现相同的变化特征,即震前DEMETER卫星1个重访周期内的平均信噪比明显减弱,而2009年相同时段和区域的平均信噪比并未出现相似变化趋势.分析认为,玉树地震孕震期信噪比衰减现象很可能是由地震-电离层耦合所致.     

玉树地震前电离层VLF信号信噪比变化

2010年4月13日23点49分38秒(世界时), 青海省玉树县发生了M_S7.1地震. 根据法国DEMETER卫星电场探测仪(ICE)实验数据, 研究震区上空半径500 km范围内电离层10—20 kHz甚低频(VLF)电场频谱信噪比发现, 玉树地震前3个不同频率的地面VLF发射站信号对应的信噪比均出现相同的变化特征, 即震前DEMETER卫星1个重访周期内的平均信噪比明显减弱, 而2009年相同时段和区域的平均信噪比并未出现相似变化趋势. 分析认为, 玉树地震孕震期信噪比衰减现象很可能是由地震-电离层耦合所致.     

电离层人工调制激发的VLF波在磁层的传播特性及应用研究

电离层人工调制可以激发甚低频（VLF）波,其中向上传播进入磁层的VLF波,不但能够用来研究磁层中的各种物理现象,且具有人工沉降高能粒子,消除辐射带等实际用途.本文使用射线追踪方法,模拟电离层调制激发的VLF波在磁层的传播路径,分析激发纬度和调制频率对传播路径和传播特性的影响;并基于低频波的色散方程和波粒共振条件,分析VLF波传播路径上与磁层高能粒子的最低共振能及其分布.研究表明,VLF波通过在磁层来回反射向更高的L-shell传播,最终稳定在某一L-shell附近.以较低的调制频率或者从较高的纬度激发的VLF波能够传播到更高的L-shell,但是,当激发纬度过高时,低频波也可能不发生磁层反射而直接进入电离层和大气层.低频波在磁层的传播过程中,在较高的纬度或者较低的L-shell能够与较高能量的电子发生共振相互作用,在较高的L-shell并且低纬地区,能够与较低能量的电子发生共振相互作用.共振谐数越高,能发生波粒共振的电子能量越高.     

Latest progress on interactions between VLF/ELF waves and energetic electrons in the inner magnetosphere

Interactions between very/extremely low frequency (VLF/ELF) waves and energetic electrons play a fundamental role in dynamics occurring in the inner magnetosphere. Here, we briefly discuss global properties of VLF/ELF waves, along with the variability of the electron radiation belts associated with wave-particle interactions and radial diffusion. We provide cases of electron loss and acceleration as a result of wave-particle interactions primarily due to such waves, and particularly some preliminary results...     

太阳活动低年电离层磁场VLF波的观测特性研究

本文利用搭载在DEMETER卫星上的感应式磁力仪(Instrument Magnetometer Search-Coil, IMSC)探测数据分析了磁场甚低频(Very Low Frequency, VLF)波功率谱的空间分布.在排除地磁扰动影响(Dst≤-30 nT,Kp≥3,AE≥200 nT)的前提下,我们给出2007年不同季节,日侧和夜侧,磁场VLF波功率谱的全球分布的背景场探测和对应的统计误差分布.利用太阳活动低年的背景场探测统计分析了磁场VLF波功率谱的频谱特性和季节变化,以及对来自空间的地磁活动的响应特征.太阳同步高度卫星探测的磁场VLF波功率谱也可以观测到源于地球特定区域里的磁场异常所产生的频谱响应.12.5~17.5 kHz频段的磁场VLF波功率谱全球分布清楚地给出了南大西洋磁场异常区的分布轮廓.夜间百慕大地磁异常影响了磁场VLF波功率谱的多个频段,并显示出类似于高纬极区磁场的地磁活动响应的观测特征.卫星磁场VLF波功率谱探测背景研究表明,源于空间的地磁扰动和源于地球的局部地磁异常都会引起电离层磁场VLF波功率谱的变化.特别是在12.5~17.5 kHz频段,我们发现了沿板块边缘分布的电离层磁场VLF波功率谱增强带,这将为我们探寻由板块运动引起的地球突发事件的电离层响应增加探测依据,为定量甄别此类事件的电离层异常信号积累实用的背景参照.     

Airborne VLF Measurements and Variations of Ground Conductivity: A Tutorial

Airborne very low frequency (VLF) data are routinely collected by national agencies and commercial companies together with other passive geophysical measurements of the static magnetic field and radiometric data. The purpose of this paper is to demonstrate that both standard three-component VLF and tensor VLF (TVLF) data contain a lot of useful quantitative and qualitative information about the electrical conductivity distribution in the upper few hundred meters of the crystalline basement. We first give a new derivation of the fundamental transfer functions (the tipper) used in the TVLF technique. We then show that the tipper can be estimated from simultaneous measurements of the wave magnetic fields from at least two transmitters with somewhat different frequencies, and present a simple model by which the maximum error introduced by the difference in frequencies can be found. Single transmitter scalar VLF maps emphasise those conductive structures that have dominant strikes in the direction of the transmitter. Multiple transmitter transfer functions are dependent only upon the underlying conductivity structure. Two dimensional structures can be quantitatively modelled by modern inversion methods developed originally for deep electromagnetic magnetotelluric (MT) soundings. In such cases three-component VLF measurements can be modelled easily upon appropriate rotation of the co-ordinate system to “strike” co-ordinates. Single frequency transfer functions (tippers) have real and imaginary parts that carry information on not only lateral contrasts in conductivity, as usually stated in text books, but, taken together, they provide a robust tool for determining the background conductivity level away from distinct conductors, and they can also be used to discriminate between deep and shallow conductors. Based upon simulations using multi-frequency data, it can be concluded that such a new development would dramatically increase the resolving power of airborne VLF measurements.     

SNR changes of VLF radio signals detected onboard the DEMETER satellite and their possible relationship to the Wenchuan earthquake

Here we used the VLF signal data received by the DEMETER satellite, transmitted from various ground VLF transmitters which are located around China, to study the changes in the signal to noise ratio (SNR) before and after the Wenchuan earthquake, which had a magnitude of 8.0. We also found that the SNRs of different frequency signals decreased significantly over the epicenter region before the earthquake, and reverted to their original levels after the earthquake. This phenomenon may be related to the earthquake. Special Foundation of Basic Research from Institute of Geophysics, China Earthquake Administration (Grant No. DQJB08Z08), National Key Technology R & D Program (Grant No. 2008BAC35B01)     

Unstructured grid based 2-D inversion of VLF data for models including topography

We present a 2-D inversion code incorporating a damped least-squares and a minimum-model approach for plane wave electromagnetic (EM) methods using an adaptive unstructured grid finite element forward operator. Unstructured triangular grids permit efficient discretization of arbitrary 2-D model geometries and, hence, allow for modeling arbitrary topography. The inversion model is parameterized on a coarse parameter grid which constitutes a subset of the forward modeling grid. The mapping from parameter to forward modeling grid is obtained by adaptive mesh refinement. Sensitivities are determined by solving a modified sensitivity equation system arising from the derivative of the finite element equations with respect to the model parameters. Firstly, we demonstrate that surface topography may induce significant effects on the EM response and in the inversion result, and that it cannot be ignored when the scale length of topographic variations is in the order of magnitude of the skin depth. Secondly, the dependency of the inversion on the starting model is discussed for VLF and VLF-R data. Thirdly, we demonstrate the inversion of a synthetic data set obtained from a model with topography. Finally, the inversion approach is applied to field data collected in a region with undulating topography.     

ZH-1卫星观测的VLF人工源信号特征分析与全波模拟

中国地震电磁监测试验卫星张衡一号（ZH-1）已于2018年2月2日成功发射,正在开展卫星数据在轨测试,并对卫星数据质量进行判定.本文对ZH-1卫星2018年5月至6月夜侧的VLF频段电场功率谱数据进行了分析.通过分析位于不同L值、具有不同发射频率的多个VLF人工源上空的卫星重访轨道观测数据,发现ZH-1卫星记录的人工源信号电场变化标准差与DEMETER卫星记录电场变化标准差几乎一致,说明ZH-1卫星观测数据具有较好的稳定性.通过重访轨道均值与全波模型计算结果对比,发现两者在数值上较为接近,在形态上较为一致,说明ZH-1卫星VLF频段电场功率谱数据具有一定的可靠性.此外,研究了VLF人工源上空及共轭区的电场分布特征和电波传播规律,并与DEMETER卫星的结果进行了对比,结果表明VLF人工源产生的电磁辐射穿透电离层后以导管或者非导管的哨声波模向共轭区传播,因为传播过程中的朗道阻尼,共轭区的电场能量比辐射源顶空更小.VLF人工源位于L < 1.5时,电磁波传播更容易发生非导管传播,VLF人工源信号导管传播模式在共轭区的电场响应相对于共轭点会发生一定程度北向偏移.

     

On different numerical inverse Laplace methods for solute transport problems

Numerical inversion is required when Laplace transform cannot be inverted analytically by manipulating tabled formulas of special cases. However, the numerical inverse Laplace transform is generally an ill-posed problem, and there is no universal method which works well for all problems. In this study, we selected seven commonly used numerical inverse Laplace transform methods to evaluate their performance for dealing with solute transport in the subsurface under uniform or radial flow condition. Such seven methods included the Stehfest, the de Hoog, the Honig–Hirdes, the Talbot, the Weeks, the Simon and the Zakian methods. We specifically investigated the optimal free parameters of each method, including the number of terms used in the summation and the numerical tolerance. This study revealed that some commonly recommended values of the free parameters in previous studies did not work very well, especially for the advection-dominated problems. Instead, we recommended new values of the free parameters for some methods after testing their robustness. For the radial dispersion, the de Hoog, the Talbot, and the Simon methods worked very well, regardless of the dispersion-dominated or advection-dominated situations. The Weeks method can be used to solve the dispersion-dominated problems, but not the advection-dominated problems. The Stehfest, the Honig–Hirdes, and the Zakian methods were recommended for the dispersion-dominated problems. The Zakian method was efficient, while the de Hoog method was time-consuming under radial flow condition. Under the uniform flow condition, all the methods could present somewhat similar results when the free parameters were given proper values for dispersion-dominated problems; while only the Simon method, the Weeks method, and the de Hoog method worked well for advection-dominated problems.     

Probing the lower ionospheric perturbations associated with earthquakes by means of subionospheric VLF/LF propagation

There have been reports for many years that the ionosphere is very sensitive to seismic effects, and the detection of ionospheric perturbations associated with earthquakes (EQs) attracts a lot of attention as a very promising candidate for short-term EQ prediction. In this review we present a possible use of VLF/LF (very low frequency (3–30 kHz)/low frequency (30–300 kHz)) radio sounding of seismo-ionospheric perturbations. In order to avoid the overlapping with my own previous reviews, we first show some pioneering results for the Kobe EQ and we try to present the latest results including the statistical evidence on the correlation between the VLF/LF propagation anomalies (ionospheric perturbations) and EQs (especially with large magnitude and with shallow depth), medium-distance (6-8 Mm) propagation anomalies, the fluctuation spectra of subionospheric VLF/LF data (the effect of atmospheric gravity waves, the effect of Earth's tides, etc.), and the mechanism of lithosphere-atmosphere-ionosphere coupling. Finally, we indicate the present situation of this kind of VLF/LF activities going on in different parts of the globe and we suggest the importance of international collaboration in this seismo-electromagnetic study.     

State of the hydrothermal activity of Soufrière of Guadeloupe volcano inferred by VLF surveys

La Soufrière (1467 m) is the active island arc volcano of Guadeloupe Island in the Lesser Antilles arc. Its historical eruptions are more or less violent phreatic outbursts the last of which, in 1976–1977, led to the evacuation of nearly 70 000 persons. The subsurface structure of the volcano consists of calderas, craters, and avalanche amphitheatres nested within the composite pile of eruptive products. Since the last magmatic eruption, dated ca. 1440 AD, the four phreatic eruptions have developed radial fractures on Soufrière dome favouring the development of a huge active hydrothermal system emphasized by a tropical environment. After the eruptions, the thermal state and the stable ground water flow are completely disorganised during several years during which the slow mineralization of rocks is becoming again preponderant. Sealing of fractures and decay of rocks permeability act as a cap for upward thermal transfers. Therefore Soufrière dome operates as a valve, resealing the hydrothermal system underlying the volcano thus providing over pressurization that could lead to the next phreatic eruption. In 1992 new small seismic swarms have appeared. Several of them are recorded every year while the emission of acid gas slowly increases.In order to recognise the superficial electrical resistive and conductive zones (less than 100 m depth) as well as the cavities on Soufrière volcano, we have made Very Low Frequency (VLF) surveys in 2000. Electrical conductive zones are clearly associated with major radial faults starting from the summit in which the hydrothermal activity takes place. In the continuation of these active hydrothermal fractures hot springs are located down slope. Conversely some of the resistive zones are associated with inactive clayed and sealed or opened faults. The distribution of the conductive zones allows detailing the state of the superficial part of the hydrothermal system of La Soufrière. The distribution of vertical clayed zones associated with major faults supposes Soufrière dome constituted of more or less consolidated blocks joined side by side and lying on the hydrothermally floor of crater Amic.