地球电磁法研究新进展——“第19届国际地球电磁感应学术研讨会”专辑

1 研讨会概况 地球电磁法是地球物理学的重要方法,它是通过观测自然的和人工产生的电场、磁场或电磁场研究自地表到地幔深处的电性结构和空间环境. 由IAGA I-2 工作组主办的两年一届的"国际地球电磁感应学术研讨会"的目的是展示和交流各种电磁方法的最新研究成果和前沿研究课题.每次研讨会的举办地点经过竞争和选举确定,1972年在Edinburgh, U.K.举行了第1届研讨会.     

Modelling the wave phenomena in acoustic and elastic media with sharp variations of physical properties using the grid‐characteristic method

This paper introduces a novel method of modelling acoustic and elastic wave propagation in inhomogeneous media with sharp variations of physical properties based on the recently developed grid‐characteristic method which considers different types of waves generated in inhomogeneous linear‐elastic media (e.g., longitudinal, transverse, Stoneley, Rayleigh, scattered PP‐, SS‐waves, and converted PS‐ and SP‐waves). In the framework of this method, the problem of solving acoustic or elastic wave equations is reduced to the interpolation of the solutions, determined at earlier time, thus avoiding a direct solution of the large systems of linear equations required by the FD or FE methods. We apply the grid‐characteristic method to compare wave phenomena computed using the acoustic and elastic wave equations in geological medium containing a hydrocarbon reservoir or a fracture zone. The results of this study demonstrate that the developed algorithm can be used as an effective technique for modelling wave phenomena in the models containing hydrocarbon reservoir and/or the fracture zones, which are important targets of seismic exploration.     

Research Note: First results of a low‐altitude unmanned aircraft system gamma survey by comparison with the terrestrial and aerial gamma survey data

In order to increase the efficiency of geological survey in difficult landscape‐morphological conditions, the authors have developed low‐altitude gamma technology surveys based on a multi‐rotor unmanned aircraft system, also known informally as Drones. These results have been compared against those from standard terrestrial and aerial gamma surveys. The successful results obtained at the reference site show that the developed methodology enables compilation of high‐quality data in difficult landscape conditions. These results yield the same level of information content as standard terrestrial gamma surveys. The exact unmanned aircraft system altitude over the earth must be maintained for carrying out high‐quality surveys. It was found that a CsI(Tl) detector with sizes 80 × 80 mm (vol. 0.4 dm³) enables obtaining of high‐quality radiometric data at altitudes of up to 70 m, at a speed of at least 20 km/h. The authors believe that low altitude ‘quasi‐terrestrial’ gamma surveys with a terrain drape may in the future replace terrestrial gamma surveys since they can offer better effectiveness at a lower cost.     

The cycling of Ni,Zn, Cu in the system “mine tailings–ground water–plants”: A case study

The comparative behaviour of Ni, Cu and Zn in the system “mine tailings–ground water–plants” has been investigated at the Ni–Cu mine site operated by INCO Ltd. Thompson Operations, Thompson, Manitoba. Oxidation of sulphide minerals causes the release of metals from exposed tailings containing Ni ∼2000 ppm, Cu ∼150 ppm and Zn ∼100 ppm to the ground water, which contains 350 mg/L Ni, 0.007 mg/L Cu, and 1.6 mg/L Zn. The metal concentration in the ground water is affected by the relative proportions of sulfide minerals, the rate of oxidation of sulphide minerals (Ni-bearing pyrrhotite > sphalerite > chalcopyrite), and the affinity of the metals for secondary Fe-phases (Ni > Zn > Cu).     

A simple predictive model of quartz solubility in water-salt-CO2 systems at temperatures up to 1000 °C and pressures up to 1000 MPa

Knowledge of the solubility of quartz over a broad spectrum of aqueous fluid compositions and T-P conditions is essential to our understanding of water-rock interaction in the Earth’s crust. We propose an equation to compute the molality of aqueous silica, m_SiO2(aq), mol·(kg H₂O)⁻¹, in equilibrium with quartz and water-salt-CO₂ fluids, as follows:

     

Properties of adjoint sea ice sensitivities to atmospheric forcing and implications for the causes of the long term trend of Arctic sea ice

