A groundwater salinity hotspot and its connection to an intermittent stream identified by environmental tracers (Mt Lofty Ranges,South Australia)

High and variable levels of salinity were investigated in an intermittent stream in a high-rainfall area (～800 mm/year) of the Mt. Lofty Ranges of South Australia. The groundwater system was found to have a local, upslope saline lens, referred to here as a groundwater salinity ‘hotspot’. Environmental tracer analyses (δ¹⁸O, δ²H, ^87/86Sr, and major elements) of water from the intermittent stream, a nearby permanent stream, shallow and deep groundwater, and soil-water/runoff demonstrate seasonal groundwater input of very saline composition into the intermittent stream. This input results in large salinity increases of the stream water because the winter wet-season stream flow decreases during spring in this Mediterranean climate. Furthermore, strontium and water isotope analyses demonstrate: (1) the upslope-saline-groundwater zone (hotspot) mixes with the dominant groundwater system, (2) the intermittent-stream water is a mixture of soil-water/runoff and the upslope saline groundwater, and (3) the upslope-saline-groundwater zone results from the flushing of unsaturated-zone salts from the thick clayey regolith and soil which overlie the metamorphosed shale bedrock. The preferred theory on the origin of the upslope-saline-groundwater hotspot is land clearing of native deep-rooted woodland, followed by flushing of accumulated salts from the unsaturated zone due to increased recharge. This cause of elevated groundwater and surface-water salinity, if correct, could be widespread in Mt. Lofty Ranges areas, as well as other climatically and geologically similar areas with comparable hydrogeologic conditions.     

Inverse modeling and uncertainty analysis of potential groundwater recharge to the confined semi-fossil Ohangwena II Aquifer,Namibia

The identification of potential recharge areas and estimation of recharge rates to the confined semi-fossil Ohangwena II Aquifer (KOH-2) is crucial for its future sustainable use. The KOH-2 is located within the endorheic transboundary Cuvelai-Etosha-Basin (CEB), shared by Angola and Namibia. The main objective was the development of a strategy to tackle the problem of data scarcity, which is a well-known problem in semi-arid regions. In a first step, conceptual geological cross sections were created to illustrate the possible geological setting of the system. Furthermore, groundwater travel times were estimated by simple hydraulic calculations. A two-dimensional numerical groundwater model was set up to analyze flow patterns and potential recharge zones. The model was optimized against local observations of hydraulic heads and groundwater age. The sensitivity of the model against different boundary conditions and internal structures was tested. Parameter uncertainty and recharge rates were estimated. Results indicate that groundwater recharge to the KOH-2 mainly occurs from the Angolan Highlands in the northeastern part of the CEB. The sensitivity of the groundwater model to different internal structures is relatively small in comparison to changing boundary conditions in the form of influent or effluent streams. Uncertainty analysis underlined previous results, indicating groundwater recharge originating from the Angolan Highlands. The estimated recharge rates are less than 1% of mean yearly precipitation, which are reasonable for semi-arid regions.     

An alternative approach to transverse and profile terrain curvature

Terrain curvature is one of the most important parameters of land surface topography. Well-established methods used in its measurement compute an index of plan or profile curvature for every single cell of a digital elevation model (DEM). The interpretation of these outputs may be delicate, especially when selected locations have to be analyzed. Furthermore, they involve a high level of simplification, contrasting with the complex and multiscalar nature of the surface curvature itself. In this paper, we present a new method to assess vertical transverse and profile curvature combining real-scale visualization and the possibility to measure these two terrain derivatives over a large range of scales. To this purpose, we implemented a GIS tool that extracts longitudinal and transverse elevation profiles from a high-resolution DEM. The performance of our approach was compared with some of the most commonly used methods (ArcMap, Redlands, CA, USA; ArcSIE, Landserf) by analyzing the terrain curvature around charcoal production sites in southern Switzerland. The different methods produced comparable results. While conventional methods quickly summarize terrain curvature in the form of a matrix of values, they involve a loss of information. The advantage of the new method lies in the possibility to measure and visualize the shape and size of the curvature, and to obtain a realistic representation of the average curvature for different subsets of spatial points. Moreover, the new method makes it possible to control the conditions in which the index of curvature is calculated.     

Fusion of multi‐resolution surface (terrestrial laser scanning) and subsurface geodata (ERT,SRT) for karst landform investigation and geomorphometric quantification

