Correction to: Leveraging HPC accelerator architectures with modern techniques — hydrologic modeling on GPUs with ParFlow

Computational Geosciences - A Correction to this paper has been published: https://doi.org/10.1007/s10596-021-10065-y     

Assessing anthropogenic impacts on limited water resources under semi-arid conditions: three-dimensional transient regional modelling in Jordan

Both increasing aridity and population growth strongly stress freshwater resources in semi-arid areas such as Jordan. The country’s second largest governorate, Irbid, with over 1 million inhabitants, is already suffering from an annual water deficit of 25 million cubic meters (MCM). The population is expected to double within the next 20 years. Even without the large number of refugees from Syria, the deficit will likely increase to more then 50 MCM per year by 2035 The Governorate’s exclusive resource is groundwater, abstracted by the extensive Al Arab and Kufr Asad well fields. This study presents the first three-dimensional transient regional groundwater flow model of the entire Wadi al Arab to answer important questions regarding the dynamic quality and availability of water within the catchment. Emphasis is given to the calculation and validation of the dynamic groundwater recharge, derived from a multi-proxy approach, including (1) a hydrological model covering a 30-years dataset, (2) groundwater level measurements and (3) information about springs. The model enables evaluation of the impact of abstraction on the flow regime and the groundwater budget of the resource. Sensitivity analyses of controlling parameters indicate that intense abstraction in the southern part of the Wadi al Arab system can result in critical water-level drops of 10 m at a distance of 16 km from the production wells. Moreover, modelling results suggest that observed head fluctuations are strongly controlled by anthropogenic abstraction rather than variable recharge rates due to climate changes.     

Kinematic analysis and analogue modelling of the Passeier- and Jaufen faults: implications for crustal indentation in the Eastern Alps

Crustal deformation in front of an indenter is often affected by the indenter’s geometry, rheology, and motion path. In this context, the kinematics of the Jaufen- and Passeier faults have been studied by carrying out paleostress analysis in combination with crustal-scale analogue modelling to infer (1) their relationship during indentation of the Adriatic plate and (2) their sensitivity in terms of fault kinematics to the geometry and motion path of Adria. The field study reveals mylonites along the Jaufen fault, which formed under lower greenschist facies conditions and is associated with top-to-the-west/northwest shear with a northern block down component. In addition, a brittle reactivation of the Jaufen shear zone under NNW–SSE to NW–SE compressional and ENE–WSW tensional stress conditions was deduced from paleostress analysis. The inferred shortening direction is consistent with fission track ages portraying Neogene exhumation of the Meran-Mauls basement south of the fault. Along the Passeier fault, deformation was only brittle to semi-ductile and paleostress tensors record that the fault was subjected to E–W extension along its northern segment varying into NW–SE compression and sinistral transpression along its southern segment. In the performed analogue experiments, a rigid, triangular shaped indenter was pushed into a sand pile resulting in the formation of a Passeier-like fault sprouting from the indenter’s tip. These kinds of north-trending tip faults formed in all experiments with shortening directions towards the NW, N, or NE. Consequently, we argue that the formation of the Passeier fault strongly corresponds to the outline of the Adriatic indenter and was only little affected by the indenter’s motion path due to induced strain partitioning in front of the different indenter segments. The associated fault kinematics along the Passeier fault including both E–W extension and NNW to NW shortening, however, is most consistent with a northward advancing Adriatic indenter.     

Textural segmentation of digital rock images into bedding units using texture energy and cluster labels

The texture of digital rock images, as recorded, for instance, with borehole imaging devices, is shown to reflect different bedding types. Textural segmentation of borehole images, therefore, subdivides the recorded sequence into bedding units. We show that a textural segmentation algorithm based on the concept of texture energy achieves good results when compared with synthetic as well as real data in which petroleum geologists have performed zonations on cores. Texture energy involves filtering of the original image with a set of texture sensitive masks. The filtering is done as a finite convolution over the size of the masks. On the resulting images the variance is computed over a relatively large sliding window, which, in its practical implementation, covers the full width of the image. The resulting nine one-dimensional curves are then clustered hierarchically into a user-determined number of image texture or lithological bedding classes. Principal component analysis previous to clustering can be used to reduce redundancy in the data. A recurring and relatively ill-defined problem in this field are macro-textures, i.e., the cyclic interbedding of two or more bedding types. We show that sliding Fourier transforms and variable mask scale can successfully address the zonation of macro-textures. In general, the method gives best results with mask sizes equivalent to 2–4 centimeters, reflecting the length scale at which the investigated geological bedding seems to have its highest variation.     

A New Equation Solver for Modeling Turbulent Flow in Coupled Matrix‐Conduit Flow Models

