水平场转换函数空间分布特征的数值模拟

应用三维有限差分法，对长方形导体的水平场转换函数分布特征进行数值模拟计算，给出了实部和虚部的分布规律，并研究了源场周期、导体导电率、导体埋深、导体长度对转换函数取值的影响．转换函数的极值分布在导体中央、4个犄角和边界附近，且随边界而移动．表明边界附近是转换函数的最佳观测位置．      

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.     

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…     

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

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

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.     

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.     

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.     

ANALYSIS AND IMPROVEMENT OF THE PRIMARY FIELD WAVEFORM OF THE GROUND TEM METHOD*

Many ground TEM systems use an on-off type of primary field waveform and measure the secondary field only in the transmitter off-time, so as to realize the pure secondary field measurement. The width of the trailing edge of the primary pulses and the duration of the transmitter on-time form a time-constant (T_c) window. Most effective excitation is achieved if the target conductor has a T_c falling into this window. To be able to measure in the transmitter on-time, two primary field waveforms are proposed. The first is a series of unipolar pulses with identical steep ramps at both leading and trailing edges. Measurements made in the on-time after the leading edges simply reduce the stacking time needed and power consumption by half. Also, long- T_c conductors are better excited. The second is a series of bipolar pulses with long ramps as the leading edges and short ramps as the trailing edges. This waveform presents two different T_c-windows at the same time and helps to improve the detection of long-T_c conductors in the presence of short- T_c ones.     

Conductivity anomalies: lower crust or asthenosphere?

The indication from surface measurements of a zone of relatively high conductivity (resistivity<200 ohm-m) at depths between 20 and 50 km has become so general over the Earth that regions without this zone can be considered anomalous. However, the depths indicated often span the lower crust and the uppermost mantle, so that before any effect can be definitely attributed to one or the other, the depth resolution in the electromagnetic measurements must be carefully considered. This paper applies the eigenvector decomposition of generalized linear inverse theory to soundings by Schlumberger resistivity, by magnetotellurics, by man-made electromagnetic fields formed by controlled current flow in grounded electric transmission lines, and by natural magnetic field variation studies to improve the bounds on depth, thickness and conductivity of a conductive layer. It is shown that many of the methods are capable of giving the depth to the top of a conductor with remarkable accuracy. Joint inversion of two or more of them offers an advantage in the separation of thickness and conductivity of both conductive and resistive layers. Natural geomagnetic field transfer functions, while accurately mapping the position of the edge of a conductor, do not provide the resolution of the other techniques, largely because the frequencies that can be practically measured at present are much too low.     

Two-Dimensional Magnetotelluric Modelling of Ore Deposits: Improvements in Model Constraints by Inclusion of Borehole Measurements

A combination of magnetotelluric (MT) measurements on the surface and in boreholes (without metal casing) can be expected to enhance resolution and reduce the ambiguity in models of electrical resistivity derived from MT surface measurements alone. In order to quantify potential improvement in inversion models and to aid design of electromagnetic (EM) borehole sensors, we considered two synthetic 2D models containing ore bodies down to 3000 m depth (the first with two dipping conductors in resistive crystalline host rock and the second with three mineralisation zones in a sedimentary succession exhibiting only moderate resistivity contrasts). We computed 2D inversion models from the forward responses based on combinations of surface impedance measurements and borehole measurements such as (1) skin-effect transfer functions relating horizontal magnetic fields at depth to those on the surface, (2) vertical magnetic transfer functions relating vertical magnetic fields at depth to horizontal magnetic fields on the surface and (3) vertical electric transfer functions relating vertical electric fields at depth to horizontal magnetic fields on the surface. Whereas skin-effect transfer functions are sensitive to the resistivity of the background medium and 2D anomalies, the vertical magnetic and electric field transfer functions have the disadvantage that they are comparatively insensitive to the resistivity of the layered background medium. This insensitivity introduces convergence problems in the inversion of data from structures with strong 2D resistivity contrasts. Hence, we adjusted the inversion approach to a three-step procedure, where (1) an initial inversion model is computed from surface impedance measurements, (2) this inversion model from surface impedances is used as the initial model for a joint inversion of surface impedances and skin-effect transfer functions and (3) the joint inversion model derived from the surface impedances and skin-effect transfer functions is used as the initial model for the inversion of the surface impedances, skin-effect transfer functions and vertical magnetic and electric transfer functions. For both synthetic examples, the inversion models resulting from surface and borehole measurements have higher similarity to the true models than models computed exclusively from surface measurements. However, the most prominent improvements were obtained for the first example, in which a deep small-sized ore body is more easily distinguished from a shallow main ore body penetrated by a borehole and the extent of the shadow zone (a conductive artefact) underneath the main conductor is strongly reduced. Formal model error and resolution analysis demonstrated that predominantly the skin-effect transfer functions improve model resolution at depth below the sensors and at distance of \(\sim \) 300–1000 m laterally off a borehole, whereas the vertical electric and magnetic transfer functions improve resolution along the borehole and in its immediate vicinity. Furthermore, we studied the signal levels at depth and provided specifications of borehole magnetic and electric field sensors to be developed in a future project. Our results suggest that three-component SQUID and fluxgate magnetometers should be developed to facilitate borehole MT measurements at signal frequencies above and below 1 Hz, respectively.     

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.     

An explanation of induced magnetic variations at Sabhawala,India

