Paleomagnetism and Detrital Zircon Geochronology of the Upper Vindhyan Sequence, Son Valley and Rajasthan, India: A ca. 1000 Ma Closure age for the Purana Basins?

The utility of paleomagnetic data gleaned from the Bhander and Rewa Groups of the “Purana-aged” Vindhyanchal Basin has been hampered by the poor age control associated with these units. Ages assigned to the Upper Vindhyan sequence range from Cambrian to the Mesoproterozoic and are derived from a variety of sources, including ⁸⁷Sr/⁸⁶Sr and δ ¹³C correlations with the global curves and Ediacara-like fossil finds in the Lakheri–Bhander limestone. New analyses of the available paleomagnetic data collected from this study and previous work on the 1073 Ma Majhgawan kimberlite, as well as detrital zircon geochronology of the Upper Bhander sandstone and sandstones from the Marwar SuperGroup suggest that the Upper Vindhyan sequence may be up to 500 Ma older than is commonly thought. Paleomagnetic analysis generated from the Bhander and Rewa Groups yields a paleomagnetic pole at 44°N, 214.0°E (A95 = 4.3°). This paleomagnetic pole closely resembles the VGP from the well-dated Majhgawan intrusion (36.8°N, 212.5°E, α₉₅ = 15.3°).Detrital zircon analysis of the Upper Bhander sandstone identifies a youngest age population at 1020 Ma. A comparison between the previously correlated Upper Bhander sandstone and the Marwar sandstone detrital suites shows virtually no similarities in the youngest detrital suite sampled. The main 840–920 Ma peak is absent in the Upper Bhander. This supports our assertion that the Upper Bhander is older than the 750–771 Ma Malani sequence, and is likely close to the age of the 1073 Ma Majhgawan kimberlite on the basis of the paleomagnetic similarities. By setting the age of the Upper Vindhyan at 1000–1070 Ma, several intriguing possibilities arise. The Bhander–Rewa paleomagnetic pole allows for a reconstruction of India at 1000–1070 Ma that overlaps with the 1073 ± 13.7 Majhgawan kimberlite VGP. Comparisons between the composite Upper Vindhyan pole (43.9°N, 210.2°E, α₉₅ = 12.2°) and the Australian 1071 ± 8 Ma Bangamall Basin sills and the 1070 Ma Alcurra dykes suggest that Australia and India were not adjacent at this time period.     

Intriguing climatic shifts in a 90 kyr old lake record from northern Russia

A 22 m long sediment core from Lake Yamozero on the Timan Ridge in northern Russia has provided evidence of intriguing climatic shifts during the last glacial cycle. An overall shallowing of the lake is reflected in the lower part of the cores, where pollen indicates a transition from glacial steppe vegetation to interstadial shrub-tundra. These beds are capped by a well-defined layer of compact clay deposited in relatively deep water, where pollen shows surrounding spruce forests and warmer-than-present summer temperatures. The most conservative interpretation is that this unit represents the last interglacial period. However, a series of Optical Stimulated Luminescence (OSL) dates suggests that it corresponds with the Early Weichselian Odderade interstadial (MIS 5a). This would imply that the Odderade interstadial was just as warm as a normal interglacial in this continental part of northern Europe. If correct, then pollen analysis, as a correlation tool, is less straightforward and the definition of an interglacial is more complex than previously thought. We discuss the validity and possible systematic errors of the OSL dates on which this age model is based, but conclude they really indicate a MIS 5a age for the warm period. Above the clay is an unconformity, most likely reflecting a period of subaerial exposure implying dry conditions. Deposition of silt under fluctuating cold climates in the Middle Weichselian continued until a second gap in the record at c . 40 kyr BP. The lake basin started to fill up again around 18 kyr BP.     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

History of crustal growth and recycling at the Pacific convergent margin of South America at latitudes 29°–36° S revealed by a U–Pb and Lu–Hf isotope study of detrital zircon from late Paleozoic accretionary systems

