Comparison between cell-centered and nodal-based discretization schemes for linear elasticity

In this paper, we study newly developed methods for linear elasticity on polyhedral meshes. Our emphasis is on applications of the methods to geological models. Models of subsurface, and in particular sedimentary rocks, naturally lead to general polyhedral meshes. Numerical methods which can directly handle such representation are highly desirable. Many of the numerical challenges in simulation of subsurface applications come from the lack of robustness and accuracy of numerical methods in the case of highly distorted grids. In this paper, we investigate and compare the Multi-Point Stress Approximation (MPSA) and the Virtual Element Method (VEM) with regard to grid features that are frequently seen in geological models and likely to lead to a lack of accuracy of the methods. In particular, we look at how the methods perform near the incompressible limit. This work shows that both methods are promising for flexible modeling of subsurface mechanics.     

Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

The Juchatengo complex: an upper-level ophiolite assemblage of late Paleozoic age in Oaxaca,southern Mexico

International Journal of Earth Sciences - The Juchatengo complex (JC) suite is located between the Proterozoic Oaxacan complex to the north and the Xolapa complex to the south, and was amalgamated...     

Contrasting sources of Late Paleozoic rhyolite magma in the Polish Lowlands: evidence from U–Pb ages and Hf and O isotope composition in zircon

International Journal of Earth Sciences - The Polish Lowlands, located southwest of the Teisseyre–Tornquist Zone, within Trans-European Suture Zone, were affected by bimodal, but dominantly...     

Rainfall thresholds for the triggering of landslides considering previous soil moisture conditions (Asturias,NW Spain)

Given its geological and climatic conditions and its rugged orography, Asturias is one of the most landslide prone areas in the North of Spain. Most of the landslides occur during intense rainfall episodes. Thus, precipitation is considered the main triggering factor in the study area, reaching average annual values of 960 mm. Two main precipitation patterns are frequent: (i) long-lasting periods of moderate rainfall during autumn and winter and (ii) heavy short rainfall episodes during spring and early summer. In the present work, soil moisture conditions in the locations of 84 landslides are analysed during two rainfall episodes, which represent the most common precipitation patterns: October–November 2008 and June 2010. Empirical data allowed the definition of available water capacity percentages of 99–100% as critical soil moisture conditions for the landslide triggering. Intensity-duration rainfall thresholds were calculated for each episode, considering the periods with sustained high soil moisture levels before the occurrence of each analysed landslide event. For this purpose, data from daily water balance models and weather stations were used. An inverse relationship between the duration of the precipitation and its intensity, consistent with published intensity-duration thresholds, was observed, showing relevant seasonal differences.     

Non-woven polypropylene fabric modified with carbon nanotubes and decorated with nanoakaganeite for arsenite removal

Due to its harmful impact on human health, the presence of heavy metals, metalloids and other toxic pollutants in drinking or irrigation water is a major concern. Recent studies have proved that nanosized adsorbents are significantly more effective than their microsized counterparts. Particular attention has been given to nanocomposites with nanoadsorbents embedded in matrixes that could provide stability to the material and contribute to eliminating problems that may appear when using conventional granular systems. This study presents the preparation of a novel hybrid filter from a commercially available polypropylene (PP) non-woven fabric matrix modified with multiwall carbon nanotubes (MWCNT) and iron oxy(hydroxide) nanoparticles, and its use in the removal of As(III). A Box–Behnken statistical experimental design has been chosen to explore relevant variables affecting the filter performance: (1) As(III) concentration, (2) pH and (3) sorbent dose. From an As(III) concentration of 10 mg L^?1, at pH 6.5 and with a sorbent dose of 5 g L^?1, the PP filter modified with MWCNT removes 10% of the initial metalloid concentration, reaching a capacity of 0.27 mg g^?1. After modification with iron oxy(hydroxide), the performance of the material is largely enhanced. The filter, under the same conditions, removes 90% of the initial As(III) concentration, reaching a capacity almost tenfold higher (2.54 mg g^?1). This work demonstrates that the developed hybrid filter is effective toward the removal of As(III) in a wide range of pHs. A cubic regression model to compute the removal of the filter as a function of pH and sorbent dose is provided.     

Toxic metals in a highly urbanized industry-impacted estuary (Bahia Blanca Estuary,Argentina): spatio-temporal analysis based on GIS

Cadmium, chromium, nickel, and lead were evaluated in the particulate fraction at one of the most industrialized estuaries at the Southwestern Atlantic Ocean, through Geographic Information System (GIS). Concentrations were analyzed at 21 stations during 2008–2010. The highest metal concentrations (Cd: 8,9; Cr: 256,49; Ni: 27,02; Pb: 78,43 µg g^??1 d. w.) were recorded at the stations located near industrial and urban discharges situated along the estuary. In addition, Pb presented a different seasonal and spatial behavior in comparison with Cd, Cr and Ni. Winter and spring presented the higher concentrations of Pb, and the inner stations presented the higher values. The estuary is considered a moderate to strongly polluted and significantly polluted according to the Index of geoaccumulation (I_geo) and the Enrichment Factor of Cd, respectively. The Multidimensional Scaling plot showed three groups of stations: the inner, associated to low levels of metals (G1), middle stations (G2) with intermediate levels and the outer (G3) with the highest ones. In addition, this work reveals the usefulness of the GIS-mapping techniques in the distribution of pollutants along an estuarine environment and the environmental quality assessment of estuarine systems.     

Evidence for the occurrence of hydrogeochemical processes in the groundwater of Khoy plain,northwestern Iran,using ionic ratios and geochemical modeling

Rapid population growth, industrialization, and agricultural expansion in the Khoy area (northwestern Iran) have led to its dependence on groundwater and degradation of groundwater quality. This study attempts to decipher the major processes and factors that degrade the groundwater quality of the Khoy plain. For this purpose, 54 groundwater samples from unconfined and confined aquifers of the plain were collected in July 2017 and analyzed for major cations and anions (Na, K, Ca, Mg, HCO₃, SO_4, and Cl), minor ions (NO₃ and F), and Al. Magnesium and bicarbonate were identified as the dominant cation and anion, respectively. Several ionic ratios and geochemical modeling using PHREEQC indicated that the most important hydrogeochemical processes to affect groundwater quality in the plain were weathering and dissolution of evaporitic and silicate minerals, mixing, and ion exchange. There were smaller effects from evaporation and anthropogenic factors (e.g., industries). Results showed that the high salinity of the groundwater in the northeast area of the plain was due to the high solubility of the evaporitic minerals, e.g., halite and gypsum. Reverse ion exchange and the contribution of mineral dissolution were more significant than ion exchange in the northeastern part of the plain. Elevated salinity of the groundwater in the southeast was attributed mostly to reverse ion exchange and somewhat to evaporation.