Induction vectors calculated for short-period events such as “bays” and “storm sudden commencements” (SSCs) indicate that the conductor causing the anomaly lies to the south of Sabhawala. The nature of vertical transfer functions obtained from spectral analysis of storm-time variations further confirms the presence of such a conductor. The nature of maximum (G_p) and minimum (G₁) response functions suggests that a two-dimensional conductor south of Sabhawala is responsible for anomalies in Z. Invoking a structural model of the Earth's crust beneath the Himalayas, in keeping with accepted orogenic theories, we propose that one structure capable of producing the desired conducting path is a rise of the asthenosphere south of the Himalayan range. Another important factor contributing to the induced variations at Sabhawala is due to a sedimentary trough running parallel to the Himalayas, again to the south. These views are confirmed by a two-dimensional model calculation. It is found that the undulation of the asthenosphere alone does not reproduce the nature of the anomalies in H and Z.     

均匀半空间导体的并矢格林函数解析式

导出了低频时谐电磁场中均匀半空间导体的并矢格林函数解析式,它由初等函数和修正贝塞尔函数的有限形状所表示．该解析式可用Kelvin函数的近似多项式计算,可显著提高并矢格林函数的计算速度．用"电磁分析方法验证模型"的第15号问题作为验证模型,结果表明,计算值和测量值吻合良好．     

Penetration of electric field from the near ground atmospheric layer to the ionosphere

A mathematical model has been proposed for describing quasi-stationary atmospheric electric fields with approximate, but fairly accurate allowance for ionospheric conductivity. It is shown that some well-known models of electric field penetration from the Earth into the ionosphere have been deemed inadequate, though they work well in the atmosphere below 50 km. In these models, the arbitrarily specified boundary condition in the upper boundary of the atmosphere omits the existing good conductor or adds not existent conductor. The maximum possible field in our model is far less than in models where ionospheric conductivity is not taken into account, but vastly larger than in models based on the approximation with infinite Pedersen conductivity in the upper ionosphere.     

New induction sounding tested in Central Europe

The hourly data of nine geomagnetic observatories situated in Central Europe have been analyzed using the generalized magnetovariation (GMV) method designed recently for induction soundings of inhomogeneous media. In this method, impedance is one of transfer functions in the differential relation between spectra of the magnetic components and their derivatives. The peculiarity of this impedance is its correspondence to the magnetotelluric one estimated from the linear relations. Three transfer functions have been estimated simultaneously for data of geomagnetic observatories, using three different routines working in the period range from three hours up to two days. Noises in the source field components have been compared with noise in the estimated plane field divergence. The multivariate errors-in-variables method was used to extract spatially and temporally coherent geomagnetic field structure from the partially incoherent geomagnetic variations. This method allows estimating reliably impedances and gradient tippers for each observatory, taking into consideration the Earth’s sphericity. The obtained responses have been used for induction soundings and for detecting a deep inhomogeneity in the region.     

An airborne tensor VLF system. From concept to realization1

Radio signals from very low frequency (VLF) transmitters distributed world-wide have been used for several decades to study the lateral variations of the electrical conductivity in the upper few hundred metres of the earth's crust. Traditionally, in airborne applications, the total magnetic fields from one or two transmitters are measured to form the basis for construction of maps that primarily show those conductive structures that are parallel or subparallel to the direction to the transmitters. The tensor VLF technique described in this paper makes use of all signals available in a predefined frequency band to construct transfer functions relating the vertical magnetic field and the two horizontal magnetic field components. These transfer functions are uniquely determined for a particular measuring site and contain information about the lateral conductivity variations in all directions. First experiences with real field data, acquired during a test survey in Sweden, show that maps of the so-called peaker, the spatial divergence of the transfer functions, give an image of the conducting structures. Most of the structures can be correlated to small valleys filled with conducting sediments or valleys underlain by conductive fracture zones in the crystalline rocks.     

地磁转换函数与汶川大地震关系的初步研究

地磁场转换函数在强震前会发生变化,这已被广泛公认。利用汶川周围9个“十五”建设的地磁台站分钟值数字记录数据,分析和研究了汶川地震前的地磁转换函数变化特征,结果表明：汶川地震的发生与周边台站地磁转换函数异常的关系是明显的,地磁短周期转换函数方法是一种物理意义明确、并有一定预报效能的好方法。     

青藏高原中南部岩石圈扩张应力场与羊八井地热异常形成机制

青藏高原中南部地区有着丰富的地热资源,其中羊八井地区存在着温度高达93℃以上,热流值高达364mW/m2的地热奇观.本文利用国际以及中国地震台网资料,详细分析了青藏高原中南部地震活动性和岩石圈脆性破裂深度范围,进而讨论了高地热异常区的岩石圈应力场背景及其动力学机制.研究结果表明,青藏高原地热异常区有高于高原内部平均活动水平的震级为6级左右的中强地震活动,其震中分布呈现出与地热活动相似的、沿近南北向断陷带的带状特征;并且在高热流区域内有高原内部罕见的中深地震活动,地震震源从100km左右深处到地表呈现柱状地震群活动的空间分布等特征.震源机制与应力场研究结果表明,虽然高原中南部应力场主压、主张应力方向与青藏高原的整体特征相符,但是地震发生类型与青藏高原周缘的挤压逆断层型地震完全不同,均属于东西向扩张力作用下的正断层型地震活动,特别是在羊八井高热流区域附近,东西向扩张应力场在岩石圈应力场中起到主导性作用,推测其控制深度可达岩石圈底部100多公里处.青藏高原地热异常区在强烈的近东西向扩张应力场作用下,岩石圈东西向扩张并发生一系列的大规模正断层活动,致使深部软流圈高温热流可以沿着活动正断层及其形成的深裂隙上涌,穿过岩石圈到达地表面,形成了高地热异常区.