A discretization and multigrid solver for a Darcy–Stokes system of three dimensional vuggy porous media

We develop a finite element discretization and multigrid solver for a Darcy–Stokes system of three-dimensional vuggy porous media, i.e., porous media with cavities. The finite element method uses low-order mixed finite elements in the Darcy and Stokes domains and special transition elements near the Darcy–Stokes interface to allow for tangential discontinuities implied by the Beavers–Joseph boundary condition. We design a multigrid method to solve the resulting saddle point linear system. The intertwining of the Darcy and Stokes subdomains makes the resulting matrix highly ill-conditioned. The velocity field is very irregular, and its discontinuous tangential component at the Darcy–Stokes interface makes it difficult to define intergrid transfer operators. Our definition is based on mass conservation and the analysis of the orders of magnitude of the solution. The coarser grid equations are defined using the Galerkin method. A new smoother of Uzawa type is developed based on taking an optimal step in a good search direction. Our algorithm has a measured convergence factor independent of the size of the system, at least when there are no disconnected vugs. We study the macroscopic effective permeability of a vuggy medium, showing that the influence of vug orientation; shape; and, most importantly, interconnectivity determine the macroscopic flow properties of the medium. This work was supported by the U.S. National Science Foundation under grants DMS-0074310 and DMS-0417431.     

The hydrochemical facies and anomalous fluoride content in the Serra Geral aquifer system,southern Brazil: a GIS approach with tectonic and principal component analysis

Groundwater with high fluoride content and water mixture patterns were studied in Serra Geral aquifer system (SGAS) using three aspects, principal component analysis (PCA), tectonic scenery and hydrochemical interpretation from 309 groundwater chemical data information from deep wells. A four-component model is suggested and explains 81% of total variance in the PCA. Six hydrochemical facies were identified. These facies suggest two different fluoride sources. Tectonic approach shows the relationship between defined hydrochemical facies and regional fracture control. The applied methodology reveals a minimum level to understand hydrochemical mixtures. The fluoride enrichment mechanisms into the groundwater are comprised in advance to guide the future uses of SGAS to the public supply.     

Dynamics of sedimentary environments in the accelerated siltation of a reservoir: the case of Alhama de Granada,southern Spain

An analysis has been made of sedimentary systems involved in the rapid silting of a reservoir constructed in 1974 in Alhama de Granada (S. Spain); in only 30 years the storage capacity of the reservoir has shrunk by 80% and its perimeter has decreased by 64%. A study of sediment lithofacies identified in a series of shallow trenches and of georadar facies identified in a series of almost 900 m lines of ground penetrating radar (GPR) images, together with a survey of surface geology, has identified 3 alluvial systems (2 transversal systems and a longitudinal system) whose deltas have filled in the reservoir. Thus, there are three phases in the evolution of the reservoir siltation: (1) an initial stage (1974–1977) typified by northward progradation of the longitudinal river delta of about 100 m year⁻¹ and an eastward progradation of the transversal system delta of about 20 m year⁻¹; (2) an intermediate stage (1977–1984) in which the longitudinal river delta progradation slowed to 25 m year⁻¹ and the axial drainage became obstructed due to the considerable eastward progradation of the transversal delta; and (3) a final phase (1984–present) in which there have been few changes in the areal distribution of the deltas apart from a southeastward expansion of the transversal delta. Generally, aggradational growth patterns (vertical accretion) have dominated in this final phase. The lithology of the source area, the slope and precipitation distribution has a significant effect not only on the sediment supply, but also indirectly on the creation of accommodation space and on the evolution of stratal growth patterns.     

Dolomitization,gypsum calcitization and silicification in carbonate–evaporite shallow lacustrine deposits

