A data exploration framework for validation and setup of hydrological models

Over the course of hydrological research projects often a large number of heterogeneous data sets are acquired from sources as diverse as boreholes, gauging stations or satellite imagery. This data then need to be integrated into models for the simulation of hydrological processes. We propose a framework for exploration of geoscientific data and visually guided preparation of such models. Data sets from a large number of sources can be imported, combined and validated to avoid potential problems due to artefacts or inconsistencies between data sets in a subsequent simulation. Boundary conditions and domain discretisations for surface and subsurface models can be created and tested regarding criteria indicating possible numerical instabilities. All data sets including simulation results can be integrated into a user-controlled 3D scene and aspects of the data can be enhanced using a number of established visualisation techniques including thresholding and user-defined transfer functions. We present the application of this framework for the preparation of a model for simulation of groundwater flow in a river catchment in southwest Germany investigated in the scope of the WESS project.     

Comments on “Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady-state” by O.A. Cirpka and A.J. Valocchi

In this comment we present a re-analysis of the analytical solution presented by Cirpka and Valocchi for steady-state concentrations of dissolved bioreactive compounds and bacterial biomass in porous media. We discuss the validity range of the analytical solution. In particular, the criterion used to determine the sustainability of biomass is revisited. This re-analysis shows that the ω$\omega$ criterion used by Cirpka and Valocchi is only a necessary but not a sufficient criterion to determine the bioreactive zones. As a consequence, the analytical solution does not provide the exact distribution of compounds throughout the domain, but can serve as upper or lower boundaries for species concentrations at a given location. These conclusions are supported by the simulation results obtained from an established reactive transport model.     

Factors affecting clarity of freshwater lakes in Brandenburg, Germany

In the German State Brandenburg, water clarity and the concentrations of the water quality components chlorophyll a, seston and gelbstoff were measured in 27 lakes. Correlation analysis showed, that spectral beam attenuation at 662 and 514 nm was mainly dependent on changes in chlorophyll a concentrations. In the UV-channel at 360 nm, beam attenuation depended mostly on gelbstoff.

Multiple linear regression provided a direct model of beam attenuation at 514 nm with the inputs of inorganic seston, chlorophyll a and gelbstoff. The specific beam attenuation coefficients were comparable to other natural waters around the world. An inverse model is presented, from which gelbstoff and chlorophyll a could be predicted with some accuracy from the inputs of beam attenuation coefficients at 514 and 360 nm. However, it became obvious that biological variability put major constraints on the predictive capacity of both the direct and the inverse model.

Furthermore, we observed a good correspondence of Secchi depth and the inverse of beam attenuation at 514 nm. The predictions of Secchi depth and chlorophyll a concentration from the inverse model were assessed in perspective of using this instrument instead of laborious chemical analysis for future trophic status classification according to LAWA (Länderarbeitsgemeinschaft Wasser). Predictions of trophic status were principally good when using calibrated models, however, quality of classification critically depended on predictions of chlorophyll a.     

Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions

Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.     

Hydrogeological investigation of shallow aquifers in an arid data-scarce coastal region (El Daba’a,northwestern Egypt)

Hydrogeological investigations in arid regions are particularly important to support sustainable development. The study area, El Daba’a in northwestern Egypt, faces scarce water resources as a result of reported climate change that particularly affects the southern Mediterranean coast and increases stress on the local groundwater reserves. This change in climate affects the area in terms of drought, over-pumping and unregulated exploration of groundwater for irrigation purposes. The hydrogeological investigation is based on a multidisciplinary data-layer analysis that includes geomorphology, geology, slope, drainage lines, soil type, structural lineaments, subsurface data, stable isotopes, and chemical analyses. The study area contains Pleistocene and middle Miocene marine limestone aquifers. Based on lithology and microfacies analysis, the middle Miocene aquifer is subdivided into two water-bearing zones. The area is affected by sets of faults and anticline folds, and these structures are associated with fractures and joints that increase permeability and facilitate the recharge of groundwater. Stable isotope data indicate that groundwater of both the Pleistocene and middle Miocene aquifers is recharged by modern precipitation. The high salinity values observed in some groundwater wells that tap both aquifers could be attributed to leaching and dissolution processes of marine salts from the aquifers’ marine limestone matrix. In addition, human activities can also contribute to an increase in groundwater salinity. A future water exploration strategy, based on the results from the multidisciplinary data-layer analysis, is proposed for the area. The derived scientific approach is transferable to other arid coastal areas with comparable conditions.     

Towards a multi-scale approach for an Earth observation-based assessment of natural resource exploitation in conflict regions

The exploitation of resources, if not properly managed, can lead to spoiling natural habitats as well as to threatening people’s health, livelihoods and security. The paper discusses a multi-scale Earth observation-based approach to provide independent information related to exploitation activities of natural resources for countries which are experiencing armed conflict. The analyses are based on medium to very high spatial resolution optical satellite data. Object-based image analysis is used for information extraction at these different scales. On a subnational level, conflict-related land cover changes as an indication of potential hot spots for exploitation activities are classified. The regional assessment provides information about potential activity areas of resource exploitation, whereas on a local scale, a site-specific assessment of exploitation areas is performed. The study demonstrates the potential of remote sensing for supporting the monitoring and documentation of natural resource exploitation in conflict regions.