Anisotropic dual-continuum representations for multiscale poroelastic materials: Development and numerical modelling

Dual-continuum (DC) models can be tractable alternatives to explicit approaches for the numerical modelling of multiscale materials with multiphysics behaviours. This work concerns the conceptual and numerical modelling of poroelastically coupled dual-scale materials such as naturally fractured rock. Apart from a few exceptions, previous poroelastic DC models have assumed isotropy of the constituents and the dual-material. Additionally, it is common to assume that only one continuum has intrinsic stiffness properties. Finally, little has been done into validating whether the DC paradigm can capture the global poroelastic behaviours of explicit numerical representations at the DC modelling scale. We address the aforementioned knowledge gaps in two steps. First, we utilise a homogenisation approach based on Levin's theorem to develop a previously derived anisotropic poroelastic constitutive model. Our development incorporates anisotropic intrinsic stiffness properties of both continua. This addition is in analogy to anisotropic fractured rock masses with stiff fractures. Second, we perform numerical modelling to test the DC model against fine-scale explicit equivalents. In doing, we present our hybrid numerical framework, as well as the conditions required for interpretation of the numerical results. The tests themselves progress from materials with isotropic to anisotropic mechanical and flow properties. The fine-scale simulations show that anisotropy can have noticeable effects on deformation and flow behaviour. However, our numerical experiments show that the DC approach can capture the global poroelastic behaviours of both isotropic and anisotropic fine-scale representations.     

Geometry and chronology of growth and drowning of Middle Triassic carbonate platforms (Cernera and Bivera/Clapsavon) in the Southern Alps (northern Italy)

The depositional architecture and the geometric relationships between platform-slope deposits and basinal sediments along with paleontological evidence indicate the time interval of the younger Anisian Reitziites reitzi ammonoid zone to largely represent the main stage of platform aggradation at the Cernera and Bivera/Clapsavon carbonate platforms. Published and new U-Pb age data of zircons from volcaniclastic layers bracketing the stratigraphic interval of platform growth constrain the duration of platform evolution to a time span shorter than 1.8±0.7m.y., probably in the order of 0.5-1m.y., reflecting fast rates of vertical platform aggradation exceeding 500 m/m.y. In the range of growth potentials for shallow-water carbonate systems estimated in relation to the time span of observation, this high rate is in agreement with values for short intervals of 10⁵-10⁶yrs (e.g., Schlager 1999). After drowning, the platforms at Cernera and Bivera/Clapsavon were blanketed by thin pelagic carbonates. On the former platform flanks the draping sediments in places comprise red nodular pelagic limestones (Clapsavon Limestone) similar in facies to the Han Bulog Limestones occurring elsewhere in Middle Triassic successions of the Mediterranean Tethys. The drowning of vast areas of former carbonate platforms possibly triggered the onset of bottom-water circulation in adjacent basins as suggested by the abrupt transition from laminated to bioturbated pelagic nodular limestones in the Buchenstein Formation which occurred close to the time of initial platform submergence. During the Late Ladinian the topographic features of the drowned platforms were onlapped by rapidly deposited, predominantly clastic successions including coarse breccias and volcanic rocks sealing and preserving the peculiar stratigraphic setting.     

Phytoreduction and volatilization of mercury by ascorbate in Arabidopsis thaliana, European beech and Norway spruce

Mercury vapor (Hg⁰) emission from plants contributes to the atmospheric Hg cycle. Young barley (Hordeum vulgare L.) plants grown on a hydroponic cultivation medium containing Hg(II) have previously been shown to increase their Hg⁰ emission significantly by reduction of Hg(II) with endogenous ascorbic acid. Regarding the potential contribution to the Hg cycle from the vast forest-covered areas, it was important to investigate this mechanism in trees. The increase in Hg⁰ emission from young European beech plants cultivated on a HgCl₂ medium exceeded that from controls by ca. tenfold and was proportional to the Hg(II) concentration. From these experiments, a flux of 12.8 μg Hg⁰/h/m² was estimated at an exposure of the roots to 20 μM Hg(II). Mercury vapor release from homogenates of Norway spruce needles exceeded that from European beech leaves by a factor of 2.3–4, i.e. in proportion to the reported AA concentrations; the reduction was maximal at alkaline pH which is typical for AA. The 8.4-fold difference in Hg⁰ release between homogenates from wild-type Arabidopsis thaliana and from its AA-deficient mutant vtc 1-1 also paralleled the reported difference in AA levels of both species. It is concluded that the phytoreduction and vaporization of Hg by AA is an important mechanism as much for Hg detoxification in trees as for Hg emission to the atmosphere. The efficiency of this process seems to result from the optimal coordination of transfer and biochemical transformation of mercuric ions and Hg vapor. There is no evidence for a relevant difference in the mechanisms of biogenic Hg(II) reduction between grass plants and trees.     

