山西临汾盆地黄土剖面重金属分布特征及其影响因素

山西临汾盆地2条黄土剖面重金属分布特征及其影响因素研究表明,Pb、Cu和Zn的垂向变化特征非常相似,而As、Cd和Hg在剖面中的分布特征与前三者不同。各重金属元素在黄土剖面中的变化特征首先均受到黄土源区物质成分的控制,其次不同程度地受到古气候变化的影响。依据Pb、Cu、Zn与各古气候替代性指标(例如磁化率,w(SiO2)/w(Al2O3)等)均呈显著相关,As、Cd和Hg与古气候替代性指标的相关性较差,推断剖面中Cu、Pb、Zn的分布特征受古气候变化的影响,其原因可能主要是这3个金属元素易以类质同象形式赋存在磁性矿物中;而As、Cd和Hg在剖面中的分布特征与古气候无关。但各重金属元素与源区指示元素的相关性表明其在剖面中的分布特征均受黄土源区物质成分的影响和控制。     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

大别山黄石-六安反射地震剖面新的地质解释

在大别山黄石至六安反射地震剖面上有很多近水平的反射体,表明大别山中心部位的中、下地壳内也有近水平的分层性。反映有流变性和动力学方面的差异,与大型薄皮构造理论对岩石圈性质的认识一致,因而在其运动过程中应服从薄皮构造的运动规律。结合以往对大别山区划分的岩石一构造组合,在前人对此反射地震剖面所作的地质解释的基础上,对反射地震剖面作了较为详细的地质解释,并建立了大别山造山带在此剖面上的两维几何结构。剖面南部为造山过程中形成的背斜构造。地表15km深度内为由碰撞混杂岩组成的扬子与中一朝大陆之间的主滑脱带。剖面中部为造山期后的侵入体。剖面北部为主滑脱带的根带（通常认为的缝合带）,根带被中生代形成的晓天磨子潭断裂带切割。剖面最北端为变质复理石中略晚于主滑脱带的反向冲断带。推溺l下地壳的断开距离在扬子大陆俯冲时（三叠纪前）规模较大,然后逐渐缩小,直至保留到今天的规模。     

On the use of gravity and magnetic anomalies for locating probable areas of metallic mineralization in South Sinai, Egypt

The present study deals with processing and interpreting the potential gravity and magnetic data in order to locate promising sites for metallic mineral occurrences in the basement rocks of South Sinai. Two promising sites were anticipated by combining the geophysical results and the available geological information. Each site was profiled and interpreted to the nearly approximated subsurface feature model, and the geometrical parameters of the expected subsurface structures may form mineralized bodies such as depth, thickness, width, dipping, density contrast, and magnetic susceptibility contrast could be delineated. Detailed ground Very Low Frequency (VLF) and Vertical Magnetic Gradient (VMG) methods were carried out in the site A2, eight VLF and VMG anomalies were determined, and their depths were calculated.     

The gravitational field of topographic-isostatic masses and the hypothesis of mass condensation II-the topographic-isostatic geoid

In Part I we focussed on a convergent representation of the gravitational potential generated bytopographic masses on top of the equipotential surface atMean Sea Level, thegeoid, and by those masses which compensate topography. Topographic masses have also been condensated, namely represented by a single layer. Part II extends the computation of the gravitational field of topographic-isostatic masses by a detailed analysis of itsforce field in terms ofvector-spherical harmonic functions. In addition, the discontinuous mass-condensated topographic gravitational force vector (head force) is given. Once we identify theMoho discontinuity asone interface of isostatically compensated topographical masses, we have computed the topographic potential and the gravitational potential which is generated by isostatically compensated masses atMean Sea Level, the geoid, and illustrated by various figures of geoidal undulations. In comparison to a data oriented global geoid computation ofJ. Engels (1991) the conclusion can be made that the assumption of aconstant crustal mass density, the basic condition for isostatic modeling, does not apply. Insteaddensity variations in the crust, e.g. betweenoceanic and continental crust densities, have to be introduced in order to match the global real geoid and its topographic-isostatic model. The performed analysis documents that thestandard isostatic models based upon aconstant crustal density areunreal.     

Insight into the structure of the Upper Rhine Graben and its basement from a new compilation of Bouguer Gravity

A compilation of gravity data from the Upper Rhine Graben (URG) is presented that includes all the main data sources from its German and French parts. This data is used to show that the URG consists of, at least, two arc-shaped and asymmetric rift units that tectonically are the basic building blocks of the graben. In this sense the URG does not differ from other continental rifts, such as the African rifts. This division should replace the now classical geomorphologic division of the URG into three segments, based on their different trends. Moreover, the gravity suggests that the faults in the central and southern segments are continuous and have the same trend, appearing to respond as a single kinematic unit. Changes in the gravity field in the graben are shown to reflect not only the structure of the graben, but also the highly variable composition of the basement. In this respect, the URG is quite different from some other Tertiary continental rifts, where possible changes in the composition of the basement are mostly masked in the gravity field by the effect of the overlying low-density sediments. This characteristic is used to study the extent of some of the main basement units that underlie the graben.     

Airborne gamma-ray spectrometry to quantify chemical erosion processes

Airborne gamma-ray spectrometry data (uranium, potassium and thorium contents) reveal geochemical heterogeneities within the monolithological Hyrôme watershed (ca. 150 km²) in the Armorican massif (western France). Our observations and computations provide important constraints on the spatial distribution and the associated magnitudes of chemical erosion processes at the scale of a small watershed. Two distinct, partially preserved, weathering profiles exhibit a strong correlation between regolith evolution and airborne-derived K/Th ratios, suggesting that the variability is linked to supergene processes. Using both airborne data and laboratory measurements on rock samples, the total net export of potassium has been estimated at 422 ± 50 kg/m² and the chemical weathering rate of potassium at 17 ± 2 kg/km²/a.     

Retrieval of stratospheric O3 and NO2 vertical profiles using zenith scattered light observations

Daily zenith scattered light intensity observations were carried out in the morning twilight hours using home-made UV-visible spectrometer over the tropical station Pune (18‡31′, 73‡51′) for the years 2000–2003. These observations are obtained in the spectral range 462–498 nm for the solar zenith angles (SZAs) varying from 87‡ to 91.5‡. An algorithm has been developed to retrieve vertical profiles of ozone (O₃) and nitrogen dioxide (NO₂) from ground-based measurements using the Chahine iteration method. This retrieval method has been checked using measured and recalculated slant column densities (SCDs) and they are found to be well matching. O₃ and NO₂ vertical profiles have been retrieved using a set of their air mass factors (AMFs) and SCDs measured over a range of 87–91.5‡ SZA during the morning. The vertical profiles obtained by this method are compared with Umkehr profiles and ozonesondes and they are found to be in good agreement. The bulk of the column density is found near layer 20–25 km. Daily total column densities (TCDs) of O₃ and NO₂ along with their stratospheric and tropospheric counterparts are derived using their vertical profiles for the period 2000–2003. The total column, stratospheric column and tropospheric column amounts of both trace gases are found to be maximum in summer and minimum in the winter season. Increasing trend is found in column density of NO₂ in stratospheric, tropospheric and surface layers, but no trend is observed in O₃ columns for above layers during the period 2000–2003     

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.