Anomalous phenomena on HF radio paths during geomagnetic disturbances

We analyze ionospheric oblique sounding data on three high-latitude and one high-latitude–midlatitude HF radio paths for February 15 and 16, 2014, when two substorms and one magnetic storm occurred. We investigate cases of anomalous propagation of signals: their reflection from sporadic layer Es, lateral reflections, type “M” or “N” modes, the presence of traveling ionospheric disturbances, and the diffusivity of signals and triplets. The most significant results are the following. In geomagnetically undisturbed times, sporadic Es-layers with reduced maximum observed frequencies (MOFEs) on three high-latitude paths were observed in both days. The values of MOFEs during disturbances are large, which leads to the screening of other oblique sounding signals reflected from the ionosphere. On all four paths, the most frequently traveling ionospheric disturbances due to the terminator were observed in quiet hours from 03:00 to 15:00 UT on the first day and from 06:00 to 13:00 UT on the second day of the experiment. In addition, both the sunset terminator and the magnetic storm on the high-latitude–mid-latitude path were found to generate traveling ionospheric disturbances jointly. No such phenomenon was found on high-latitude paths.     

Numerical analysis of one-dimensional nonlinear acoustic wave

Numerical investigations on one-dimensional nonlinear acoustic wave with third and fourth order nonlinearities are presented using high-order finite-difference (HFD) operators with a simple flux-limiter (SFL) algorithm. As shown by our numerical tests, the HFDSFL method is able to produce more stable, accurate and conservative solutions to the nonlinear acoustic waves than those computed by finite-difference combined with the flux-corrected-transport algorithm. Unlike the linear acoustic waves, the nonlinear acoustic waves have variable phase velocity and waveform both in time-space (t-x) domain and frequency-wavenumber (f-k) domain; of our special interest is the behaviour during the propagation of nonlinear acoustic waves: the waveforms are strongly linked to the type of medium nonlinearities, generation of harmonics, frequency and wavenumber peak shifts. In seismic sense, these characteristics of nonlinear wave will introduce new issues during such seismic processing as Normal Moveout and f-k filter. Moreover, as shown by our numerical experiment for a four-layer model, the nonlinearities of media will introduce extra velocity errors in seismic velocity inversion.     

Device for the physical modeling of nonlinear phenomena in seismic exploration

The instrument required for investigating the acoustic nonlinearity of artificial and natural rock specimens should combine high sensitivity and a high noise immunity of measurements. The methods of the induction and measurement of the acoustic field engaged in the device under consideration make it possible to investigate weak nonlinear phenomena. The device was tested on a number of models.     

Mapping of tidal current and associated nonlinear currents in the Xiangshan Bay by coastal acoustic tomography

Ocean Dynamics - Herein, tidal current and associated nonlinear currents were measured using nine coastal acoustic tomography (CAT) systems surrounding the tomography domain in Xiangshan Bay from...     

Constraints imposed on crustal heat production by the linear heat flow relation

The class of three-dimensional radiogenic heat source functions, which satisfy a linear relation between heat flow and surface heat production, is constructed by expanding the source in terms of Fourier transforms in the horizontal directions and Laguerre functions in the vertical. The linear relation then becomes a set of linear equations among the expansion coefficients. In particular, it is shown that the preferred exponential decay model can only hold if heat flow is constant at the Earth's surface. Furthermore, if a heat source function has any lateral variation, then the linear relation cannot survive erosion of any uniformly thick slab. For a given heat flow data set, construction of a source solution, closest in the mean-square sense to the exponential decay model, is demonstrated. When applied to a heat flow data set from central England, the construction gives vertical profiles which locally show significant departures from the exponential model.     

Laboratory and field investigations of rockburst phenomena using concurrent geotomographic imaging and acoustic emission/microseismic techniques

The progress made on three phases of a research project, started in 1986 to investigate mining induced seismicity/rockburst phenomena using concurrent geotomographic imaging and microseismic monitoring techniques, is described. Phase I is the geotomographic software development and laboratory calibration trials. Phase II is the enhancement of traditional microseismic monitoring instrumentation with a waveform acquisition system, so that source mechanism studies can be carried out on mining induced seismic events. Phase III is the field trials of the hybrid technique which will be used to monitor changing rock mass physical properties, in response to mining. Preliminary results from all three phases are given, together with an outline of current and future research planned.Presented at the Fred Leighton Memorial Workshop on Mining Induced Seismicity, Montreal 1987.     

