Simulation of DNAPL infiltration and spreading behaviour in the saturated zone at varying flow velocities and alternating subsurface geometries

The influence of varying groundwater flow velocities on DNAPL infiltration and spreading behaviour was investigated by multiphase modelling using TMVOC and PetraSim. The multiphase models were calibrated by results of previously conducted laboratory experiments for the complete spatio-temporal range of the experiments. The small scale 2D scenario modelling was applied to qualify and quantify changes in position, architecture, geometry and dissolution of a TCE body in a fully saturated homogeneous sandy medium. The applied flow velocities ranging from 0.05 up to 40.00 m/day exhibited that the DNAPL TCE is affected even at the lowest flow velocity in its position, its size and its architecture. Additionally, several impermeable lenses with simple geometry were assumed in the model, to investigate the influence of stratified subsoil. In the experimental set-ups, the DNAPL body reacts more sensitive to the applied groundwater flow velocities than to the geometrical set-up of the scenarios. A possible consequence can be the transportation and displacement of a DNAPL pool due to natural or anthropogenic induced high groundwater flow velocities, as by Pump and Treat facilities, complicating site investigation process and planning of remediation activities.     

Investigating sediments and rock structures beneath a river using underwater ERT

This article presents hydrogeophysical investigations performed in a well-developed, long-term hydrogeological gypsum karst research site where subsurface evaporite dissolution has led to the subsidence of a river dam and an adjacent highway; both constructed on gypsum-containing rock, southeast of Basel, Switzerland. An observation system was set up to improve the protection of surface and subsurface water resources during remedial construction measures of the highway and in order to understand the processes, as well as the temporal evolution, of rock water interaction (flow and dissolution). However, no detailed hydrogeological information beneath the river could be derived from the previous investigations. To supplement the basic knowledge on this area, underwater Electrical Resistivity Tomography (ERT) measurements were conducted in the river bed upstream of the dam. The ERT-data are interpreted together with drill-core information and a conceptual 3D-Model of the area behind the dam and beneath the river. Results help to delineate weathered zones, associated faults and the thickness of sediment deposits behind the dam, as well as to locate voids within the local karst system. The combination of the ERT and modeling allows the optimization of future site-specific remedial construction measures.     

Reversed North Atlantic gyre dynamics in present and glacial climates

The dynamics of the North Atlantic subpolar gyre (SPG) are assessed under present and glacial boundary conditions by investigating the SPG sensitivity to surface wind-stress changes in a coupled climate model. To this end, the gyre transport is decomposed in Ekman, thermohaline, and bottom transports. Surface wind-stress variations are found to play an important indirect role in SPG dynamics through their effect on water-mass densities. Our results suggest the existence of two dynamically distinct regimes of the SPG, depending on the absence or presence of deep water formation (DWF) in the Nordic Seas and a vigorous Greenland?CScotland ridge (GSR) overflow. In the first regime, the GSR overflow is weak and the SPG strength increases with wind-stress as a result of enhanced outcropping of isopycnals in the centre of the SPG. As soon as a vigorous GSR overflow is established, its associated positive density anomalies on the southern GSR slope reduce the SPG strength. This has implications for past glacial abrupt climate changes, insofar as these can be explained through latitudinal shifts in North Atlantic DWF sites and strengthening of the North Atlantic current. Regardless of the ultimate trigger, an abrupt shift of DWF into the Nordic Seas could result both in a drastic reduction of the SPG strength and a sudden reversal in its sensitivity to wind-stress variations. Our results could provide insight into changes in the horizontal ocean circulation during abrupt glacial climate changes, which have been largely neglected up to now in model studies.     

Comparison of models calculating the sky view factor used for urban climate investigations

