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%.     

班公湖—怒江构造带西段三叠纪—侏罗纪构造—沉积演化

班公湖－怒江构造带西段在大地构造位置上处于特提斯构造域东端，横跨班公湖－怒江断裂带。三叠纪－株罗纪期间，其构造－沉积演化经历了大陆初始裂谷（T）、原洋裂谷（J1）、残余弧后盆地（J2－J3）阶段。初始裂谷阶段的拉张是呈南断北超的半地堑式由东向西进行的，逐渐形成地堑式原洋裂谷盆地。中晚侏罗世，南部新特提斯洋壳开始北各俯冲,产生的区域挤压应力使原洋裂谷逐渐封闭，裂谷盆地的小洋壳表现出以南向俯冲为主的双向式腑冲，同时伴生区域热沉降，盆地具残余弧后盆地的性质。该阶段,羌南地区发育碳酸盐岩为主的稳定陆缘沉积，冈度斯－念青唐古拉板片北部则形成广泛南超的近源碎屑沉积。     

Shallow seismicity, stress distribution and crustal deformation pattern in the Andaman-West Sunda arc and Andaman Sea, northeastern Indian Ocean

Shallow seismicity and available source mechanisms in the Andaman–westSunda arc and Andaman sea region suggest distinct variation in stressdistribution pattern both along and across the arc in the overriding plate.Seismotectonic regionalisation indicates that the region could be dividedinto eight broad seismogenic sources of relatively homogeneousdeformation. Crustal deformation rates have been determined for each oneof these sources based on the summation of moment tensors. The analysisshowed that the entire fore arc region is dominated by compressive stresseswith compression in a mean direction of N23^°, and the rates ofseismic deformation velocities in this belt decrease northward from 5.2± 0.65 mm/yr near Nias island off Sumatra and 1.12 ±0.13 mm/yr near Great Nicobar islands to as much as 0.4 ±0.04 mm/yr north of 8^°N along Andaman–Nicobar islandsregion. The deformation velocities indicate, extension of 0.83 ±0.05 mm/yr along N343^° and compression of 0.19 ±0.01 mm/yr along N73^° in the Andaman back arc spreadingregion, extension of 0.18 ± 0.01 mm/yr along N125^° andcompression of 0.16 ± 0.01 mm/yr along N35^° in NicobarDeep and west Andaman fault zone, compression of 0.84 ±0.12 mm/yr N341^° and extension of 0.77 ± 0.11 mm/yralong N72^° within the transverse tectonic zone in the Andamantrench, N-S compression of 3.19 ± 0.29 mm/yr and an E-Wextension of 1.24 ± 0.11 mm/yr in the Semangko fault zone ofnorth Sumatra. The vertical deformation suggests crustal thinning in theAndaman sea and crustal thickening in the fore arc and Semangko faultzones. The apparent stresses calculated for all major events range between0.1–10 bars and the values increase with increasing seismic moment.However, the apparent stress estimates neither indicate any significantvariation with faulting type nor display any variation across the arc, incontrast to the general observation that the fore arc thrust events showhigher stress levels in the shallow subduction zones. It is inferred that theoblique plate convergence, partial subduction of 90^°E Ridge innorth below the Andaman trench and the active back arc spreading are themain contributing factors for the observed stress field within the overridingplate in this region.     

Enrichissement en uranium authigène dans les sédiments glaciaires de l'océan AustralEnrichments in authigenic uranium in glacial sediments of the Southern Ocean

Four sediment cores from the Polar frontal zone and the Antarctic zone in the Indian sector of the Southern Ocean present an increase of authigenic uranium during glacial periods. We show that this increase in uranium is due to a combination of (i) an increase in the lateral transport of organic matter, (ii) a decrease in the oxygen in deep waters, and (iii) a process of diagenesis. It appears that uranium concentration cannot be used as a proxy of palaeoproductivity in the Southern Ocean, as previously suggested by Kumar et al. in 1995. To cite this article: L. Dezileau et al., C. R. Geoscience 334 (2002) 1039–1046.     

从尼玛地区地质新资料看中特提斯洋的构造演化

班公湖—怒江结合带南缘存在一套中、晚侏罗世稳定浅海碎屑岩沉积，属残余海盆地沉积，表明尼玛地区的俯冲消减机制在中侏罗世以后已结束。结合带南侧三叠系确哈拉群为一套半深水—深水沉积，是陆架边缘沉积序列，代表结合带打开之初的较早期沉积。确哈拉群之上不整合覆盖了一套中侏罗世钙碱性岛弧火山岩系，是班公湖—怒江结合带在早侏罗世向南俯冲对应的滞后弧火山岩。在结合带南侧80～100km范围内分布着一条东西长超过100km的中晚侏罗世后碰撞强过铝花岗岩带，属班公湖—怒江结合带向南俯冲碰撞作用的后碰撞阶段产物。综合认为，尼玛地区中特提斯洋是在三叠纪打开、中侏罗世以前向南俯冲闭合的。结合区域上该结合带闭合时间有早有晚、俯冲方向有南有北的事实，提出中特提斯是一个具有众多互不相通的、时代早晚各不相同的小洋盆共同组成的多岛洋，其间存在许多大小不一、运动方向和性质各不相同的地体。不同时期、不同方向的弧—弧碰撞、弧—陆碰撞造山(造陆)机制是解释中特提斯洋发展演化诸多问题的理想模式。     

New structural and petrologic data on Mesozoic schists in the Rhodope (Bulgaria): geodynamic implications

Structural analysis of low-grade rocks highlights the allochthonous character of Mesozoic schists in southeastern Rhodope, Bulgaria. The deformation can be related to the Late Jurassic–Early Cretaceous thrusting and Tertiary detachment faulting. Petrologic and geochemical data show a volcanic arc origin of the greenschists and basaltic rocks. These results are interpreted as representing an island arc-accretionary complex related to the southward subduction of the Meliata–Maliac Ocean under the supra-subduction back-arc Vardar ocean/island arc system. This arc-trench system collided with the Rhodope in Late Jurassic times. To cite this article: N.G. Bonev, G.M. Stampfli, C. R. Geoscience 335 (2003).