A multi‐method research design based on terrestrial laser scanning, GIS, geophysical prospecting (electrical resistivity tomography, refraction seismics) and sedimentology is applied for the first time to investigate enclosed karst depressions in an integrated way. Fusing multi‐resolution surface and subsurface geodata provides profound insights into the formation, geometry and geomorphologic processes of dolines. The studied landforms, which are located in the Dikti Mountains of East Crete, are shown to be filled by loose sediments of thicknesses of up to 30 m that mainly consist of fine‐grained material overlying solid bedrock at depths below 35 to 45 m. By combining subsurface observations with geomorphometric calculations, local doline genesis can be traced back to initial collapse of fractured bedrock followed by subsequent infilling with colluvials. In order to define crucial methodological requirements and guidelines for data fusion, both the impact of different elevation models and the influence of data resolution are assessed. Surface volumes of depressions derived by the digital surface model are 7–21% higher than the results obtained from the terrain model due to vegetation. Similarly, estimates of infill volume calculated on the basis of geophysical outcomes and elevation data differ by up to 13%. Calculations of the landforms' current volumes (i.e. total surface and subsurface volume), however, are fairly insensitive to raster resolution. Hence, the distinct geomorphologic properties of landforms (e.g. shape, terrain roughness, slope inclination) substantially determine the geomorphometric analysis of both surface and subsurface data. As shown by the findings, data fusion to integrate digital terrain, geophysical and sedimentological datasets of varied resolutions benefits geomorphologic studies and helps provide a comprehensive image of landforms. Copyright © 2013 John Wiley & Sons, Ltd.     

Uncertainties of geochemical modeling during CO<Subscript>2</Subscript> sequestration applying batch equilibrium calculations

One of the uncertainties in the field of carbon dioxide capture and storage (CCS) is caused by the parameterization of geochemical models. The application of geochemical models contributes significantly to calculate the fate of the CO₂ after its injection. The choice of the thermodynamic database used, the selection of the secondary mineral assemblage as well as the option to calculate pressure dependent equilibrium constants influence the CO₂ trapping potential and trapping mechanism. Scenario analyses were conducted applying a geochemical batch equilibrium model for a virtual CO₂ injection into a saline Keuper aquifer. The amount of CO₂ which could be trapped in the formation water and in the form of carbonates was calculated using the model code PHREEQC. Thereby, four thermodynamic datasets were used to calculate the thermodynamic equilibria. Furthermore, the equilibrium constants were re-calculated with the code SUPCRT92, which also applied a pressure correction to the equilibrium constants. Varying the thermodynamic database caused a range of 61% in the amount of trapped CO₂ calculated. Simultaneously, the assemblage of secondary minerals was varied, and the potential secondary minerals dawsonite and K-mica were included in several scenarios. The selection of the secondary mineral assemblage caused a range of 74% in the calculated amount of trapped CO₂. Correcting the equilibrium constants with respect to a pressure of 125 bars had an influence of 11% on the amount of trapped CO₂. This illustrates the need for incorporating sensitivity analyses into reaction pathway modeling.     

Using horizontal and vertical building structure to constrain indirect sensor orientation

This paper presents a method to integrate linear horizontal, vertical and right-angled scene structures into the bundle adjustment of image sequences. An increasing number of airborne image acquisition systems is available and equipped with non-metric small- or medium-frame cameras and no or insufficiently accurate INS devices. In cases where the data is to be used for the production of geo-spatial data, where a certain accuracy and precision is required, an indirect sensor orientation, possibly including self-calibration, needs to be performed. The idea which led to the presented approach is to reduce the number of GCPs necessary for this task by applying the mentioned scene structures. The method directly uses the linear structures, visible at man-made objects as fictive observations within the adjustment, while self-calibration of intrinsic camera parameters and lens distortion is included as well.Experiments with two datasets demonstrate that, through this method, only limited GCP information is required to obtain satisfactory results. In fact, in one experiment using oblique images, several scene constraints were provided and only the datum was defined by ground control. The residuals at check points from this setup were similar to the traditional case where several well-distributed GCPs were available in the scene. In the second experiment the ability of this approach to support the bundle adjustment was shown for a UAV dataset. Although no GCP and camera calibration information was available, the visual inspection of adjusted object points and the residuals at horizontal structures confirmed the ability of the method to align an image block with the structure, as embodied in the defined scene constraints. Despite the convincing outcome of the experiments, it needs to be mentioned that some manual work is still involved in defining the constraints. In future work the issue of automation will be addressed.     

Elimination of Carbamazepine,Diclofenac and Naproxen from Treated Wastewater by Nanofiltration

Most conventional wastewater treatment plants remove very small amounts of micropollutants, such as pharmaceuticals. Here, the ability of two different types of submerged nanofiltration flat sheet modules to remove pharmaceuticals from wastewater is analyzed. The two nanofiltration membranes were used at relatively low pressures of only 0.3 and 0.7 bar. At such low pressures, the membranes did not retain salts to a great extent. This is advantageous in wastewater treatment because no salt concentrate is produced. Carbamazepine was retained only slightly by the nanofiltration membranes, whereas approximately 60% of diclofenac and naproxen were retained by both membranes. This level of effectiveness might not be enough to justify the use of such a system as an additional treatment step in wastewater treatment plants.