The electric properties of eclogites from the Bohemian Massif under high temperatures and pressures

Summary The D.C. and A.C. electric conductivity and the dielectric constant of a set of 8 eclogites from the Bohemian Massif is studied under high temperatures (200–900°C) and pressures (up to 20 kb). The electric conductivity of the studied eclogites is strongly affected by the content of symplectite minerals: their conductivity increases with their concentration. The average increase of conductivity with pressure between 1 kb and 20 kb is about 50% and is manifested particularly below 4 kb. The correlation between the electric conductivity and the dielectric constant is very high (r=0.975) and allows for the values of one parameter to be estimated from those of another.     

Laboratory measurements of electrical properties of rocks and minerals

In this paper all important results on laboratory measurements of electrical properties of rocks and minerals for the last four years are presented. Although basic results from all over world are reviewed, East European and U.S.S.R. works are described more extensively. Only D.C. conductivity results are considered, however all principal conclusions hold in similar form for both D.C. and A.C. conductivities. At the end the most important areas of further study and measurement are depicted and the main problems of the future are outlined.     

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.     

利用深井观测速度型记录的最大振幅测定近震震级的一种方法

根据上海台网观测记录到的上海及邻近地区1985年至1990年期间179次地震事件的资料,以虹桥台、南汇台的两口深井观测地震记录为例,探讨如何利用深井观测速度型记录的最大振幅来测定近震震级。对于400公里以内的地震,本文根据统计规律,提出:在已知系统速度灵敏度K,地震记录的最大振幅A和地震震中距Δ的前提下,速度型记录的近震震级可以表示为M_v=logA/K+blog(Δ)+c误差分析表明,M_v与M_L是比较接近的,结果是可用的。     

Electrical conductivity of magmatic liquids: effects of temperature,oxygen fugacity and composition

The effects of temperature, f_O2 and composition on the electrical conductivity of silicate liquids have been experimentally determined from 1200 to 1550°C under a range of f_O2 conditions sufficient to change the oxidation state of Fe from predominantly Fe²⁺ to Fe³⁺. Oxidation of ferrous to ferric iron in the melt has no measurable effect on the conductivity of melts with relatively low ratios of divalent to univalent cations. Under strongly oxidizing conditions a minor decrease of conductivity is detected inth highΣM²/ΣM⁺ ratios. It is concluded that for purposes of estimating the conductivity of magmatic liquids, f_O2 may be ignored to a first approximation. Both univalent and divalent cation transport is involved in electrical conduction. Melts relying heavily on divalent cations for conduction, i.e. melts with relatively large ΣM²⁺/ΣM⁺ ratios, show strong departures from Arrheenius temperature dependence with the apparent activation energies decreasing steadily as the temperature increases. Conductivities dominated by the univalent cations, in melts with relatively small ΣM²⁺/ΣM⁺ ratios, show classical Arrhenius temperature dependence. These observations are discussed in terms of the general characteristics of the melt structure.Compositional variations within the magmatic range account for much less than an order of magnitude variation in electrical conductivity at a fixed temperature. This observation, combined with previous measurements of the conductivity of olivine (A. Duba, H.C. Heard and R. Schock, 1974) make it possible to state with reasonable confidence that melts occurring within the mantle will be more conductive by 3–4 orders of magnitude than their refractory residues. Potential applications to geothermometry are discussed.     

A 3D ERT study of solute transport in a large experimental tank