Detrital zircon provides a powerful archive of continental growth and recycling processes. We have tested this by a combined laser ablation ICP-MS U–Pb and Lu–Hf analysis of homogeneous growth domains in detrital zircon from late Paleozoic coastal accretionary systems in central Chile and the collisional Guarguaráz Complex in W Argentina. Because detritus from a large part of W Gondwana is present here, the data delineate the crustal evolution of southern South America at its Paleopacific margin, consistent with known data in the source regions.Zircon in the Guarguaráz Complex mainly displays an U–Pb age cluster at 0.93–1.46 Ga, similar to zircon in sediments of the adjacent allochthonous Cuyania Terrane. By contrast, zircon from the coastal accretionary systems shows a mixed provenance: Age clusters at 363–722 Ma are typical for zircon grown during the Braziliano, Pampean, Famatinian and post-Famatinian orogenic episodes east of Cuyania. An age spectrum at 1.00–1.39 Ga is interpreted as a mixture of zircon from Cuyania and several sources further east. Minor age clusters between 1.46 and 3.20 Ga suggest recycling of material from cratons within W Gondwana.The youngest age cluster (294–346 Ma) in the coastal accretionary prisms reflects a so far unknown local magmatic event, also represented by rhyolite and leucogranite pebbles. It sets time marks for the accretion history: Maximum depositional ages of most accreted metasediments are Middle to Upper Carboniferous. A change of the accretion mode occurred before 308 Ma, when also a concomitant retrowedge basin formed.Initial Hf-isotope compositions reveal at least three juvenile crust-forming periods in southern South America characterised by three major periods of juvenile magma production at 2.7–3.4 Ga, 1.9–2.3 Ga and 0.8–1.5 Ga. The ¹⁷⁶Hf/¹⁷⁷Hf of Mesoproterozoic zircon from the coastal accretionary systems is consistent with extensive crustal recycling and addition of some juvenile, mantle-derived magma, while that of zircon from the Guarguaráz Complex has a largely juvenile crustal signature. Zircon with Pampean, Famatinian and Braziliano ages (< 660 Ma) originated from recycled crust of variable age, which is, however, mainly Mesoproterozoic. By contrast, the Carboniferous magmatic event shows less variable and more radiogenic ¹⁷⁶Hf/¹⁷⁷Hf, pointing to a mean early Neoproterozoic crustal residence. This zircon is unlikely to have crystallized from melts of metasediments of the accretionary systems, but probably derived from a more juvenile crust in their backstop system.     

Sulfur isotope measurement of sulfate and sulfide by high-resolution MC-ICP-MS

We have developed a technique for the accurate and precise determination of ³⁴S/³²S isotope ratios (δ³⁴S) in sulfur-bearing minerals using solution and laser ablation multiple-collector inductively coupled plasma mass spectrometry (MC-ICP-MS). We have examined and determined rigorous corrections for analytical difficulties such as instrumental mass bias, unresolved isobaric interferences, blanks, and laser ablation- and matrix-induced isotopic fractionation. Use of high resolution sector-field mass spectrometry removes major isobaric interferences from O₂⁺. Standard-sample bracketing is used to correct for the instrumental mass bias of unknown samples. Background on sulfur masses arising from memory effects and residual oxygen-tailing are typically minor (< 0.2‰, within analytical error), and are mathematically removed by on-peak zero subtraction and by bracketing of samples with standards determined at the same signal intensity (within 20%). Matrix effects are significant (up to 0.7‰) for matrix compositions relevant to many natural sulfur-bearing minerals. For solution analysis, sulfur isotope compositions are best determined using purified (matrix-clean) sulfur standards and sample solutions using the chemical purification protocol we present. For in situ analysis, where the complex matrix cannot be removed prior to analysis, appropriately matrix-matching standards and samples removes matrix artifacts and yields sulfur isotope ratios consistent with conventional techniques using matrix-clean analytes. Our method enables solid samples to be calibrated against aqueous standards; a consideration that is important when certified, isotopically-homogeneous and appropriately matrix-matched solid standards do not exist. Further, bulk and in situ analyses can be performed interchangeably in a single analytical session because the instrumental setup is identical for both. We validated the robustness of our analytical method through multiple isotope analyses of a range of reference materials and have compared these with isotope ratios determined using independent techniques. Long-term reproducibility of S isotope compositions is typically 0.20‰ and 0.45‰ (2σ) for solution and laser analysis, respectively. Our method affords the opportunity to make accurate and relatively precise S isotope measurement for a wide range of sulfur-bearing materials, and is particularly appropriate for geologic samples with complex matrix and for which high-resolution in situ analysis is critical.     

Delineation of groundwater protection zones by the backward particle tracking method: theoretical background and GIS-based stochastic analysis

The backward particle tracking method, an effective and powerful tool that can be used to delineate groundwater protection zones, is presented. The theoretical background and insights on the applicability of this method are provided. Moreover, the present work enriches the backward particle tracking method with an uncertainty analysis concerning the porosity values, applying a Monte Carlo (MC) approach, coupled with the use of geographical information systems (GIS). As an application example, a wellfield in the Komotini area, Greece, is investigated. The present study may serve as a potential guideline for wellfield delineation, particularly in areas like Greece where lack of data related to the hydrogeological system is often a problem.     

Geomorphic changes in Ras Al-Subiyah area,Kuwait

