Computation of Large Anisotropic Seismic Heterogeneities (CLASH)

A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.     

A traveltime reciprocity discrepancy in the Podvin & Lecomte time3d finite difference algorithm

We have found that the extensively used finite difference scheme time3d produces time fields which are 'asymmetric' in the sense that traveltimes computed to the right of the source are faster than traveltimes computed to the left. All finite difference schemes will, as they are approximations to the wave equation, to some extent fail to obey reciprocity perfectly. We show, however, that the errors in time3d may be significant—and unnecessarily large. An asymmetry in the point source initialization has been identified, and after correction time3d produces time fields with an improved reciprocity.     

Electromagnetic methods to characterize the Savoia di Lucania waste dump (Southern Italy)

The aim of this work is the joint application and integration of non-invasive geoelectrical methods for studying the landfill of Savoia di Lucania (Southern Italy). This landfill for its engineering features and small dimensions (70 m × 30 m × 6 m) represents an optimal test site to assess a geophysical survey protocol for municipal solid waste landfills investigation and monitoring. The landfill of Savoia di Lucania has been built with a reinforced concrete material and coated with a high-density polyethylene (HDPE) liner. Three electrical resistivity tomographies (ERT), two self-potential (SP) map surveys and one induced polarization (IP) section have been performed, both in the surrounding area and inside the waste landfill. The geophysical investigations have well defined some buried boundaries of the landfill basin and localized the leachate accumulation zones inside the dumpsite. Comparison of our results with other engineering and geological investigations could be the key for evaluating the integrity of the HDPE liner. Finally, the joint use of the ERT, IP and SP methods seems to be a promising tool for studying and designing new monitoring systems able to perform a time-lapse analysis of waste landfill geometry and integrity.     

地面地震层析技术在金属矿勘查中的试验研究

由于金属矿区地质构造的复杂性和复杂的地表条件, 适合于探测层状介质的地震反射方法技术在金属矿勘查中的应用受到了限制。为提高地震记录的横向分辨率, 探测复杂的地质构造及隐伏岩体等, 需研究相应的地震方法技术, 地面地震层析技术在金属矿勘查中的试验研究就是为解决上述问题而开展的。本文就非线性地面地震层析技术的基本原理进行了阐述, 并就该方法在金属矿勘查中的试验效果进行了讨论。     

2004年印度洋海啸印尼亚齐省灾情遥感监测与分析

2004年12月底印尼苏门答腊岛西北海岸发生里氏9.0级地震,引发印度洋沿岸的印尼、斯里兰卡、泰国、印度等国家先后遭受强烈海啸。印度尼西亚亚齐省是遭受海啸灾害最严重的地区之一。为了评估海啸灾害造成的严重影响,应用遥感和地理信息系统技术,以2005-01-07英国DMC小卫星遥感影像为基础数据源,结合收集到的该地区灾前(2001年)LandsatETM 遥感影像,对亚齐省进行了监测和快速评估。监测结果表明距离震中位置较近的亚齐省苏门答腊岛西海岸遭受影响较严重,亚齐省苏门答腊岛西海岸由于距离震中位置较近遭受严重影响,而苏门答腊岛东海岸则影响较轻,受灾最严重的地区是亚齐省首府班达亚齐市附近。经统计,整个亚齐省遭受海啸严重影响的地区面积约为483km2,其中城市建成区受淹面积约为26·2km2,乡村(镇)地区受淹面积约为430·5km2。文章的最后还结合数字高程模型对海啸灾情的结果进行了分析。结果表明,绝大多数因海啸而受淹地区高程低于20m。     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

World-wide synthetic tide parameters for gravity and vertical and horizontal displacements

The response of the Earth’s crust to the direct effect of lunisolar gravitational forcing is known as the body tide. The body tide is superimposed by surface-loading forces due to the pressure of the periodically varying ocean tide acting on the Earth, called ocean tide loading (OTL). Both body tide and OTL can be decomposed into components of the same frequency known as tidal parameters. However, OTL is more complicated than body tides because of the dynamic effects of the ocean. Estimating OTL requires a model of the ocean tides and knowledge of the elastic properties of the solid Earth. Thus, synthetic tide parameters (amplitude factors and phase leads) have been developed here on a world-wide grid for gravity and positional displacements. The body tide contributions were added to the oceanic contribution to provide the Earth tide response. The accuracy and reliability of the synthetic tidal parameters have been estimated by comparing observed gravity and vertical-displacement tide parameters with those interpolated from our synthetic model, which shows good agreement. Tests also indicate that the synthetic tide parameters provide realistic gravimetric and displacements for practical use in tidal prediction.     

