Upper-mantle structure of the Baltic Shield below the Swedish National Seismological Network (SNSN) resolved by teleseismic tomography

Upper-mantle structure under the Baltic Shield is studied using non-linear high resolution teleseismic P -phase tomography. Observed relative arrival-time residuals from 52 teleseismic earthquakes recorded by the Swedish National Seismological Network (SNSN) are inverted to delineate the structure of the upper mantle. The network consists of 47 (currently working) three-component broad-band stations located in an area about 450 km wide and 1450 km long. In order to reduce complications due to possible significant three-dimensionality of Earth structure, events chosen for this study lay close to in-line with the long-axis of the array  (±30°)  . Results indicate P -wave velocity perturbations of ±3 per cent down to at least 470 km below the network. The size of the array allows inversion for structures even at greater depths, and lateral variations of velocity at depths of up to 680 km appear to be resolved. Below the central part of the array (60°–64° N), where ray coverage is best, the data reveals a large region of relatively low velocity at depths of over about 300 km. At depths less than about 250–300 km, the models include a number of features, including an apparent slab-like structure dipping gently towards the north.     

2D magnetic resonance tomography applied to karstic conduit imaging

Karstic conduits play a crucial role for water supply in many parts of the world. However, the imaging of such targets is generally a difficult task for most geophysical methods. Magnetic Resonance Sounding (MRS) is a geophysical method designed for imaging of water bearing structures. Initially, MRS was developed for characterizing horizontally stratified aquifers. However, when applying a 1D MRS measuring setup to the imaging of 2D–3D targets, the size of which may be much smaller than the loop, the accuracy and the lateral resolution may not be sufficient. We have studied the possibility of simultaneously processing several MRS aligned along a profile to perform a Magnetic Resonance Tomography (MRT). This work emphasizes the gain of resolution for 2D–3D imagery of MRT versus the interpolation of 1D inversion results of MRS along the same profile. Numerical modelling results show that the MRT response is sensitive to the size and location of the 2D target in the subsurface. Sensitivity studies reveal that by using the coincident transmitting/receiving (TX/RX) setup and shifting the loop around the anomaly area, the depth, section and position of a single karstic conduit with a size smaller than the MRS loop size can be resolved. The accuracy of the results depends on the noise level and signal level, the latter parameter being linked to the depth and volume of the karstic conduit and the water content in the limestone matrix. It was shown that when applying MRT to the localization of 2D anomalies such as karstic conduits, the inclination of the geomagnetic field, the orientation of the MRT profile and the angle of crossover of the conduit by the MRT profile must be taken into account. Otherwise additional errors in interpretation should be expected. A 2D inversion scheme was developed and tested. Both numerical and experimental results confirm the efficiency of the developed approach.     

基于线弹性位错模型反演1997年西藏玛尼Mw7.5级地震的干涉雷达同震形变场——Ⅱ滑动分布反演

1997年11月8日西藏M_w7.5级玛尼地震是干涉雷达技术应用于地震观测以来的一次重要事件.在第一部分中，我们应用广泛使用的Okada线弹性位错模型，假设断层的各个分段滑动量均匀，反演得到断层各个分段的几何参数和均匀滑动量.本部分的反演进一步去除滑动均匀假设，并利用更能反映断层真实状态的角形元位错模型（线弹性），在第一部分反演得到断层几何的基础上，反演断层面的静态位错分布.反演结果表明，线弹性滑动分布模型能够更好地解释观测数据，进一步提高反演的数据拟合程度.最终得到了断层面上的走滑和倾滑位错分布.首次得到的断层面滑动分布显示断层面滑动在浅部（0～12 km）比较集中，地震破裂长度约170 km，最大左旋走滑位移达4.8 m；反演结果还表明局部段落存在较大倾滑位移，量值达到1.9 m，这在断层模型中是不能忽略的，它可能是断层两侧形变不对称的原因之一；反演得到的标量地震矩为2.18×10²⁰ N·m，相当于矩震级M_w7.5，与Velasco等利用地震波形反演得到的结果一致.     

Crust and Upper Mantle Structure in the Caribbean Region by Group Velocity Tomography and Regionalization

An overview of the crust and upper mantle structure of Central America and the Caribbean region is presented as a result of the processing of more than 200 seismograms recorded by digital broadband stations from SSSN and GSN seismic networks. Group velocity dispersion curves are obtained in the period range from 10s to 40s by FTAN analysis of the fundamental mode of the Rayleigh waves; the error of these measurements varies from 0.06 and 0.09 km/s. From the dispersion curve, seven tomographic maps at different periods and with average spatial resolution of 500 km are obtained. Using the logical combinatorial classification techniques, eight main groups of dispersion curves are determined from the tomographic maps and eleven main regions, each one characterized by one kind of dispersion curves, are identified. The average dispersion curves obtained for each region are extended to 150s by adding data from a larger-scale tomographic study (Vdovin et al., 1999) and inverted using a nonlinear procedure. A set of models of the S-wave velocity vs. depth in the crust and upper mantle is found as a result of the inversion process. In six regions we identify a typically oceanic crust and upper mantle structure, while in the other two the models are consistent with the presence of a continental structure. Two regions, located over the major geological zones of the accretionary crust of the Caribbean region, are characterized by a peculiar crust and upper mantle structure, indicating the presence of lithospheric roots reaching, at least, about 200 km of depth.     

