垂直分量磁称记录台资料在电磁感应法中的应用

一、引言用电磁感应方法研究局部地球电导率结构时经常引用转换函数的概念。转换函数反映了地磁变化的源场与感应场之间的线性关系。由于地磁短周期变化的源场在中低纬度地区是近似均匀的,而其感应场又是由地下感应电流产生,因此转换函数反映了地下电导率分布的某些特证。所以,分析转换函数的空间分布及随频率的变化,可做为研究局部地球电性结构的     

地磁水平场转换函数的短期前兆

利用昌黎、白家疃1972～1997年的磁变仪资料，计算了水平场台际转换函数和单台垂直场转换函数．通过对转换函数变化全过程的分析，注意到仅仅水平场转换函数在唐山地震之前有明显的短期前兆．这给地震预报提供了一个新的可能捕捉到短期前兆的手段．      

地磁水平场转换函数在唐山地震前的异常反应

搜集了昌黎和白家疃1972～1984年的磁变仪资料．19组资料的计算结果表明:① 白家疃的单台垂直场转换函数A，B及昌黎的B没有明显异常，但昌黎的Au，Av有明显异常变化．这一结果与用量图方法求得的a，b基本一致．表明异常局限在余震区及其附近；② 台际转换函数Cu，Cv，Fu，Fv在唐山地震前都有一些变化，但虚部Cv，Fv的变化较为显著．这点与对台湾花莲地震的研究结果一致．但Beamish(1982)对Carlisle 5级地震的研究结果，是实部Fu，Cu的变化较为显著．      

走滑断层对洋中脊热结构与流动场的影响规律及对太平洋南部边界RRF型三联点的解释

选取太平洋板块南部边界的板块相对运动速度不同的两个洋脊－洋脊－转换断层（RRF）型三联点,即麦夸里（Macquarie）三联点和南太平洋三联点,为研究对象,通过数值模拟的方法,研究该类型三联点走滑断层边界两侧的板块相对运动速度对三联点附近地区地幔流动场和温度结构的影响。模拟结果表明:太平洋南部边界RRF三联点走滑断层边界两侧的板块相对运动速度控制着三联点附近的温度分布和地幔流动;随着走滑断层边界两侧板块相对运动速度的增加,转换断层相对滑动速度增加,温度上升,距洋脊边界100 km范围内的地幔流体速度变大;麦夸里三联点和南太平洋三联点处3个板块的相对运动,使得三联点的转换断层边界浅部产生剪应力集中,导致震源深度集中在15—25 km;同时相对运动产生的地幔流动引起温度结构变化,该变化控制着地形变化。      

三维大地电磁感应参量研究

本文讨论了大地电磁感应转换函数在使用上有争议的问题.用几个不同的三维(3D)模型计算了感应矢量和转换函数,特别对它们的虚部进行了分析.方形和矩形导电体的感应具有二维(2D)特征,其感应曲线可用来估计导体的埋深和范围.由二个埋深不同导体组成的结构,则表现出非二维特性,其感应函数曲线变化复杂,实虚矢量夹角大,两个转换函数的值即使在主轴方向也可能大小相当.      

水平磁场台际转换函数Gu的大小--有限差分法模型计算研究

计算的几种模型表明,用改进的三维有限差分法程序计算的异常体中心点上的水平磁场实部可达1．53,这就意味着水平场转换函数Cu也可达1．53,与单台垂直磁场的转换函数有相同量级的大小。     

中国地磁基本场模式建立方法探讨

利用７０年代的地磁资料，讨论了测点和密度和测点坐标变换对建立中国地磁基本场模式的影响，结果表明，只选用分布较为均匀的２００个左右的测点建立基本场模式是可行的，对测点的球面坐标进行Ｍｅｒｃａｔｏｒ投影变换是必要的，可以显著减小模式的场值分布在边界附近出现的畸变现象；３阶到５阶模式均主要反映正常的基本磁场分布，而５阶模式效果最好。     

基于尺度空间技术的归一化Facet模型位场边界识别

边界识别是位场数据处理解释中的重要环节,传统边界识别方法通常不能均衡深、浅部地质体边界.基于尺度空间技术和归一化的Facet模型检测算子,本文开发了一种带通空间滤波和边缘检测相结合的边界识别方法,有效地提高位场数据边界识别的精度和可靠性.为了验证本文算法的有效性和稳定性,分析了不同尺度空间函数和检测算子对算法的影响,并且对比了传统边界识别方法的效果.理论模拟和实际数据分析表明,利用位场垂向二阶导数进行的基于尺度空间技术的归一化Facet模型边界识别方法不仅算法的稳定性强,而且可以避免高阶导数对噪声干扰放大作用,同时均衡深部和浅部地质体边界,从而可以更精确地识别地质体的形态.     

