The 2007 Bengkulu earthquake, its rupture model and implications for seismic hazard

The 12 September 2007 great Bengkulu earthquake (M _w 8.4) occurred on the west coast of Sumatra about 130 km SW of Bengkulu. The earthquake was followed by two strong aftershocks of M _w 7.9 and 7.0. We estimate coseismic offsets due to the mainshock, derived from near-field Global Positioning System (GPS) measurements from nine continuous SuGAr sites operated by the California Institute of Technology (Caltech) group. Using a forward modelling approach, we estimated slip distribution on the causative rupture of the 2007 Bengkulu earthquake and found two patches of large slip, one located north of the mainshock epicenter and the other, under the Pagai Islands. Both patches of large slip on the rupture occurred under the island belt and shallow water. Thus, despite its great magnitude, this earthquake did not generate a major tsunami. Further, we suggest that the occurrence of great earthquakes in the subduction zone on either side of the Siberut Island region, might have led to the increase in static stress in the region, where the last great earthquake occurred in 1797 and where there is evidence of strain accumulation.     

Abrasion-derived sediments under intensified winds at the latest Pleistocene leading edge of the advancing Sinai-Negev erg

Quaternary desert loess and sandstone-loessite relationships in the geological record raise questions regarding causes and mechanisms of silt formation and accretion. In the northern Sinai-Negev desert carbonate terrain, only sand abrasion in active erg could have produced the large quantities of quartzo-feldspathic silts constituting the late Quaternary northwestern Negev loess. In the continuum of source (medium to fine sand of dunes) to sink (silts in loess) the very fine sand is unaccounted for in the record. This weakens the sand abrasion model of silt formation as a global process. Here, we demonstrate that, as predicted by experiments, abrasion by advancing dunes generated large quantities of very fine sand (60-110 μm) deposited within the dune field and in close proximity downwind. This very fine sand was generated 13-11 ka, possibly synchronous with the Younger Dryas under gusty sand/dust storms in the southeastern Mediterranean and specifically in the northern Sinai-Negev erg. These very fine sands were washed down slope and filled small basins blocked by the advancing dunes; outside these sampling basins it is difficult to identify these sands as a distinct product. We conclude that ergs are mega-grinders of sand into very fine sand and silt under windy Quaternary and ancient aeolian desert environments.     

Arsenic mobility in fluvial environment of the Ganga Plain,northern India

In the northern part of the Indian sub-continent, the Gomati River (a tributary of the Ganga River) was selected to study the dynamics of Arsenic (As) mobilization in fluvial environment of the Ganga Plain. It is a 900-km-long, groundwater-fed, low-gradient, alluvial river characterized by monsoon-controlled peaked discharge. Thirty-six water samples were collected from the river and its tributaries at low discharge during winter and summer seasons and were analysed by ICP-MS. Dissolved As and Fe concentrations were found in the range of 1.29–9.62 and 47.84–431.92 μg/L, respectively. Arsenic concentration in the Gomati River water has been detected higher than in its tributaries water and characteristically increases in downstream, attributed to the downstream increasing of Fe₂O₃ content, sedimentary organic carbon and silt-clay content in the river sediments. Significant correlation of determination (r ² = 0.68) was also observed between As and Fe concentrations in the river water. Arsenic concentrations in the river water are likely to follow the seasonal temperature variation and reach the level of World Health Organization’s permissible limit (10 μg/L) for drinking water in summer season. The Gomati River longitudinally develops reducing conditions after the monsoon season that mobilize As into the river water. First, dissolved As enters into pore-water of the river bed sediments by the reductive dissolution of Fe-oxides/hydroxides due to microbial degradation of sedimentary organic matter. Thereafter, it moves upward as well as down slope into the river water column. Anthropogenically induced biogeochemical processes and tropical climatic condition have been considered the responsible factors that favour the release of As in the fluvial environment of the Ganga Plain. The present study can be considered as an environmental alarm for future as groundwater resources of the Ganga–Brahmaputra Delta are seriously affecting the human–environment relationship at present.     

Assessment of tunnel instability—a numerical approach

This paper outlines the application of numerical modeling to predict deformation and stability of tunnel to be excavated in Bansagar, M.P., India. To meet the ever-increasing demand of transportation, energy, and other infrastructure projects, a large volume of rock tunneling is being carried out throughout the world. The geotechnical properties along the route of the 1,800-m long tunnel in the Bansagar region of India have been studied. The rock mass rating and rock mass quality systems were employed for empirical rock mass quality determination. Numerical analysis for the stress–strain distribution of the tunnel excavation and support systems was also carried out. In order to simulate the excavation of tunnel (NATM) at a depth of 150 m below the ground , a series of finite element analyses using Mohr-coulomb elasto-plastic constitutive model has been carried out using PLAXIS 2D. The stability of tunnel has been analyzed, and stress pattern have been discussed.     

