3-D magnetotelluric inversion and model validation with gravity data for the investigation of flood basalts and associated volcanic rifted margins

20 magnetotelluric (MT) soundings were collected on the Isle of Skye, Scotland to provide a high-resolution three-dimensional (3-D) electrical resistivity model of a volcanic province within the framework of a project jointly interpreting gravity, seismic, geological and MT data. The full 3-D inversion of the MT data jointly interpreted with gravity data reveals upper crustal structure. The main features of the model are interpreted in conjunction with previous geological mapping and borehole data. Our model extends to 13 km depth, several kilometres below the top of the Lewisian basement. The top of the Lewisian basement is at approximately 7–8 km depth and the topography of its surface was controlled by Precambrian rifting, during which a 4.5 km thick sequence of Torridonian sediments was deposited. The Mesozoic sediments above, which can reach up to 2.2 km thick, have small-scale depocentres and are covered by up to 600 m of Tertiary lava flows. The interpretation of the resistivity model shows that 3-D MT inversion is an appropriate tool to image sedimentary structures beneath extrusive basalt units, where conventional seismic reflection methods may fail.     

Numerical inversion of deformation caused by pressure sources: application to Mount Etna (Italy)

The interpretation of geodetic data in volcanic areas is usually based on analytical deformation models. Although numerical finite element (FE) modelling allows realistic features such as topography and crustal heterogeneities to be included, the technique is not computationally convenient for solving inverse problems using classical methods. In this paper, we develop a general tool to perform inversions of geodetic data by means of 3-D FE models. The forward model is a library of numerical displacement solutions, where each entry of the library is the surface displacement due to a single stress component applied to an element of the grid. The final solution is a weighted combination of the six stress components applied to a single element-source. The pre-computed forward models are implemented in a global search algorithm, followed by an appraisal of the sampled solutions. After providing extended testing, we apply the method to model the 1993–1997 inflation phase at Mt Etna, documented by GPS and EDM measurements. We consider four different forward libraries, computed in models characterized by homogeneous/heterogeneous medium and flat/topographic free surface. Our results suggest that the elastic heterogeneities of the medium can significantly alter the position of the inferred source, while the topography has minor effect.     

Effects of galvanic distortion on magnetotelluric data over a three-dimensional regional structure

Although the galvanic distortion due to local, near-surface inhomogeneities is frequency-independent, its effect on the magnetotelluric data becomes, in a 3-D structure, frequency-dependent. Therefore, both the apparent resistivity and the phase responses are disturbed, and a correction should be carried out prior to the 3-D interpretation in order to retrieve the 3-D regional impedance tensor. In many cases, the structure is 2-D for depths corresponding to a first range of periods and 3-D for longer periods (called 2-D/3-D). For these cases, a simple method which allows us to retrieve the 3-D regional impedance tensor (except the static shift) is presented. The method proposed uses the Groom & Bailey decomposition of the distortion matrix for the short periods. Three examples are presented: two using synthetic data and one employing real data. These examples show the effect of the galvanic distortion over a regional 2-D/3-D model and the retrieval of the regional transfer functions from the distorted ones.     

基于遥感数据的元江县植被/土地利用类型研究

元江县山区面积占全县总面积的96．8％，垂直高差达2253m，地形复杂，植被类型存在明显的垂直分布。只用一个解译体系来解译这一有较大差异地区的影像是有困难的。因此本研究为提高影像解译的准确性，在对2002年2月9日覆盖元江全景的TM卫星影像进行数据校正、影像融合、影像分层、计算机自动监督分类的基础上，利用全球卫星定位系统(GPS)，结合野外踏察及地形图、土地利用相关信息对机读结果进行目视检验和GIS分析校正，最终划分出18类植被类型、9类土地利用类型。判读结果经检验，土地利用的判读精度为84．2％，植被类型的判读精度为80．6％。     

Non-linear crustal corrections in high-resolution regional waveform seismic tomography

