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.     

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.     

The stratigraphical signature of the Anthropocene in England and its wider context

The Anthropocene deposits of England, here regarded as those formed after ～1950 CE, are now extensive, take various forms, and may be characterized and recognized by a number of stratigraphic signals, such as artificial radionuclides, pesticide residues, microplastics, enhanced fly ash levels, concrete fragments and a novel variety of ‘technofossils’ and neobiotic species. They include the uppermost parts of both ‘natural’ deposits such as the sediment layers formed in lakes and estuaries, and more directly human-made or human-influenced ones such as landfill deposits and the ‘artificial ground’ beneath urban areas and around major constructions. ‘Negative deposits’ include the worked areas of quarries and regions such as the English Fenland, where thick peat deposits have ablated to leave a strongly modified underlying landscape, and extend beneath into the subterranean realm as mine workings, metro systems and boreholes. The production of these is still rapidly increasing and evolving in character, while the early signs of global change, such as warming, sea level rise, and modifications to biotic assemblages, are beginning to further modify the emerging geology of this new phase of Earth history.     

Surface Geophysical Methods used to Verify the Karst Geological Structure in the Built-up Area: A Case Study of Specific Engineering-Geological Conditions

This article presents a research study of complex limestone karst engineering-geological conditions in the municipality Valaská near Banská Bystrica in Slovakia. The aim of the study is to demonstrate the impossibility of spatial identification of cave spaces using surface geophysical methods due to the specific engineering-geological conditions of a thick surface layer of anthropogenic fill containing highly heterogeneous anthropogenic material. Its maximum thickness is 3 m. Another specificific condition of the study area is its location in the built-up area, due to which the applicability of geophysical methods was limited. The article contains methodological recommendations to be used in analogous geological conditions with karst structures topped with anthropogenic fill, which complicates the identification of cave spaces. The recommended solution herein is the identification of the cave system using underground mapping of the karst and its projection onto the surface for which surface geophysical methods have been combined.     

Small-strain behaviour of unsaturated silty clay: experiments and model interpretation

Acta Geotechnica - The experimental data dealing with the so-called small strain stiffness of soils are indispensable in developing and calibrating advanced numerical models. A literature review...     

Effects of cadmium, zinc and nitrogen status on non-protein thiols in the macroalgae Enteromorpha spp. from the Scheldt Estuary (SW Netherlands, Belgium) and Thermaikos Gulf (N Aegean Sea, Greece)

Enteromorpha prolifera (Scheldt Estuary) and E. linza (Thermaikos Gulf) were incubated at three salinities with 100 and 200microgL(-1)Cd and Zn. The objective was to measure effects of Cd, Zn and nitrogen (N) status on the pools of metal-binding non-protein thiols: glutathione and phytochelatins, (gamma-glutamyl-cysteinyl)(n)-glycine (PC). In E. linza, ammonium pools were higher, but amino acid pools, total N and protein contents were lower than in E. prolifera. Reduced glutathione (GSH) pools were positively correlated with free glutamate and protein contents. In E. linza GSH pools increased and the ratio of reduced to oxidized glutathione (GSH:(GSH+0.5GSSG)), an indicator of oxidative stress, decreased with Cd contents, indicating Cd-induced glutathione oxidation. Total glutathione pools (reduced plus oxidized) ranged from 16nmolSgdwt(-1) in controls (at 0.5micromolCdgdwt(-1)) to 179nmolSgdwt(-1) (at 1.9micromolCdgdwt(-1)) at the highest cadmium dosage. Cadmium stimulated PC synthesis in E. prolifera which suggests that in N-rich algae, glutathione pools were high enough for PC synthesis. In both species GSH and protein increased with Zn contents, whereas GSH:(GSH+0.5GSSG) decreased, which would indicate Zn-induced oxidative stress; in E. linza, at the highest salinity the glutathione redox ratio decreased from 0.61 (at 2.9micromolZngdwt(-1)) to 0.26 (at 4.9nmolSgdwt(-1)) (at 0.5molCdgdwt(-1)). PCs were not synthesized in response to Zn, which may have resulted in Zn-induced GSH oxidation. The presence of both oxidative effects (Cd, Zn) and detoxification (Cd) could be identified by observing the responses of glutathione and PC pools to metal stress.     

The presolar grain inventory of fine‐grained chondrule rims in the Mighei‐type (CM) chondrites

We investigated the inventory of presolar silicate, oxide, and silicon carbide (SiC) grains of fine‐grained chondrule rims in six Mighei‐type (CM) carbonaceous chondrites (Banten, Jbilet Winselwan, Maribo, Murchison, Murray and Yamato 791198), and the CM‐related carbonaceous chondrite Sutter's Mill. Sixteen O‐anomalous grains (nine silicates, six oxides) were detected, corresponding to a combined matrix‐normalized abundance of ~18 ppm, together with 21 presolar SiC grains (~42 ppm). Twelve of the O‐rich grains are enriched in ¹⁷O, and could originate from low‐mass asymptotic giant branch stars. One grain is enriched in ¹⁷O and significantly depleted in ¹⁸O, indicative of additional cool bottom processing or hot bottom burning in its stellar parent, and three grains are of likely core‐collapse supernova origin showing enhanced ¹⁸O/¹⁶O ratios relative to the solar system ratio. We find a presolar silicate/oxide ratio of 1.5, significantly lower than the ratios typically observed for chondritic meteorites. This may indicate a higher degree of aqueous alteration in the studied meteorites, or hint at a heterogeneous distribution of presolar silicates and oxides in the solar nebula. Nevertheless, the low O‐anomalous grain abundance is consistent with aqueous alteration occurring in the protosolar nebula and/or on the respective parent bodies. Six O‐rich presolar grains were studied by Auger Electron Spectroscopy, revealing two Fe‐rich silicates, one forsterite‐like Mg‐rich silicate, two Al‐oxides with spinel‐like compositions, and one Fe‐(Mg‐)oxide. Scanning electron and transmission electron microscopic investigation of a relatively large silicate grain (490 nm × 735 nm) revealed that it was crystalline åkermanite (Ca₂Mg[Si₂O₇]) or a an åkermanite‐diopside (MgCaSi₂O₆) intergrowth.     

Evaluation of the Presence of Methane in Złoczew Lignite: Comparison with Other Lignite Deposits in Poland

Natural Resources Research - Catastrophic cases of methane explosion during exploratory drilling within the Be?chatów ortholignite deposit have led to testing for methane in other Polish...