Seismic,petrologic, and geodetic analyses of the 1999 dome-forming eruption of Guagua Pichincha volcano,Ecuador

Guagua Pichincha, located 14 km west of Quito, Ecuador, is a stratovolcano bisected by a horseshoe-shaped caldera. In 1999, after some months of phreatic activity, Guagua Pichincha entered into an eruptive period characterized by the extrusion of several dacitic domes, vulcanian eruptions, and pyroclastic flows. We estimated the three-dimensional (3-D) P-wave velocity structure beneath Guagua Pichincha using a tomographic inversion method based on finite-difference calculations of first-arrival times. Hypocenters of volcano-tectonic (VT) earthquakes and long-period (LP) events were relocated using the 3-D P-wave velocity model. A low-velocity anomaly exists beneath the caldera and may represent an active volcanic conduit. Petrologic analysis of eruptive products indicates a magma storage region beneath the caldera, having a vertical extent of 7–8 km with the upper boundary at about sea level. This zone coincides with the source region of deeper VT earthquakes, indicating that a primary magma body exists in this region. LP swarms occurred in a cyclic pattern synchronous with ground deformation during magma extrusions. The correlation between seismicity and ground deformation suggests that both respond to pressure changes caused by the cyclic eruptive behavior of lava domes.     

星载极化SAR数据模拟

介绍了自行开发的星载极化SAR数据计算机模拟系统，此系统采用的模拟方法遵循了星载极化SAR系统实际运行的情况，按发射脉冲时间序列，分别依次计算各发射脉冲相应的回波数据，与常用的以点目标为基础逐点模拟方法相比，计算效率高，现贴近实际系统，便于模拟各种误差的影响和保持其对所有目标回波影响的一致性和相关性。给出了星载极化SAR数据模拟的实例。     

SAR图像斑点噪声抑制的本质

合成孔径雷达（SAR）图像的斑点噪声阻碍SAR的应用。在过去20多年中，提出了许多SAR图像去除斑点噪声方法，虽然每种去除斑点噪声方法都能有效去除斑点噪声，但都不同程度地损失了图像的边缘信息，本文首先以数学物理的观点描述了SAR图像斑点噪声抑制问题。其次分析了斑点噪声的统计特性和几种典型的去除斑点噪声方法并进行了简单的比较，在最后研究了SAR图像斑点噪声抑制的本质，这对开发既能有效去除斑点噪声，又能有效保持边缘信息的去除噪声方法十分重要。     

极化SAR干涉测量模拟研究

极化SAR干涉测量将极化测量技术和干涉测量技术巧妙地结合，能够提高干涉测量精度，尤其是在植被覆盖的情况下，不仅能够测量地表地形，而且使得估计植被高度成为可能。因此，近几年来国际上对极化SAR干涉测量技术的研究越来越重视，极化SAR干涉测量已成为SAR应用技术发展的重要方向。通过模拟技术深入研究极化SAR干涉测量理论和方法，并提出了一种新的极化SAR干涉测量改进算法，模拟实验结果表明，该算法可以进一步提高测量精度。     

一种层析成像的正演和反演方法

作者在本文中提出了一种层析成像的正演和反演方法，首先对求解区作网络划分，由给定的正演速度模型在网格点上双线性插值，获得整个网格节点的速度分布，然后，对波阵面作一近似，在此基础上，用最小走时最短路径射线追踪法完成正演，反演则是通过正演的方法实现的，理论模型算例表明，本方法效果较好。     

星载SAR复数图像的配准

星载干涉雷达测量已经被证实用于生成大范围（乃至全球）的DEM或监测地表位移具有广阔的应用前景，其中SAR复数图像序列的精确配准是保证生成有效干涉相位图的重要环节之一。本文分析了卫星InSAR系统中所涉及的几种坐标系统变换问题，并提出了基于此的几何配准方案，即基于卫星轨道状态矢量、SAR成像几何进行配准的算法；同时，还设计了基于几何配准结果和能量影像相关技术的配准算法，最后使用ERS－1／2、JERS－1和RADARSAT卫星获取的3种SAR复数图像对分别做几何配准和相关配准试验和比较，并对卫星SAR图像质量和轨道数据质量作了定性分析。