Based on adjoint sensitivities of the coupled Massachusetts Institute of Technology ocean–sea ice circulation model, the potential influence of thermodynamic atmospheric forcing on the interannual variability of the September sea ice area (AREA) and volume (VOLUME) in the Arctic is investigated for the three periods 1980–1989, 1990–1999 and 2000–2009. Sensitivities suggest that only large forcing anomalies prior to the spring melting onset in May can influence the September sea ice characteristics while even small changes in the atmospheric variables during subsequent months can significantly influence September sea ice state. Specifically, AREA close to the ice edge in the Arctic seas is highly sensitive to thermodynamic atmospheric forcing changes from June to July. In contrast, VOLUME is highly sensitive to atmospheric temperature changes occurring during the same period over the central parts of the Arctic Ocean. A comparison of the sea ice conditions and sensitivities during three different periods reveals that, due to the strong decline of sea ice concentration and sea ice thickness, sea ice area became substantially more sensitive to the same amplitude thermodynamic atmospheric forcing anomalies during 2000–2009 relative to the earlier periods. To obtain a quantitative estimate of changes that can be expected from existing atmospheric trends, adjoint sensitivities are multiplied by monthly temperature differences between 1980s and two following decades. Strongest contributions of surface atmospheric temperature differences to AREA and VOLUME changes are observed during May and September. The strongest contribution from the downward long-wave heat flux to AREA changes occurs in September and to VOLUME changes in July–August. About 62 % of the AREA decrease simulated by the model can be explained by summing all contributions to the thermodynamic atmospheric forcing. The changing sea ice state (sensitivity) is found to enhance the decline and accounts for about one third of the explained reduction. For the VOLUME decrease, the explained fraction of the decrease is only about 37 %.     

Alpine tectonic evolution of a Jurassic subduction-accretionary complex: Deformation,kinematics and Ar/Ar age constraints on the Mesozoic low-grade schists of the Circum-Rhodope Belt in the eastern Rhodope-Thrace region,Bulgaria-Greece

Deformation of the Circum-Rhodope Belt Mesozoic (Middle Triassic to earliest Lower Cretaceous) low-grade schists underneath an arc-related ophiolitic magmatic suite and associated sedimentary successions in the eastern Rhodope-Thrace region occurred as a two-episode tectonic process: (i) Late Jurassic deformation of arc to margin units resulting from the eastern Rhodope-Evros arc–Rhodope terrane continental margin collision and accretion to that margin, and (ii) Middle Eocene deformation related to the Tertiary crustal extension and final collision resulting in the closure of the Vardar ocean south of the Rhodope terrane. The first deformational event D₁ is expressed by Late Jurassic NW-N vergent fold generations and the main and subsidiary planar-linear structures. Although overprinting, these structural elements depict uniform bulk north-directed thrust kinematics and are geometrically compatible with the increments of progressive deformation that develops in same greenschist-facies metamorphic grade. It followed the Early-Middle Jurassic magmatic evolution of the eastern Rhodope-Evros arc established on the upper plate of the southward subducting Maliac-Meliata oceanic lithosphere that established the Vardar Ocean in a supra-subduction back-arc setting. This first event resulted in the thrust-related tectonic emplacement of the Mesozoic schists in a supra-crustal level onto the Rhodope continental margin. This Late Jurassic-Early Cretaceous tectonic event related to N-vergent Balkan orogeny is well-constrained by geochronological data and traced at a regional-scale within distinct units of the Carpatho-Balkan Belt. Following subduction reversal towards the north whereby the Vardar Ocean was subducted beneath the Rhodope margin by latest Cretaceous times, the low-grade schists aquired a new position in the upper plate, and hence, the Mesozoic schists are lacking the Cretaceous S-directed tectono-metamorphic episode whose effects are widespread in the underlying high-grade basement. The subduction of the remnant Vardar Ocean located behind the colliding arc since the middle Cretaceous was responsible for its ultimate closure, Early Tertiary collision with the Pelagonian block and extension in the region caused the extensional collapse related to the second deformational event D₂. This extensional episode was experienced passively by the Mesozoic schists located in the hanging wall of the extensional detachments in Eocene times. It resulted in NE-SW oriented open folds representing corrugation antiforms of the extensional detachment surfaces, brittle faulting and burial history beneath thick Eocene sediments as indicated by 42.1–39.7 Ma ⁴⁰Ar/³⁹Ar mica plateau ages obtained in the study. The results provide structural constraints for the involvement components of Jurassic paleo-subduction zone in a Late Jurassic arc-continental margin collisional history that contributed to accretion-related crustal growth of the Rhodope terrane.     