This paper describes and interprets the mineral and facies assemblages that occur in carbonate–evaporite shallow lacustrine deposits, considering the importance of the processes pathway (i.e. dolomitization, gypsum calcitization and silicification). The Palaeogene deposits of the Deza Formation (Almazán Basin, central‐northern Spain) are selected as a case study to determine the variety of physicochemical processes taking place in carbonate–evaporite shallow lakes and their resulting diagenetic features. Dolostones are the predominant lithology and are composed mainly of dolomite with variable amounts of secondary calcite (5 to 50%), which mainly mimic lenticular gypsum (pseudomorphs). Five morphological types of dolomite crystal were identified as follows: dolomite tubes, dolomite cylinders, rhombohedral dolomite, spheroidal and quasi‐rhombohedral dolomite, and cocoon‐shaped dolomite. The dolomite cylinders and tubes are interpreted as the dolomitized cells of a widespread microbial community. The sequence of diagenetic processes started with growth of microlenticular interstitial gypsum in a calcareous mud deposited on the playa margin mudflats, and that sometimes included microbial sediments. Immediately following growth of gypsum, dolomite replaced the original calcite (or possibly aragonite) muds, the microbial community and the gypsum. Partial or total replacement of gypsum by dolomite was related mainly to the biomineralization of endolithic microbial communities on gypsum crystals. Later calcitization took place under vadose, subaerial exposure conditions. The development of calcrete in distal alluvial settings favoured the release of silica and subsequent silicification on the playa margin mudflats. Stable isotope compositions of calcite range from ?9·02 to ?5·83‰ δ¹³C_PDB and ?7·10 to 1·22‰ δ¹⁸O_PDB; for the dolomite, these values vary from ?8·93 to ?3·96‰ δ¹³C_PDB and ?5·53 to 2·4‰ δ¹⁸O_PDB. Quartz from the cherts has δ¹⁸O_SMOW values ranging from 27·1 to 31·1‰. Wide variation and relatively high δ¹⁸O_SMOW values for dolomite indicate evaporitic and closed hydrological conditions; increased influx of meteoric waters reigned during the formation of secondary calcite spar.     

River tributaries as sediment sinks: Processes operating where the Tapajós and Xingu rivers meet the Amazon tidal river

The Amazon River is the largest fluvial source of fresh water and sediment to the global ocean and has the longest tidally influenced reach in the world. Two major rivers, the Tapajós and Xingu, enter the Amazon along its tidal reach. However, unlike most fluvial confluences, these are not one‐way conduits through which water and sediment flow downstream towards the sea. The drowned‐river valleys (rias) at the confluences of the Tapajós and Xingu with the Amazon River experience water‐level fluctuations associated not only with the seasonal rise and fall of the river network, but also with semidiurnal tides that propagate as far as 800 km up the Amazon River. Superimposed seasonal and tidal forcing, distinct sediment and temperature signatures of Amazon and tributary waters, and antecedent geomorphology combine to create mainstem–tributary confluences that act as sediment traps rather than sources of sediment. Hydrodynamic measurements are combined with data from sediment cores to determine the distribution of tributary‐derived and Amazon‐derived sediment within the ria basins, characterize the sediment‐transport mechanisms within the confluence areas and estimate rates of sediment accumulation within both rias. The Tapajós and Xingu ria basins trap the majority of the sediment carried by the tributaries themselves in addition to ca 20 Mt year^?1 of sediment sourced from the Amazon River. These findings have implications for the interpretation of stratigraphy associated with incised‐valley systems, such as those that dominated the transfer of sediment to the oceans during lowstands in sea level.     

Asbestiform and non-asbestiform morphologies in a talc and vermiculite mine from the province of Córdoba (Argentina): a case study

In this work, a talc and vermiculite mine from the province of Córdoba (Argentina) was investigated with special emphasis on the occurrence of asbestiform and non-asbestiform phases. The meta-ultramafic rock was studied by a multimethodological approach, complementing field studies with petrographic-mineralogical, compositional and morphological analyses. Samples were examined by stereomicroscopy, polarizing light microscopy, SEM–EDS, XRD, DSC-TGA and FTIR. Complementary, compositional and textural analyses were performed with FE-SEM–EDS and EPMA. Talc-rich veins with a laminar and fibrous appearance were at first recognized. However, the fibrous morphology observed both in the field and by microscopy is due to an apparent habit because of the sample orientation. To avoid erroneous interpretations, studies by secondary electron images (SEM) are fundamental to carrying out this type of analysis. Tremolite was identified in different zones of the outcrop; however, only ~40% of the crystals located in the vermiculite zone have dimensions to be considered as asbestiform fibres in the range of respirable particles. In these types of complex deposits affected by superimposed metamorphic, igneous and deformational events, multimethodological approaches are necessary to develop models of occurrence of asbestiform morphologies that may be applicable to other with similar characteristics.     

Estuarine flooding in urban areas: enhancing vulnerability assessment

The assessment of vulnerability provides valuable knowledge in the risk assessment steps of a risk governance process. Given the multiscale, multilevel, and multisectorial aspects of flood risk, the diversified entities that directly and indirectly intervene in risk management require specific outputs from the assessment studies. Urban areas in estuarine margins are particularly exposed and vulnerable to flooding. Such interface conditions are found in the Old City Centre of the Seixal, located in the Tagus estuary, Portugal. Here, two distinct methodologies were applied for the assessment of territorial vulnerability. A regional, lower-scale, methodology explores the application of the statistical procedure based on the SoVI^® at the statistical block level. A second, local and higher-scale, methodology is based in data collected through field matrices at the building and statistical sub-block level. Comparison of results revealed that the lower-scale assessment provides information on the vulnerability drivers at the regional and municipal level. Nevertheless, only at a higher-scale, it is possible to characterize and differentiate the smaller geographical units of analysis that compose the Old City Centre of Seixal. The lower-scale vulnerability assessment allows a strategic response, based on adaptation measures such as spatial planning, institutional capacity building and public awareness. The local level assessment provides more accurate knowledge to support local emergency planning and the allocation of operational and material resources at the urban level. Nevertheless, rather than antagonistic, both models can be considered as complementary, having in mind the requirements of an holistic flood risk governance model.

     

