Emplacement ages and sources of kimberlites and related rocks in southern Africa: U–Pb ages and Sr–Nd isotopes of groundmass perovskite

Groundmass perovskite has been dated by LA-ICPMS in 135 kimberlites and related rocks from 110 localities across southern Africa. Sr and/or Nd isotopes have been analysed by LA-MC-ICPMS in a subset of these and integrated with published data. The age distribution shows peaks at 1,600–1,800, 1,000–1,200, 500–800 and 50–130 Ma. The major “bloom” of Group I kimberlites at ca 90 ± 10 Ma was preceded by a slow build-up in magmatic activity from ca 180 Ma. The main pulse of Group II kimberlites at 120–130 Ma was a distinct episode within this build-up. Comparison of the isotopic data with seismic tomography images suggests that metasomatized subcontinental lithospheric mantle (SCLM) with very low ε _Nd and high ⁸⁷Sr/⁸⁶Sr, (the isotopic signature of Group II kimberlites) was focused in low-Vs zones along translithospheric structures. Such metasomatized zones existed as early as 1,800 Ma, but were only sporadically tapped until the magmatic build-up began at ca 180 Ma, and contributed little to the kimberlitic magmas after ca 110 Ma. We suggest that these metasomatized volumes resided in the deep SCLM and that their low-melting point components were “burned off” by rising temperatures, presumably during an asthenospheric upwelling that led to SCLM thinning and a rise in the ambient geotherm between 120 and 90 Ma. The younger Group I kimberlites therefore rarely interacted with such SCLM, but had improved access to shallower volumes of differently metasomatized, ancient SCLM with low ⁸⁷Sr/⁸⁶Sr and intermediate ε _Nd (0–5). The kimberlite compositions therefore reflect the evolution of the SCLM of southern Africa, with metasomatic-enrichment events from as early as 1.8 Ga, through a major thermal and compositional change at ca 110 Ma, and the major kimberlite “bloom” around 90 Ma.     

Efficient block preconditioners for the coupled equations of pressure and deformation in highly discontinuous media

Large‐scale simulations of flow in deformable porous media require efficient iterative methods for solving the involved systems of linear algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps in material properties, which is often the case in geological applications, such as basin evolution at regional scales. The success of iterative methods for this type of problems depends strongly on finding effective preconditioners. This paper investigates how the block‐structured matrix system arising from single‐phase flow in elastic porous media should be preconditioned, in particular for highly discontinuous permeability and significant jumps in elastic properties. The most promising preconditioner combines algebraic multigrid with a Schur complement‐based exact block decomposition. The paper compares numerous block preconditioners with the aim of providing guidelines on how to formulate efficient preconditioners. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulating stochastically excited oscillations: The mode lifetime of ξ Hya

The discovery of solar-like oscillations in the giant star ξ Hya (G7 III) was reported by Frandsen et al. (2002). Their frequency analysis was very limited due to alias problems in the data set (caused by single-site observations). The extent to which the aliasing affected their analysis was unclear due to the unknown damping time of the stellar oscillation modes. In this paper we describe a simulator created to generate time series of stochastically excited oscillations, which takes as input an arbitrary window function and includes both white and non-white noise. We also outline a new method to compare a large number of simulated time series with an observed time series to determine the damping time, amplitude, and limited information on the degree of the stochastically excited modes. For ξ Hya we find the most likely amplitude to be ∼ 2 m s⁻¹, in good agreement with theory (Houdek and Gough, 2002), and the most likely damping time to be ∼ 2 days, which is much shorter than the theoretical value of 15–20 days calculated by Houdek and Gough (2002).     

The influence of central black holes on gravitational lenses

Recent observations indicate that many if not all galaxies host massive central black holes. In this paper we explore the influence of black holes on the lensing properties. We model the lens as an isothermal ellipsoid with a finite core radius plus a central black hole. We show that the presence of the black hole substantially changes the critical curves and caustics. If the black hole mass is above a critical value, then it will completely suppress the central images for all source positions. Realistic central black holes are likely to have masses below this critical value. Even in such subcritical cases, the black hole can suppress the central image when the source is inside a zone of influence, which depends on the core radius and black hole mass. In the subcritical cases, an additional image may be created by the black hole in some regions, which for some radio lenses may be detectable with high-resolution and large dynamic range VLBI maps. The presence of central black holes should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.     

Quaternary volcanic ash zones on the Iceland Plateau,southern Norwegian Sea

Four Quaternary volcanic ash zones in the southern Norwegian Sea have been investigated in core P57-7 from the Iceland Plateau. Both the geochemical composition and morphological variation of each ash layer have been studied. The four volcanic ash zones appear in the light oxygen isotope stages 1, 5, 7 and 11. The ash zones are composed of transparent platy grains, light brown transparent grains, brown blocky and black blocky grains and white/transparent pumice, and each zone shows a distinct stratigraphic evolution. The geochemical results show a mixture of basaltic and silicic grains in each ash zone, and that each zone contains grains from more than one eruption. The geochemical investigations strongly suggest that all the ashes are derived from Iceland. The youngest ash zone includes two layers, which based on their geochemical composition and stratigraphic position are correlated with the ¹⁴C-dated Vedde Ash (10 600 yr BP) and Saksunarvatn Ash (9 100 yr BP). Possible sources on Iceland for these layers are discussed.     

