Optical tomography of the aurora and EISCAT

Tomographic reconstruction of the three-dimensional auroral are emission is used to obtain vertical and horizontal distributions of the optical auroral emission. Under the given experimental conditions with a very limited angular range and a small number of observers, algebraic reconstruction methods generally yield better results than transform techniques. Different algebraic reconstruction methods are tested with an auroral are model and the best results are obtained with an iterative least-square method adapted from emission-computed tomography. The observation geometry used during a campaign in Norway in 1995 is tested with the are model and root-mean-square errors, to be expected under the given geometrical conditions, are calculated. Although optimum geometry was not used, root-mean-square errors of less than 2% for the images and of the order of 30% for the distribution could be obtained. The method is applied to images from real observations. The correspondence of original pictures and projections of the reconstructed volume is discussed, and emission profiles along magnetic field lines through the three-dimensionally reconstructed arc are calibrated into electron density profiles with additional EISCAT measurements. Including a background profile and the temporal changes of the electron density due to recombination, good agreement can be obtained between measured profiles and the time-sequence of calculated profiles. These profiles are used to estimate the conductivity distribution in the vicinity of the EISCAT site. While the radar can only probe the ionosphere along the radar beam, the three-dimensional tomography enables conductivity estimates in a large area around the radar site.Former address: MPE Garching     

Depth determination from a non-stationary magnetic profile for scaling geology

The conventional spectral analysis method for interpretation of magnetic data assumes stationary spatial series and a white‐noise source distribution. However, long magnetic profiles may not be stationary in nature and source distributions are not white. Long non‐stationary magnetic profiles can be divided into stationary subprofiles following Wiener filter theory. A least‐squares inverse method is used to calculate the scaling exponents and depth values of magnetic interfaces from the power spectrum. The applicability of this approach is demonstrated on non‐stationary synthetic and field magnetic data collected along the Nagaur–Jhalawar transect, western India. The stationarity of the whole profile and the subprofiles of the synthetic and field data is tested. The variation of the mean and standard deviations of the subprofiles is significantly reduced compared with the whole profile. The depth values found from the synthetic model are in close agreement with the assumed depth values, whereas for the field data these are in close agreement with estimates from seismic, magnetotelluric and gravity data.     

High angle prestack depth migration with absorption compensation

The absorption effect of actual subsurface media can weaken wavefield energy, decrease the dominating frequency, and further lead to reduced resolution. In migration, some actions can be taken to compensate for the absorption effect and enhance the resolution. In this paper, we derive a one-way wave equation with an attenuation term based on the timespace domain high angle one-way wave equation. A complicated geological model is then designed and synthetic shot gathers are simulated with acoustic wave equations without and with an absorbing term. The derived one-way wave equation is applied to the migration of the synthetic gathers without and with attenuation compensation for the simulated shot gathers. Three migration profiles are obtained. The first and second profiles are from the shot gathers without and with attenuation using the migration method without compensation, the third one is from the shot gathers with attenuation using the migration method with compensation. The first and third profiles are almost the same, and the second profile is different from the others below the absorptive layers. The amplitudes of the interfaces below the absorptive layers are weak because of their absorption. This method is also applied to field data. It is concluded from the migration examples that the migration method discussed in this paper is feasible.     

PSEUDOREFLECTION PROFILING METHOD: AN EFFICIENT COMPLEMENT TO CDP METHOD1

A method for pseudoreflection profiling (PRP) is presented. By using transmission seismograms, this method produces time sections similar to those obtained using the common depth point (CDP) method. Pseudoreflection profiles are obtained using the following procedure. Equivalent transmission responses at sites on soft deposits are derived by Wiener filter estimation using seismograms acquired on the sites and on the basement rock. These responses are substituted into the Claerbout-Kunetz equation to obtain the pseudoreflection profiles. A 3D subsurface model is produced when the PRP method is applied to explosion seismograms. The time sections thus obtained are consistent with those obtained using the CDP method in a neighbouring area. The maximum station interval for obtaining coherent profiles with this method is estimated to be one quarter of the wavelength of the transmitted waves. This technique should be used in combination with another method such as the CDP method because it has two weak points: the velocity structure cannot be obtained directly and undesirable waves cannot be eliminated through processing.     