Karst aquifers represent dual flow systems consisting of a highly conductive conduit system embedded in a less permeable rock matrix. Hybrid models iteratively coupling both flow systems generally consume much time, especially because of the nonlinearity of turbulent conduit flow. To reduce calculation times compared to those of existing approaches, a new iterative equation solver for the conduit system is developed based on an approximated Newton–Raphson expression and a Gauß–Seidel or successive over‐relaxation scheme with a single iteration step at the innermost level. It is implemented and tested in the research code CAVE but should be easily adaptable to similar models such as the Conduit Flow Process for MODFLOW‐2005. It substantially reduces the computational effort as demonstrated by steady‐state benchmark scenarios as well as by transient karst genesis simulations. Water balance errors are found to be acceptable in most of the test cases. However, the performance and accuracy may deteriorate under unfavorable conditions such as sudden, strong changes of the flow field at some stages of the karst genesis simulations.     

Emerging Technologies and Synergies for Airborne and Space-Based Measurements of Water Vapor Profiles

A deeper understanding of how clouds will respond to a warming climate is one of the outstanding challenges in climate science. Uncertainties in the response of clouds, and particularly shallow clouds, have been identified as the dominant source of the discrepancy in model estimates of equilibrium climate sensitivity. As the community gains a deeper understanding of the many processes involved, there is a growing appreciation of the critical role played by fluctuations in water vapor and the coupling of water vapor and atmospheric circulations. Reduction of uncertainties in cloud-climate feedbacks and convection initiation as well as improved understanding of processes governing these effects will result from profiling of water vapor in the lower troposphere with improved accuracy and vertical resolution compared to existing airborne and space-based measurements. This paper highlights new technologies and improved measurement approaches for measuring lower tropospheric water vapor and their expected added value to current observations. Those include differential absorption lidar and radar, microwave occultation between low-Earth orbiters, and hyperspectral microwave remote sensing. Each methodology is briefly explained, and measurement capabilities as well as the current technological readiness for aircraft and satellite implementation are specified. Potential synergies between the technologies are discussed, actual examples hereof are given, and future perspectives are explored. Based on technical maturity and the foreseen near-mid-term development path of the various discussed measurement approaches, we find that improved measurements of water vapor throughout the troposphere would greatly benefit from the combination of differential absorption lidar focusing on the lower troposphere with passive remote sensors constraining the upper-tropospheric humidity.     

3-D assessment of peak-metamorphic conditions by Raman spectroscopy of carbonaceous material: an example from the margin of the Lepontine dome (Swiss Central Alps)

This study monitors regional changes in the crystallinity of carbonaceous matter (CM) by applying Micro-Raman spectroscopy to a total of 214 metasediment samples (largely so-called Bündnerschiefer) dominantly metamorphosed under blueschist- to amphibolite-facies conditions. They were collected within the northeastern margin of the Lepontine dome and easterly adjacent areas of the Swiss Central Alps. Three-dimensional mapping of isotemperature contours in map and profile views shows that the isotemperature contours associated with the Miocene Barrow-type Lepontine metamorphic event cut across refolded nappe contacts, both along and across strike within the northeastern margin of the Lepontine dome and adjacent areas. Further to the northeast, the isotemperature contours reflect temperatures reached during the Late Eocene subduction-related blueschist-facies event and/or during subsequent near-isothermal decompression; these contours appear folded by younger, large-scale post-nappe-stacking folds. A substantial jump in the recorded maximum temperatures across the tectonic contact between the frontal Adula nappe complex and surrounding metasediments indicates that this contact accommodated differential tectonic movement of the Adula nappe with respect to the enveloping Bündnerschiefer after maximum temperatures were reached within the northern Adula nappe, i.e. after Late Eocene time.     

Effects of the peat acid digestion protocol on geochemically and morphologically diverse tephra deposits

Tephra shards for electron probe microanalysis are most efficiently extracted from peat using acid digestion, which removes organic material that hinders density separation methods. However, strong acids are known to alter glass chemical compositions, and several studies have examined how acid digestion affects rhyolitic volcanic glass. The focus on rhyolitic tephra in these studies leaves considerable uncertainty, as the dissolution rates of natural glasses (including tephra) are determined by the chemical composition and surface area/volume ratio, both of which vary in tephra deposits. Here, we use duplicate samples of basaltic, trachydacitic and rhyolitic tephra to examine physical and geochemical alteration following acid digestion. Scanning electron microscope imagery reveals no discernible degradation of glass surfaces, and electron probe microanalysis results from duplicate samples are statistically indistinguishable. These findings suggest the acid digestion protocol for organic peats does not significantly alter glass geochemistry regardless of shard morphologies or geochemical compositions.