Retrograde zoning in garnets of Elat-Wadi Magrish metamorphic rocks

Garnets of the progressively metamorphosed Elat-Wadi Watir metasediments show an inverse zonation; MnO concentration is relatively higher at the margins than in the center of the crystals while MgO shows the opposite trend. The zoning is probably caused by a retrograde metamorphism during which the garnet was partly replaced by chlorite. The control of the temperature from which garnets began to cool on zoning effectiveness is discussed.     

Isotopic and REE studies of lunar basalt 12038: implications for petrogenesis of aluminous mare basalts

We report Sr, Nd, and Sm isotopic studies of lunar basalt 12038, one of the so-called aluminous mare basalts. A precise internal Rb-Sr isochron yields a crystallization age of 3.35±0.09 AE and initial⁸⁷Sr/⁸⁶Sr=0.69922?2 (2σ error limits, 1AE=10⁹ years, λ(⁸⁷Rb)=0.0139AE^?1). An internal Sm-Nd isochron yields an age of 3.28±0.23AE and initial¹⁴³Nd/¹⁴⁴Nd=0.50764?28. Present-day¹⁴³Nd/¹⁴⁴Nd is less than the “chondritic” value, i.e. ?(Nd, 0)=?2.3±0.4 where ?(Nd) is the deviation of¹⁴³Nd/¹⁴⁴Nd from chondritic evolution, expressed as parts in 10⁴. At the time of crystallization ?(Nd, 3.2AE)=1.5±0.6.We have successfully modeled the evolution of the Sr and Nd isotopic compositions and the REE abundances within the framework of our earlier model for Apollo 12 olivine-pigeonite and ilmenite basalts. The isotopic and trace element features of 12038 can be modeled as produced by partial melting of a cumulate mantle source which crystallized from a lunar magma ocean with a chondrite-normalized REE pattern of constant negative slope. Chondrite-normalized La/Yb=2.2 for this hypothetical magma ocean pattern. A plot of I(Sr) versus ?(Nd) for the Apollo 12 basalts clearly shows the influence of varying proportions of olivine, clinopyroxene, orthopyroxene, and plagioclase in the basalt source regions. A small percentage of plagioclase (～5%) in the 12038 source apparently is responsible for low I(Sr) and ?(Nd) in this basalt. Aluminous mare basalts from Mare Crisium (Luna 24) and by inference Mare Fecunditatis (Luna 16) occupy locations on the I(Sr)-?(Nd) plot similar to that of 12038, implying that some basalts from three widely separated lunar regions came from plagioclase-bearing source regions. A summary of model calculations for mare basalts shows a record of lunar mantle solidification during the period when REE abundances in the lunar magma ocean increased from ～20× chondritic to >100× chondritic. Although there is a general trend from olivine to clinopyroxene-dominated source regions with progressive magma ocean evolution, significant mineralogical heterogeneities in mantle composition apparently formed at any given stage of evolution, as evidenced in particular by the three Apollo 12 magma types.     

Surface ruptures induced by the devastating 1068 AD earthquake in the southern Arava valley, Dead Sea Rift, Israel

The Elat fault (a segment of the Dead Sea Transform) runs along the southern Arava valley (part of the Dead Sea Rift, Israel) forming a complex fault zone that displays a time-dependent seismic behaviour. Paleoseismic evidence shows that this fault zone has generated at least 15 earthquakes of magnitude larger than M 6 during the late Pleistocene and the Holocene. However, at present the Elat fault is one of the quietest segments of the Dead Sea Transform, lacking even microsesimicity. The last event detected in the southern Arava valley occurred in the Avrona playa and was strong enough to have deformed the playa and to change it from a closed basin with internal drainage into an open basin draining to the south.Paleoseismological, geophysical and archaeological evidences indicate that this event was the historical devastating earthquake, which occurred in 1068 AD in the eastern Mediterranean region. According to the present study this event was strong enough to rupture the surface, reactivate at least two fault branches of the Elat fault and vertically displace the surface and an early Islamic irrigation system by at least 1 m. In addition, the playa area was uplifted between 2.5 and 3 m along the eastern part of the Elat fault shear zone. Such values are compatible with an earthquake magnitude ranging between M 6.6 and 7. Since the average recurrence interval of strong earthquakes during the Holocene along the Elat fault is about 1.2 ± 0.3 ky and the last earthquake occurred more about 1000 years ago, the possibility of a very strong earthquake in this area in the future should be seriously considered in assessing seismic hazards.     