We compare 3-D upper mantle anisotropic structures beneath the North American continent obtained using standard and improved crustal corrections in the framework of Non-linear Asymptotic Coupling Theory (NACT) applied to long period three component fundamental and higher mode surface waveform data. Our improved approach to correct for crustal structure in high-resolution regional waveform tomographic models goes beyond the linear perturbation approximation, and is therefore more accurate in accounting for large variations in Moho topography within short distances as observed, for instance, at ocean–continent margins. This improved methodology decomposes the shallow-layer correction into a linear and non-linear part and makes use of 1-D sensitivity kernels defined according to local tectonic structure, both for the forward computation and for the computation of sensitivity kernels for inversion. The comparison of the 3-D upper mantle anisotropic structures derived using the standard and improved crustal correction approaches shows that the model norm is not strongly affected. However, significant variations are observed in the retrieved 3-D perturbations. The largest differences in the velocity models are present below 250 km depth and not in the uppermost mantle, as would be expected. We suggest that inaccurate crustal corrections preferentially map into the least constrained part of the model and therefore accurate corrections for shallow-layer structure are essential to improve our knowledge of parts of the upper mantle where our data have the smallest sensitivity.     

Reconstructing pre-erosion topography using spatial interpolation techniques:A validation-based approach

Understanding the topographic context preceding the development of erosive landforms is of major relevance in geomorphic research, as topography is an important factor on both water and mass movement-related erosion, and knowledge of the original surface is a condition for quantifying the volume of eroded material. Although any reconstruction implies assuming that the resulting surface reflects the original topography, past works have been dominated by linear interpolation methods, incapable of generating curved surfaces in areas with no data or values outside the range of variation of inputs. In spite of these limitations, impossibility of validation has led to the assumption of surface representativity never being challenged. In this paper, a validation-based method is applied in order to define the optimal interpolation technique for reconstructing pre-erosion topography in a given study area. In spite of the absence of the original surface, different techniques can be nonetheless evaluated by quantifying their capacity to reproduce known topography in unincised locations within the same geomorphic contexts of existing erosive landforms. A linear method (Triangulated Irregular Network, TIN) and 23 parameterizations of three distinct Spline interpolation techniques were compared using 50 test areas in a context of research on large gully dynamics in the South of Portugal. Results show that almost all Spline methods produced smaller errors than the TIN, and that the latter produced a mean absolute error 61.4% higher than the best Spline method, clearly establishing both the better adjustment of Splines to the geomorphic context considered and the limitations of linear approaches. The proposed method can easily be applied to different interpolation techniques and topographic contexts, enabling better calculations of eroded volumes and denudation rates as well as the investigation of controls by antecedent topographic form over erosive processes.     

Ground deformation in a viscoelastic medium composed of a layer overlying a half-space: a comparison between point and extended sources

We obtain and compare analytical and numerical solutions for ground displacement caused by an overpressurized magma chamber placed in a linear viscoelastic medium composed of a layer over a half-space. Different parameters such as size, depth and shape of the chamber, crustal rheology and topography are considered and discussed. Numerical solutions for an axisymmetric extended source are computed using a finite element method (FEM). Analytical solutions for a point source are obtained using the dislocation theory and the propagator matrix technique. In both cases, the elastic solutions are used together with the correspondence principle of linear viscoelasticity to obtain the solution in the Laplace transform domain. Viscoelastic solutions in the time domain are derived inverting the Laplace transform using the Prony series method. The differences between the results allow us to constrain the applicability of the point source and the flat surface hypothesis, which are usually implicitly assumed when analytical solutions are derived. The effect of the topography is also considered. The results obtained show that neglecting the topographic effects may, in some cases, introduce an error greater than that implicit in the point-source hypothesis. Therefore, for an adequate modelling and interpretation of the time-dependent displacements, topography must be considered.     

基于DEM的山区遥感图像地形校正方法

基于地表反射为各向异性的假定,可以将遥感图像的地形校正分为光照和反射率校正两部分。基于6S大气校正模型和数字高程模型(DEM)对山区遥感图像进行光照和反射率校正,通过消除遥感图像中大气和地形的影响,反演得到地表反射率。通过对江西兴国县TM图像的实验,实现了对山区遥感图像的地形校正。对模型中的大气参数进行了敏感性分析,发现大气参数的误差对地表反射率的反演精度影响较小。     

Seismic imaging of the laterally varying D" region beneath the Cocos Plate