Selenium and other trace elements in phosphate rock of Bayovar–Sechura (Peru)

Selenium is an essential trace element in regard to human and other animals' health, but it presents a narrow range of dietary deficiency and toxic excess. Some phosphate rocks are rich in selenium. Furthermore it is necessary to control selenium's and other trace elements’ compositions in the phosphate rocks, because phosphorite ores of different geographical sites are varied, and even ores from the same outcrop can vary substantially in trace element levels.     

Morphological and geochemical analysis of the Laguna Blanca/Zapaleri obsidian source in the Atacama Puna

This paper describes a large obsidian deposit located along the southern banks of Laguna Blanca, on the eastern slope of the Jarellón volcanic caldera near the Chilean–Bolivian border. The obsidian at this site occurs in flows or sheets of deflated black or reddish to brown pebbles, redeposited on the shores of a lake. Blocks of obsidian are only found around the caldera rim. Here we analyze the shape of the obsidian pebbles and their geochemistry, comparing them with previously published data. The results indicate that the geochemical composition of the samples strongly matches previous analyses of obsidian from cultural contexts. This obsidian source was one of the most important sites of obsidian procurement since at least the Formative Period in the Atacama Puna region. © 2010 Wiley Periodicals, Inc.     

Sulfur K-edge XANES analysis of natural and synthetic basaltic glasses: Implications for S speciation and S content as function of oxygen fugacity

XANES analyses at the sulfur K-edge were used to determine the oxidation state of S species in natural and synthetic basaltic glasses and to constrain the fO₂ conditions for the transition from sulfide (S²⁻) to sulfate (S⁶⁺) in silicate melts. XANES spectra of basaltic samples from the Galapagos spreading center, the Juan de Fuca ridge and the Lau Basin showed a dominant broad peak at 2476.8 eV, similar to the spectra obtained from synthetic sulfide-saturated basalts and pyrrhotite. An additional sharp peak at 2469.8 eV, similar to that of crystalline sulfides, was present in synthetic glasses quenched from hydrous melts but absent in anhydrous glasses and may indicate differences in sulfide species with hydration or presence of minute sulfide inclusions exsolved during quenching. The XANES spectra of a basalt from the 1991 eruption of Mount Pinatubo, Philippines, and absarokitic basalts from the Cascades Range, Oregon, USA, showed a sharp peak at 2482.8 eV, characteristic of synthetic sulfate-saturated basaltic glasses and crystalline sulfate-bearing minerals such as hauyne. Basaltic samples from the Lamont Seamount, the early submarine phase of Kilauea volcano and the Loihi Seamount showed unequivocal evidence of the coexistence of S²⁻ and S⁶⁺ species, emphasizing the relevance of S⁶⁺ to these systems. XANES spectra of basaltic glasses synthesized in internally-heated pressure vessels and equilibrated at fO₂ ranging from FMQ − 1.4 to FMQ + 2.7 showed systematic changes in the features related to S²⁻ and S⁶⁺ with changes in fO₂. No significant features related to sulfite (S⁴⁺) species were observed. These results were used to construct a function that allows estimates of S⁶⁺/ΣS from XANES data. Comparison of S⁶⁺/ΣS data obtained by S Kα shifts measured with electron probe microanalysis (EPMA), S⁶⁺/ΣS obtained from XANES spectra, and theoretical considerations show that data obtained from EPMA measurements underestimate S⁶⁺/ΣS in samples that are sulfate-dominated (most likely because of photo-reduction effects during analysis) whereas S⁶⁺/ΣS from XANES provide a close match to the expected theoretical values. The XANES-derived relationship for S⁶⁺/ΣS as a function of fO₂ indicates that the transition from S²⁻ to S⁶⁻ with increasing fO₂ occurs over a narrower interval than what is predicted by the EPMA-derived relationship. The implications for natural systems is that small variation of fO₂ above FMQ + 1 will have a large effect on S behavior in basaltic systems, in particular regarding the amount of S that can be transported by basaltic melts before sulfide saturation can occur.