青藏高原北部地区构造变形特征及与强震关系

通过对1993～1999、1999～2001年青藏高原北部地区GPS水平运动资料的非震负位错模型反演和形变应变场时空演化分析,结合地质构造和有资料积累以来的强震活动,研究块体及边界带的构造变形特征,以及与强震孕育的关系．结果表明：(1)5级以上,特别是6级左右及6级以上地震多发生在区域应变场剪切应变率或面膨胀率高值区、边缘或其附近,尤其是与区域主干断裂构造运动一致的剪切应变率高值区;形变差异显著的块体边界附近;沿块体边界断裂带的高应变积累的闭锁区段及其附近．(2)1999～2001年较1993～1999年大范围的水平运动速率明显减弱,但应变率的减弱幅度不大,可能揭示昆仑山口西8．1大震孕育对北东向的构造应力传递起了某种“阻隔和调制作用”及太平洋板块的西进作用加强．(3)目前应变场分布高值区、块体间差异运动显著的边界带及高应变积累闭锁段有阿尼玛卿断裂中东段与鄂拉山断裂交汇区、祁连山构造带中东段及与海原断裂西段交汇区、日月山一拉脊山断裂与西秦岭北缘断裂交汇区、庄浪河断裂与西秦岭北缘断裂交汇区、六盘山断裂附近．     

平面SV波在饱和半空间中沉积谷地周围的散射

采用一种特殊的间接边界积分方程法,求解了平面SV波在饱和半空间中任意形状沉积谷地周围的二维散射问题。结合饱和半空间中膨胀波源和剪切波源格林函数,由分布在沉积和半空间交界面附近两虚拟波源面上的波源分别构造沉积内外的散射波场,由交界面连续条件建立方程并求解确定虚拟波源密度,总波场反应即可由自由波场和散射波场叠加而得。然后通过边界条件验算、退化解答与现有结果的比较以及稳定性检验,验证了方法的计算精度。通过一组典型算例,研究了平面SV波在饱和半空间中沉积谷地周围散射的基本规律,详细给出了不同参数情况沉积谷地附近地表位移幅值和孔隙水压,着重分析了入射SV波频率和角度、边界渗透条件、沉积孔隙率等因素对场地反应的影响,得出了一些有益的结论。     

Numerical modeling calculation for the spatial distribution characteristics of horizontal field transfer functions

Introduction There are two kinds of methods for the mathematic study on electro-magnetic induction in geo-field: one is analytic method, such as integral equation method; the other is numerical model-ing method, such as finite difference method and finite element method. The analytic method can only be applied to the conductor with very simple shape, such as sphere, circular cylinder, etc. With the increasing of computing speed and popularity of computer, the numerical modeling methods are use…     

A heating process of Kuchi-erabu-jima volcano,Japan, as inferred from geomagnetic field variations and electrical structure

Since August 2000, we have recorded the total intensity of the geomagnetic field at the summit area of Kuchi-erabu-jima volcano, where phreatic eruptions have repeatedly occurred. A time series analysis has shown that the variations in the geomagnetic field since 2001 have a strong relationship to an increase in volcanic activity. These variations indicate thermal demagnetization of the subsurface around the presently active crater. The demagnetization source for the early variations, until summer 2002, was estimated at about 200 m below sea level. For the variations since 2003, the source was modeled on the basis of the expansion of a uniformly magnetized ellipsoid. The modeling result showed that the source is located at 300 m above sea level beneath the crater. We carried out an audio-frequency magnetotelluric survey with the aim of obtaining a relation between the demagnetization source and the shallow structure of the volcano. A two-dimensional inversion applied to the data detected two good conductors, a shallow thin one which is restricted to a region around the summit area, while the other extends over the edifice at depths between 200 and 800 m. These conductors are regarded as clay-rich layers with low permeability, which were assumed to be generated through hydrothermal alteration. The demagnetization source for the early variations was possibly located at the lower part of the deep conductor and the source after 2003 lies between the two conductors, where groundwater is considered to be abundant. Based on these results, as well as on seismological, geodetic, and geochemical information, we propose a heating process of the Kuchi-erabu-jima volcano. In the initial stage, high-temperature volcanic gases supplied from the deep-seated magma remained temporarily at the level around the lower part of the less permeable deep conductor since the ascent path had not yet been established. Then, when the pathway developed as a result of repeated earthquakes, it became possible for a massive flux of volcanic gases to ascend through the conductor. The high temperature gases reached the aquifer located above the conductor and the heat was efficiently transported to the surrounding rocks through the groundwater. As a consequence, an abrupt increase of the gas flux and diffusion of the heat through the aquifer occurred and the high-temperature zone expanded. Since the high-temperature zone is located beneath another conductor, which acts as caprock, we assume that the energy of the phreatic explosion is accumulated there.     