We use an axisymmetric, spherical Earth finite difference algorithm to model SH -wave propagation through cross-sections of laterally varying lower mantle models beneath the Cocos Plate derived from recent data analyses. Synthetic seismograms with dominant periods as short as 4 s are computed for several models: (1) a D" reflector 264 km above the core–mantle boundary with laterally varying S -wave velocity increases of 0.9–2.6 per cent, based on localized structures from a 1-D double-array stacking method; (2) an undulating D" reflector with large topography and uniform velocity increase obtained using a 3-D migration method and (3) cross-sections through the 3-D mantle S -wave velocity tomography model TXBW. We apply double-array stacking to assess model predictions of data. Of the models explored, the S -wave tomography model TXBW displays the best overall agreement with data. The undulating reflector produces a double Scd arrival that may be useful in future studies for distinguishing between D" volumetric heterogeneity and D" discontinuity topography. Synthetics for the laterally varying models show waveform variability not observed in 1-D model predictions. It is challenging to predict 3-D structure based on localized 1-D models when lateral structural variations are on the order of a few wavelengths of the energy used, particularly for the grazing geometry of our data. Iterative approaches of computing synthetic seismograms and adjusting model characteristics by considering path integral effects are necessary to accurately model fine-scale D" structure.     

A quantitative methodology to extract regional magnetotelluric impedances and determine the dimension of the conductivity structure

A methodology to determine quantitatively the dimensionality of the dominant conducting structures and the resolution of the structural parameters in magnetotelluric data is presented. In addition, the method recovers the regional impedance responses when the regional structure can be characterized, at least approximately, as 1- or 2-D. The methodology is based upon three general models of the MT tensor, each of which has a distinct parameterization and physical interpretation. A weighted statistical residual describes quantitatively the fit of the model response to the data within the scatter of the measured data and hence permits: (1) tests of dimensionality, (2) determination of the appropriate strike angle, and (3) recovery of the regional responses.
The method has been tested extensively with synthetic data and proven to be successful. These synthetic studies give insight into the different physical parameterizations and the stability of the parameters determined. We describe and illustrate some of these synthetic studies. With field data, the methodology is not always as straightforward, but its application to a great many sites has proven valuable. Data from two closely spaced sites, which are both affected by strong but very different 3-D effects, are analysed to illustrate the geological significance of the results. The analyses reveal and recover regional responses within the data which indicate the presence of electrical anisotropy located deep in the crust and upper mantle. Analyses of the entire data set, of which these two form a part, confirm this finding.     

月表形貌特征研究进展及趋势分析

在月球探测及科学研究过程中,对月表形貌特征的认识和分析是月球探测计划的关键任务之一,有助于理解和揭示月球的形态特征及空间分异规律,对于分析月球岩石构造、估算月表年龄、反演月壤厚度、恢复月球的起源和演化历史等都具有重要意义。从月表形貌特征的塑造及表达、宏观特征的分析与定量刻画、月表撞击坑的识别与分类及空间分异特征等几方面进行了归纳和分析,提出了未来月表形貌的研究方向,包括全月球形貌特征的多级分区划分、撞击坑类型的划分、全月球撞击坑空间分异特征及规律的研究、基于撞击坑空间分布探讨月球演化、基于形貌特征的比较行星学研究等。     

Modelling of ground deformation and gravity fields using finite element method: an application to Etna volcano

Elastic finite element models are applied to investigate the effects of topography and medium heterogeneities on the surface deformation and the gravity field produced by volcanic pressure sources. Changes in the gravity field cannot be interpreted only in terms of gain of mass disregarding the ground deformation of the rocks surrounding the source. Contributions to gravity changes depend also on surface and subsurface mass redistribution driven by dilation of the volcanic source. Both ground deformation and gravity changes were firstly evaluated by solving a coupled axisymmetric problem to estimate the effects of topography and medium heterogeneities. Numerical results show significant discrepancies in the ground deformation and gravity field compared to those predicted by analytical solutions, which disregard topography, elastic heterogeneities and density subsurface structures. With this in mind, we reviewed the expected gravity changes accompanying the 1993–1997 inflation phase on Mt Etna by setting up a fully 3-D finite element model in which we used the real topography, to include the geometry, and seismic tomography, to infer the crustal heterogeneities. The inflation phase was clearly detected by different geodetic techniques (EDM, GPS, SAR and levelling data) that showed a uniform expansion of the overall volcano edifice. When the gravity data are integrated with ground deformation data and a coupled FEM modelling was solved, a mass intrusion could have occurred at depth to justify both ground deformation and gravity observations.     