The sky view factor (SVF) describes the surface geometry and is a commonly used and important measure in urban climate investigations whose aim is the exploration of effects of a complex urban surface on climatological processes in built-up areas. A selection of methods and models for calculating the SVF was compared. For this purpose, fish eye images were taken at several locations in the city of Szeged, southern Hungary. The fish eye images equidistantly follow linear transects to cover a range of SVF values and to analyze the reaction of the methods to a continuously changing environment. The fish eye pictures were evaluated by three methods: the method according to Steyn (Atmos-Ocean 18(3):245?C258, 1980) implemented in a GIS-Script, the ??Edit free sky view factor?? tool of the RayMan model and BMSkyView. The SVF values at the coordinates of the fish eye pictures were calculated with three numerical models (SkyHelios, ArcView SVF extension, and SOLWEIG) with a 3D building data base as input. After comparing the results of the first run, a deviation occurs. The deviation disappears after implementing an option to include a weighting factor in some of the models.     

Total least-squares adjustment of condition equations

The usual least-squares adjustment within an Errors-in-Variables (EIV) model is often described as Total Least-Squares Solution (TLSS), just as the usual least-squares adjustment within a Random Effects Model (REM) has become popular under the name of Least-Squares Collocation (without trend). In comparison to the standard Gauss-Markov Model (GMM), the EIV-Model is less informative whereas the REM is more informative. It is known under which conditions exactly the GMM or the REM can be equivalently replaced by a model of condition equations or, more generally, by a Gauss-Helmert Model. Similar equivalency conditions are, however, still unknown for the EIV-Model once it is transformed into such a model of condition equations. In a first step, it is shown in this contribution how the respective residual vector and residual matrix look like if the TLSS is applied to condition equations with a random coefficient matrix to describe the transformation of the random error vector. The results are demonstrated using a numeric example which shows that this approach may be valuable in its own right.     

The chemical composition including the Rare Earth Elements of the three major glass types of Europe and the Orient used in late antiquity and the Middle Ages

Sets of 20 soda ash glasses, 16 soda lime glasses and 23 wood ash glasses mainly from excavations in Europe (additional soda ash glasses from Egypt) were analysed on 61 chemical elements. Average SiO₂ is about 62% in soda glasses and 50% in wood ash glasses. The three groups of glasses contain on average 13% Na₂O, 18% Na₂O and 13% K₂O as fluxes to lower the melting temperature of quartz at their production. The starting materials beside quartz were halophytic plant ash for soda ash glass, trona (Na₃H(CO₃)₂·2H₂O) and lime (clamshells) for soda lime glass and beech ash for wood ash glass. Each of the three major glass types contains specific Rare Earth Element (REE) concentrations mainly contained in quartz and its intergrown minerals. 50 Paleozoic and Mesozoic sandstones from Central Europe represent the quartz composition. The REE pattern of these glasses apparently indicates major compositional stages of the Continental Earth's Crust. The boron to lithium and sodium to potassium ratios as in seawater suggest reactions of materials for soda glass with seawater. Negative Ce anomalies in the three glasses are caused by reactions of quartz with seawater.     

Dimensional analysis of the earthquake-induced pounding between inelastic structures

In this paper the seismic response of inelastic structures with unilateral contact is revisited with dimensional analysis. All physically realizable contact types are captured via a non-smooth complementarity approach. The implementation of formal dimensional analysis leads to a condensed presentation of the response and unveils remarkable order even though two different types of non-linearity coexist in the response: the boundary non-linearity of unilateral contact and the inelastic behaviour of the structure itself. It is shown that regardless the intensity and frequency content of the excitation, all response spectra become self-similar when expressed in the appropriate dimensionless terms. The proposed approach hinges upon the notion of the energetic length scale of an excitation which measures the persistence of ground shaking to impose deformation demands. Using the concept of persistency which is defined for excitations with or without distinct pulses, the response is scaled via meaningful novel intensity measures: the dimensionless gap and the dimensionless yield displacement. The study confirms that contact may have a different effect on the response displacements of inelastic structures depending on the spectral region. In adjacent inelastic structures, such as colliding buildings or interacting bridge segments, contact is likely to alter drastically the excitation frequencies’ at which the system is most vulnerable. Finally, it is shown that the proposed approach yields maximum response displacements which correlate very well with the persistency of real earthquakes for a bridge system with considerably complex behaviour.     

