Revised age of proximal deposits in the Zagros foreland basin and implications for Cenozoic evolution of the High Zagros

The regionally extensive, coarse-grained Bakhtiyari Formation represents the youngest synorogenic fill in the Zagros foreland basin of Iran. The Bakhtiyari is present throughout the Zagros fold-thrust belt and consists of conglomerate with subordinate sandstone and marl. The formation is up to 3000 m thick and was deposited in foredeep and wedge-top depocenters flanked by fold-thrust structures. Although the Bakhtiyari concordantly overlies Miocene deposits in foreland regions, an angular unconformity above tilted Paleozoic to Miocene rocks is expressed in the hinterland (High Zagros).

The Bakhtiyari Formation has been widely considered to be a regional sheet of Pliocene–Pleistocene conglomerate deposited during and after major late Miocene–Pliocene shortening. It is further believed that rapid fold growth and Bakhtiyari deposition commenced simultaneously across the fold-thrust belt, with limited migration from hinterland (NE) to foreland (SW). Thus, the Bakhtiyari is generally interpreted as an unmistakable time indicator for shortening and surface uplift across the Zagros. However, new structural and stratigraphic data show that the most-proximal Bakhtiyari exposures, in the High Zagros south of Shahr-kord, were deposited during the early Miocene and probably Oligocene. In this locality, a coarse-grained Bakhtiyari succession several hundred meters thick contains gray marl, limestone, and sandstone with diagnostic marine pelecypod, gastropod, coral, and coralline algae fossils. Foraminiferal and palynological species indicate deposition during early Miocene time. However, the lower Miocene marine interval lies in angular unconformity above ~ 150 m of Bakhtiyari conglomerate that, in turn, unconformably caps an Oligocene marine sequence. These relationships attest to syndepositional deformation and suggest that the oldest Bakhtiyari conglomerate could be Oligocene in age.

The new age information constrains the timing of initial foreland-basin development and proximal Bakhtiyari deposition in the Zagros hinterland. These findings reveal that structural evolution of the High Zagros was underway by early Miocene and probably Oligocene time, earlier than commonly envisioned. The age of the Bakhtiyari Formation in the High Zagros contrasts significantly with the Pliocene–Quaternary Bakhtiyari deposits near the modern deformation front, suggesting a long-term (> 20 Myr) advance of deformation toward the foreland.     

Observed daily large-scale rainfall patterns during BOBMEX-1999

A daily rainfall dataset and the corresponding rainfall maps have been produced by objective analysis of rainfall data. The satellite estimate of rainfall and the raingauge values are merged to form the final analysis. Associated with epochs of monsoon these rainfall maps are able to show the rainfall activities over India and the Bay of Bengal region during the BOBMEX period. The intra-seasonal variations of rainfall during BOBMEX are also seen using these data. This dataset over the oceanic region compares well with other available popular datasets like GPCP and CMAP. Over land this dataset brings out the features of monsoon in more detail due to the availability of more local raingauge stations.     

Geochemistry and mineralogy of shallow alluvial aquifers in Daudkandi upazila in the Meghna flood plain, Bangladesh

The shallow alluvial aquifers of the delta plains and flood plains of Bangladesh, comprises about 70% of total land area are mostly affected by elevated concentrations of arsenic (As) in groundwater exposing a population of more than 35 million to As toxicity. Geochemical studies of shallow alluvial aquifer in the Meghna flood plain show that the uppermost yellowish grey sediment is low in As (1.03 mg/kg) compared to the lower dark grey to black sediment (5.24 mg/kg) rich in mica and organic matter. Sequential extraction data show that solid phase As bound to poorly crystalline and amorphous metal (Fe, Mn, Al)-oxyhydroxides is dominant in the grey to dark grey sediment and reaches its maximum level (3.05 mg/kg) in the mica rich layers. Amount of As bound to sulphides and organic matter also peaks in the dark grey to black sediment. Vertical distributions of major elements determined by X-ray fluorescence (XRF) show that iron (Fe₂O₃), aluminum (Al₂O₃) and manganese (MnO) follow the general trend of distribution of As in the sediments. Concentrations of As, Mn, Fe, HCO₃ ⁻, SO₄ ²⁻ and NO₃ ⁻ in groundwater reflect the redox status of the aquifer and are consistent with solid phase geochemistry. Mineralogical analysis by X-ray diffraction (XRD) and scanning electron microscopy (SEM) fitted with energy dispersive X-ray spectrometer (EDS) revealed dominance of crystalline iron oxides and hydroxides like magnetite, hematite and goethite in the oxidised yellowish grey sediment. Amorphous Fe-oxyhydroxides identified as grain coating in the mica and organic matter rich sediment suggests weathering of biotite is playing a critical role as the source of Fe(III)-oxyhydroxides which in turn act as sink for As. Presence of authigenic pyrite in the dark grey sediment indicates active reduction in the aquifer.     

