裂缝性岩石裂缝扩展的光弹构造物理模拟——以南华北盆地户部寨沙河阶组为例

以南华北盆地户部寨第三系沙河阶组粉砂岩低渗透储层为例,总结了其中发育的5种类型构造裂缝组合,即:①南北向雁形张节理组成的共轭剪节理带;②东西向雁形张节理组成的共轭剪节理带;③缠结状多组剪切构造裂缝;④组合状交切的多组剪切构造裂缝及张裂缝;⑤散状无序排列的构造裂缝。在此基础上,建立光弹性试验构造物理模型,沿现应力场96°方向,均匀地对该模型逐步增加外载力。结果表明,沿各组合裂缝尖端分别形成对称于裂缝走向、且主轴垂直于裂缝的Ⅰ型等色线条纹;对称于裂缝走向、且主轴平行于裂缝走向的Ⅱ型条纹和斜交于裂缝走向的复合型(Ⅲ型)条纹。第①类组合裂缝尖端以Ⅰ型等色线条纹的拉张应力和张剪应力为主;第②类组合裂缝尖端以复合型(Ⅲ型)等色线条纹的挤压应力为主;第③、④和⑤类组合裂缝尖端的应力分布较为复杂。当模型外载荷接近1000kg时,裂缝尖端起裂,并开始发育显微裂缝,而当外载荷达1100-1200 kg时,沿裂缝尖端部位形成宏观裂缝,并开始扩展。5种类型的裂缝组合扩展发育程度由强至弱依次为:①,③,⑤,②,④。其中第①类共轭裂缝组合起裂最早,扩展最快,裂缝间的连通趋势最好,说明规则分布的构造裂缝组合,在垂直于外施加载衙时更易达到连通的目的。     

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

Constraints on the application of carbon isotope thermometry in high- to ultrahigh-temperature metamorphic terranes

Nine marble horizons from the granulite facies terrane of southern India were examined in detail for stable carbon and oxygen isotopes in calcite and carbon isotopes in graphite. The marbles in Trivandrum Block show coupled lowering of δ¹³C and δ¹⁸O values in calcite and heterogeneous single crystal δ¹³C values (? 1 to ? 10‰) for graphite indicating varying carbon isotope fractionation between calcite and graphite, despite the granulite facies regional metamorphic conditions. The stable isotope patterns suggest alteration of δ¹³C and δ¹⁸O values in marbles by infiltration of low δ¹³C–δ¹⁸O‐bearing fluids, the extent of alteration being a direct function of the fluid‐rock ratio. The carbon isotope zonation preserved in graphite suggests that the graphite crystals precipitated/recrystallized in the presence of an externally derived CO₂‐rich fluid, and that the infiltration had occurred under high temperature and low f_O2 conditions during metamorphism. The onset of graphite precipitation resulted in a depletion of the carbon isotope values of the remaining fluid+calcite carbon reservoir, following a Rayleigh‐type distillation process within fluid‐rich pockets/pathways in marbles resulting in the observed zonation. The results suggest that calcite–graphite thermometry cannot be applied in marbles that are affected by external carbonic fluid infiltration. However, marble horizons in the Madurai Block, where the effect of fluid infiltration is not detected, record clear imprints of ultrahigh temperature metamorphism (800–1000 °C), with fractionations reaching <2‰. Zonation studies on graphite show a nominal rimward lowering δ¹³C on the order of 1 to 2‰. The zonation carries the imprint of fluid deficient/absent UHT metamorphism. Commonly, calculated core temperatures are > 1000 °C and would be consistent with UHT metamorphism.     

Momentum transport of wave zero during March: A possible predictor for the Indian summer monsoon

Analysis of monthly momentum transport of zonal waves at 850 hPa for the period 1979 to 1993, between ‡S and ‡N for January to April, using zonal (u) and meridional (v) components of wind taken from the ECMWF reanalysis field, shows a positive correlation (.1% level of significance) between the Indian summer monsoon rainfall (June through September) and the momentum transport of wave zero TM(0) over latitudinal belt between 25‡S and 5‡N (LB) during March. Northward (Southward) TM(0) observed in March over LB subsequently leads to a good (drought) monsoon season over India which is found to be true even when the year is marked with the El-Nino event. Similarly a strong westerly zone in the Indian Ocean during March, indicates a good monsoon season for the country, even if the year is marked with El-Nino. The study thus suggests two predictors, TM(0) over LB and the strength of westerly zone in the Indian Ocean during March.     

Timing of Fracture Formation and Duration of Mineralization at the Hishikari Deposit, Southern Kyushu, Japan

Abstract. The Hishikari epithermal gold deposit consists of the Honko-Sanjin zone and Yamada zone. On the basis of observed crosscutting relationships of veins, vein system in the Hishikari deposit was classified into the early and late veins. They differ from orientation, width, control by unconformity, the amount of displacement and structural features, but have similarity in gold grade and Au/Ag ratios. K-Ar ages are presented for the classified veins to clarify timing of fracturing and duration of mineralization.
Age determination revealed that the early veins in the Honko-Sanjin zone range from 0.86 to 1.11 Ma and that the late veins range from 0.73 to 0. 84 Ma. The economically most important mineralization in the Hishikari deposit occurred around 0.90 Ma in a very short period (about 0.05 million years).
The distribution of mineralization age in the whole Hishikari deposit, occurrences of slickensides and ground water temperatures imply that mineralization lasted longer in the southern part of the deposit. As compared the early veins with late veins, structural changes including width, strike, dip and vein characteristics occurred at around 0.85 Ma.     

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.