A high resolution, cross-borehole, 3D electrical resistivity tomography (ERT) study of solute transport was conducted in a large experimental tank. ERT voxels comprising the time sequence of electrical images were converted into a 3D array of ERT estimated fluid conductivity breakthrough curves and compared with direct measurements of fluid conductivity breakthrough made in wells. The 3D ERT images of solute transport behaviour were also compared with predictions based on a 3D finite-element, coupled flow and transport model, accounting for gravity induced flow caused by concentration differences.The tank (dimensions 185×245×186 cm) was filled with medium sand, with a gravel channel and a fine sand layer installed. This heterogeneous system was designed to complicate solute transport behaviour relative to a homogeneous sand tank, and to thus provide a challenging but insightful analysis of the ability of 3D ERT to resolve transport phenomena. Four ERT arrays and 20 piezometers were installed during filling. A NaCl tracer (conductivity 1.34 S/m) was injected and intensively monitored with 3D ERT and direct sampling of fluid chemistry in piezometers.We converted the bulk conductivity estimate for 250 voxels in the ERT imaged volume into ERT estimated voxel fluid conductivity by assuming that matrix conduction in the tank is negligible. In general, the ERT voxel response is in reasonable agreement with the shape of fluid conductivity breakthrough observed in six wells in which direct measurements of fluid conductivity were made. However, discrepancies occur, particularly at early times, which we attribute to differences between the scale of the image voxels and the fluid conductivity measurement, measurement errors mapped into the electrical inversion and artificial image roughness resulting from the inversion.ERT images revealed the 3D tracer distribution at 15 times after tracer injection. The general pattern and timing of solute breakthrough observed with ERT agreed with that predicted from the flow/transport modelling. However, the ERT images indicate a vertical component of tracer transport and preferential flow paths in the medium sand. We attribute this to transient vertical gradients established during tracer injection, and heterogeneity caused by sorting of the sand resulting from the filling procedure. In this study, ERT provided a unique dataset of 250 voxel breakthrough curves in 1.04 m³. The use of 3D ERT to generate an array of densely sampled estimated fluid conductivity breakthrough curves is a potentially powerful tool for quantifying solute transport processes.     

基于二维图像的数字岩心电导率计算方法研究

电导率是表征岩石电学性质的重要物理参数,在地质资源勘查和测井解释等领域发挥着巨大作用.快速、准确地确定岩石电导率具有重要的理论和实践意义.作为近年来发展的一种岩石物理数值模拟工具,数字岩心技术在定量计算电导率等物性参数方面应用广泛.三维微观结构的准确获取是数字岩心技术计算岩石电导率的关键,但传统获取岩石三维微观结构的方法较为复杂费时.为了方便快速地通过数字岩心技术计算岩石的电导率,本文研究了岩石二维与三维数字岩心的电导率联系.我们基于微米级X射线CT扫描得到的三个砂岩样品的微观结构信息建立了三维数字岩心,并通过有限元法计算的三维数字岩心电导率与实验数据的对比验证数值计算方法的有效性.随后我们数字地扩展了岩石的孔隙,产生了较大孔隙度的三维数字岩心样本,在此基础上,计算了三维数字岩心和相应二维数字岩心的电导率,并通过Archie公式分别拟合了电导率与孔隙度之间的关系,得到了相应的胶结系数.结果表明,三维数字岩心的胶结系数小于二维数字岩心的胶结系数,且二者的比值与岩石实测孔隙度呈线性负相关关系.以该联系为纽带,通过二维图像快速计算得到的电导率与孔隙度关系,确定了三维数字岩心的电导率与孔隙度关系,并进一步通过三维数字岩心的孔隙度计算其电导率.该方法计算得到的人工砂岩样品的电导率与其三维数字岩心电导率相关系数高于96%,验证了基于二维图像的数字岩心电导率计算方法的有效性.本文的研究结果为快速、准确地计算岩石电导率提供了新的思路,在油气勘探开发中有广阔的应用前景.     

The influence of oxidation state on the electrical conductivity of olivine

The electrical conductivity of a single crystal of San Carlos olivine (Fo₉₂, 0.16 wt.% Fe₂O₃) has been measured as a function of temperature and oxygen fugacity (fO₂). After heating to 1338°C at fO₂ = 10^?12 atm., the conductivity at 950°C was 10^?5 (ohm-m)^?1, almost 3 orders of magnitude less than that measured in air. This decrease is due to the reduction of Fe³⁺ to Fe²⁺. Further heating to 1500°C at fO₂ = 10^?14 atm., decreased the electrical conductivity at 950°C to 10^?6 (ohm-m)^?1. When recovered at room temperature, the speciment had Fo₉₆ composition and contained small, opaque blebs distributed throughout the crystal. Derivations of temperature profiles for the earth's mantle from conductivity-depth models must take account of the important role played by iron oxidation state in the electrical conductivity of olivine.     

Bulk density effects on soil hydrologic and thermal characteristics: A numerical investigation