平面海底界面上球体目标的声散射建模研究

根据目标声散射的T矩阵方法、将界面附近点源激发的声场利用复像方法展开、并考虑到界面对散射声场的影响,建立了平面海底界面上方目标的声散射模型.与已有模型相比,该模型能够计算更一般情况点源激发的三维声散射场.通过计算考察基于复像方法的有关系数证明了该模型的准确性,并给出了不同入射方向的界面上刚性球体和弹性球壳的散射声场算例,对目标散射幅度的分布以及随角度的变化规律做了分析.最后,设计实验对该模型做了验证.     

Doppler selection of HF radiosignals on long paths

The long-term registration of the Doppler spectra of HF radiosignals has been performed on the Moscow-Akademik Vernadsky Ukrainian Antarctic station path. It has been revealed that the spectra are split when the solar terminator crosses direct and return radio lines. The spectral and energy characteristics of direct and return signals have been calculated within the scope of the asymptotic theory of long-range propagation of decametric radiowaves.     

Effects of pore fluids on quasi-static shear modulus caused by pore-scale interfacial phenomena

It is evident from the laboratory experiments that shear moduli of different porous isotropic rocks may show softening behaviour upon saturation. The shear softening means that the shear modulus of dry samples is higher than of saturated samples. Shear softening was observed both at low (seismic) frequencies and high (ultrasonic) frequencies. Shear softening is stronger at seismic frequencies than at ultrasonic frequencies, where the softening is compensated by hardening due to unrelaxed squirt flow. It contradicts to Gassmann's theory suggesting that the relaxed shear modulus of isotropic rock should not depend upon fluid saturation, provided that no chemical reaction between the solid frame and the pore fluid. Several researchers demonstrated that the shear softening effect is reversible during re-saturation of rock samples, suggesting no permanent chemical reaction between the solid frame and the pore fluid. Therefore, it is extremely difficult to explain this fluid–rock interaction mechanism theoretically, because it does not contradict to the assumptions of Gassmann's theory, but contradicts to its conclusions. We argue that the observed shear softening of partially saturated rocks by different pore fluids is related to pore-scale interfacial phenomena effects, typically neglected by the rock physics models. These interface phenomena effects are dependent on surface tension between immiscible fluids, rock wettability, aperture distribution of microcracks, compressibility of microcracks, porosity of microcracks, elastic properties of rock mineral, fluid saturation, effective stress and wave amplitude. Derived equations allow to estimate effects of pore fluids and saturation on the shear modulus and mechanical strength of rocks.     

非对称圆弧片状声源在井旁地层中产生的声场

提出了一种可用于随钻三维反射声波成像测井的圆弧片状声源,推导了该声源声学性质在波数-频率域内的数学描述,并利用实轴积分的方法对该声源在井旁地层中产生的声场进行求解,考察了声源的线度、频率等因素对该声源向地层中辐射的声场的影响.研究结果表明,特定尺寸和特定频率下的圆弧片状声源产生的纵波场的水平和垂直指向性图中仅具有一个明显的主瓣,且主瓣三分贝角宽较窄,方位分辨率较高,指向性良好,适用于随钻三维反射成像测井;SV波场和SH波场的指向性图中存在两个或者多个角瓣,SV波场在垂直于井轴方向上辐射的能量为0,另外在井中发射和接收的反射SH波场和SV波场会相互干扰,所以难以利用SV波场和SH波场进行三维反射声波成像;声源的线度和频率等因素对其辐射声场影响较大,考虑到激励效率和方位分辨率等因素,对于本文描述的井孔模型,选择圆周角在75°和90°之间, 主频为12 kHz左右的声源是合适的.     

Heating of beam ions by ion acoustic waves

Satellite measurements show that ion beams above the auroral acceleration region are heated to hundreds of eV in a direction perpendicular to the magnetic field. We show that ion acoustic waves may be responsible for much of this heating. Even in the absence of a positive slope in the velocity distribution of the beam ions, ion acoustic waves can be generated by a fan instability. We present analytical estimates of the wave growth rate and ion beam heating rate. These estimates, which are confirmed by particle simulations, indicate that the perpendicular temperature of the beam ions will increase by 30 eV/s, or by 1 eV in 20–25 km. From the simulations we also conclude that the heating saturates at a perpendicular temperature around 200 eV, which is consistent with observations.     