GPS无线电掩星后向传播方法

在大气多路径传播条件下,讨论了两种无线电信号的反演方法:几何光学方法和后向传播方法.当大气存在多路径效应时,采用几何光学方法反演大气参数会引入较大误差.后向传播方法将无线电信号从多路径区域反推至单路径区域,减少了大气多路径效应的影响.为了比较不同的反演方法,利用多相位屏模型,数值模拟了大气多路径条件下无线电信号在大气中...     

Magnetic mapping for the detection and characterization of UXO: Use of multi-sensor fluxgate 3-axis magnetometers and methods of interpretation

Fluxgate vector magnetometers can be powerful instruments in magnetic mapping and in the detection of unexploded ordnance (UXO). Being lightweight and having low power requirements, they allow the development of multi-sensor use. The main problems in using fluxgate magnetometers arise from calibration errors and drift but these can be overcome using a quick and simple method of calibration in the field. This method also has the advantage of compensating permanent and induced magnetic fields generated by magnetized objects carried with the sensors. Measurement accuracy is similar to that obtained with scalar magnetometers. Multi-magnetometer profiling allows direct inversion of raw magnetic data along profiles in order to locate and characterize dipoles typically generated by UXO. The method used is a non-linear inverse procedure that estimates the three coordinates and magnetization vector of the dipole. Application of the technique to both a synthetic case and a field example illustrates the advantages of the method compared with the use of the analytic signal.     

高分辨折射地震资料处理方法及其应用

高分辨折射探测剖面近几年在我国的一些重要地学基础研究项目，特别是在中国地震局“十五” 期间开展的城市活断层探测及复杂上部地壳精细结构研究中得到了越来越广泛的应用.本文对在用传统有限差分方法处理高分辨折射地震资料时的缺陷及其改进算法的研究进展进行了较全面的讨论，对一些缺陷我们已给出了解决的方案，并已成功地应用于实际人工地震高分辨折射测量资料的分析处理，取得了较好的效果，还给出了将改进方法应用于深反射剖面初至波处理解释的一个应用实例，对于另一些缺陷我们在对国内外该领域研究现状调研的基础上，给出了目前我们正在进行中的解决方案.     

Reconstruction of the global lightning distribution from ELF electromagnetic background signals

The global lightning distribution has been obtained for the first time as the result of an inverse problem for the natural extremely low-frequency (ELF) electromagnetic background signals measured simultaneously at three stations: Lekhta (Karelia, Russia), Moshiri (Hokkaido, Japan), and West Greenwich (Rhode Island, USA). We employed a two-step technique consisting of: (1) inversion of ELF field spectra to a source distance distribution (distance profile) for each station of the network and (2) treating the obtained distance profiles as projections in the tomographic method. The structure of the source distance profiles reconstructed from experimental field spectra, varies during a day in a manner specific for each station, which can be explained generally by the temporal distribution of the sources in the tropical continental areas, so-called world thunderstorm centres (WTCs). The comparison of the obtained source distributions with the satellite data on lightning collected by Optical Transient Detector (OTD) and Lightning Imaging Sensor (LIS) shows satisfactory agreement, which proves that our newly developed inversion technique is a suitable method for mapping global lightning.     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。     

An objective algorithm for reconstructing the three-dimensional ocean temperature field based on Argo profiles and SST data

While global oceanic surface information with large-scale, real-time, high-resolution data is collected by satellite remote sensing instrumentation, three-dimensional (3D) observations are usually obtained from in situ measurements, but with minimal coverage and spatial resolution. To meet the needs of 3D ocean investigations, we have developed a new algorithm to reconstruct the 3D ocean temperature field based on the Array for Real-time Geostrophic Oceanography (Argo) profiles and sea surface temperature (SST) data. The Argo temperature profiles are first optimally fitted to generate a series of temperature functions of depth, with the vertical temperature structure represented continuously. By calculating the derivatives of the fitted functions, the calculation of the vertical temperature gradient of the Argo profiles at an arbitrary depth is accomplished. A gridded 3D temperature gradient field is then found by applying inverse distance weighting interpolation in the horizontal direction. Combined with the processed SST, the 3D temperature field reconstruction is realized below the surface using the gridded temperature gradient. Finally, to confirm the effectiveness of the algorithm, an experiment in the Pacific Ocean south of Japan is conducted, for which a 3D temperature field is generated. Compared with other similar gridded products, the reconstructed 3D temperature field derived by the proposed algorithm achieves satisfactory accuracy, with correlation coefficients of 0.99 obtained, including a higher spatial resolution (0.25° × 0.25°), resulting in the capture of smaller-scale characteristics. Finally, both the accuracy and the superiority of the algorithm are validated.     