Calcium isotope (δCa) fractionation along hydrothermal pathways, Logatchev field (Mid-Atlantic Ridge, 14°45′N)

We investigate the Logatchev Hydrothermal Field at the Mid-Atlantic Ridge, 14°45′N to constrain the calcium isotope hydrothermal flux into the ocean. During the transformation of seawater to a hydrothermal solution, the Ca concentration of pristine seawater ([Ca]_SW) increases from about 10 mM to about 32 mM in the hydrothermal fluid endmember ([Ca]_HydEnd) and thereby adopts a δ^44/40Ca_HydEnd of −0.95 ± 0.07‰ relative to seawater (SW) and a ⁸⁷Sr/⁸⁶Sr isotope ratio of 0.7034(4). We demonstrate that δ^44/40Ca_HydEnd is higher than that of the bedrock at the Logatchev field. From mass balance calculations, we deduce a δ^44/40Ca of −1.17 ± 0.04‰ (SW) for the host-rocks in the reaction zone and −1.45 ± 0.05‰ (SW) for the isotopic composition of the entire hydrothermal cell of the Logatchev field. The values are isotopically lighter than the currently assumed δ^44/40Ca for Bulk Earth of −0.92 ± 0.18‰ (SW) [Skulan J., DePaolo D. J. and Owens T. L. (1997) Biological control of calcium isotopic abundances in the global calcium cycle. Geochim. Cosmochim. Acta61,(12) 2505-2510] and challenge previous assumptions of no Ca isotope fractionation between hydrothermal fluid and the oceanic crust [Zhu P. and Macdougall J. D. (1998) Calcium isotopes in the marine environment and the oceanic calcium cycle. Geochim. Cosmochim. Acta62,(10) 1691-1698; Schmitt A. -D., Chabeaux F. and Stille P. (2003) The calcium riverine and hydrothermal isotopic fluxes and the oceanic calcium mass balance. Earth Planet. Sci. Lett. 6731, 1-16]. Here we propose that Ca isotope fractionation along the fluid flow pathway of the Logatchev field occurs during the precipitation of anhydrite. Two anhydrite samples from the Logatchev Hydrothermal Field show an average fractionation of about Δ^44/40Ca = −0.5‰ relative to their assumed parental solutions. Ca isotope ratios in aragonites from carbonate veins from ODP drill cores indicate aragonite precipitation directly from seawater at low temperatures with an average δ^44/40Ca of −1.54 ± 0.08‰ (SW). The relatively large fractionation between the aragonite precipitates and seawater in combination with their frequent abundance in weathered mafic and ultramafic rocks suggest a reconsideration of the marine Ca isotope budget, in particular with regard to ocean crust alteration.     

Groundwater Evolution and Mineral Alteration Reactions in the Basaltic Rock Sequence of Mt. Wasserkuppe, Germany: A Case Study

Groundwater sampling was accomplished in the basaltic sequence of the Rh?n mountain range, Germany, in order to investigate hydrochemical groundwater evolution and to delineate mineral alteration reactions involved in natural weathering. The hydrochemical compositions of near-surface groundwaters indicate a Ca/Mg–HCO₃ type with near-neutral pH and evolve to a Na–HCO₃ type with high pH at greater depth. Column experiments were performed with basaltic and phonolitic rock samples to determine individual mineral alteration reactions. The basic reactions could be related to the alteration of olivine, Ca-pyroxene, plagioclase, pyrrhotite, and feldspathoids under formation of secondary clay minerals (smectites, illite) and goethite. The mineral alteration reactions deduced from the leaching experiments by inverse modelling were found to be consistent with the mineral reactions associated with the natural groundwaters. The reactions calculated for groundwater evolution involve the alteration of primary and secondary minerals to produce low-T mineral phase. The conversion of secondary Na-beidellite to illite occurs at a later stage of groundwater evolution, reducing the concentrations of K⁺ and Mg²⁺. Near-surface groundwaters do not indicate significant cation exchange. Initial cation exchange requires elevated pH values, with Mg²⁺ removed from solution preferred to Ca²⁺. Na-alkalisation of the groundwaters at greater depth suggests the exchange of Na⁺ for Mg²⁺ and Ca²⁺ on Na-beidellite, supported by cation exchange on coatings of iron hydroxides as alteration products. Among the mature high-pH groundwater at greater depth, the dissolution of anorthite and albite has significant effect on groundwater composition.     

