基于球面重力反演苏拉威西地区莫霍面结构

基于球坐标系的地球物理反演能有效避免地球曲率的影响, 适用于大尺度构造研究.本文基于重力异常数据在球坐标系下反演莫霍面深度, 结合数据误差及光滑正则化项建立反演目标函数并求解, 同时将该方法应用于苏拉威西地区.苏拉威西地区具有复杂的断裂系统、年轻的俯冲带, 是研究俯冲起始机制等前沿科学问题的理想场所.目前研究区的地球物理观测尚不充分, 缺乏对莫霍面形成有效约束的地震数据, 对研究区莫霍面的整体认识较少.本文基于卫星重力观测数据, 通过匹配滤波方法提取与研究区莫霍面结构相关的重力异常, 并结合频谱分析确定该地区的莫霍面深度参考值.在反演中, 通过两次随机子抽样交叉验证选择最优的超参数, 包括正则化因子、莫霍面密度差以及参考莫霍面深度, 迭代反演获得莫霍面深度.反演结果表明: 研究区莫霍面平均深度为20.0 km, 深度变化范围为9.2~33.3 km.总体上, 海区莫霍面浅, 约10.0~20.0 km, 陆区莫霍面深, 约25.0~33.0 km, 该结果与Crust1.0全球模型、前人重力反演结果以及地震数据基本相符, 总体上反映了苏拉威西地区的莫霍面变化特征.

     

Numerical modeling of lock-exchange gravity/turbidity currents by a high-order upwinding combined compact difference scheme

This study presents two-dimensional direct numerical simulations for sediment-laden current with higher density propagating forward through a lighter ambient water.The incompressible NavierStokes equations including the buoyancy force for the density difference between the light and heavy fluids are solved by a finite difference scheme based on a structured mesh.The concentration transport equations are used to explore such rich transport phenomena as gravity and turbidity currents.Within the framework of an Upwinding Combined Compact finite Difference(UCCD)scheme,rigorous determination of weighting coefficients underlies the modified equation analysis and the minimization of the numerical modified wavenumber.This sixth-order UCCD scheme is implemented in a four-point grid stencil to approximate advection and diffusion terms in the concentration transport equations and the first-order derivative terms in the Navier-Stokes equations,which can greatly enhance convective stability and increase dispersive accuracy at the same time.The initial discontinuous concentration field is smoothed by solving a newly proposed Heaviside function to prevent numerical instabilities and unreasonable concentration values.A two-step projection method is then applied to obtain the velocity field.The numerical algorithm shows a satisfying ability to capture the generation,development,and dissipation of the Kelvin-Helmholz instabilities and turbulent billows at the interface between the current and the ambient fluid.The simulation results also are compared with the data in published literatures and good agreements are found to prove that the present numerical model can well reproduce the propagation,particle deposition,and mixing processes of lock-exchange gravity and turbidity currents.     

Systematic Methodology for Velocity-dependent Gravity Modelling of Density Crustal Cross-Sections, using an Optimization Procedure

The paper describes a new method of gravity modelling of the density cross section along a DSS profile. We propose a systematic procedure of gravity interpretation using structural information contained both in the velocity distribution of the seismic model and the geometry of its layers. The procedure is based on the optimisation of the density parameters of the model. The values of the parameters are limited by relationships between the seismic velocity and the density (optimisation with inequalities). The formulation of the method allows for discussion of the presence of isostasy and the distribution of compensating masses. The discussion is very important for the considered crust-lithosphere system. The presented calculations are done for one Carpathian profile (CEL01). The analysis of the properties of the fitting procedure provides the characteristic depth of subcrustal compensating density structures, the position of crucial tectonic boundaries along the profiles, locations of the crustal areas having exceptional lithology, the regional density differences of subcrustal matter and the estimation of regional deviations from the isostatic equilibrium.     

An Upper Pliocene coarse pumice breccia generated by a shallow submarine explosive eruption,Milos, Greece