Ionospheric composition measurement by EISCAT using a global fit procedure

The aim of this paper is to show the gain of the global fit procedure compared to the classical gate-by-gate analysis in the simultaneous determination of the temperatures and composition profiles in the lower F region. A method for the evaluation of the transition altitude and the width of the transition region between molecular and atomic ions is described in the case of the EISCAT CP1 experiment; it is validated using a simulated experiment. It is shown that alternating-code gates have to be processed together with the long-pulse gates. Our new method is applied on six days of actual data. The results obtained during quiet periods are compared with independent estimations.     

Evaluation of shear wave velocity profile using SPT based uphole method

The uphole method is a field seismic test which uses receivers on the ground surface and an underground source. A modified form of the uphole method is introduced in order to obtain efficiently the shear wave velocity (V_S) profile of a site. This method is called the standard penetration test (SPT)-uphole method because it uses the impact energy of the split spoon sampler in the SPT test as a source. Since the SPT-uphole method can be performed simultaneously with the SPT test it is economical and not labor intensive compared to the original uphole methods which use small explosives or a mechanical source. Field testing and interpretation procedures for the proposed method are described. To obtain reliable travel time information of the shear wave, the first peak point of the shear wave using two component geophones is recommended. Through a numerical study using the finite element method (FEM), the procedure of the proposed method was verified. Finally, the SPT-uphole method was performed at several sites, and the field applicability of the proposed method was verified by comparing the V_S profiles determined by the SPT-uphole method with the profiles determined by the downhole, the spectral analysis of surfaces waves (SASW) method and from the SPT-N values.     

UNE METHODE D'ANALYSE FINE DES PROFILS AEROMAGNETIQUES *

A method to analyse aeromagnetic profiles is proposed: several horizontal (width-) and vertical (amplitude-) parameters are deduced from the anomalies. These parameters are then compared with diagrams drawn on a logarithmic scale. The comparison yields depth and magnetization of the causative bodies. The method lends itself equally well to a quick estimate and to a thorough interpretation, depending on the number of parameters deduced and on the number of diagrams. It can be applied also to the anomalies of the total field and its vertical gradient.     

Comparative study of electron density from incoherent scatter measurements at Arecibo with the IRI-95 model during solar maximum

Arecibo (18.4 N, 66.7 W) incoherent scatter (IS) observations of electron density N(h) are compared with the International Reference Ionosphere (IRI-95) during midday (10/14 h), for summer, winter and equinox, at solar maximum (1981). The N(h) profiles below the F2 peak, are normalized to the peak density NmF2 of the F region and are then compared with the IRI-95 model using both the standard B0 (old option) and the Gulyaeva-B0 thickness (new option). The thickness parameter B0 is obtained from the observed electron density profiles and compared with those obtained from the IRI-95 using both the options. Our studies indicate that during summer and equinox, in general, the values of electron densities at all the heights given by the IRI model (new option), are generally larger than those obtained from IS measurements. However, during winter, the agreement between the IRI and the observed values is reasonably good in the bottom part of the F2 layer but IRI underestimates electron density at F1 layer heights. The IRI profiles obtained with the old option gives much better results than those generated with the new option. Compared to the observations, the IRI profiles are found to be much thicker using Gulyaeva-B0 option than using standard B0.     

二维介质柱的电磁成像研究