遥感影像地形校正研究进展及其比较实验

地形校正作为复杂地形区遥感影像预处理的重要步骤,对提高地表参数遥感定量化精度具有重要意义。为此,在简述地形校正含义与目标的基础上,回顾并总结了国内外各种地形校正方法并将其划分为基于波段比、DEM和超球面3类方法,以期为地形校正及相关研究提供参考。在DEM的支持下,采用11种地形校正方法对ETM+影像进行了校正比较实验,研究表明:(1)VECA、b、C、Teillet-回归、SCS+C、Minnaert和Minnaert-SCS校正7种地形校正效果较好,可用于遥感影像的地形校正;而Cosine-T、Cosine-C、SCS和Cosine-b校正存在过度校正现象,不宜选择。(2)VECA与b校正模型校正效果最好,且VECA校正比b校正可操作性更强。在此基础上,从地形效应的理论诠释与数学表达、DEM数据、地形校正应用研究3个方面探讨了目前该领域存在的一些问题和难点,并对今后可能的工作重点和研究方向提出了建议。     

Effects of volcano topography on seismic broad‐band waveforms

Volcano seismology often deals with rather shallow seismic sources and seismic stations deployed in their near field. The complex stratigraphy on volcanoes and near‐field source effects have a strong impact on the seismic wavefield, complicating the interpretation techniques that are usually employed in earthquake seismology. In addition, as most volcanoes have a pronounced topography, the interference of the seismic wavefield with the stress‐free surface results in severe waveform perturbations that affect seismic interpretation methods. In this study we deal predominantly with the surface effects, but take into account the impact of a typical volcano stratigraphy as well as near‐field source effects. We derive a correction term for plane seismic waves and a plane‐free surface such that for smooth topographies the effect of the free surface can be totally removed. Seismo‐volcanic sources radiate energy in a broad frequency range with a correspondingly wide range of different Fresnel zones. A 2‐D boundary element method is employed to study how the size of the Fresnel zone is dependent on source depth, dominant wavelength and topography in order to estimate the limits of the plane wave approximation. This approximation remains valid if the dominant wavelength does not exceed twice the source depth. Further aspects of this study concern particle motion analysis to locate point sources and the influence of the stratigraphy on particle motions. Furthermore, the deployment strategy of seismic instruments on volcanoes, as well as the direct interpretation of the broad‐band waveforms in terms of pressure fluctuations in the volcanic plumbing system, are discussed.     

基于矩匹配算法的山区影像地形辐射校正方法研究

针对山区遥感影像地形辐射校正问题,应用现有模型(如C校正)进行地形辐射校正很难达到理想效果。引入矩匹配算法,利用DEM数据计算坡度、坡向等地形信息,以特定坡度和坡向数据为参考依据,对影像进行地形辐射校正。通过北京房山区SPOT5影像进行试验,表明该方法能在很大程度上消除地形阴影,更好地反映阴影区域的细节信息,同时光谱特性保真程度较好,原模糊的影像区域通过处理基本上达到有效识别地物的要求。     

Earthquake location in strongly heterogeneous media

Traveltime computation methods for strongly heterogeneous 3-D media developed during recent years are well suited for earthquake location. We present here a new method based on the traveltime algorithm of Podvin-Lecomte, related to the inverse problem formulation of Tarantola & Valette. The Podvin-Lecomte method, based on the Huygens principle, is very robust and allows arbitrary surface topography and station placement even for borehole instruments. First arrival traveltimes are computed for each of the recording stations using a fine 3-D velocity mesh (up to 10⁶ cells on a workstation). The traveltime grid allows the use of the Tarantola & Valette formulation, which enables a full non-linear approach. The solution is given as a 3-D probability density function of hypocentre coordinates, which accounts for the arrival time measurements as well as a priori information for the location, the accuracy of both the arrival time readings and the computation of the theoretical traveltimes. This powerful method called 3DGRIDLOC gives the location of the induced seismicity of the gas field of Lacq (France) using 443 520 cells of a 3-D velocity mesh and the observations from nine recording stations, one of which is located at the bottom of a 3880 m deep borehole. Location of synthetic foci as well as more than 500 actual earthquakes shows the real advantages of this new method versus the classical HYPO71. A new insight into the induced seismicity is now possible: induced seismicity may occur as far away as 10 km from the gas reservoir and involve a much greater volume of rock than expected using earlier locations.     