Trace element diagenesis in pyrite-rich sediments of the Achterwasser lagoon, SW Baltic Sea

The early diagenesis of trace elements (V, Cr, Co, Cu, Zn, As, Cd, Ba, U) in anoxic sediments of the Achterwasser, a shallow lagoon in the non-tidal Oder estuary in the Baltic Sea, was investigated in the context of pyrite formation. The dissolved major redox parameters show a two-tier distribution with transient signals in the occasionally re-suspended fluid mud layer (FM) and a permanently established diagenetic sequence in the sediment below. Intense microbial respiration leads to rapid depletion of O₂ within the uppermost mm of the FM. The reduction zones of Mn, Fe and sulfate overlap in the FM and in the permanently anoxic sediment section which appears to be a typical feature of estuarine sediments, under low-sulfate conditions. Degrees of pyritization (DOP) range from 50% in the FM to remarkably high values > 90% at 50 cm depth. Pyrite formation at the sediment surface is attributed to the reaction of Fe-monosulfides with intermediate sulfur species via the polysulfide pathway. By contrast, intense pyritization in the permanently anoxic sediment below is attributed to mineral growth via adsorption of aqueous Fe-sulfide complexes onto pyrite crystals which had originally formed in the surface layer.The studied trace elements show differential behavior patterns which are closely coupled to the diagenetic processes described above: (i) Zn, Cu and Cd are liberated from organic matter in the thin oxic layer of the sediment and diffuse both upwards across the sediment/water boundary and downwards to be trapped as monosulfides, (ii) V, Cr, Co and As are released during reductive dissolution of Mn- and Fe-oxyhydroxides, (iii) U removal from pore water occurs concomitantly to Fe reduction in the FM and is attributed to reduction of U(VI) to U(IV), (iv) the Ba distribution is controlled by reductive dissolution of authigenic barite in the sulfate reduction zone coupled with upward diffusion and re-precipitation. The incorporation of trace elements into pyrite is most intense for Co, Mn and As, intermediate for Cu and Cr and little to negligible for U, Zn, Cd, V and Ba. The observed trend is largely in agreement with previous studies and may be explained with differing rates for ligand exchange. Slow and fast ligand exchange and thus precipitation kinetics are also displayed by downcore increasing (Mn, Cr, Co and As) or constantly low (Zn, Cu, Cd) pore water concentrations. The downward increasing degrees of trace metal pyritization (DTMP) for Co, Cu, Zn and As are, in analogy to pyrite growth, assigned to adsorption of sulfide complexes or As oxyanions onto preexisting pyrite minerals.     

Altimetry-derived tide model for improved tide and water level forecasting along the European continental shelf

Ocean Dynamics - With the continued rise in global mean sea level, operational predictions of tidal height and total water levels have become crucial for accurate estimations and understanding of...     

Causes and consequences of long-term groundwater overabstraction in Jordan

Hydrogeology Journal - In 2017, a comprehensive review of groundwater resources in Jordan was carried out for the first time since 1995. The change in groundwater levels between 1995 and 2017 was...     

Cartograms for Use in Forecasting Weather-Driven Natural Hazards

This paper evaluates the potential of using cartograms for visualizing and interpreting forecasts of weather-driven natural hazards in the context of global weather forecasting and early warning systems. The use of cartograms is intended to supplement traditional cartographic representations of the hazards in order to highlight the severity of an upcoming event. Cartogrammetric transformations are applied to forecasts of floods, heatwaves, windstorms and snowstorms taken from the European Centre for Medium-range Weather Forecasts (ECMWF) forecast archive. Key cartogram design principles in standard weather forecast visualization are tested. Optimal cartogram transformation is found to be dependent on geographical features (such as coastlines) and forecast features (such as snowstorm intensity). For highly spatially autocorrelated weather variables used in analysing several upcoming natural hazards such as 2m temperature anomaly, the visualization of the distortion provides a promising addition to standard forecast visualizations for highlighting upcoming weather-driven natural hazards.