A preparation zone for volcanic explosions beneath Naka-dake crater,Aso volcano,as inferred from magnetotelluric surveys

The 1st crater of Naka-dake, Aso volcano, is one of the most active craters in Japan, and known to have a characteristic cycle of activity that consists of the formation of a crater lake, drying-up of the lake water, and finally a Strombolian-type eruption. Recent observations indicate an increase in eruptive activity including a decrease in the level of the lake water, mud eruptions, and red hot glows on the crater wall. Temporal variations in the geomagnetic field observed around the craters of Naka-dake also indicate that thermal demagnetization of the subsurface rocks has been occurring in shallow subsurface areas around the 1st crater. Volcanic explosions act to release the energy transferred from magma or volcanic fluids. Measurement of the subsurface electrical resistivity is a promising method in investigating the shallow structure of the volcanic edifices, where energy from various sources accumulates, and in investigating the behaviors of magma and volcanic fluids. We carried out audio-frequency magnetotelluric surveys around the craters of Naka-dake in 2004 and 2005 to determine the detailed electrical structure down to a depth of around 1 km. The main objective of this study is to identify the specific subsurface structure that acts to store energy as a preparation zone for volcanic eruption. Two-dimensional inversions were applied to four profiles across the craters, revealing a strongly conductive zone at several hundred meters depth beneath the 1st crater and surrounding area. In contrast, we found no such remarkable conductor at shallow depths beneath the 4th crater, which has been inactive for 70 years, finding instead a relatively resistive body. The distribution of the rotational invariant of the magnetotelluric impedance tensor is consistent with the inversion results. This unusual shallow structure probably reflects the existence of a supply path of high-temperature volcanic gases to the crater bottom. We propose that the upper part of the conductor identified beneath the 1st crater is mainly composed of hydrothermally altered zone that acts both as a cap to upwelling fluids supplied from deep-level magma and as a floor to infiltrating fluid from the crater lake. The relatively resistive body found beneath the 4th crater represents consolidated magma. These results suggest that the shallow conductor beneath the active crater is closely related to a component of the mechanism that controls volcanic activity within Naka-dake.     

A PHYSICAL SCALE MODEL STUDY OF THE COMPARATIVE PERFORMANCE OF TWO MODES OF OPERATION FOR FIXED-LOOP TURAM-TYPE EM SYSTEMS1

The results of a physical scale model study of the conventional mode of operation of fixed loop electromagnetic systems and an alternative mode called the tx-parallel mode in which traverses are run parallel to the long axis of the rectangular transmitter loop are presented. The results show that over thick or dipping conductors, the tx-parallel configuration provides coupling with the target which is comparable with that provided by the conventional configuration. In addition, the tx-parallel configuration is shown to provide more consistent indications of the direction and magnitude of conductor dip. Over wide conductors, where separate conventional surveys are needed to define the opposite edges of such conductors, it is shown that only a single tx-parallel survey is needed to locate both edges of the conductor. The tx-parallel results were found to allow better resolution of the individual anomalies due to closely spaced parallel conductors. The tx-parallel response of identical sheet conductors of opposite dip indicated that the response of the separate sheets could be recognized even when the two sheets were placed at zero separation. This was found to be due in part to spatial displacement of the individual current vortices within each conductor owing to their mutual repulsion.     

Short-term precursor of the geomagnetic horizontal field transfer functions

IntroductionThe history of studying the changes in deep geoelectric conductivity of an earthquake preparation using the parameters of geomagnetic short-period variations has been about 30 years. At first, Japanese calculated the real transfer functions as well as the a and b coefficients of Parkinson vector by taking the measurement of graph and using empirical formula. There are some examples, such as the 1923 Kanto M=7.9 earthquake (Yanagihara, 1972), the 1966 Tashkent M=5.5 earthquake, t…     

Distribution of the aphid Hyalopterus pruni Geoffr. within and between habitats of Common Reed Phragmites australis (Cav.) Trin. ex Steudel as a result of migration and population growth