The onset of extension during lithospheric shortening:a two-dimensional thermomechanical model for lithospheric unrooting

We model the evolution of the lithosphere during its shortening and consequent gravitational collapse with special emphasis on the induced variations in the surface stress regime and dynamic topography. In particular, we analyse the conditions leading, immediately after lithospheric failure, to local extension, eventually coeval with compression. Different crustal rheologies and kinematic conditions as well as thermally imposed mechanical rupture are considered. Numerical calculations have been performed by using a 2-D finite element code that couples the thermal and mechanical equations for a Newtonian rheology with a temperature-dependent viscosity. The results show that, after the failure of a gravitationally unstable lithospheric root, the replacement of lithospheric mantle by warmer asthenospheric material induces a considerable variation in the dynamic topography and in the surface stress regime. The occurrence of local extension, its intensity and its spatial distribution depend mainly on whether convergence continues throughout the process or ceases after or before the lithospheric failure. Similarly, uplift/subsidence and topographic inversion are controlled by kinematic conditions and crustal rheology. Mechanical rupture produces drastic changes in the surface stress regime and dynamic topography but only for a short time period, after which the system tends to evolve like a continuous model.     

From ACH tomographic models to absolute velocity models

The ACH method, a widely used tomographic inverse method, is characterized by the use of relative residuals in order to avoid possible biases coming from outside the target volume. The ACH method thus does not really retrieve the 3-D structure of the target volume, but instead leads to velocity contrasts relative to the layer average of the velocity, this average value remaining unknown ( Aki et al. 1977 ). Two artefacts derive from this particularity: (1) velocity contrasts are known only in the horizontal direction and it is not possible, in a strict mathematical sense, to estimate the contrasts in the vertical direction with ACH alone; (2) negative anomalies are often interpreted as low velocities, whereas negative anomalies may correspond to high velocities if the average value of the corresponding layer is sufficiently high. The converse is true of positive anomalies. We show with synthetic data how these artefacts can affect the interpretation of tomographic images. We propose to correct the artefacts by reintroducing the 1-D regional average model, and show in synthetic experiments how effective this correction can be.
  The application of this procedure to data recorded in the Kunlun region shows that the retrieval of the absolute values of the 3-D velocity model is helpful for interpreting the tomographic images and better defining which features are anomalous.     

贵州省乌江流域人口分布与地形的关系

研究流域人口分布与地形的关系,有助于了解地理环境对人口分布的影响。以贵州省乌江流域为研究区,基于DEM数据和人口数据,提取海拔、坡度、起伏度等地形因子,研究人口分布与地形因子的关系。研究结论：① 人口数量在海拔800~1400 m的地区超过60%,人口密度在1000~1200 m的地区最高。② 人口数量与人口密度随着坡度上升总体呈下降趋势。③ 人口密度在起伏度小于50 m的地区超过1000人/km ²,人口密度随着起伏度上升总体呈下降趋势。     

Adaptive unstructured grid finite element simulation of two-dimensional magnetotelluric fields for arbitrary surface and seafloor topography

We present an adaptive unstructured triangular grid finite element approach for effectively simulating plane-wave diffusive electromagnetic fields in 2-D conductivity structures.
The most striking advantage of irregular grids is their potential to incorporate arbitrary geometries including surface and seafloor topography. Adaptive mesh refinement strategies using an a posteriori error estimator yield most efficient numerical solutions since meshes are only refined where required.
We demonstrate the robustness of this approach by comparison with analytical solutions and previously published numerical simulations. Maximum errors may systematically be reduced to, for example, 0.8 per cent for the apparent resistivity and 0.2° in the phase.
An additional accuracy study of the thickness of the air layer in E-polarization suggests to keep a minimum thickness depending on lateral conductivity contrasts within the earth.
Furthermore, we point out the new quality and flexibility of our simulation technique by addressing two marine magnetotelluric applications. In the first case, we discuss topographic effects associated with a synthetic sinusoidal sea bottom model and in the second case, we show a close-to-reality scenario using real bathymetry data from the East Pacific Rise at 17°S.