Experimental study of the carbonation of partially serpentinized and weathered peridotites

Carbonation of partially serpentinized and weathered peridotites was studied experimentally under hydrothermal conditions (T: 200 °C, P_CO2: 130-180 bars). Experiments were performed in a closed system using whole rock drill core samples (height: 1 cm, diameter: 1 cm) as starting material. The initial samples were composed mainly of meshwork serpentine, relicts of primary olivine and an olivine weathering product (deweylite assemblage). Two types of solutions, each with a total salt content corresponding to that of average seawater (35 g/L dissolved salts), were used: (1) a Na-Ca-Cl solution (12.5 g/L CaCl₂ + 22.5 g/L NaCl) and (2) a NaCl solution (35 g/L NaCl). After 15-25 days of experimental treatment, the samples were partly covered with carbonates. In addition, noticeable carbonation reactions had occurred below the sample surfaces within zones with thicknesses up to 250 μm. In the Na-Ca-Cl solution, both the olivine relicts and the deweylite assemblage were partly replaced by calcite along the surrounding serpentine veins. However, the extent of calcitization was found to be considerably larger for the deweylite assemblage than for the olivine. Bulk fluid analyses show an increase in the Mg and Si concentrations with reaction time. In the NaCl solution, the deweylite assemblage was partly dissolved resulting in large voids within the reaction zone. In contrast, the olivine was replaced by magnesite. Under the conditions of our experiments, the meshwork serpentine was not reactive, but aided fluid infiltration into the rock samples. The experimentally produced microtextures closely resemble those found in natural examples. Our study elucidates the mechanisms by which carbonates form in ultramafic rocks under relatively high P_CO2-T conditions and particularly in the presence of Ca-bearing aqueous solutions. The existence of a serpentine meshtexture and the presence of weathering products formed from primary Mg-silicates may have significant beneficial effects on in situ CO₂ mineral sequestration in ultramafic rocks.     

Mass fluxes of xenobiotics below cities: challenges in urban hydrogeology

Urban areas are the focus of major ecological, social and economical activity. They are thus also prime locations of increasing conflict with regard to water use and water protection. As a direct and/or indirect consequence of urban land use and human activity, urban water systems are frequently polluted with organic contaminants including waste water-born xenobiotics such as pharmaceuticals, personal care products (collectively known as PPCPs) and endocrine-active substances. This study reviews new integrated methodologies including flux calculations as well as chemical investigations for determining the impact of human activities on urban water systems and on processes within the urban watershed. The use of indicator substances, representing different contaminant sources and pathways, integral pumping tests and mass balance approaches are suitable alternatives within these environments. The issues are explored using contaminant mass balance examples from Halle/Saale and Leipzig, Germany.     

The sinkhole enigma in the Alpine Foreland, Southeast Germany: Evidence of impact-induced rock liquefaction processes

Sudden collapse of the Quaternary soil to form sinkholes on the order of meters and tens of meters has been a geologic phenomenon within living memory in a localized area north of Lake Chiemsee in Southeast Germany. Failing a satisfying explanation, a relation with an undefined glaciation process has always been proposed. Excavations and geophysical measurements at three newly affected sites show underground features such as prominent sandy-gravelly intrusions and extrusions typical of rock liquefaction processes well known to occur during strong earthquakes. Since strong earthquakes can reasonably be excluded to have affected the area under discussion, it has been suggested that the observed widespread liquefaction is related with the recently proposed Holocene Chiemgau meteorite impact event. Except for one earlier proposed but unassertive relation between impact and liquefaction, the obviously direct association of both processes in the Chiemgau area emphasizes that observed paleoliquefaction features need not necessarily have originated solely from paleoseismicity but can provide a recognizable regional impact signature.