The aphid Hyalopterus pruni migrates to its secondary host Common Reed in late spring and early summer. Reed fields are known to be more heavily infested at the edges compared to the centre, as the migrants are attracted to the colour of the reed. There is also evidence that reeds growing at wet sites are more heavily infested than at dry sites. This paper shows that both proceses interact in determining the distribution of aphids in different habitats.The observations were carried out in the vast reed field of the Zuidelijk Flevoland polder (The Netherlands) during a number of years after its reclamation in 1968. In the period of migration flying aphids were caught with sticky traps situated in and around the reed field. The distribution of aphids on the reed was determined at different times during the season along transects through the centre of the polder. The existence early in the season (June) of a higher infestation at the edges than in the centre of the field was confirmed. The density of winged aphids in this early period was higher in wet than in dry habitats, which may have been a result of a preference of the settling migrants for reed standing in water. Later in the season (August) the dry habitat still had a relatively low infestation in the centre, but at the wet sites the infestation was high everywhere. This could be attributed to the much quicker population growth that was found in the wet habitat and subsequent short range movement of newly born winged aphids.     

Development of a methodology for assessing the environmental impact of radioactivity in Northern Marine environments

The requirement to assess the impacts of radioactivity in the environment explicitly and transparently is now generally accepted by the scientific community. A recently developed methodology for achieving this end for marine ecosystems is presented within this paper. With its clear relationship to an overarching system, the marine impact assessment is built around components of environmental transfer, ecodosimetry and radiobiological effects appraisal relying on the use of "reference organisms". Concentration factors (CFs), dynamic models and, in cases where parameters are missing, allometry have been employed in the consideration of radionuclide transfer. Dose conversion coefficients (DCCs) have been derived for selected flora and fauna using, inter alia, dose attenuation and chord distribution functions. The calculated dose-rates can be contextualised through comparison with dose-rates arising from natural background and chronic dose-rates at which biological effects have been observed in selected "umbrella" endpoints.     

Regional induction studies: A review of methods and results

The geomagnetic skin-effect is specified by setting three length scales in relation to each other: L₁ for the overhead source. L₂ for the lateral non-uniformity of the subsurface conductor, L₃ for the depth of penetration of a quasi-uniform transient field into this conductor. Relations for the skin-effect of a quasi-uniform source in layered conductors are generalized to include sources of any given geometry by introducing response kernels as functions of frequency and distance. They show that only those non-uniformities of the source which occur within a distance comparable to L₃ from the point of observation are significant. The skin-effect of a quasi-uniform source in a laterally non-uniform earth is expressed by linear transfer functions for the surface impedance and the surface ratio of vertical/horizontal magnetic variations. In the case of elongated structures and E-polarisation of the source, a modified apparent resistivity is defined which as a function of depth and distance gives a first orientation about the internal distribution of conductivity. The skin-effect of a non-uniform source in a non-uniform earth is considered for stationary and “running” sources. Recent observations on the sea floor and on islands indicate a deep-seated change of conductivity at the continent—ocean transition, bringing high conductivity close to the surface, a feature which may not prevail, however, over the full width of the ocean. There is increasingly reliable evidence for high conductivities (0.02 to 0.1 micro ^?1 m^?1) at subcrustal or even at crustal depth beneath certain parts of the continents, in some cases without obvious correlation to geological structure.     

Insights into the Structure and Surface Geology of Isla Socorro,Mexico, from Airborne Magnetic and Gamma-Ray Surveys

The island of Socorro is located in the eastern Pacific Ocean, 650 km off the coast of Mexico. It is a rare example of an oceanic volcanic island whose above sea level volume is made up mostly of peralkaline trachytes and rhyolites, with subordinate mafic rocks. Subaerial volcanism started several hundred thousand years ago and continues until recent times. We present an investigation of surface and subsurface geology of the island, based on the first detailed extensive geophysical survey on the island. Acquired airborne magnetic and gamma-ray data were compared to existing geological information and supplemented with field investigations and satellite imagery. Magnetic data show a wide minimum in the central part of the island, possibly connected to a high-temperature zone in the deeper central portion of the volcano, likely to be due to a still hot magma body. The data also depict two parallel edges possibly suggesting the existence of a nested caldera. Analysis on upward continued magnetic data by recent imaging techniques highlighted two deep sources located around 5 km b.s.l., interpreted as feeding structures that are now filled with crystalline rocks. Gamma-ray data have been interpreted through integration with the geological survey results. Several previously known volcanic deposits have been identified based on radioelement distribution, and others have been redefined based on field evidence. A new succession of volcanic members is proposed, to be verified through more detailed geological mapping, geochemical analyses of rock samples and radiometric dating.