Reactive transport modelling of Cr(VI) treatment by cast iron under fast flow conditions

The Cr(VI) reduction behaviour of five different types of grey cast iron shavings meant for construction of a permeable reactive barrier (PRB) in Thun, Switzerland, was investigated by performing batch and column experiments. Closed system batch experiments clearly demonstrated that the largest shavings also containing the largest spherical C inclusions, were best qualified for the PRB due to their fast Cr(VI) reduction rate. However, the column experiment performed with this type of material revealed that a complete and long term Cr(VI) treatment by the planned PRB was questionable due to the site-specific high groundwater flow velocities and the almost O₂ and CaCO₃ saturated aquifer conditions. The experimentally observed Cr(VI) breakthrough is explained as a result of a decline of reactivity provoked by a strong passivation by observed Fe hydroxides (FeOOH-polymorphs) and carbonates (calcite and aragonite).The column experiment was simulated using the reactive transport modelling code CrunchFlow. Iron cycling and intra-aqueous reactions were incorporated into a previously described reaction network in order to model the strong decline in reactivity of grey cast iron. All key parameters like aqueous species concentrations and mineralogical evolution of the column were successfully modelled. The modelling results confirmed that the observed Cr(VI) breakthrough was caused by surface passivation of the Fe shavings. Complete oxidation of the initially present mass of shavings is not predicted to occur during the expected PRB lifetime of several decades. The model seems to be robust, and it is expected that an application of the calibrated model in 2D to field sites will allow a quantitative evaluation for the performance of planned PRB’s in such environments.     

Passive image interferometry for monitoring crustal changes with ambient seismic noise

The technique to retrieve Green’s functions by cross-correlation of seismic noise that is now often referred to as seismic interferometry opened the way to a completely new type of seismological observation – the continuous monitoring of changes in the subsurface medium. In this paper we review Passive image interferometry (PII) as the technique that makes use of the possibility to repeatedly retrieve Green’s functions in order to precisely infer medium changes by comparing Green’s functions that are retrieved at different times. PII is capable of identifying velocity changes below 0.1% and was applied successfully in different tectonic situations. We review three investigations of PII that demonstrate the wide range of possible applications. The first application that introduced the concept of PII in 2006 revealed hydrology related changes of subsurface velocity at Merapi volcano. The second application we present was the first noise based observation of co-seismic changes. And we finally review a recent investigation that documented velocity variations in the shallow subsurface of the moon.     

Development and assessment of protective winter covers for marble statuaries of the Schlossbrücke, Berlin (Germany)

The present study documents the results of an inter-disciplinary model project that was planned with the aim of developing an innovative winter covering system for marble statuaries located on the Schlossbrücke (Berlin). Such a system would need to fulfil the various requirements for structural stability, aesthetics, climate and practical use. This applied research represents the first complex scientific study of the sustainability of a winter covering system. The study is characterised by the use of complex scientific instruments such as special laboratory analysis and numerical simulation tools. The interaction between the environment and the artefacts in connection with the innovative winter covering structures were studied by extensive climatic monitoring.     

Sulfate reduction below the sulfate-methane transition in Black Sea sediments

A sudden increase in salinity about 7000 years ago caused a shift in the deposited sediments of the Black Sea from limnic to brackish-marine. Due to the development of an anoxic deep water basin and a relatively high sulfate concentration, organic matter is mineralized primarily through sulfate reduction in modern Black Sea sediments. Earlier studies showed that sulfate-reducing bacteria are abundant within the limnic sub-surface sediment in spite of extremely low concentrations of sulfate and organic carbon. A main objective of the present study was therefore to understand the depth distribution of sulfate reduction across the different sediment units, even deep below the sulfate-methane transition. Our study combined experimental measurements of sulfate reduction using ³⁵S radiotracer with analyses of sulfur and iron geochemistry in pore water and sediment. Potential sulfate reduction rates were measured with ³⁵S in sediment samples that were amended with sulfate and organic substrates and incubated in time-series up to 48 h. Sulfate reduction could thereby be detected and quantified at depths where concentrations of sulfate were otherwise too low to enable calculation of the rates. The results demonstrate that sulfate-reducing bacteria are active several meters below the sulfate-methane transition in Black Sea sediments. The cryptic sulfate reduction below the sulfate-methane transition may be driven by sulfate produced from re-oxidation of reduced sulfur species with oxidized iron minerals buried in the deep limnic sediment.     

A High Order Perturbation Analysis of the Sitnikov Problem

The Sitnikov problem is one of the most simple cases of the elliptic restricted three body system. A massless body oscillates along a line (z) perpendicular to a plane (x,y) in which two equally massive bodies, called primary masses, perform Keplerian orbits around their common barycentre with a given eccentricity e. The crossing point of the line of motion of the third mass with the plane is equal to the centre of gravity of the entire system. In spite of its simple geometrical structure, the system is nonlinear and explicitly time dependent. It is globally non integrable and therefore represents an interesting application for advanced perturbative methods. In the present work a high order perturbation approach to the problem was performed, by using symbolic algorithms written in Mathematica. Floquet theory was used to derive solutions of the linearized equation up to 17th order in e. In this way precise analytical expressions for the stability of the system were obtained. Then, applying the Courant and Snyder transformation to the nonlinear equation, algebraic solutions of seventh order in z and e were derived using the method of Poincaré–Lindstedt. The enormous amount of necessary computations were performed by extensive use of symbolic programming. We developed automated and highly modularized algorithms in order to master the problem of ordering an increasing number of algebraic terms originating from high order perturbation theory.     

