A general correspondence principle for time-domain electromagnetic wave and diffusion fields

A general correspondence principle is presented that relates any time-domain electromagnetic diffusion field to an electromagnetic wavefield in a 'corresponding' configuration. The principle applies to arbitrarily inhomogeneous and anisotropic media and arbitrary transmitters and receivers. For the correspondence between the two types of electromagnetic fields to hold, the electric conductivity in the diffusive case and the permittivity in the wavefield case should have the same spatial variation, while the permeability distributions in space in the two cases are to be identical. Essential steps in the derivation of the correspondence principle are the use of the time Laplace transformation of causal signals, taken at real, positive values of the transform parameter, the Schouten-Van der Pol theorem in the theory of the Laplace transformation, and the reliance upon Lerch's theorem of the uniqueness of the interrelation between causal field quantities and their time-Laplace-transform representations at real, positive values of the transform parameter. Correspondence is then established between the tensorial Green's functions in the two cases, where the Green's functions are the point-receiver responses (either electric or magnetic field) to point-transmitter excitations (either electric- or magnetic-current source).
Through the correspondence principle, all transient electromagnetic wavefields (where losses are neglected) have as a counterpart a transient diffusive electromagnetic field (where the electric displacement current is neglected). The interrelation yields the tool to compare quantitatively the potentialities of the two types of fields in transient electromagnetic geophysical prospecting.
Finally, a general medium-parameter scaling law for time-domain electromagnetic wavefields is presented.     

Calculating the two-dimensional magnetotelluric Jacobian in finite elements using reciprocity

To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p _a), the impedance phase (φ) and the vertical magnetic field ( K _zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude.     

Genesis and distribution of mineral waters as a consequence of recent lithospheric dynamics: the Rhenish Massif,Central Europe

Three major, interdependent processes control the genesis and distribution of mineral and thermal waters in the Rhenish Massif, Central Europe: (a) Magmatic processes in the upper mantle provide most of the CO₂ to produce bicarbonate waters in shallow aquifers. (b) Extension of the brittle upper crust enables the ascent of sodium chloride waters. (c) Uplift and erosion shape the massif's relief, which determines the extent of flow systems and the distribution of thermal springs. The chemistry of mineral waters further depends on the aquifers' mineral composition. A comprehensive set of hydrological, chemical, tectonic and geophysical data on the Rhenish Massif has been compiled. It was used to classify the mineral waters and to map the spatial distribution of water properties. The composition of cuttings from several representative wells producing different water types shows that the hydrothermal alteration of the aquifer rocks consists mainly of kaolinization of chlorite and dissolution of feldspar. Numerical transport simulations favour two modes of groundwater flow: topography-driven flow and the pressure-driven ascent of basement brines along active faults. Thermal convection is less important.     

变质反应的代数识别方法

在回顾变质矿物共生组合和变质反应图解分析法的基础上，论述了录求变质矿物组合之间、同一组合内矿物之间的平衡变质反应的代数方法，即矩阵高斯变换，线性回归分析、矩阵奇值分解。     

内蒙大水清金矿带交代作用及其与金矿化关系

大水清金矿带是华北克拉通中有代表性的改造型绿岩金矿之一。它产于前东武纪绿岩地体中安家营子石英二长岩体内外接触带附近的剪切带中。金矿化交代主要受后阶段张性脆裂控制，而剪切带早阶段韧性剪切作用无明显矿化交代，后者主要为成矿交代提供有益的江游交代构造环，境即成矿交代蚀变作用和剪切带同步的两阶段演化模式。大水清金矿带的交代作用有四个阶段：即用长石化阶段，绿泥石化阶段，黄铁绢英岩化和石英－硫化物阶段以及石英－碳酸盐化阶段。第三阶段由五个亚阶段组成：即黄铁绢英岩化和强硅化－强黄铁矿化，石英－多金属硫化物交代，石英－贫扶矿充填交代和细粒黄铁矿充填交代。其中前四个亚阶段是本区金的主要成矿期．当韧性剪切转化为张性脆裂时产生的骤然扩容和伴随发生的交代体系转变为富硅体系的钾交代和石英－硫化物交代是金矿成矿的关键。本文还研究和探讨了该金矿成矿交代流体的性质及交代的物理化学环境。     

蚀变带元素地球化学成矿预测法──以金龙山微细浸染型金矿为例

金矿化作用伴围岩蚀变，且多受断裂控制，反映在地球化学上出现某些元素的异常。在已发现一害矿化体并掌握一般控制矿特征的地区，利用已有工程在蚀变带内连续捡块采样，分析有关批示元素，总结有关的地球化学指标，建立地球化学预测模式，对已知矿体及外围异常进行了地球化学评价。     

判定地下水二维水动力弥散参数的直线图解法

弥散参数的求得，是正确认识地下水污染运移及对其控制的前提。本文根据极值理论和非线性方程线性化手段，导出了浓度和时间变化的线性关系，得出了确定弥散参数的直线图解法。它可以弥补其它方法的不足。借用实际弥散试验资料对其进行了检验，表明所建数学模型与实际模型之间吻合良好，即所求得的参数结果可信。     

Seismic Ground Motion Expected for the Eastern District of Naples

The seismic ground motion of a test area in the eastern district of Naples is computed with a hybrid technique based on the mode summation and the finite difference methods. This technique allows us the realistic modelling of source and propagation effects, including local soil conditions. In the modelling, we consider the 1980 Irpinia earthquake, a good example of strong shaking for the area of Naples, which is located about 90 km from the epicenter.The detailed geological setting is reconstructed from a large number of drillings. The sub-soil is mainly formed by alluvial (ash, stratified sand and peat) and pyroclastic materials overlying a pyroclastic rock (yellow neapolitan tuff), representing the neapolitan bedrock. The detailed information available on mechanical properties of the sub-soil and its geometry warrants the application of the sophisticated hybrid technique.As expected, the sedimentary cover causes an increase of the signal's amplitudes and duration. If thin peat layers are present, the amplification effects are reduced, and the peak ground accelerations are similar to those observed for the bedrock model. This can be explained by the backscattering of wave energy at such layers, that tend to seismically decouple the upper from the lower part of the structure.For SH-waves, the influence of the variations of the S-wave velocities on the spectral amplification is studied, by considering locally measured velocities and values determined from near-by down-hole measurements. The comparison between the computed spectral amplifications confirms the key role of an accurate determination of the seismic velocities of the different layers.The comparison performed between a realistic 2-D seismic response and a standard 1-D response, based on the vertical propagation of waves in a plane layered structure, shows considerable difference, from which it is evident that serious caution must be taken in the modelling of expected ground motion at a specific site.     

Regionalization of unit hydrograph Parameters: 2. Uncertainty analysis

Hydrologic model parameters obtained from regional regression equations are subject to uncertainty. Consequently, hydrologic model outputs based on the stochastic parameters are random. This paper presents a systematic analysis of uncertainty associated with the two parameters, N and K, in Nash's IUH model from different regional regression equations. The uncertainty features associated with N and K are further incorporated to assess the uncertainty of the resulting IUH. Numerical results indicate that uncertainty of N and K from the regional regression equations are too significant to be ignored.