Geochemistry and mineralogy of arsenic in (natural) anaerobic groundwaters

Here new data from field bioremediation experiments and geochemical modeling are reported to illustrate the principal geochemical behavior of As in anaerobic groundwaters. In the field bioremediation experiments, groundwater in Holocene alluvial aquifers in Bangladesh was amended with labile water-soluble organic C (molasses) and MgSO₄ to stimulate metabolism of indigenous SO₄-reducing bacteria (SRB). In the USA, the groundwater was contaminated by Zn, Cd and SO₄, and contained <10 μg/L As under oxidized conditions, and a mixture of sucrose and methanol were injected to stimulate SRB metabolism. In Bangladesh, groundwater was under moderately reducing conditions and contained ∼10 mg/L Fe and ∼100 μg/L As. In the USA experiment, groundwater rapidly became anaerobic, and dissolved Fe and As increased dramatically (As > 1000 μg/L) under geochemical conditions consistent with bacterial Fe-reducing conditions. With time, groundwater became more reducing and biogenic SO₄ reduction began, and Cd and Zn were virtually completely removed due to precipitation of sphalerite (ZnS) and other metal sulfide mineral(s). Following precipitation of chalcophile elements Zn and Cd, the concentrations of Fe and As both began to decrease in groundwater, presumably due to formation of As-bearing FeS/FeS₂. By the end of the six-month experiment, dissolved As had returned to below background levels. In the initial Bangladesh experiment, As decreased to virtually zero once biogenic SO₄ reduction commenced but increased to pre-experiment level once SO₄ reduction ended. In the ongoing experiment, both SO₄ and Fe(II) were amended to groundwater to evaluate if FeS/FeS₂ formation causes longer-lived As removal. Because As-bearing pyrite is the common product of SRB metabolism in Holocene alluvial aquifers in both the USA and Southeast Asia, it was endeavored to derive thermodynamic data for arsenian pyrite to better predict geochemical processes in naturally reducing groundwaters. Including the new data for arsenian pyrite into Geochemist’s Workbench, its stability field completely dominates in reducing Eh–pH space and “displaces” other As-sulfides (orpiment, realgar) that have been implied to be important in previous modeling exercises and reported in rare field conditions.     

Considerations for conducting incubations to study the mechanisms of As release in reducing groundwater aquifers

Microbial Fe reduction is widely believed to be the primary mechanism of As release from aquifer sands in Bangladesh, but alternative explanations have been proposed. Long-term incubation studies using natural aquifer material are one way to address such divergent views. This study addresses two issues related to this approach: (1) the need for suitable abiotic controls and (2) the spatial variability of the composition of aquifer sands. Four sterilization techniques were examined using orange-colored Pleistocene sediment from Bangladesh and artificial groundwater over 8 months. Acetate (10 mM) was added to sacrificial vials before sterilization using either (1) 25 kGy of gamma irradiation, (2) three 1-h autoclave cycles, (3) a single addition of an antibiotic mixture at 1× or (4) 10× the typical dose, and (5) a 10 mM addition of azide. The effectiveness of sterilization was evaluated using two indicators of microbial Fe reduction, changes in diffuse spectral reflectance and leachable Fe(II)/Fe ratios, as well as changes in P-extractable As concentrations in the solid phase. A low dose of antibiotics was ineffective after 70 days, whereas autoclaving significantly altered groundwater composition. Gamma irradiation, a high dose of antibiotics, and azide were effective for the duration of the experiment.     

Coupled hydrogeological and reactive transport modelling of the Simpevarp area (Sweden)

The Simpevarp area is one of the alternative sites being considered for the deep geological disposal of high level radioactive waste in Sweden. In this paper, a coupled regional groundwater flow and reactive solute transport model of the Simpevarp area is presented that integrates current hydrogeological and hydrochemical data of the area. The model simulates the current hydrochemical pattern of the groundwater system in the area. To that aim, a conceptual hydrochemical model was developed in order to represent the dominant chemical processes. Groundwater flow conditions were reproduced by taking into account fluid-density-dependent groundwater flow and regional hydrogeologic boundary conditions. Reactive solute transport calculations were performed on the basis of the velocity field so obtained. The model was calibrated and sensitivity analyses were carried out in order to investigate the effects of heterogeneities of hydraulic conductivity in the subsurface medium. Results provided by the reactive transport model are in good agreement with much of the measured hydrochemical data. This paper emphasizes the appropriateness of the use of reactive solute transport models when water-rock interaction reactions are involved, and demonstrates what powerful tools they are for the interpretation of hydrogeological and hydrochemical data from site geological repository characterization programs, by providing a qualitative framework for data analysis and testing of conceptual assumptions in a process-oriented approach.     