The Filakopi Pumice Breccia (FPB) is a very well exposed, Pliocene volcaniclastic unit on Milos, Greece, and has a minimum bulk volume of 1 km³. It consists of three main units: (A) basal lithic breccia (4–8 m) mainly composed of angular to subangular, andesitic and dacitic clasts up to 2.6 m in diameter; (B) very thickly bedded, poorly sorted pumice breccia (16–17 m); and (C) very thick, reversely graded, grain-supported, coarse pumice breccia (6.5–20 m), at the top. The depositional setting is well constrained as shallow marine (up to a few hundred metres) by overlying fossiliferous and bioturbated mudstone. This large volume of fine pumice clasts is interpreted to be the product of an explosive eruption from a submarine vent because: (1) pumice clasts are the dominant component; (2) the coarse pumice clasts (>64 mm) have complete quenched margins; (3) very large (>1 m) pumice clasts are common; (4) overall, the formation shows good hydraulic sorting; and (5) a significant volume of ash was deposited together with the coarsest pyroclasts.The bed forms in units A and B suggest deposition from lithic-rich and pumiceous, respectively, submarine gravity currents. In unit C, the coarse (up to 6.5 m) pumice clasts are set in matrix that grades upwards from diffusely stratified, fine (1–2 cm) pumice clasts at the base to laminated shard rich mud at the top. The coarse pumice clasts in unit C were settled from suspension and the framework was progressively infilled by fine pumice clasts from waning traction currents and then by water-settled ash. The FPB displays important features of the products of submarine explosive eruptions that result from the ambient fluid being seawater, rather than volcanic gas or air. In particular, submarine pyroclastic deposits are characterised by the presence of very coarse juvenile pumice clasts, pumice clasts with complete quenched rims, and good hydraulic sorting.Electronic Supplementary Material Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.Editorial responsibility: J. Donelly-Nolan     

Long term variations in the mesospheric mean flow observed at Grahamstown (South Africa) and Adelaide (Australia)

The seasonal and interannual behaviour of monthly mean winds at a height of 90 km recorded at Grahamstown (33.3°S, 26.5°E) and Adelaide (34.5°S, 138.5°E) between 1987 and 1994 are compared. The zonal wind is found to be consistently stronger at Grahamstown and is always eastward, whereas at Adelaide it sometimes reverses. Maxima tend to occur near the solstices, the primary maximum during summer at Grahamstown, in agreement with satellite results, and during winter at Adelaide. The meridional wind also tends to be stronger at Grahamstown, but at both stations is predominantly northward with a maximum in summer and generally not as strong as the zonal component. This seasonal behaviour is reasonably well understood in terms of the interaction of the mean flow with gravity waves propagating up from below, with coriolis forces also playing an important role in the case of the meridional wind. Satellite observations do not generally support the idea that longitudinal differences between the stations could be attributed to the presence of a tropospheric/stratospheric stationary wave. It is suggested that these differences are more probably associated with local effects. Interannual zonal wind patterns at the two sites are similar over the summer months but are less well correlated during the rest of the year. The underlying causes of this variability are not well understood but are most probably global in nature, at least during the summer.     

一种卫星重力梯度数据的ARMA最优滤波模型构建方法

由于卫星重力梯度观测的有色噪声特性和海量观测特征,在利用直接法进行重力场模型的最小二乘求解时,观测值的协方差阵为超大型的非对角阵,这给数值求解带来了极大困难.本文提出了一种基于先验误差功率谱密度的最优ARMA滤波模型构建方法,结合法方程的分块求解策略,可实现对卫星重力梯度观测值的高效滤波处理.数值仿真结果表明,利用最优ARMA滤波器进行时域滤波后,法方程的态性得到了明显改善,重力梯度观测值中的有色噪声得到了有效的"白化"处理,大地水准面精度得到了显著提升.     