从体等效原理出发得到用介质柱介电参数表示的体等效电流,并得到一个积分方程组求解散射问题．利用变分原理,对此方程组求变分,并化为矩阵方程组,消去含介质柱内总场变化的项,得到散射场变化与介质柱介电参数变化之间关系的矩阵方程．即反演方程．用算例检验此法的反演精度、收敛性能、多目标反演能力、抗噪声性能及利用不完备测量散射信息反演目标的能力．结果表明,此法是一种实用性很强的二维介质柱电磁成像方法．     

CSAMT和重力方法在狮子湖温泉深部地球物理勘查中的应用

作为一种新型开发的绿色能源,地热资源被越来越多的人所重视.开采埋藏较深的地热资源风险大,因此开发前的地质和地球物理勘查是十分必要的.采用单一的地球物理方法勘探具有很大的风险,因此采用多种方法进行综合调查可以降低单一方法的风险,取得较好的效果.本文以狮子湖温泉为例,研究了地球物理方法在温泉勘探的应用.本次勘探采用可控源音频大地电磁法(CSAMT-Controlled Source Audio-frequency Magnetotellurics)和微重力测量,CSAMT是针对大地电磁测深法场源的随机性和信号微弱,提出的一种采用可以控制的人工场源改进方案.仪器采用美国Zonge公司生产的GDP-32Ⅱ.该方法由人工向地下供入音频谐变电流建立电磁场,通过仪器在地面接收从地下反馈来的信息,根据不同时代、岩性地层电性特征达到勘查目的.为此我们在测区做了两条剖面,从CSAMT反演图推断,自上而下可分为3个电阻率层,该剖面视电阻率具有很好的层理特征,反映了厚大的中新生代地层（Edn、K）覆盖.剖面西部有明显的泥盆系地层（D）存在.根据视电阻率的变化特点,可以推断这条剖面的5条断层.在剖面中部距地表400～800米深处存在一明显低阻区,推测应为含水破碎带或低阻泥岩.对数据进行二维反演,可以清楚的看到利于储水的盆地构造.CSAMT方法受静态效应影响很大.静态效应位移可能是由地形和电阻率的浅部的横向变化引起的,既是不可避免的,也是不可预测的.因此对数据进行预处理是十分必要的.减少静态效应的影响办法有以下三种：（1）、对效应进行理论计算;（2）、采用空间滤波和相位积分等处理方法;（3）、使用独立的、无静态效应的测量方法.计算静态效应理论值在理论上是简捷的,但在实际的野外条件下,由于无法预测引起静态效应的物体的几何尺寸和电性参数,因此这种方法无法得到可靠的校正值.空间滤波处理是目前广泛采用的一类方法,Bostick（1986）提出了消除MT数据中静态效应的电磁列阵剖面法（EMAP）.EMAP法由于采用连续的剖面测量,可采用窗口可变的自适应空间滤波器-汉宁窗（Hanning window）或叫余弦钟形滤波器消除静态效应.但是这种方法提供了静态效应的要求的数据密度,这就增大了大量的额外测量,提高了获得数据的代价.因此,我们做重力勘探与CSAMT相对比.微重力测量使用美国产LCR-D型重力仪.实测的微重力异常是地下由浅到深各类地质体的物性差异在地面综合叠加的效应,其中包括界面起伏、岩性不均匀、地壳与壳下物质的厚度变化等诸多地质因素在内.实测的重力异常值经过固体潮改正、零点漂移改正、布格改正、正常场改正之后,得到改正后的重力异常值.从微重力反演结果看来,自西往东重力异常逐渐减小,程台阶下降趋势,并趋于平缓,验证了CSAMT的异常结果.从而弥补了CSAMT法的不足.通过这两种方法,我们大致查明了新生代红层盆地的产出形态.该红层盆地西侧边缘位于青山小学一带,自西往东变深;并反演计算出新生代红层与晚古生代泥盆系地层的分界面.深部地球物理勘探的方法有很多,各有各自的优缺点,我们不能从单一的一种方法得到的结论来判断地下地质构造.CSAMT法勘探深度大,但由于本身的物理特性,导致静态效应、近场效应等影响甚大.严重影响我们对地下地质目标体深度的判断,重力和CSAMT法的相互验证,很好的说明了多种地球物理方法综合测定的优势.     

Doppler-shifted auroral H <beta> emission: a comparison between observations and calculations