Iron deficiency in pyrrhotite of suevites from the Chesapeake Bay impact crater,USA—A consequence of shock metamorphism?

Abstract– Pyrrhotite from suevite of the 35 Ma Chesapeake Bay impact structure (CBIS) shows a shock metamorphism and we report on several mineralogical and magnetic features. Pyrrhotite shows strong brittle deformation with a high density of stacking faults, twinning parallel to the hexagonal (001) planes and average fault distances in the order of 10 nm. Although the determination of a superstructure was not possible due to the lattice defects, the reflections of the NiAs subcell, which is typical of all pyrrhotite modifications, were clearly detected. This phase is ferrimagnetic with a Curie temperature (T_C) between 350 and 365 °C, and suevite with this phase does not show the 34 K transition. The most peculiar feature is the low metal/sulfur ratio of 0.81, which indicates a distinctly higher vacancy concentration than for 4C pyrrhotite and a composition close to smythite (Fe₉S₁₁). This phase carries a stable natural remanent magnetization and is relatively hard magnetic. Steep inclinations of the natural remanent magnetization vector, however, suggest that this phase has been remagnetized by the drilling process. A possible explanation is the magnetic domain size of faultless areas of about 10 nm in diameter, which is at the lower limit of the single domain size near the threshold, below which superparamagnetic behavior occurs. The low thermal stability of this phase excludes postshock heating above 300 °C for the suevite of the CBIS. Our results imply that the iron‐deficient pyrrhotite is produced by shock metamorphism, although an iron loss due to shock has never been reported before for pyrrhotite.     

Advancing tsunami risk assessment by improving spatio-temporal population exposure and evacuation modeling

Tsunamis are among the most destructive and lethal of coastal hazards. These are time-specific events, and despite directly affecting a narrow strip of coastline, a single occurrence can have devastating effects and cause massive loss of life, especially in urbanized coastal areas. In this work, in order to consider the time dependence of population exposure to tsunami threat, the variation of spatio-temporal population distribution in the daily cycle is mapped and analyzed in the Lisbon Metropolitan Area. High-resolution daytime and nighttime population distribution maps are developed using ‘intelligent dasymetric mapping,’ that is, applying areal interpolation to combine best-available census data and statistics with land use and land cover data. Workplace information and mobility statistics are considered for mapping daytime distribution. In combination with a tsunami hazard map, information on infrastructure, land use and terrain slope, the modeled population distribution is used to assess people’s evacuation speed, applying a geospatial evacuation modeling approach to the city of Lisbon. The detailed dynamic population exposure assessment allows producing both daytime and nighttime evacuation time maps, which provide valuable input for evacuation planning and management. Results show that a significant amount of population is at risk, and its numbers increase dramatically from nighttime to daytime, especially in the zones of high tsunami flooding susceptibility. Also, full evacuation can be problematic in the daytime period, even if initiated immediately after a major tsunami-triggering earthquake. The presented approach greatly improves tsunami risk assessment and can benefit all phases of the disaster management process.     

Impact of climate change on groundwater recharge in a small catchment in the Black Forest, Germany

Temporal and spatial changes of the hydrological cycle are the consequences of climate variations. In addition to changes in surface runoff with possible floods and droughts, climate variations may affect groundwater through alteration of groundwater recharge with consequences for future water management. This study investigates the impact of climate change, according to the Special Report on Emission Scenarios (SRES) A1B, A2 and B1, on groundwater recharge in the catchment area of a fissured aquifer in the Black Forest, Germany, which has sparse groundwater data. The study uses a water-balance model considering a conceptual approach for groundwater-surface water exchange. River discharge data are used for model calibration and validation. The results show temporal and spatial changes in groundwater recharge. Groundwater recharge is progressively reduced for summer during the twenty-first century. The annual sum of groundwater recharge is affected negatively for scenarios A1B and A2. On average, groundwater recharge during the twenty-first century is reduced mainly for the lower parts of the valley and increased for the upper parts of the valley and the crests. The reduced storage of water as snow during winter due to projected higher air temperatures causes an important relative increase in rainfall and, therefore, higher groundwater recharge and river discharge.     

The IWAS-ToolBox: Software coupling for an integrated water resources management

Numerical modeling of interacting flow and transport processes between different hydrological compartments, such as the atmosphere/land surface/vegetation/soil/groundwater systems, is essential for understanding the comprehensive processes, especially if quantity and quality of water resources are in acute danger, like e.g. in semi-arid areas and regions with environmental contaminations. The computational models used for system and scenario analysis in the framework of an integrated water resources management are rapidly developing instruments. In particular, advances in computational mathematics have revolutionized the variety and the nature of the problems that can be addressed by environmental scientists and engineers. It is certainly true that for each hydro-compartment, there exists many excellent simulation codes, but traditionally their development has been isolated within the different disciplines. A new generation of coupled tools based on the profound scientific background is needed for integrated modeling of hydrosystems. The objective of the IWAS-ToolBox is to develop innovative methods to combine and extend existing modeling software to address coupled processes in the hydrosphere, especially for the analysis of hydrological systems in sensitive regions. This involves, e.g. the provision of models for the prediction of water availability, water quality and/or the ecological situation under changing natural and socio-economic boundary conditions such as climate change, land use or population growth in the future.