The Ras Al-Subiyah area is considered one of the most promising areas in Kuwait for future development. This development will include a new town called Subiyah and its associated infrastructure. This area is also being considered as the location for connection between Boubyan Island, which is now undergoing major development and the Kuwait mainland. The present study investigates the geomorphology of the Ras Al-Sabiyah area in the northern sector of Kuwait. The study area is generally flat, and it is located west of the Jal Az-Zor escarpment. It is bordered on the east by the Khor Al-Sabiyah tidal channel and on the south by Kuwait Bay. The area receives sediments from several sources; currently the most important are aeolian sediments and the deposition of mud delivered through the Khor Al-Sabiyah from the Iraqi marshes. The study area has been subjected to severe environmental changes due to the Gulf wars and the drainage of Iraqi marshes and the associated artificial changes in fluvial system. Twenty-two surface sediments were collected from the Ras Al-Subiyah area. Samples were collected to include the main geomorphologic characteristic features of the study area. Field observations and remote sensing images from 1990 and 2001 were used to produce an updated geomorphologic map for the Ras Al-Subiyah and a map showing geomorphic changes between 1990 and 2001. Grain size of the surface sediment ranges from gravel to medium sand. In general, grain size statistical analysis indicates that most of the areas are composed of two or more classes of sands transported and deposited from different sources including aeolian, sabkhas, river and the bays. The variability in the grain size statistical parameters may be attributed to the complexity of surface morphology as well as the diversity in the type of depositional environment in the Ras Al-Subiyah area. The total area subjected to change during the 12-year period (1990–2001) is about 32 km² as calculated using GIS techniques.     

Impact evaluation on groundwater of the infiltration of an urban wastewater geopurification system placed over a detrital aquifer in a semiarid region (Spain)

Over a period of 4 years and 4 months, the geopurification installations at Dehesas de Guadix (Granada, Spain) were monitored to determine the impact on soil and groundwater of the controlled discharge of urban wastewater, and also to identify the best indicators of the entry of the recharged water into the aquifer. The installations are located in an area where the climate is Mediterranean sub-arid, with an average precipitation of less than 287 mm/year, and a rate of evapotranspiration that is almost three times greater. The system was controlled by determining the balance of majority nutrients and boron in the soil and in the groundwater, both at the points affected directly by the wastewater discharge and at others. The quantity of mass discharged was relatively large (COD 14,656 g/m², NO₃ 85 g/m², NO₂ 4 g/m², NH₄ 2,425 g/m², PO₄ 1,143 g/m², K 1,531 g/m², B 63 g/m²). It was observed that the elimination of nutrients within the soil (COD 97.5%, PO₄ 94.4%, K 59.17%, N _total 18.8%, B 12.69%) was very efficient except for the nitrogen, which nevertheless did not reach the groundwater, as it was eliminated at deep levels of the unsaturated zone. Only 12.69% of the boron was removed, and appreciable, increasing amounts of this element did reach the groundwater. Unexpectedly, none of the majority nutrients behaved as a reliable indicator of the impact on groundwater; despite this, the boron and the bicarbonate did clearly reflect the arrival of the recharged water, and are proposed as the best indicators.     

Doline probability map using logistic regression and GIS technology in the central Ebro Basin (Spain)

In the surroundings of Zaragoza, karstification processes are especially intense in covered karst areas where fluvial terraces lie directly on Tertiary evaporites. Since the beginning of Quaternary, these processes have lead to the development of collapse and subsidence dolines with a wide range of sizes, which have significant economic impacts. To reduce economic impact and increase safety, a regional analysis of this phenomenon is needed for spatial management. Therefore, a probability map of dolines was developed using logistic regression and geographic information system (GIS) techniques. This paper covers the selection of input data, manipulation of data using the GIS technology, and the use of logistic regression to generate a doline probability map. The primary variable in the doline development in this area is geomorphology, represented by the location of endorheic areas and different terrace levels. Secondary variables are the presence of irrigation and the water table gradient.     

Determining effective wellbore permeability from a field pressure test: a numerical analysis of detection limits

We propose a simple pressure test that can be used in the field to determine the effective permeability of existing wellbores. Such tests are motivated by the need to understand and quantify leakage risks associated with geological storage of CO₂ in mature sedimentary basins. If CO₂ is injected into a deep geological formation, and the resulting CO₂ plume encounters a wellbore, leakage may occur through various pathways in the “disturbed zone” surrounding the well casing. The effective permeability of this composite zone, on the outside of the well casing, is an important parameter for models of leakage. However, the data that exist on this key parameter do not exist in the open literature, and therefore specific field tests need to be done in order to reduce the uncertainty inherent in the leakage estimates. The test designed and analyzed herein is designed to measure effective wellbore permeability within a low-permeability caprock, bounded above and below by permeable reservoirs, by pressurizing the reservoir below and measuring the response in the reservoir above. Alternatively, a modified test can be performed within the caprock without directly contacting the reservoirs above and below. We use numerical simulation to relate pressure response to effective well permeability and then evaluate the range of detection of the effective permeability based on instrument measurement error and limits on fracture pressure. These results can guide field experiments associated with site characterization and leakage analysis.