Two sounding rockets equipped with photometers and particle detectors have been flown into proton auroras. The measured altitude dependence of the proton flux is compared with calculations based upon known energy-range relations for protons in air. Expressions suitable for numerical calculations of Doppler profiles at arbitrary angles to the geomagnetic field and at different heights within an aurora are developed. Profiles due to some typical proton spectra have been calculated and it is shown that altitude profiles at some wavelengths are more sensitive to the shape of the proton spectrum than are profiles at other wavelengths. Variations in the H <beta> Doppler profile versus height for several angles with the magnetic field is studied. Profiles, as generated by the actually measured protons in the energy range 1 keV to 1 MeV, have been calculated and are compared with direct optical measurements made by ground and rocket photometers. The rocket photometers took measurements at different wavelengths within the Doppler profile. The correspondence between calculations and measurements is generally good. The total H <beta> is calculated and fair agreement with the measured intensity is found.     

Inversion of SASW dispersion curves based on maximum flexibility coefficients in the wave number domain

Spectral analysis of surface waves (SASW) is a nondestructive in-situ testing method that is used to determine stiffness profiles of soil and pavement sites based on dispersion characteristics of Rayleigh-type surface waves.Inversion of the Rayleigh wave dispersion curve of a site provides information on the variation of shear-wave velocity with depth. In the inversion procedures currently used for SASW tests, the field dispersion curve is matched with a theoretical dispersion curve obtained for the fundamental mode of surface wave propagation.In order to overcome difficulties associated with the presence of multi-modes in SASW signals, a new inversion method based on the maximum vertical flexibility coefficient is introduced in this paper. Unlike root-searching methods, the new method easily identifies the predominant propagation modes. In this new approach, the simplex method is used to match field and theoretical dispersion curves automatically. The purpose of this paper is to present the details of the new method and to demonstrate its advantages.     

Accuracy Estimation of the Modern Core Magnetic Field Models Using DMA-Methods for Recognition of the Decreased Geomagnetic Activity in Magnetic Observatory Data

A new approach to identify the signals of the Earth’s main magnetic field (core field) based on the magnetic observatory data processing is suggested. The algorithms implemented in the approach are based on the Discrete mathematical analysis (DMA). The developed method is used to analyze the data from 49 midand low-latitude observatories of the INTERMAGNET network collected in 2011–2015. The results are compared with the classical method for determining the periods of low magnetic activity of external origin which is adopted by the International Association for Geomagnetism and Aeronomy (IAGA). The advantages of the suggested new approach are demonstrated. Based on the data records for the selected time intervals, the time series of the core field components and their secular variations are obtained for each observatory. These data are compared to the values predicted by the most accurate core field models: SIFM, CHAOS-6, and EMM-2015. The accuracy of the models is estimated using a set of statistical parameters: Pearson’s coefficient of linear correlation, Spearman’s and Kendall’s coefficients of rank correlation, the mean and median values over the data sets, and the mean difference between the data obtained by the suggested method from observatory measurements and the model predictions.     

Statistical equilibrium distributions of baroclinic vortices in a rotating two-layer model at low Froude numbers

The point-vortex equilibrium statistical model of two-layer baroclinic quasigeostrophic vortices in an unbounded f-plane is examined. A key conserved quantity, angular momentum, serves to confine the vortices to a compact domain, thereby justifying the statistical mechanics model, and also eliminating the need for boundary conditions in a practical method for its resolution. The Metropolis method provides a fast and efficient algorithm for solving the mean field non-linear elliptic PDEs of the equilibrium statistical theory. A verification of the method is done by comparison with the exact Gaussian solution at the no interaction limit of zero inverse temperature. The numerical results include a geophysically and computationally relevant power law for the radii at which the most probable vortex distribution is non-vanishing: For fixed total circulation, and fixed average angular momentum, the radii of both layers are proportional to the square root of the inverse temperature β. By changing the chemical potentials μ of the runs, one is able to model the most probable vorticity distributions for a wide range of total circulation and energy. The most probable vorticity distribution obtained at low positive temperatures are consistently close to a radially symmetric flat-top profiles. At high temperatures, the radially symmetric vorticity profiles are close to the Gaussian distribution.