Geometric modelling and object-oriented software concepts applied to a heterogeneous fractured network from the Grimsel rock laboratory

Discrete fracture network simulations are computationally intensive and usually time-consuming to construct and configure. This paper presents a case study with techniques for building a 3D finite element model of an inhomogeneous fracture network for modelling flow and tracer transport, combining deterministic and stochastic information on fracture aperture distributions. The complex intersected fractures represent a challenge for geometrical model design, mesh quality requirements and property allocations. For the integrated and holistic modelling approach, including the application of numerical and analytical simulation techniques, new object-oriented concepts in software engineering are implemented to ensure a resourceful and practicable software environment.     

The electric field gradient in natural iron-doped chrysoberyl Al<Subscript>2</Subscript>BeO<Subscript>4</Subscript> and sinhalite MgAlBO<Subscript>4</Subscript> single crystals

This work reports on the evaluation of the electric field gradient (EFG) in natural chrysoberyl Al₂BeO₄ and sinhalite MgAlBO₄ using two different procedures: (1) experimental, with single crystal Mössbauer spectroscopy (SCMBS) on the three principal sections of each sample and (2) a “fully quantitative” method with cluster molecular orbital calculations based on the density functional theory. Whereas the experimental and theoretical results for the EFG tensor are in quantitative agreement, the calculated isomer shifts and optical d–d-transitions exhibit systematic deviations from the measured values. These deviations indicate that the substitution of Al and Mg with iron should be accompanied by considerable local expansion of the coordination octahedra.     

Sauerstoff-Isotopenuntersuchungen an Mineralen eines metamorphen Profils der Hohen Tauern, Österreich

Oxygen isotope analyses on minerals of a metamorphic profile from Hohe Tauern, Austria indicate that the oxygen isotopes of the coexisting minerals quartz, garnet, and biotite have been reequilibrated during alpidic metamorphism. The K/Ar-ages of these biotites are in the range of 20–30 my. Metamorphic temperatures can be calculated from the quartz-biotite and quartz-garnet isotopic fractionations. The metamorphic temperature in the northern part of the region under investigation is 440° C and increases in direction to Tauernhauptkamm to about 600° C. Indicated by petrographic investigations, two samples from this area show no remarkable metamorphic influence. K/Ar-determinations on hornblendes and biotites yielded variscan ages.The isotopic temperatures of these samples, calculated by the oxygen isotope fractionations of quartz-biotite, quartz-magnetite, and feldspar-magnetite are equal within the limits of instrumental error and represent the crystallisation temperature of a granodioritic magma of 660–670° C.     

A Critical Analysis of Transverse Dispersivity Field Data

Transverse dispersion, or tracer spreading orthogonal to the mean flow direction, which is relevant e.g, for quantifying bio-degradation of contaminant plumes or mixing of reactive solutes, has been studied in the literature less than the longitudinal one. Inferring transverse dispersion coefficients from field experiments is a difficult and error-prone task, requiring a spatial resolution of solute plumes which is not easily achievable in applications. In absence of field data, it is a questionable common practice to set transverse dispersivities as a fraction of the longitudinal one, with the ratio 1/10 being the most prevalent. We collected estimates of field-scale transverse dispersivities from existing publications and explored possible scale relationships as guidance criteria for applications. Our investigation showed that a large number of estimates available in the literature are of low reliability and should be discarded from further analysis. The remaining reliable estimates are formation-specific, span three orders of magnitude and do not show any clear scale-dependence on the plume traveled distance. The ratios with the longitudinal dispersivity are also site specific and vary widely. The reliability of transverse dispersivities depends significantly on the type of field experiment and method of data analysis. In applications where transverse dispersion plays a significant role, inference of transverse dispersivities should be part of site characterization with the transverse dispersivity estimated as an independent parameter rather than related heuristically to longitudinal dispersivity.     

Validation of turbulence closure parameterisations for stably stratified flows using the PROVESS turbulence measurements in the North Sea

A number of parameterisations for the simulation of mixing processes in the thermocline are compared and tested against the microstructure data of the PROVESS campaigns, conducted in the northern part of the North Sea during the autumn of 1998. The transport term in the turbulent kinetic energy equation is parameterised via the introduction of a third stability function S_k for turbulent energy diffusion. The formulations are compared with a simpler scheme based upon limiting conditions for turbulence variables. Improved results are obtained with a new form of S_k. The best agreement is, however, found with the simpler limiting scheme. This is explained in terms of a turbulence length scale theory for stably stratified turbulence. In agreement with previous laboratory and ocean data it is found that the ratios of the Thorpe and Kolmogorov scales to the Ozmidov length scale approach critical limiting values in the thermocline. The first of these conditions is satisfied when limiting conditions are implemented into the scheme, providing the necessary minimum value for the dissipation rate, whereas the schemes without limiting conditions fail to produce this critical ratio. The basic reason for this failure is that the Thorpe scale is overestimated, which is shown to be connected to an even larger overprediction of the dissipation rate of temperature variance. To investigate the impact of non-resolved advective processes and salinity stratification on the turbulence predictions, additional numerical experiments were conducted using a simple scheme for data assimilation. The best agreement is found again with the limiting scheme, which is able to make reasonable predictions for the dissipation rate without knowing the detailed shape of the mean stratification profile. It is shown that advective transport due to tidally and wind-driven motions has a non-negligible impact on vertical mixing. This is seen in the data and the models by periodic enhancements of turbulent mixing inside the thermocline.     

Protracted river response to medieval earthquakes