Soil bulk density (ρ_b) is commonly treated as static in studies of land surface dynamics. Magnitudes of errors associated with this assumption are largely unknown. Our objectives were to (a) quantify ρ_b effects on soil hydrologic and thermal properties and (b) evaluate effects of ρ_b on surface energy balance and heat and water transfer. We evaluated 6 soil properties, volumetric heat capacity, thermal conductivity, soil thermal diffusivity, water retention characteristics, hydraulic conductivity, and vapour diffusivity, over a range of ρ_b, using a combination of 6 models. Thermal conductivity, water retention, hydraulic conductivity, and vapour diffusivity were most sensitive to ρ_b, each changing by fractions greater than the associated fractional changes in ρ_b. A 10% change in ρ_b led to 10–11% change in thermal conductivity, 6–11% change in saturated and residual water content, 49–54% change in saturated hydraulic conductivity, and 80% change in vapour diffusivity. Subsequently, 3 field seasons were simulated with a numerical model (HYDRUS‐1D) for a range of ρ_b values. When ρ_b increased 25% (from 1.2 to 1.5 Mg m^?3), soil temperature variation decreased by 2.1 °C in shallow layers and increased by 1 °C in subsurface layers. Surface water content differed by 0.02 m³ m^?3 for various ρ_b values during drying events but differences mostly disappeared in the subsurface. Matric potential varied by >100 m of water. Surface energy balance showed clear trends with ρ_b. Latent heat flux decreased 6%, sensible heat flux increased 9%, and magnitude of ground heat flux varied by 18% (with a 25% ρ_b increase). Transient ρ_b impacted surface conditions and fluxes, and clearly, it warrants consideration in field and modelling investigations.     

Thermal properties of vesicular rhyolite

Thermal diffusivity of rhyolite melt and rhyolite foam (70–80% porosity) has been measured using the radial heat transfer method. Cylindrical samples (length 50–55 mm, diameter 22 mm) of rhyolite melt and foam have been derived by heating samples of Little Glass Mountain obsidian. Using available data on heat capacity and density of rhyolite melt, the thermal conductivity of samples has been determined. The difference in thermal conductivity between rhyolite melt and foam at igneous temperatures ( 1000°C) is about one order of magnitude. The effect of thermal insulation of magmas due to vesiculation and foaming of the top layer is discussed in terms of the data obtained using a simple illustrative model of magma chamber convection.     

Semianalytical models of sprite formation from plasma inhomogeneities

A spherical plasma inhomogeneity located at mesospheric altitudes in a thundercloud quasi-electrostatic field is considered as a possible cause of sprite formation. A simple semianalytical model of ionization instability in a quasi-electrostatic field, the value of which is larger than the air breakdown value, is developed on the assumption that plasma ball conductivity is controlled by impact ionization and electron attachment to neutrals. After several simplifications, the problem is reduced to a system of ordinary differential equations for the average conductivity and plasma ball radius. The analytical estimates and numerical simulation indicate that the predicted expansion rate and acceleration of the plasma inhomogeneity boundary are close in magnitude to the values observed during high-speed imaging of sprite development.     

In-Situ Electrical Conductivity and Permeability of Mid-Crustal Rocks from the Ktb Drilling: Consequences for High Conductive Layers in the Earth Crust

Freshly cored samples from a microprofile (7011–7013m in depth) of the German Continental Deep Drilling Project (KTB) were taken to measure the complex electrical conductivity (1 kHz up to 1 MHz), porosity, BET-surface, permeability and density. The porosity ranged about 1 vol%, while the permeability k varied from 16.05 µD to > 0.01 µD for in-situ pressure conditions. The permeability decreased about 2 orders in magnitude up to pressures of 200 MPa. Conductivity was measured in the same pressure range on 1 M NaCl saturated samples. Thin sections and SEM analysis revealed an enrichment of carbon and ilmenite (about 1 vol%) on inner cleavage cracks of mica, thus causing an unusual high (ranging from 4.2 × 10^-3 S/m to 67 × 10^-3 S/m) being orders of magnitude higher than normally measured on such types of rocks (about 300 × 10^-6 S/m). An inverse pressure dependence of was detected on some of the samples. Electronic conduction was confirmed by least-squares-fits of model data to the frequency dispersion of the conductivity and by measuring the time dependence of the volume conductivity and its frequency dispersion. Thus the dominating role of the reconnected network of carbon and ilmenite on the enhanced volume conductivity was proved. An increase of the conductivity due to hydrofracturing by high pore fluid pressures plays a less important role.     

3D modelling of near-surface, environmental effects on AEM data