First U–Pb SHRIMP age of the Hauterivian stage, Neuquén Basin, Argentina

A high-resolution ion-microprobe (SHRIMP) U–Pb zircon age from a tuff layer intercalated in the ammonoid bearing sedimentary succession of the Neuquén Basin in Argentina provides a robust geochronologic date to add to the absolute ages and to improve the relative chronology of the Early Cretaceous Hauterivian stage. The tuff layer appears interbedded between shales of the upper member (Agua de la Mula) of the Agrio Formation within the Spitidiscus riccardii ammonoid zone (base of the Late Hauterivian) yielding a date of 132.5 ± 1.3 Ma. This date confirms and supports an accurate correlation between the ammonoid biostratigraphy of the Neuquén Basin with the Western Mediterranean Province of the Tethys during the Early Cretaceous and matches with the most recently published time scale. It also casts doubts on the validity of K–Ar ages on glauconite-grains recently reported from the Lower Cretaceous of the Vocontian Basin of France.     

An operational multiscale system for hazards prediction, mapping, and response

By definition, a crisis is a situation that requires assistance to be managed. Hence, response to a crisis involves the merging of local and non-local emergency response personnel. In this situation, it is critical that each participant: (1) know the roles and responsibilities of each of the other participants; (2) know the capabilities of each of the participants; and (3) have a common basis for action. For many types of natural disasters, this entails having a common operational picture of the unfolding events, including detailed information on the weather, both current and forecasted, that may impact on either the emergency itself or on response activities. The Consequences Assessment Tool Set (CATS) is a comprehensive package of hazard prediction models and casualty and damage assessment tools that provides a linkage between a modeled or observed effect and the attendant consequences for populations, infrastructure, and resources, and, hence, provides the common operational picture for emergency response. The Operational Multiscale Environment model with Grid Adaptivity (OMEGA) is an atmospheric simulation system that links the latest methods in computational fluid dynamics and high-resolution gridding technologies with numerical weather prediction to provide specific weather analysis and forecast capability that can be merged into the geographic information system framework of CATS. This paper documents the problem of emergency response as an end-to-end system and presents the integrated CATS–OMEGA system as a prototype of such a system that has been used successfully in a number of different situations.     

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

Drought and institutional adaptation in the Great Plains of Alberta and Saskatchewan, 1914–1939

Agriculture in the southern Great Plains of Canada has been particularly vulnerable to prolonged episodes of drought. Using climate data and a precipitation minus potential evapotranspiration index, the extent of the region’s exposure to drought is examined. Between 1914 and 1917, the Dry Belt was particularly vulnerable to drought, whereas after 1928, a much larger region known as the Palliser Triangle covering most of southern Alberta and Saskatchewan was much more exposed to drought. These droughts provoked major institutional adaptation, in particular the establishment of the Special Areas Board by the Government of Alberta, and the creation of the Prairie Farm Rehabilitation Administration by the Government of Canada. Both organizations have proved to be relatively permanent public adaptations to the natural hazard of drought in the region. Moreover, these earlier experiences with prolonged drought as well as institution-building may be of value in helping the residents of the Palliser Triangle adapt to predicted climate changes in the future as well as anticipate some of the barriers to effective institutional adaptation.     

Tidal inundation mapping under enhanced land subsidence in Semarang,Central Java Indonesia

Tidal inundation by high tide under enhanced land subsidence is a damaging phenomenon and a major threat to the Semarang urban area in Indonesia. It impacts on economic activities, as well as the cost of an emergency program and causes interruption of pubic services, danger of infectious diseases and injury to human lives. This study examines a spatial analysis tool on the GIS-raster system for the tidal inundation mapping based on the subsidence-benchmark data and modified detail digital elevation model. Neighborhood operation and iteration model as a spatial analysis tool have been applied in order to calculate the encroachment of the tidal inundation on the coastal area. The resulting map shows that the tidal flood spreads to the lowland area and causes the inundation of coastal settlement, infrastructure, as well as productive agricultural land, i.e., the fish-pond area. The monitoring of the vulnerable area due to the tidal inundation under the scenario of extended land subsidence plays an important role in long-term coastal zone management in Semarang.