A BME solution of the inverse problem for saturated groundwater flow

In most real-world hydrogeologic situations, natural heterogeneity and measurement errors introduce major sources of uncertainty in the solution of the inverse problem. The Bayesian Maximum Entropy (BME) method of modern geostatistics offers an efficient solution to the inverse problem by first assimilating various physical knowledge bases (hydrologic laws, water table elevation data, uncertain hydraulic resistivity measurements, etc.) and then producing robust estimates of the subsurface variables across space. We present specific methods for implementing the BME conceptual framework to solve an inverse problem involving Darcys law for subsurface flow. We illustrate one of these methods in the case of a synthetic one-dimensional case study concerned with the estimation of hydraulic resistivity conditioned on soft data and hydraulic head measurements. The BME framework processes the physical knowledge contained in Darcys law and generates accurate estimates of hydraulic resistivity across space. The optimal distribution of hard and soft data needed to minimize the associated estimation error at a specified sampling cost is determined. This work was supported by grants from the National Institute of Environmental Health Sciences (Grant no. 5 P42 ES05948 and P30ES10126), the National Aeronautics and Space Administration (Grant no. 60-00RFQ041), the Army Research Office (Grant no. DAAG55-98-1-0289), and the National Science Foundation under Agreement No. DMS-0112069.     

Tectonics of a fast spreading center: A Deep-Tow and sea beam survey on the East Pacific rise at 19°30′ S

Sea Beam and Deep-Tow were used in a tectonic investigation of the fast-spreading (151 mm yr^-1) East Pacific Rise (EPR) at 19°30 S. Detailed surveys were conducted at the EPR axis and at the Brunhes/Matuyama magnetic reversal boundary, while four long traverses (the longest 96 km) surveyed the rise flanks. Faulting accounts for the vast majority of the relief. Both inward and outward facing fault scarps appear in almost equal numbers, and they form the horsts and grabens which compose the abyssal hills. This mechanism for abyssal hill formation differs from that observed at slow and intermediate spreading rates where abyssal hills are formed by back-tilted inward facing normal faults or by volcanic bow-forms. At 19°30 S, systematic back tilting of fault blocks is not observed, and volcanic constructional relief is a short wavelength signal (less than a few hundred meters) superimposed upon the dominant faulted structure (wavelength 2–8 km). Active faulting is confined to within approximately 5–8 km of the rise axis. In terms of frequency, more faulting occurs at fast spreading rates than at slow. The half extension rate due to faulting is 4.1 mm yr^-1 at 19°30 S versus 1.6 mm yr^-1 in the FAMOUS area on the Mid-Atlantic Ridge (MAR). Both spreading and horizontal extension are asymmetric at 19°30 S, and both are greater on the east flank of the rise axis. The fault density observed at 19°30 S is not constant, and zones with very high fault density follow zones with very little faulting. Three mechanisms are proposed which might account for these observations. In the first, faults are buried episodically by massive eruptions which flow more than 5–8 km from the spreading axis, beyond the outer boundary of the active fault zone. This is the least favored mechanism as there is no evidence that lavas which flow that far off axis are sufficiently thick to bury 50–150 m high fault scarps. In the second mechanism, the rate of faulting is reduced during major episodes of volcanism due to changes in the near axis thermal structure associated with swelling of the axial magma chamber. Thus the variation in fault spacing is caused by alternate episodes of faulting and volcanism. In the third mechanism, the rate of faulting may be constant (down to a time scale of decades), but the locus of faulting shifts relative to the axis. A master fault forms near the axis and takes up most of the strain release until the fault or fault set is transported into lithosphere which is sufficiently thick so that the faults become locked. At this point, the locus of faulting shifts to the thinnest, weakest lithosphere near the axis, and the cycle repeats.     