Depth Estimation from Gravity Data Using the Maximum Entropy Method (MEM) and the Multi Taper Method (MTM)

The depth determination from the gravity data in frequency domain is carried out using the classical fast Fourier transform (FFT) method utilizing scaling properties of ensemble of anomalous source. The problem of calculating power spectrum from the FFT is well described in the literature. Here, the application of other high-resolution methods of power spectrum calculation, such as maximum entropy method (MEM) and multi-taper method (MTM) are explored to estimate depth to anomalous sources. At the outset, the FFT, the MEM and the MTM are tested on synthetic gravity data, generated for different types of synthetic models and then all these methods are applied to the field gravity data of the Bengal basin. The MTM with scaling is found to be superior for providing the detailed subsurface information rather than the MEM and the FFT methods in the case of synthetic as well as field examples.     

Detection of model voids by identifying reverberation phenomena in GPR records

In order to improve the reliability of the ground penetrating radar (GPR) method in identifying subsurface sinkholes and karst cavities, laboratory investigations have been performed. The main objective of this work was to examine the relationship between horizontal/vertical voids dimensions and wavelengths of various antennas, and the corresponding GPR responses. Emphasis was given to the investigation of the factors that cause the appearance of reverberation phenomena in the signal pattern.The tests were conducted in 5 m × 10 m area by 2-m-deep trench filled with homogenous, dry sand. The voids models (empty fiberglass cylinders in diameters of 0.6 m, 1.0 m, 1.5 m and 2.4 m, with various heights) were buried vertically with their tops at depths of between 0.7 and 1.5 m. Investigations were performed for the various model conditions by towing 500, 300 and 100 MHz antennas along a pre-established grid, for the various model conditions.The GPR data collected using the 500 MHz bistatic antenna above the 1.0-m- and the 1.5-m-diameter cylinders, and using the 300 MHz bistatic antenna above the 1.5-m-diameter cylinder, confirmed the presence of a reverberation phenomenon, i.e. a strong convex signal pattern, containing a series of high amplitude extending oscillations with reduced frequency.Based on past practical GPR experience of void detection and presently obtained experimental data, two rules of thumbs may be adopted for the prediction of the appearance of resonant radar pictures:

1. The void diameter larger than the wavelength in air of the antenna used.

2. The vertical size of the empty void not significantly smaller than its horizontal dimension.

The strong reverberations generated by the inner surface of the void targets were found to approximate standing waves generated in cylindrical waveguides and waveguide resonators. The theoretical, experimental and practical results obtained concur.     

Complex exposure histories of chert clasts in the late Pleistocene shorelines of Lake Lisan,southern Israel

Activities of ²⁶Al and ¹⁰Be in five chert clasts sampled from two beach ridges of late Pleistocene Lake Lisan, precursor of the Dead Sea in southern Israel, indicate low rates of chert bedrock erosion and complex exposure, burial, and by inference, transport histories. The chert clasts were derived from the Senonian Mishash Formation, a chert‐bearing chalk, which is widely exposed in the Nahal Zin drainage basin, the drainage system that supplied most of the material to the beach ridges. Simple exposure ages, assuming only exposure at the beach ridge sampling sites, range from 35 to 354 ky; using the ratio ²⁶Al/¹⁰Be, total clast histories range from 0·46 to 4·3 My, unrelated to the clasts' current position and exposure period on the late Pleistocene beach ridges, 160–177 m below sea level. Optically stimulated luminescence dating of fine sediments from the same and nearby beach ridges yielded ages of 20·0 ± 1·4 ka and 36·1 ± 3·3 ka. These ages are supported by the degree of soil development on the beach ridges and correspond well with previously determined ages of Lake Lisan, which suggest that the lake reached its highest stand around 27 000 cal. years BP . If the clasts were exposed only once and than buried beyond the range of significant cosmogenic nuclide production, then the minimum initial exposure and the total burial times before delivery to the beach ridge are in the ranges 50–1300 ky and 390–3130 ky respectively. Alternatively, the initial cosmogenic dosing could have occurred during steady erosion of the source bedrock. Back calculating such rates of rock erosion suggests values between 0·4 and 12 m My^?1. The relatively long burial periods indicate extended sediment storage as colluvium on slopes and/or as alluvial deposits in river terraces. Some clasts may have been stored for long periods in abandoned Pliocene and early Pleistocene routes of Nahal Zin to the Mediterranean before being transported again back into the Nahal Zin drainage system and washed on to the shores of Lake Lisan during the late Pleistocene. Copyright © 2003 John Wiley & Sons, Ltd.