LA-ICP-MS analysis of single fluid inclusions in a quartz crystal (Madan ore district,Bulgaria)

A “long-living” crystal of barren quartz from Kroushev Dol Pb-Zn deposit (Madan district, Rhodope Mountains, Bulgaria) was studied. The semitransparent base part (the “root”) of the crystal contains abundant inclusions, predominantly along healed cracks, while the upper half or third of the crystal is clear and poor in inclusions. In order to analyze fluid inclusions in the quartz crystal, it was cut into 4 pieces across and along the c-axis and doubly-polished sections were prepared. Fluid inclusions trapped in this quartz supply information about the temporal evolution of paleofluids depositing ore minerals.     

Analytical Approach to Calculate Quasi-thermal Noise Spectrum in Irregular Plasma

Basic properties of quasi-thermal noise spectrum in irregular plasma have been investigated using an analytical point of view. A simple formula for the plasma frequency splitting effect has been obtained for ionospheric conditions. A passive electric antenna, immersed in a stable plasma, detects the fluctuations of the electric potential due to thermal motion of the ambient particles. Properties of this quasi-thermal noise spectrum in homogeneous plasma are relatively well known and are effectively used for diagnostics of space plasma (Aksenov et al., 1978; Trakhtengerts and Chugunov, 1978; Kellog, 1981; Meyer-Vernet and Perche, 1989). Especially, in the Earth's ionosphere or solar wind plasma, random irregularities of electron density are always present. These irregularities may substantially change properties of these media through electromagnetic radiation and may also modify quasi-thermal noise spectrum, which depends on the effective dielectric permittivity tensor. This tensor is defined as the dielectric permittivity tensor of some imaginary `effective' regular medium in which the field of point source is the same as the mean field in the medium with random irregularities (Ryzhov, Tamoikin and Tatarskii, 1965; Ryzhov 1968). Since the correlation function of electric field fluctuations in random medium may be expressed through the effective dielectric constant tensor (Ryzhov, 1968), it may be used for direct calculation of quasi-thermal noise spectrum. In Zabotin et al. (2000), the Born approximation was used to calculate numerically the effective dielectric permittivity tensor and the modified noise spectrum while we analytically estimate herein the modified noise spectrum. This revised version was published online in July 2006 with corrections to the Cover Date.     

Eddy energy sources and mesoscale eddies in the Sea of Okhotsk

Based on eddy-permitting ocean circulation model outputs, the mesoscale variability is studied in the Sea of Okhotsk. We confirmed that the simulated circulation reproduces the main features of the general circulation in the Sea of Okhotsk. In particular, it reproduced a complex structure of the East-Sakhalin current and the pronounced seasonal variability of this current. We established that the maximum of mean kinetic energy was associated with the East-Sakhalin Current. In order to uncover causes and mechanisms of the mesoscale variability, we studied the budget of eddy kinetic energy (EKE) in the Sea of Okhotsk. Spatial distribution of the EKE showed that intensive mesoscale variability occurs along the western boundary of the Sea of Okhotsk, where the East-Sakhalin Current extends. We revealed a pronounced seasonal variability of EKE with its maximum intensity in winter and its minimum intensity in summer. Analysis of EKE sources and rates of energy conversion revealed a leading role of time-varying (turbulent) wind stress in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk in winter and spring. We established that a contribution of baroclinic instability predominates over that of barotropic instability in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk. To demonstrate the mechanism of baroclinic instability, the simulated circulation was considered along the western boundary of the Sea of Okhotsk from January to April 2005. In April, the mesoscale anticyclonic eddies are observed along the western boundary of the Sea of Okhotsk. The role of the sea ice cover in the intensification of the mesoscale variability in the Sea of Okhotsk was discussed.