Mountain rivers respond to strong earthquakes by rapidly aggrading to accommodate excess sediment delivered by co-seismic landslides. Detailed sediment budgets indicate that rivers need several years to decades to recover from seismic disturbances, depending on how recovery is defined. We examine three principal proxies of river recovery after earthquake-induced sediment pulses around Pokhara, Nepal's second largest city. Freshly exhumed cohorts of floodplain trees in growth position indicate rapid and pulsed sedimentation that formed a fan covering 150 km² in a Lesser Himalayan basin with tens of metres of debris between the 11th and 15th centuries AD. Radiocarbon dates of buried trees are consistent with those of nearby valley deposits linked to major medieval earthquakes, such that we can estimate average rates of re-incision since. We combine high-resolution digital elevation data, geodetic field surveys, aerial photos, and dated tree trunks to reconstruct geomorphic marker surfaces. The volumes of sediment relative to these surfaces require average net sediment yields of up to 4200 t km^–2 yr^–1 for the 650 years since the last inferred earthquake-triggered sediment pulse. The lithological composition of channel bedload differs from that of local bedrock, confirming that rivers are still mostly evacuating medieval valley fills, locally incising at rates of up to 0.2 m yr^–1. Pronounced knickpoints and epigenetic gorges at tributary junctions further illustrate the protracted fluvial response; only the distal portions of the earthquake-derived sediment wedges have been cut to near their base. Our results challenge the notion that mountain rivers recover speedily from earthquakes within years to decades. The valley fills around Pokhara show that even highly erosive Himalayan rivers may need more than several centuries to adjust to catastrophic perturbations. Our results motivate some rethinking of post-seismic hazard appraisals and infrastructural planning in active mountain regions. © 2018 John Wiley & Sons, Ltd.     

Solid–liquid–air coupling in multiphase porous media

This paper addresses various issues concerning the modelling of solid–liquid–air coupling in multiphase porous media with an application to unsaturated soils. General considerations based on thermodynamics permit the derivation and discussion of the general form of field equations; two cases are considered: a three phase porous material with solid, liquid and gas, and a two phase porous material with solid, liquid and empty space. Emphasis is placed on the presentation of differences in the formulation and on the role of the gas phase. The finite element method is used for the discrete approximation of the partial differential equations governing the problem. The two formulations are then analysed with respect to a documented drainage experiment carried out by the authors. The merits and shortcomings of the two approaches are shown. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling karst denudation on a synthetic landscape

This contribution presents the results of a numerical study of karst denudation on limestone plateaux. The landscape evolution model used incorporates not only long-range fluvial and short-range hill-slope processes, but also large-scale chemical dissolution of limestone surfaces. The relative efficiencies of fluvial and chemical processes are of equal importance to the landscape evolution of a plateau dropping to sea level along an escarpment. While fluvial processes have an impact confined mostly to river channels, the karst denudation process is more uniform, removing material also from the plateau surface. The combined effect of both processes results in a landscape evolution almost twice as effective as the purely erosional evolution of an insoluble landscape.     

In situ cone penetration tests at the active Dashgil mud volcano,Azerbaijan: Evidence for excess fluid pressure,updoming, and possible future violent eruption

Active mud volcanism is a global phenomenon that represents a natural hazard by self-igniting eruptions and the continuous emission of methane gas in both marine and continental settings. Mud domes are often found in compressional tectonic settings such as the Caucasus orogenic wedge. Dashgil mud volcano, the most prominent of >200 features in Azerbaijan, has erupted vigorously in historic times. For several years, we have observed variations in the activity of Dashgil dome, including transients in methane flux, build-up of extrusive mud cones on the main feature, and flexural polygonal cracks adjacent to the main crater lake and new mud cones. In spring 2007, we carried out in situ CPTU (Cone Penetration Testing with Pore Pressure measurement) experiments in the crestal area of Dashgil. Our data suggest that the central portion of the crater lake, which hosts the conduit for gas (and possible mud) ascent, shows both low sediment shear strength (<5–20 kPa) and excess pore fluid pressures between 15 and 30 kPa supra-hydrostatic at 1 m sub-bottom depth. In situ cone resistance as a measure for undrained shear strength is as low as 150 kPa in the conduit, whereas the mud is found rather stiff in all other testing locations (300–700 kPa, probably a result of deeply buried shales of the Maikop formation parts of which now liquefy and ascend). Pore pressure is low in the centre of the conduit, probably because of rapidly migrating gas. It increases to 30 kPa at the lake bottom and deep flank, then decreases upslope on the lake flank, and reaches hydrostatic values at the crater rim. From the overpressured region beneath the fluid-filled crest of Dashgil dome, combined with the other observations, we suspect to currently witness an ongoing period of updoming. The presence of sintered mudstones from explosive eruptions in 1908 and 1928 (and most likely before) suggests that a similar violent activity may occur in the near future.     

Hydraulic zonation of the lagoons of Marano and Grado,Italy. A modelling approach

The hydraulic regime-based zonation scheme of the Lagoons of Marano and Grado (Italy) has been derived by means of numerical models. A finite element modelling system has been used to describe the water circulation taking in account different forces such as tide, wind and rivers. The model has been validated by comparing the simulation results against measured water levels, salinity and water temperature data collected in several stations inside the lagoons. The analysis of water circulation, salinity and spatial distribution of passive tracers released at the inlets, led to a physically-based division of the lagoons system into six subbasins. The derived classification scheme is of crucial value for understanding the renewal capacity and pollutants distribution patterns in the lagoon.