Evaluation of the accuracy of the EIGEN-1S and -2 CHAMP-derived gravity field models by satellite crossover altimetry

Since the advent of CHAMP, the first in a series of low-altitude satellites being almost continuously and precisely tracked by GPS, a new generation of long-wavelength gravitational geopotential models can be derived. The accuracy evaluation of these models depends to a large extent on the comparison with external data of comparable quality. Here, two CHAMP-derived models, EIGEN-1S and EIGEN-2, are tested with independent long-term-averaged single satellite crossover (SSC) sea heights from three altimetric satellites (ERS-1, ERS-2 and Geosat). The analyses show that long-term averages of crossover residuals still are powerful data to test CHAMP gravity field models. The new models are tested in the spatial domain with the aid of ERS-1/-2 and Geosat SSCs, and in the spectral domain with latitude-lumped coefficient (LLC) corrections derived from the SSCs. The LLC corrections allow a representation of the satellite-orbit-specific error spectra per order of the models spherical harmonic coefficients. These observed LLC corrections are compared to the LLC projections from the models variance–covariance matrix. The excessively large LLC errors at order 2 found in the case of EIGEN-2 with the ERS data are discussed. The degree-dependent scaling factors for the variance-covariance matrices of EIGEN-1S and –2, applied to obtain more realistic error estimates of the solved-for coefficients, are compatible with the results found here.     

Geophysical Expression of Low Sulphidation Epithermal Au-Ag Deposits and Exploration Implications –Examples from the Hokusatsu Region of SW Kyushu, Japan–

Abstract: Low sulphidation epithermal Au-Ag deposits constitute a major exploration target for many mining companies worldwide. These deposits are commonly characterized by quartz vein systems which formed at temperatures often in the range 150–250°C giving rise to major alteration halos. Accompanying mineralization and associated alteration is the destruction and redistribution of magnetite in the host rocks and the formation of various K-bearing minerals such as adularia, illite and alunite. These alteration halos can often cover an area of at least several square kilometres and as such have the potential to be detected during reconnaissance geological surveys utilizing high resolution airborne aeromagnetic and radiometric sensors.
Analysis of airborne magnetic and radiometric data recently acquired over the Hokusatsu Region of SW Kyushu indicates that major low sulphidation epithermal Au-Ag deposits are characterized by several distinct geophysical features varying with the degree of post-mineralization cover and level of erosion. Several deposits coincide with extensive magnetic "flat" and high radiometric K–counts. In addition, several workers have shown that a very close association exists between Bouguer gravity highs and epithermal Au mineralization in the Hokusatsu Region. The combined use of areomagnetics, radiometrics and gravity data therefore provides an exceptionally effective approach to exploration for low sulphidation epithermal Au–Ag deposits in the Hokusatsu Region of Japan.     

吉林东部重力异常的分形特征及其地质意义

分形是研究复杂地质地球物理问题的手段之一。本文利用分数布朗噪声模型,采用能谱法,计算了吉林东部重力异常的分维数,分维数等值图为地质体的划分提供了新的证据。     

陕渝黔桂1800 km超长探测剖面重力异常场特征及深部地壳结构探榷

通过对跨越华北克拉通、秦岭造山带、扬子克拉通几大地质构造单元的,北起陕西榆林经秦岭过重庆鱼泉、贵州贵阳并向南到广西凭祥全长1810 km超长重力探测剖面的数据进行处理分析和解释,构建了沿剖面的二维地壳密度结构模型,并详细分析了沿剖面壳内各界面与Moho界面展布的深部结构和构造特征、构划出了该剖面的深部断裂分布,探讨了剖面辖区跨越的克拉通、造山带、接触带或耦合带等一系列的区域构造的差异,同时对其可能的地质构造含义进行初步了解释,以期能对深化认识该剖面跨越地区的地壳结构、构造单元划分及动力学等研究,提供相关重力场的依据.