Interpretation of aeromagnetic data across the central crystalline shield area of Nigeria

Summary. A residual map of the total magnetic field (above 25 000 nT base) is presented for a portion of the central crystalline shield area of Nigeria and overlapping small portions of the Chad basin and the Benue rift (8°30'−12° 00'lat, and 7°−10°30' long). The map (based on a dataset digitized from recently released aeromagnetic sheets of Nigeria) leads to four results. (1) A magnetic boundary, evident on the map, separates the Younger Granite complexes into two groups. The groups are petrologically different, and the boundary may be a fault line with uplift to the south. (2) South of the boundary the map is dominated by a system of sub-parallel anomalies striking NE–SW, possibly representing major tectonic trends, and a set of fractures through which the Younger Granite complexes were intruded. The trend of the system parallels the Benue rift and lineaments in the oceanic crust off West Africa. (3) Negative magnetic anomalies lie over most of the known ring complexes, and over some suspected buried ring complexes and other intrusions. (4) 2½-and 3-D modelling shows that the larger complexes extend to 12 km depth, and the smaller ones to 6 km. They have nearly vertical sides, and magnetization contrasts range from 0.3 to 0.5 A m⁻¹.     

INVERSE DDA OF A SYSTEM OF BLOCKS IN SPHERE COORDINATES

1　IntroductionDiscontinuousDeformationAnalysis(DDA)developedbyShiin 1 980’s[1 ,2 ] maybeusedtocalculatethedeformationanddisplacementsinamulti blocksystemandtoanalyzetherelationbe tweenforcesanddisplacementsinthesystem .Ifthesurveydisplacementofeverypointordirectionisknown,alldisplacementsanddeformationsofblocks,eventheglide ,theclosureandthestretchoftheboundary planeofblocks ,canbecalculated .TheresultisoptimallyfittedaccordingtotheLeastSquarePrinciple.DDAismainlyusedinrockblockssystems…     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Annual variation of atmospheric carbonyl sulfide in the marine atmosphere in the Southern Indian Ocean

Measurements of the concentrations of carbonyl sulfide (COS) in the marine atmosphere were made over a period of two years in the southern Indian Ocean (Amsterdam Island, 37°50 S–77°31 E; March 1987–February 1988 and April 1989–February 1990). The mean atmospheric COS concentration for the whole period was 475±48 pptv (n=544). Atmospheric COS concentrations show no significant seasonal variation with a summer to winter ratio of 1.05. Taking into account the observed variability of the atmospheric COS concentration (10%), a value of 1.4 yr is estimated as a lower limit for the atmospheric COS lifetime. A comparison of the COS data at Amsterdam Island with those obtained in the Southern Hemisphere in the past 12 yr does not reveal any significant trend in the tropospheric background COS mixing ratio.     

Crustal and upper mantle structures of the brazilian coast

The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.     

The role of heterogeneous strain in the development and preservation of a polymetamorphic record in high-P granulites, western Canadian Shield

Mafic rocks in the Chipman domain of the Athabasca granulite terrane, western Canadian Shield, provide the first well‐documented record of two distinct high‐P granulite facies events in the same domain in this region. Textural relations and the results of petrological modelling (NCFMASHT system) of mafic granulites are interpreted in terms of a three‐stage tectonometamorphic history. Stage 1 involved development of the assemblage Grt + Cpx + Qtz ± Pl (M₁) from a primary Opx‐bearing igneous precursor at conditions of 1.3 GPa, 850–900 °C. Field and microstructural observations suggest that M₁ developed synchronously with an early S₁ gneissic fabric. Stage 2 is characterized by heterogeneous deformation (D₂) and synkinematic partial retrogression of the peak assemblage to an amphibole‐bearing assemblage (M₂). Stage 3 involved a third phase of deformation and a return to granulite facies conditions marked by the prograde breakdown of amphibole (Amph₂) to produce matrix garnet (Grt_3a) and the coronitic assemblage Cpx_3b + Opx_3b + Ilm_3b + Pl_3b (M_3b) at 1.0 GPa, 800–900 °C. M₁ and M_3b are correlated with 2.55 and 1.9 Ga metamorphic generations of zircon, respectively, which were dated in a separate study. Heterogeneous strain played a crucial role in both the development and preservation of these rare examples of multiple granulite facies events within single samples. Without this fortuitous set of circumstances, the apparent reaction history could have incorrectly led to an interpretation involving a single‐cycle high‐grade event. The detailed P–T–t–D history constructed for these rocks provides the best evidence to date that much of the east Lake Athabasca region experienced long‐term lower crustal residence from 2.55 to 1.9 Ga, and thus the region represents a rare window into the reactivation and ultimate stabilization processes of cratonic lithosphere.     

饱和软黏土动力学特性循环扭剪试验研究

以饱和软黏土的循环扭剪试验研究为基础,阐明不固结不排水条件下饱和软黏土的动力学特性。在循环扭剪和循环三轴两种不同试验应力状态下,通过研究不同围压和不同静、动应力组合下饱和软黏土的应力等效破坏关系和应变等效破坏关系,提出在循环荷载作用下饱和软黏土的循环破坏同样遵循Mises屈服准则,且在Mises屈服准则下,饱和软黏土的循环强度趋于不变量。循环破坏的过程可以等效为一种拟静力弹塑性循环蠕变,建立了饱和软黏土循环累积变形随静荷载、循环荷载以及循环破坏振次之间的关系式。上述结论分别从应力和变形两个方面阐述饱和软黏土的动力学特性,与试验应力状态和围压无关,可以推广到一般应力状态下。     

东菲律宾海深水铁锰结壳发育站位沉积物的粒度及黏土矿物学特征

大洋铁锰结壳是潜在的重要海洋资源,多年来世界各国对此已作了大量调查.传统意义上的铁锰结壳均分布于3000m水深以浅的海山基岩上.中国科学院海洋研究所在对东菲律宾海的地质调查中采集到大量褐黑色深水铁锰结壳.在300多个站位的取样中,在44个站位采集到结壳样品,可见结壳在研究区的分布较为广泛.该结壳发育于4000~5000m水深的致密至半固结沉积物的表层,厚度一般为2~3cm,最厚可达7cm.结壳为板(层)状.     

伶仃洋沉积动力特点的研究

位于珠江三角洲东侧的伶仃洋，因径流下泄与潮流进退的流向不一，使它各分流口门的出口水道都有主槽和支槽之分，即都有主干水道和分汊水道。陆架高盐海水入侵又使伶仃洋内沉积动力过程在空间分布上发生差异，如沉积物分布有粗－细－稍粗之分；而水体中的密度、速度差异，常常产生锋带，对水下地形的发展有不可忽视的影响。因而全面认识发生在伶仃洋内的沉积动力作用，对深水航道的选线极为重要。     

渤海三维海洋温度和海流数值预报

"十·五"期间,我国开展了三维海洋温度和海流数值预报的业务化研究工作.经过3年的努力,渤海三维海洋温度和海流数值预报系统研制完成,并于2003年10月,开始试预报.本文对该预报系统以及运行情况进行了介绍,并分析了所存在的问题和发展方向.