The evolution of lithospheric mantle beneath the Kalahari Craton and its margins

The compositional structure and thermal state of the subcontinental lithospheric mantle (SCLM) beneath the Kalahari Craton and the surrounding mobile belts have been mapped in space and time using >3400 garnet xenocrysts from >50 kimberlites intruded over the period 520–80 Ma. The trace-element patterns of many garnets reflect the metasomatic refertilisation of originally highly depleted harzburgites and lherzolites, and much of the lateral and vertical heterogeneity observed in the SCLM within the craton is the product of such metasomatism. The most depleted, and possibly least modified, SCLM was sampled beneath the Limpopo Belt by early Paleozoic kimberlites; the SCLM beneath other parts of the craton may represent similar material modified by metasomatism during Phanerozoic time. In the SW part of the craton, the SCLM sampled by “Group 2” kimberlites (>110 Ma) is thicker, cooler and less metasomatised than that sampled by “Group 1” kimberlites (mostly ≤95 Ma) in the same area. Therefore, the extensively studied xenolith suite from the Group 1 kimberlites probably is not representative of primary Archean SCLM compositions. The relatively fertile SCLM beneath the mobile belts surrounding the craton is interpreted as largely Archean SCLM, metasomatised and mixed with younger material during Paleoproterozoic to Mesoproterozoic rifting and compression. This implies that at least some of the observed secular evolution in SCLM composition worldwide may reflect the reworking of Archean SCLM. There are strong correlations between mantle composition and the lateral variations in seismic velocity shown by detailed tomographic studies. Areas of relatively low Vp within the craton largely reflect the progressive refertilisation of the Archean root during episodes of intraplate magmatism, including the Bushveld (2 Ga) and Karroo (ca. 180 Ma) events; areas of high Vp map out the distribution of relatively less metasomatised Archean SCLM. The relatively low Vp of the SCLM beneath the mobile belts around the craton is consistent with its fertile composition. The seismic data may be used to map the lateral extent of different types of SCLM, taking into account the small lateral variations in the geotherm identified using the techniques described here.     

Seismic tomography of Yunnan region using short-period surface wave phase velocity

Introduction The Yunnan region is located on the east margin of the collision zone between the Indian and Eurasian plates; it belongs to the south section of the N-S Seismic Belt of China and is the junc-ture of the Yangtze metaplatform, Songpan-Garz?fold system, Sanjiang fold system and South China fold system. This region has complex tectonic movements, crisscross faults and frequently occurring strong earthquakes, and hence it is one of the regions with the strongest earthquake ac-tivi…     

Comparison of Large-Eddy Simulation Data with Spatially Averaged Measurements Obtained by Acoustic Tomography – Presuppositions and First Results

An attempt is made to compare results oflarge-eddy simulation (LES) in a convective boundarylayer using the model PALM with experimental data obtained from acoustic travel time tomography.This method provides two-dimensional data arrays, which are considered as more suitable forLES-validation than classical local orline-integrated measurements, because the tomographic data are area- or volume-averaged.For a quantitative comparison with experimental data in general, some prerequisites have to be considered: First of all, the initial and boundary conditions of the LES model have to be provided correctly by the experiment. Considering measurement errors, a sensitivity study was performed to investigate the influence of inaccurate initial and boundary conditions on the simulation results.This showed that for determining some boundary conditions, such as the surface temperature and the roughness length, high measurement accuracies are necessary, which are difficult to reach or which at least require considerable extra measurement efforts.The initial and boundary conditions provided by the Lindenberg experiment in 1999 turned out to be of insufficient accuracy to allow quantitative comparisons.However, a qualitative comparison was performed instead to investigate if the acoustic tomography method is a proper method for comparisons with LES models in general.It showed a good qualitative agreement with some quantitative differences. These differences can partly be explained by the sensitivity of the LES to initial and boundary conditions and by the limitations of the acoustic tomography.     