Site effects at Euroseistest—I. Determination of the valley structure and confrontation of observations with 1D analysis

This paper describes the process of construction of the 2D model of Volvi's geological structure and results of empirical and theoretical approaches to the evaluation of site response at Euroseistest. The construction of the 2D model is based on a re-interpretation of the available geophysical and geotechnical data in an effort to improve the definition of the subsoil structure at Euroseistest in terms of the most important parameters needed to model site response. The results of this re-interpretation are compared with a previous published 2D model of the same alluvial valley. Different analysis of the measurements and different criteria in the synthesis of data have led to a different model, even if both studies had access to the same field measurements. This underscores the fact that a model results of an interpretation and is not uniquely determined by the data, no matter how detailed they are. The well known subsoil structure opened the possibility to correlate the geometry and the dynamic properties of the 2D model with the results of site response determined from a detailed analysis of two events in frequency and time domains and 1D numerical modeling. The study of site response shows the important effect of the lateral variations on the ground motion and suggests that the contribution of locally generated surface waves to the resonant peak may be important. In the case of Volvi's graben, the limitations of the 1D approximation to simulate ground motion under complex soil conditions in both frequency and time domains are also shown. This paper lays the ground for a companion article dealing with 2D site effects in this basin.     

THE SCATTERING OF SH-CHANNEL WAVES BY A FAULT IN A COAL SEAM*

The imaging of faults in coal seams by the in-seam seismic method has now become standard practice. In the UK over 300 surveys have been undertaken and the technique is now part of the exploration arsenal of colliery planners. From these users comes the pressure for two major improvements, namely an increase in range and target identification. This paper is directed towards the latter problem. It has long been recognized that the reflected channel waves must contain information on the fault structure that caused the reflection, and model experiments have been undertaken to investigate the reflection process. Only recently, however, have attempts been made to quantify the reflection process. Calculations using both the finite-difference and finite-element techniques have been carried out, and estimates of the reflection coefficient as a function of frequency have been obtained. The object of this paper is to extend these considerations by calculating analytically the scattering matrix of an SH-channel wave after interaction at a fault plane. The scattering matrix is calculated as a function of frequency, hade angle, and fault throw. The method employed is based on the decomposition of the incident SH-channel wave into Fourier components, the calculation of plane wave reflection and transmission coefficients within the constraints of geometrical acoustics, and finally the synthesis of the scattering matrix by application of the Helmholtz-Kirchhoff integral. The calculation throughout is restricted to normal modes.     

EISCAT/CRRES observations: nightside ionospheric ion outflow and oxygen-rich substorm injections

We present combined observations made near midnight by the EISCAT radar, all-sky cameras and the combined released and radiation efects satellite (CRRES) shortly before and during a substorm. In particular, we study a discrete, equatorward-drifting auroral arc, seen several degrees poleward of the onset region. The arc passes through the field-aligned beam of the EISCAT radar and is seen to be associated with a considerable upflow of ionospheric plasma. During the substorm, the CRRES satellite observed two major injections, 17 min apart, the second of which was dominated by O⁺ ions. We show that the observed are was in a suitable location in both latitude and MLT to have fed O⁺ ions into the second injection and that the upward flux of ions associated with it was sufficient to explain the observed injection. We interpret these data as showing that arcs in the nightside plasma-sheet boundary layer could be the source of O⁺ ions energised by a dipolarisation of the mid- and near-Earth tail, as opposed to ions ejected from the dayside ionosphere in the cleft ion fountain.