An Analytical Solution for Mechanical Responses Induced by Temperature and Air Pressure in a Lined Rock Cavern for Underground Compressed Air Energy Storage

Mechanical responses induced by temperature and air pressure significantly affect the stability and durability of underground compressed air energy storage (CAES) in a lined rock cavern. An analytical solution for evaluating such responses is, thus, proposed in this paper. The lined cavern of interest consists of three layers, namely, a sealing layer, a concrete lining and the host rock. Governing equations for cavern temperature and air pressure, which involve heat transfer between the air and surrounding layers, are established first. Then, Laplace transform and superposition principle are applied to obtain the temperature around the lined cavern and the air pressure during the operational period. Afterwards, a thermo-elastic axisymmetrical model is used to analytically determine the stress and displacement variations induced by temperature and air pressure. The developments of temperature, displacement and stress during a typical operational cycle are discussed on the basis of the proposed approach. The approach is subsequently verified with a coupled compressed air and thermo-mechanical numerical simulation and by a previous study on temperature. Finally, the influence of temperature on total stress and displacement and the impact of the heat transfer coefficient are discussed. This paper shows that the temperature sharply fluctuates only on the sealing layer and the concrete lining. The resulting tensile hoop stresses on the sealing layer and concrete lining are considerably large in comparison with the initial air pressure. Moreover, temperature has a non-negligible effect on the lined cavern for underground compressed air storage. Meanwhile, temperature has a greater effect on hoop and longitudinal stress than on radial stress and displacement. In addition, the heat transfer coefficient affects the cavern stress to a higher degree than the displacement.     

Spalling Experiments on Large Hard Rock Specimens

Specimens of coarse-grained Äspö diorite were axially compressed to observe stress-induced spalling. The specimens had a novel design characterized by two manufactured large radius notches on opposite sides. The tangential stress occurring in the notches aimed to represent the tangential loading around a circular opening. Fracture stages were monitored by acoustic emission measurements. Rock chips were formed similar to those found in situ, which indicates a similar fracture process. Slabs were cut out from the specimens and impregnated using a fluorescent material to visualize the cracks. The cracks were subsequently examined by the naked eye and by means of microscopy images, from which fracture paths could be identified and related to different minerals and their crystallographic orientations. The microscopy analyses showed how the stress field and the microstructure interact. Parallel cracks were formed 2–4 mm below the surface, sub-parallel to the direction of the maximum principal stress. The crack initiation, the roles of minerals such as feldspar, biotite and quartz and their grain boundaries and crystallographic directions are thoroughly studied and discussed in this paper. Scale effects, which relate to the stress gradient and microstructure, are discussed.     

Strenulagus (Mammalia: Lagomorpha) from the Middle Eocene Irdin Manha Formation of the Erlian Basin, Nei Mongol, China

New material of Strenulagus solaris, a stem lagomorph, from the lower part of the Middle Eocene Irdin Manha Formation at the Irdin Manha and Huheboerhe localities, Erlian Basin, Nei Mongol, China includes remains of isolated cheek teeth, fragmentary upper incisors(d I2) and postcranial elements, reported from the species for the first time. The tarsal bones display advanced features, such as a wide calcaneofibular articulation surface and, along with an increase in hypsodonty, the new data document the early stage of lagomorph evolution in the Middle Eocene of the Erlian Basin. This finding is stratigraphically well constrained and probably predates that of the type material of S. solaris from the Khaychin Formation, Bugin-Tsav Basin, Mongolia.     

Geostatistics as a groundwater exploration planning tool: case of a brackish-saline aquifer

This paper discusses a geostatistical approach to model a groundwater aquifer in 3-D. The study aims at utilizing geostatistics as a tool for characterizing zones of better-water quality in a brackish-saline aquifer. In particular, the geostatistical model was constructed to characterize the aquifer’s salinity, represented by total dissolved solids (TDS), using logs of porosity and resistivity. Quality-checked estimated TDS vertical profiles were employed to construct and model horizontal and vertical semivariograms. Parameters of semivariogram models were used to develop both the kriging plan and the generated model. Results of this modeling process are shown in the form of horizontal salinity distribution maps. The aquifer was sliced into 20 layers, each 20 m thick, to represent its overall thickness. Salinity layers maps reflect vertical stratification of TDS concentrations in the aquifer and show that water quality deteriorates with depth and toward the northern part of the aquifer. Relatively better-quality water (TDS ≤10,000 mg/l) can be found at depths between 100 and 250 m below the aquifer’s top in both eastern and southeastern parts. Water in the same interval to the western and southwestern parts reflects the presence of higher TDS concentration. From a planning point of view, it is more feasible to target the eastern part of the aquifer for pumping and desalination purposes. In addition, the generated model could be utilized as an initial condition for flow simulation.