Crustal heterogeneity and seismotectonics of the region around Beijing, China

A detailed three-dimensional (3-D) P-wave velocity model of the crust and uppermost mantle under the Chinese capital (Beijing) region is determined with a spatial resolution of 25 km in the horizontal direction and 4–17 km in depth. We used 48,750 precise P-wave arrival times from 2973 events of local crustal earthquakes, controlled seismic explosions and quarry blasts. These events were recorded by a new digital seismic network consisting of 101 seismic stations equipped with high-sensitivity seismometers. The data are analyzed by using a 3-D seismic tomography method. Our tomographic model provides new insights into the geological structure and tectonics of the region, such as the lithological variations and large fault zones across the major geological terranes like the North China Basin, the Taihangshan and the Yanshan mountainous areas. The velocity images of the upper crust reflect well the surface geological and topographic features. In the North China Basin, the depression and uplift areas are imaged as slow and fast velocities, respectively. The Taihangshan and Yanshan mountainous regions are generally imaged as broad high-velocity zones, while the Quaternary intermountain basins show up as small low-velocity anomalies. Velocity changes are visible across some of the large fault zones. Large crustal earthquakes, such as the 1976 Tangshan earthquake (M=7.8) and the 1679 Sanhe earthquake (M=8.0), generally occurred in high-velocity areas in the upper to middle crust. In the lower crust to the uppermost mantle under the source zones of the large earthquakes, however, low-velocity and high-conductivity anomalies exist, which are considered to be associated with fluids. The fluids in the lower crust may cause the weakening of the seismogenic layer in the upper and middle crust and thus contribute to the initiation of the large crustal earthquakes.     

Imaging low-velocity anomalies with the aid of seismic tomography