Modelling the wave phenomena in acoustic and elastic media with sharp variations of physical properties using the grid‐characteristic method

This paper introduces a novel method of modelling acoustic and elastic wave propagation in inhomogeneous media with sharp variations of physical properties based on the recently developed grid‐characteristic method which considers different types of waves generated in inhomogeneous linear‐elastic media (e.g., longitudinal, transverse, Stoneley, Rayleigh, scattered PP‐, SS‐waves, and converted PS‐ and SP‐waves). In the framework of this method, the problem of solving acoustic or elastic wave equations is reduced to the interpolation of the solutions, determined at earlier time, thus avoiding a direct solution of the large systems of linear equations required by the FD or FE methods. We apply the grid‐characteristic method to compare wave phenomena computed using the acoustic and elastic wave equations in geological medium containing a hydrocarbon reservoir or a fracture zone. The results of this study demonstrate that the developed algorithm can be used as an effective technique for modelling wave phenomena in the models containing hydrocarbon reservoir and/or the fracture zones, which are important targets of seismic exploration.     

南京地磁台地铁干扰特征分析与抑制处理

地铁运行时产生的漏电流而形成的附加磁场会对附近的地磁观测造成严重干扰。通过对地铁供电系统的模拟分析得出了地铁干扰的一些基本特征。选取了2008年5月1日到30日南京地磁台受地铁干扰的日变数据和泉州参考地磁台日变数据,利用Welch功率谱估计法进行功率谱对比分析,发现地铁干扰主要对地磁场的Z分量、D分量的功率谱的高频部分造成影响,并且影响的频带范围几乎一致,主要集中在0．004-0．03Hz范围,而对H分量的影响比较微弱。在此基础上,通过小波变换对受干扰数据的Z分量和D分量进行干扰抑制处理,处理结果表明该方法能够有效抑制地铁造成的地磁场干扰。因此本文研究也提供了地磁台站干扰抑制的一种方法。     

Detection of model voids by identifying reverberation phenomena in GPR records

In order to improve the reliability of the ground penetrating radar (GPR) method in identifying subsurface sinkholes and karst cavities, laboratory investigations have been performed. The main objective of this work was to examine the relationship between horizontal/vertical voids dimensions and wavelengths of various antennas, and the corresponding GPR responses. Emphasis was given to the investigation of the factors that cause the appearance of reverberation phenomena in the signal pattern.The tests were conducted in 5 m × 10 m area by 2-m-deep trench filled with homogenous, dry sand. The voids models (empty fiberglass cylinders in diameters of 0.6 m, 1.0 m, 1.5 m and 2.4 m, with various heights) were buried vertically with their tops at depths of between 0.7 and 1.5 m. Investigations were performed for the various model conditions by towing 500, 300 and 100 MHz antennas along a pre-established grid, for the various model conditions.The GPR data collected using the 500 MHz bistatic antenna above the 1.0-m- and the 1.5-m-diameter cylinders, and using the 300 MHz bistatic antenna above the 1.5-m-diameter cylinder, confirmed the presence of a reverberation phenomenon, i.e. a strong convex signal pattern, containing a series of high amplitude extending oscillations with reduced frequency.Based on past practical GPR experience of void detection and presently obtained experimental data, two rules of thumbs may be adopted for the prediction of the appearance of resonant radar pictures:

1. The void diameter larger than the wavelength in air of the antenna used.

2. The vertical size of the empty void not significantly smaller than its horizontal dimension.

The strong reverberations generated by the inner surface of the void targets were found to approximate standing waves generated in cylindrical waveguides and waveguide resonators. The theoretical, experimental and practical results obtained concur.     

Effects of experimental uncertainty on the calculation of hillslope flow paths