东亚南亚夏季降水对海温异常响应的数值模拟

利用IAP AGCM-Ⅱ大气环流模式，模拟了赤道中太平洋地区海温异常对东亚和南亚夏季（6～8月）降水的影响。结果表明：当1～4月该海温呈正异常时，中南半岛和我国东部降水将明显偏多，孟加拉湾附近降水偏少；而该海温呈负异常时，则我国华中和孟加拉湾附近降水偏多；对以上海温异常最敏感处为中南半岛和孟加拉湾附近地区，其中前者降水与该海温异常呈正相关，后者则呈负相关；该海温异常还会在中高纬地区激发出强迫波列。     

DMS and SO2 at Baring Head, New Zealand: Implications for the Yield of SO2 from DMS

Atmospheric dimethyl sulfide (DMS) and sulfur dioxide (SO₂) concentrations were measured at Baring Head, New Zealandduring February and March 2000. Anti-correlated DMS and SO₂ diurnalcycles, consistent with the photochemical production of SO₂ from DMS, were observed in clean southerly air off the ocean. The data is used to infer a yield of SO₂ from DMS oxidation. The estimated yields are highly dependent on assumptions about the DMS oxidation rate. Fitting the measured data in a photochemical box model using model-generated OH levels and the Hynes et al. (1986) DMS + OH rate constant suggests that theSO₂ yield is 50–100%, similar to current estimates for the tropical Pacific.However, the observed amplitude of the DMS diurnal cycle suggests that the oxidation rate is higher than that used by the model, and therefore, that theSO₂ yield is lower in the range of 20–40%.     

3D shear-wave velocity structure of the Eifel plume, Germany

The Quaternary Eifel volcanic fields, situated on the Rhenish Massif in Germany, are the focus of a major interdisciplinary project. The aim is a detailed study of the crustal and mantle structure of the intraplate volcanic fields and their deep origin. Recent results from a teleseismic P-wave tomography study reveal a deep low-velocity structure which we infer to be a plume in the upper mantle underneath the volcanic area [J.R.R. Ritter et al., Earth Planet. Sci. Lett. 186 (2001) 7-14]. Here we present a travel-time investigation of 5038 teleseismic shear-wave arrivals in the same region. First, the transverse (T) and radial (R) component travel-time residuals are treated separately to identify possible effects of seismic anisotropy. A comparison of 2044 T- and 2994 R-component residuals demonstrates that anisotropy does not cause any first-order travel-time effects. The data sets reveal a deep-seated low-velocity anomaly beneath the volcanic region, causing a delay for teleseismic shear waves of about 3 s. Using 3773 combined R- and T-component residuals, an isotropic non-linear inversion is calculated. The tomographic images reveal a prominent S-wave velocity reduction in the upper mantle underneath the Eifel region. The anomaly extends down to at least 400 km depth. The velocity contrast to the surrounding mantle is depth-dependent (from −5% at 31-100 km depth to at least −1% at 400 km depth). At about 170-240 km depth the anomaly is nearly absent. The resolution of the data is sufficient to recover the described features, however the anomaly in the lower asthenosphere is underestimated due to smearing and damping. The main anomaly is similar to the P-wave model except the latter lacks the ‘hole’ near 200 km depth, and both are consistent with an upper mantle plume structure. For plausible anhydrous plume material in the uppermost 100 km of the mantle, an excess temperature as great as 200-300 K is estimated from the seismic anomaly. However, 1% partial melt reduces the required temperature anomaly to about 100 K. The temperature anomaly associated with the deeper part of the plume (250 to about 450 km depth) is at least 70 K. However, this estimate is quite uncertain, because the amplitude of the shear-wave anomaly may be larger than the modelled one. Another possibility is water in the upwelling material. The gap at 170-240 km depth could arise from an increase of the shear modulus caused by dehydration processes which would not affect P-wave velocities as much. An interaction of temperature and compositional variations, including melt and possibly water, makes it difficult to differentiate quantitatively between the causes of the deep-seated low-velocity anomaly.