This study considers the three-dimensional (3D) modelling of compact, at-surface conductive bodies on frequency domain airborne electromagnetic (AEM) survey data. The context is the use of AEM data for environmental and land quality applications. The 3D structures encountered are typically conductive, of limited thickness (<20 m) and form ‘point’ source locations carrying potential environmental risk. The scale of such bodies may generate single-profile, ‘bulls-eye’ anomalies. In attempts to recover geological information, such anomalies may be considered to represent noise. In environmental AEM, the correct interpretation of such features is important. The study uses a combination of theoretical models and trial-fixed-wing survey data obtained in populated areas of the UK. Scale issues are discussed in terms of the volumetric footprints of the induced electric field generated by systems flown at both low and high elevation. One of the primary uses of AEM survey data lies in the assessment of conductivity maps. These are typically obtained using one-dimensional (1D) conductivity models at individual measurement points. In order to investigate the limitations of this approach, 3D modelling of conductive structures with dimensions less than 350×350 m and thicknesses extending to 20 m has been carried out. A 1D half space inversion of the data obtained at each frequency is then used to assess the behaviour of the spatial information. The results demonstrate that half space conductivity values obtained over compact 3D targets generally provide only apparent conductivity results. For thin, at-surface bodies, conductivity values are biased to lower values than the true conductivity except at high frequency. The spatial perturbation to both coupling ratios and 1D conductivity models can be laterally extensive. The results from 3D modelling indicate that the use of horizontal derivatives applied to the conductivity models offers enhanced edge detection. The practical application of such derivatives to both regional- and local-scale survey data is presented.. The special case of a near-surface, metallic pipeline has been modelled. The problem constitutes an inductive limit (current gathering) response in which the perturbation is largely confined to the in-phase coupling ratios. The main perturbations, in data and conductivity models, are within about 40 m of each side of the pipeline. The maximum perturbation to the conductivity model is only a factor of 1.5 above background. Detailed survey data across a former compact landfill (about 100×100 m) are used to compare the model behaviour predicted by the 3D modelling with survey results. The survey, conducted at two separate altitudes, provides a demonstration of 3D effects on 1D survey models as a function of frequency and elevation. Although the nature of the landfill materials and their location are not known precisely, the mapping information appears realistic.     

Upscaling of flow in heterogeneous porous formations: Critical examination and issues of principle

We consider heterogeneous media whose properties vary in space and particularly aquifers whose hydraulic conductivity K may change by orders of magnitude in the same formation. Upscaling of conductivity in models of aquifer flow is needed in order to reduce the numerical burden, especially when modeling flow in heterogeneous aquifers of 3D random structure. Also, in many applications the interest is in average values of the dependent variables over scales larger or comparable to the conductivity length scales. Assigning values of the conductivity K_b to averaging domains, or computational blocks, is the topic of a large body of literature, the problem being of wide interest in various branches of physics and engineering. It is clear that upscaling causes loss of information and at best it can render a good approximation of the fine scale solution after averaging it over the blocks.The present article focuses on upscaling approaches dealing with random media. It is not meant to be a review paper, its main scope being to elucidate a few issues of principle and to briefly discuss open questions. We show that upscaling can be usually achieved only approximately, and the result may depend on the particular upscaling scheme adopted. The typically scarce information on the statistical structure of the fine-scale conductivity imposes a strong limitation to the upscaling problem. Also, local upscaling is not possible in nonuniform mean flows, for which the upscaled conductivity tensor is generally nonlocal and it depends on the domain geometry and the boundary conditions. These and other limitations are discussed, as well as other open topics deserving further investigation.     

二维大地电磁资料频域逆散射反演

逆散射反演是由已知散射波场反演地下介质中散射体各参数了分布的方法，本文从这一思想出发，针对二维大地电磁资料，在BORN近似的前提下，提出了一种新的反地下电导率参数的方法，经模型试算与实际资料的反演验证正确，这一方法特别适于地下无强电性分界面的地区。     

Quantifying the effects of three-dimensional subsurface heterogeneity on Hortonian runoff processes using a coupled numerical,stochastic approach