Theoretical considerations (Snell's law) suggest that low-velocity fanomalies are undersampled and therefore should be poorly resolved by inversion schemes based on ray-tracing methods. A synthetic study considering a 40×20 m low-velocity anomaly (300 m/s) placed at the center of a 400×50 m block with gradient background velocity model (from 3000 m/s at the surface to 4000 m/s at the base) indicates that the low ray density in ray-tracing coverage diagrams of tomographic inversion studies can be used as evidence for the existence of low-velocity anomalies. Combined normal incidence seismic reflection images and the velocity models obtained by tomographic inversions of first-arrival travel times form an efficient scheme to resolve low-velocity anomalies such as fracture zones. Furthermore, the velocity models derived from tomographic inversions are used in a wave equation datuming algorithm to account for statics caused by a strongly laterally variable shallow surface (weathering) layer and provide seismic reflection images of fracture zones (low-velocity anomaly) within a granitic pluton.     

自电位层析成像的理论与实验研究

为研究自电位异常场的物理机制，本文根据Patella半无限空间自电位层析成像的理论方法，对有限边界条件下的图像重建进行了研究，根据四边形界面处自电位的测量数据，重建出异常电荷的二维概率分布.文中对图像重建进行了数值模拟，分析了改进成像精度的措施，并对受压岩石的破裂过程开展了实验模拟研究，成像结果表明，实验中发现的岩石破裂的裂隙同正电荷出现的最大概率区域的位置基本一致，反映了裂隙尖端处正负电荷分离的微观过程.     

Seismic tomography imaging of an unstable embankment

Seismic tomography imaging was employed to make a diagnosis and choose a remedy for an embankment supported by a retaining wall showing clear evidence of structural instability. The geometry and structural characteristics of the wall, the inside geometry of cracks, and the physical parameters of the underground materials were the primary objectives of the geophysical survey. Seismic data were acquired along two vertical sections each one delimited by the lines of the sources within an inclined borehole and of the receivers lying on the ground. For each section, a total of 744 travel times were inverted to obtain compressional-wave velocities on a regular rectangular grid of squared cells (1 m × 1 m) using an inversion algorithm which is based on the simultaneous iterative reconstruction technique (SIRT). Ultrasonic tests carried out in laboratory on intact specimens, together with other supplemental site information, improved tomography resolution, allowing global and node constraints that forced the solution of the inverse problem to match known boundary values. Tomography imaged two high-velocity zones separated by a curved-shape low-velocity zone. The first ones were found to be compatible with the concrete retaining wall and with the schist materials, ranging from highly weathered to intact schist, while the low-velocity zone was interpreted as filling materials and/or completely decomposed rock schist. The combination of seismic tomography and laboratory measurements allowed extrapolation of important parameters over a large volume of rock mass, otherwise only representative of small rock samples near boreholes. It also enabled engineering characterization of subsurface rock mass, providing useful and accurate information to design a remedy for the embankment.     

Natural gradient experiment on transport of jet fuel derived hydrocarbons in an unconfined sandy aquifer

This paper deals with a field experiment, combining the push–pull and tracer tests, conducted under natural gradient conditions at the international Oslo airport. The studied aquifer, showing very complex hydrogeological settings, has been contaminated by a jet fuel spill. The tracer solutes—bromide, toluene, o-xylene, 1,2,4-trimethylbenzene, 1,3,5-trimethylbenzene and naphthalene, have been injected into the plume. Their migration and changes in concentration of the electron acceptors and metabolic by-products have been monitored. Fast removal of both the non-reactive tracer as well as the aromatic organics has been observed. The tracer pulse could only be detected 2 m downgradient from the injection points. At this point, toluene and o-xylene have been completely removed, however, trimethylbenzenes and naphthalene have been detected. Their depletion, based on calculations of available electron acceptors, can, to a large extent, be accounted for intrinsic biodegradation, with Fe(III) and sulphate reductions as the major controlling processes.     

Rayleigh wave tomography in the North Atlantic: high resolution images of the Iceland, Azores and Eifel mantle plumes

Presented in this paper is a high resolution S_v-wave velocity and azimuthal anisotropy model for the upper mantle beneath the North Atlantic and surrounding region derived from the analysis of about 9000 fundamental and higher-mode Rayleigh waveforms. Much of the dataset comes from global and national digital seismic networks, but to improve the path coverage a number of instruments at coastal sites in northwest Europe, Iceland and eastern Greenland was deployed by us and a number of collaborators. The dense path coverage, the wide azimuthal distribution and the substantial higher-mode content of the dataset, as well as the relatively short path-lengths in the dataset have enabled us to build an upper mantle model with a horizontal resolution of a few hundred kilometers extending to 400 km depth. Low upper mantle velocities exist beneath three major hotspots: Iceland, the Azores and Eifel. The best depth resolution in the model occurs in NW Europe and in this area low S_v-velocities in the vicinity of the Eifel hotspot extend to about 400 km depth. Major negative velocity anomalies exist in the North Atlantic upper mantle beneath both Iceland and the Azores hotspots. Both anomalies are, above 200 km depth, 4–7% slow with respect to PREM and elongated along the mid-Atlantic Ridge. Low velocities extend to the south of Iceland beneath the Reykjanes Ridge where other geophysical and geochemical observations indicate the presence of hot plume material. The low velocities also extend beneath the Kolbeinsey Ridge north of Iceland, where there is also supporting geochemical evidence for the presence of hot plume material. The low-velocity upper mantle beneath the Kolbeinsey Ridge may also be associated with a plume beneath Jan Mayen. The anomaly associated with the Azores extends from about 25°N to 45°N along the ridge axis, which is in agreement with the area influenced by the Azores Plume, predicted from geophysical and geochemical observations. Compared to the anomaly associated with Iceland, the Azores anomaly is elongated further along the ridge, is shallower and decays more rapidly with depth. The fast propagation direction of horizontally propagating S_v-waves in the Atlantic south of Iceland correlates well with the east–west ridge-spreading direction at all depths and changes to a direction close to NS in the vicinity of Iceland.