Measurement uncertainty is a key hindrance to the quantification of water fluxes at all scales of investigation. Predictions of soil‐water flux rely on accurate or representative measurements of hydraulic gradients and field‐state hydraulic conductivity. We quantified the potential magnitude of errors associated with the parameters and variables used directly and indirectly within the Darcy – Buckingham soil‐water‐flux equation. These potential errors were applied to a field hydrometric data set collected from a forested hillslope in central Singapore, and their effect on flow pathway predictions was assessed. Potential errors in the hydraulic gradient calculations were small, approximately one order of magnitude less than the absolute magnitude of the hydraulic gradients. However, errors associated with field‐state hydraulic conductivity derivation were very large. Borehole (Guelph permeameter) and core‐based (Talsma ring permeameter) techniques were used to measure field‐saturated hydraulic conductivity. Measurements using these two approaches differed by up to 3\9 orders of magnitude, with the difference becoming increasingly marked within the B horizon. The sensitivity of the shape of the predicted unsaturated hydraulic conductivity curve to ±5% moisture content error on the moisture release curve was also assessed. Applied moisture release curve error resulted in hydraulic conductivity predictions of less than ±0\2 orders of magnitude deviation from the apparent conductivity. The flow pathways derived from the borehole saturated hydraulic conductivity approach suggested a dominant near‐surface flow pathway, whereas pathways calculated from the core‐based measurements indicated vertical percolation to depth. Direct tracer evidence supported the latter flow pathway, although tracer velocities were approximately two orders of magnitude smaller than the Darcy predictions. We conclude that saturated hydraulic conductivity is the critical hillslope hydrological parameter, and there is an urgent need to address the issues regarding its measurement further. Copyright © 2000 John Wiley & Sons, Ltd.     

Limitations on gyration constants for main symmetry systems of rocks

In the case of propagation of plane elastic waves in anisotropic gyrotropic media, Christoffel tensor is complex; its real part contains stiffnesses and an imaginary part includes components of the fifth-rank gyration tensor. Inequalities relating stiffnesses and gyration constants are derived from the conditions for potential energy to be positive. The necessary and sufficient conditions for the positive definiteness of the complex matrix of stiffnesses and gyration constants are used. Sets of inequalities are obtained for two types of rocks belonging to acentric limit groups ∞∞ and ∞. These inequalities provide a possibility to carry out modelling of elastic wave propagation in the media considered, setting the values of gyration constants not arbitrarily but in accordance with physical laws.     

Seismic diagnosis of railway substructures by using secondary acoustic emission

One difficult task for the seismic diagnosis of existing structures is how to nondestructively evaluate the damage degree of invisible substructures, such as embedded foundations. To diagnose substructures efficiently, a method for nondestructive inspection is developed by applying acoustic emission (AE) technique. As a newly proposed method, characteristics of secondary AE induced by train operations were investigated, and experiments using model piles and in-situ AE monitoring of in-service railway bridges conducted under railroad traffic, from which it was demonstrated that the proposed method is practicable enough to detect invisible defects in structures. A new index, known as RTRI (ratio of Repeated Train load at the onset of AE activity to Relative maximum load for Inspection period) is proposed for structural damage qualification based on the results of in-situ AE monitoring.     

Time domain scattering of acoustic plane waves by vertical faults

Acoustic plane wave scattering at a vertical fault structure represents the simplest two-dimensional model of geophysical exploration that can be investigated by analytical techniques. The exact and complete solution, in the time domain, for the scattering of the pressure field of an acoustic plane wave normally incident on a vertical fault structure is determined adapting previous results given for the frequency domain. The wave form of the pressure field of the incident plane wave is expressed by a causal time function that decays exponentially with time at every point above the fault (z<0). The zero-order term of the scattered pressure field has been computed above the fault. This zero-order term consists of an inverse Fourier transform which reduces to a closed expression forx=0, and contains an integral of a Hankel function forx#0. The high frequency part of the inverse Fourier transform forx#0 is computed employing asymptotic expressions for the Hankel function. The integral of the asymptotic expression of the Hankel function reduces to: (i) a Fresnel integral which contains a plane wave term for |x||z|; and (ii) a stationary point plane wave term plus an upper limit term for |x|=O(|z|). For the latter case the plane wave term cancels, leaving a cylindrical wave emanated from the edge of the fault. The wave front is well defined in shape, in phase and in amplitude. The amplitude of the scattered field is discontinuous atx=0, presents a jump and is well defined for |x| small and is rather smooth for |x| large.