The impact of three-dimensional subsurface heterogeneity in the saturated hydraulic conductivity on hillslope runoff generated by excess infiltration (so-called Hortonian runoff) is examined. A fully coupled, parallel subsurface–overland flow model is used to simulate runoff from an idealized hillslope. Ensembles of correlated, Gaussian random fields of saturated hydraulic conductivity are used to create uncertainty in spatial structure. A large number of cases are simulated in a parametric manner with the variance of the hydraulic conductivity varied over orders of magnitude. These cases include rainfall rates above, equal and below the geometric mean of the hydraulic conductivity distribution. These cases are also compared to theoretical representations of runoff production based on simple assumptions regarding (1) the rainfall rate and the value of hydraulic conductivity in the surface cell using a spatially-indiscriminant approach; and (2) a percolation-theory type approach to incorporate so-called runon. Simulations to test the ergodicity of hydraulic conductivity on hillslope runoff are also performed. Results show that three-dimensional stochastic representations of the subsurface hydraulic conductivity can create shallow perching, which has an important effect on runoff behavior that is different than previous two-dimensional analyses. The simple theories are shown to be very poor predictors of the fraction of saturated area that might runoff due to excess infiltration. It is also shown that ergodicity is reached only for a large number of integral scales (∼30) and not achieved for cases where the rainfall rate is less than the geometric mean of the saturated hydraulic conductivity.     

历史有感地震目录编制原则方法的研究——以《中国历史有感地震目录》编制为例

在吸收历史地震研究、历史政区地理研究等多方面成果基础上,基于历史有感地震记载特点,探索了历史有感地震目录的编制原则,给出了地震日期、震中位置、震级的确定方法。以《中国地震历史资料汇编》等为依据,编制了公元前618年-公元1949年期间中国境内3≤M<443/4级历史有感地震目录9 121条。目录数量为中国同期文字记载的43/4级以上破坏性地震的7倍多,比全球其他国家或地区最早的同类目录早1 695年,约占全球同期M≤4.5级有感地震目录的2/3。     

Estimation of hydraulic conductivity in an alluvial system using temperatures

Well water temperatures are often collected simultaneously with water levels; however, temperature data are generally considered only as a water quality parameter and are not utilized as an environmental tracer. In this paper, water levels and seasonal temperatures are used to estimate hydraulic conductivities in a stream-aquifer system. To demonstrate this method, temperatures and water levels are analyzed from six observation wells along an example study site, the Russian River in Sonoma County, California. The range in seasonal ground water temperatures in these wells varied from <0.2 degrees C in two wells to approximately 8 degrees C in the other four wells from June to October 2000. The temperature probes in the six wells are located at depths between 3.5 and 7.1 m relative to the river channel. Hydraulic conductivities are estimated by matching simulated ground water temperatures to the observed ground water temperatures. An anisotropy of 5 (horizontal to vertical hydraulic conductivity) generally gives the best fit to the observed temperatures. Estimated conductivities vary over an order of magnitude in the six locations analyzed. In some locations, a change in the observed temperature profile occurred during the study, most likely due to deposition of fine-grained sediment and organic matter plugging the streambed. A reasonable fit to this change in the temperature profile is obtained by decreasing the hydraulic conductivity in the simulations. This study demonstrates that seasonal ground water temperatures monitored in observation wells provide an effective means of estimating hydraulic conductivities in alluvial aquifers.     

Evidence and effects of a wave-driven nonlinear current in the equatorial electrojet

Ionospheric two-stream waves and gradientdrift waves nonlinearly drive a large-scale (D.C.) current in the E-region ionosphere. This current flows parallel to, and with a comparable magnitude to, the fundamental Pedersen current. Evidence for the existence and magnitude of wave-driven currents derives from a theoretical understanding of E-region waves, supported by a series of nonlinear 2D simulations of two-stream waves and by data collected by rocket instruments in the equatorial electrojet. Wave-driven currents will modify the large-scale dynamics of the equatorial electrojet during highly active periods. A simple model shows how a wave-driven current appreciably reduces the horizontally flowing electron current of the electrojet. This reduction may account for the observation that type-I radar echoes almost always have a Doppler velocity close to the acoustic speed, and also for the rocket observation that electrojet regions containing gradientdrift waves do not appear also to contain horizontally propagating two-stream waves. Additionally, a simple model of a gradient-drift instability shows that wavedriven currents can cause nonsinusoidal electric fields similar to those measured in situ.     

Electrical conductivity of a spinel lherzolite and a garnet peridotite to 1550°C: relevance to the effects of partial melting

Electrical conductivity σ of two ultramafic rocks (a spinel lherzolite and a garnet peridotite) has been investigated to melting temperature at 1 bar under known oxygen fugacity environment. The electrical conductivity of the two rocks is found to increase with degree of partial melting and an ～ 15% melt fraction is necessary for the electrical conductivity to increase by ～ 1 order of magnitude. For a given melt fraction electrical conductivity of a spinel lherzolite is lower than that of a garnet peridotite and may be attributed to the differences in the composition of the melts formed.