Abstract: When Fe(II) bearing groundwaters surface in streams, particulate authigenic Fe-rich material is produced by oxidation. Such freshly precipitated Fe minerals may be transported as suspended sediment and have a profound impact on the fate of trace metals and nutrients in rivers. The objective of this study was to monitor changes in mineralogy and composition of authigenic material from its source to streams of increasing order. Groundwaters, surface waters, and suspended sediment in streams of different order were sampled in the Kleine Nete catchment (Belgium), a lowland with Fe-rich groundwaters (3.5-53.8 mg Fe/L; pH 6.3-6.9). Fresh authigenic material (>0.45 mu m) was produced by oxidising filtered (<0.45 mu m) groundwater and surface water. This material contained, on average, 44% Fe, and smaller concentrations of C, P, and Ca. Iron EXAFS (Extended X-ray Absorption Fine Structure) spectroscopy showed that the Fe was present as poorly crystalline hydrous ferric oxides with a structure similar to that of ferrihydrite. The Fe concentration in the suspended sediment samples decreased to 36-40% (stream order 2), and further to 18-26% (stream order 4 and 5). Conversely, the concentrations of organic C, Ca, Si, and trace metals increased with increasing stream order, suggesting mixing of authigenic material with suspended sediment from a different source. The Fe speciation in the suspended sediment was similar to that in fresh authigenic material, but more Fe-Fe interactions were observed, i.e. it was increasingly hydrolysed, suggesting ageing reactions. The suspended sediment in the streams of order 4 and 5 is estimated to contain between 31% and 59% of authigenic material, but more data are needed to refine this estimate. The authigenic material is an important sink for P in these streams which may alleviate the eutrophication risk in this catchment. (C) 2013 Elsevier Ltd. All rights reserved.

Abstract: Global water quality (WQ) modeling is an emerging field. In this article, we identify the missing linkages between global and basin/local-scale WQ models, and discuss the possibilities to fill these gaps. We argue that WQ models need stronger linkages across spatial scales. This would help to identify effective scale-specific WQ management options and contribute to future development of global WQ models. Two directions are proposed to improve the linkages: nested multiscale WQ modeling towards enhanced water management, and development of next-generation global WQ models based-on basin/local-scale mechanistic understanding. We highlight the need for better collaboration among WQ modelers and policy-makers in order to deliver responsive water policies and management strategies across scales.

Abstract: Urban environments provide opportunities for greater resource efficiency and the fostering of urban ecosystems. Brownfield areas are a typical example of underused land resources. Brownfield redevelopment projects that include green infrastructure allow for further ecosystems to be accommodated in urban environments. Green infrastructure also deliver important urban ecosystem services (UES) to local residents, which can greatly contribute to improving quality of life in cities. In this case study, we quantify and assess the economic value of five UES for a brownfield redevelopment project in Antwerp, Belgium. The assessment is carried out using the “Nature Value Explorer” modelling tool. The case includes three types of green infrastructure (green corridor, infiltration gullies and green roofs) primarily intended to connect nature reserves on the urban periphery and to avoid surface runoff. The green infrastructure also provides air filtration, climate regulation, carbon sequestration and recreation ecosystem services. The value of recreation far exceeds other values, including the value of avoided runoff. The case study raises crucial questions as to whether existing UES valuation approaches adequately account for the range of UES provided and whether such approaches can be improved to achieve more accurate and reliable value estimates in future analyses.

Abstract: A quasi 3D modelling approach was developed by integrating a crop growth (LINGRA-N) and a hydrological model (Hydrus-1D) to simulate and visualize water flow, soil-water storage, water stress and crop yield over a heterogeneous sandy field. We assessed computational efficiency and uncertainty with low to high-spatial resolution input factors (soil-hydraulic properties, soil-layer thickness and groundwater level) and evaluated four irrigation scenarios (no, current, optimized and triggered) to find the optimal and cost-effective irrigation scheduling. Numerical results showed that the simulation uncertainty was reduced when using the high-resolution information while a fast performance was maintained. The approach accurately determined the field scale irrigation requirements, taking into account spatial variations of input information. Optimal irrigation scheduling is obtained by triggered-irrigation resulting in saving up to similar to 300% water as compared to the current-irrigation, while yield increased similar to 1%. Overall, the approach can be useful to help decision makers and applicants in precision farming. (C) 2017 Published by Elsevier Ltd.

Abstract: In environmental assessments, the characterization of urban soils relies heavily on invasive investigation, which is often insufficient to capture their full spatial heterogeneity. Non-invasive geophysical techniques enable rapid collection of high-resolution data and provide a cost-effective alternative to investigate soil in a spatially comprehensive way. This paper presents the results of combining multi-receiver electromagnetic induction and stepped-frequency ground penetrating radar to characterize a former garage site contaminated with petroleum hydrocarbons. The sensor combination showed the ability to identify and accurately locate building remains and a high-density soil layer, thus demonstrating the high potential to investigate anthropogenic disturbances of physical nature. In addition, a correspondence was found between an area of lower electrical conductivity and elevated concentrations of petroleum hydrocarbons, suggesting the potential to detect specific chemical disturbances. We conclude that the sensor combination provides valuable information for preliminary assessment of urban soils.

Abstract: The fate of iron (Fe) may affect that of phosphorus (P) and arsenic (As) in natural waters. This study addresses the removal of Fe, P, and As from streams in lowland catchments fed by reduced, Fe-rich groundwater (average: 20 mg Fe L-1). The concentrations of dissolved Fe (<0.45 mu m) in streams gradually decrease with increasing hydraulic residence time (travel time) of the water in the catchment. The removal of Fe from streamwater is governed by chemical reactions and hydrological processes: the oxidation of ferrous iron (Fe(II)) and the subsequent formation of particulate Fe oxyhydroxides proceeds as the water flows through the catchment into increasingly larger streams. The Fe removal exhibits first-order kinetics with a mean half-life of 12 h, a value in line with predictions by a kinetic model for Fe(II) oxidation. The Fe concentrations in streams vary seasonally: they are higher in winter than in summer, due to shorter hydraulic residence time and lower temperature in winter. The removal of P and As is much faster than that of Fe. The average concentrations of P and As in streams (42 mu g P L-1) and 1.4 mu g As L-1) are 1 order of magnitude below those in groundwater (393 mu g P L-1 and 17 mu g As L-1). This removal is attributed to fast sequestration by oxidizing Fe when the water enters oxic environments, possibly by adsorption on Fe oxyhydroxides or by formation of ferric phosphates. The average P and As concentrations in groundwater largely exceed local environmental limits for freshwater (140 mu g P L-1 and 3 mu g As L((-1)), but in streams, they are below these limits. Naturally occurring Fe in groundwater may alleviate the environmental risk associated with P and As in the receiving streams.

Abstract: The soil at industrial sites is frequently characterized by very heterogeneous properties, which are often related to physical disturbance and contamination. A conventional approach to characterize the soil, with only a limited number of invasive observations, fails to capture the full extent of soil heterogeneity. Proximal soil sensing provides efficient tools to record spatially dense soil information. Nevertheless, because the output of most sensors is affected by more than one soil property, the simultaneous characterization of different soil properties requires the use of multiple sensors. Here, we apply multi-receiver electromagnetic induction (EMI) and stepped frequency ground penetrating radar (GPR) to survey a former gasworks site in a seaport area of Belgium. We used the EMI and GPR sensors in a motorized system to obtain densely sampled measurements of apparent electrical conductivity, apparent magnetic susceptibility and contrasts in relative dielectric permittivity. Our study shows that the sensors give detailed information on the variation in these electromagnetic soil properties. Interpretation of the variation in terms of the stratification of the soil was hampered by localized anthropogenic disturbances. However, the sensors provided complementary information that enabled the identification, discrimination and accurate location of several of these localized disturbances, including underground utility services such as electric cables, buried structures such as the remains of foundations and contamination by salts. Because these represent typical targets in industrial site investigation, we conclude that multi-receiver EMI and stepped frequency GPR provide a useful set of tools to expedite the investigation of industrial sites.

Abstract: Knowledge on groundwater-surface water interaction and especially on exchange fluxes between streams and aquifers is an important prerequisite for the study of transport and fate of contaminants and nutrients in the hyporheic zone. One possibility to quantify groundwater-surface water exchange fluxes is by using heat as an environmlental tracer. Modern field equipment including multilevel temperature sticks and the novel open-source analysis tool LPML make this technique ever more attractive. The recently developed LPML method solves the one-dimensional fluid flow and heat transport equation by combining a local polynomial method with a maximum likelihood estimator. In this study, we apply the LPML method on field data to quantify the spatial and temporal variability of vertical fluxes and their uncertainties from temperature-time series measured in a Belgian lowland stream. Over several months, temperature data were collected with multilevel temperature sticks at the streambed top and at six depths for a small stream section. Long-term estimates show a range from gaining fluxes of -291 mm day(-1) to loosing fluxes of 12 mm day(-1); average seasonal fluxes ranged from -138 mm day(-1) in winter to -16 mm day(-1) in summer. With our analyses, we could determine a high spatial and temporal variability of vertical exchange fluxes for the investigated stream section. Such spatial and temporal variability should be taken into account in biogeochemical cycling of carbon, nutrients and metals and in fate analysis of contaminant plumes. In general, the stream section was gaining during most of the observation period. Two short-term high stream stage events, seemingly caused by blockage of the stream outlet, led to a change in flow direction from gaining to losing conditions. We also found more discharge occurring at the outer stream bank than at the inner one indicating a local flow-through system. With the conducted analyses, we were able to advance our understanding of the regional groundwater flow system. Copyright (C) 2015 John Wiley & Sons, Ltd.

Abstract: Finding a correspondence between soil hydraulic properties, such as saturated hydraulic conductivity (Ks) and apparent electrical conductivity (ECa) as an easily measurable parameter, may be a way forward to estimate the spatial distribution of hydraulic properties at the field scale. In this study, the spatial distributions of Ks, of soil ECa measured by a DUALEM-21S sensor and of soil physical properties were investigated in a sandy grassland. To predict field scale Ks, the statistical relationship between co-located soil Ks, and EMI-ECa was evaluated. Results demonstrated the large spatial variability of all studied properties with Ks being the most variable one (CV = 86.21%) followed by ECa (CV >= 53.77%). A significant negative correlation was found between In-transformed Ks and ECa (r = 0.83; P <= 0.01) at two depths of exploration (0-50 and 0-100 cm). This site specific relation between In Ks and ECa was used to predict saturated hydraulic conductivity over 0-50 cm depth for the whole field. The empirical relation was validated using an independent dataset of measured Ks. The statistical results demonstrate the robustness of this empirical relation with mean estimation error MEE = 0.46 (cm h(-1)), root-mean-square estimation errors RMSEE = 0.74 (cm h(-1)), coefficient of determination r(2) = 0.67 and coefficient of model efficiency Ce = 0.64. The relationship was then used to produce a detailed map of Ks for the whole field. The result will allow model predictions of spatially distributed water content in view of irrigation management. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: Organic carbon introduction in the soil to initiate remedial measures, nitrate infiltration due to agricultural practices or sulphate intrusion owing to industrial usage can influence the redox conditions and pH, thus affecting the mobility of heavy metals in soil and groundwater. This study reports the fate of Zn and Cd in sandy aquifers under a variety of plausible in-situ redox conditions that were induced by introduction of carbon and various electron acceptors in column experiments. Up to 100% Zn and Cd removal (from the liquid phase) was observed in all the four columns, however the mechanisms were different. Metal removal in column K1 (containing sulphate), was attributed to biological sulphate reduction and subsequent metal precipitation (as sulphides). In the presence of both nitrate and sulphate (K2), the former dominated the process, precipitating the heavy metals as hydroxides and/or carbonates. In the presence of sulphate, nitrate and supplemental iron (Fe(OH)(3)) (K3), metal removal was also due to precipitation as hydroxides and/or carbonates. In abiotic column, K4, (with supplemental iron (Fe(OH)(3)), but no nitrate), cation exchange with soil led to metal removal. The results obtained were modeled using the reactive transport model PHREEQC-2 to elucidate governing processes and to evaluate scenarios of organic carbon, sulphate and nitrate inputs. (C) 2010 Elsevier B.V. All rights reserved.

Abstract: Stimulated anaerobic dechlorination is generally considered a valuable step for the remediation of aquifers polluted with chlorinated ethenes (CEs). Correct simulation and prediction of this process in situ, however, require good knowledge of the associated biological reactions. The aim of this study was to evaluate the dechlorination reaction in an aquifer contaminated with trichloroethene (TCE) and its daughter products, discharging into the Zenne River. Different carbon sources were used in batch cultures and these were related to the dechlorination reaction, together with the monitored biomarkers. Appropriate kinetic formulations were assessed. Reductive dechlorination of TCE took place only when external carbon sources were added to microcosms, and occurred concomitant with a pronounced increase in the Dehalococcoides mccartyi cell count as determined by 16S rRNA gene-targeted qPCR. This indicates that native dechlorinating bacteria are present in the aquifer of the Zenne site and that the oligotrophic nature of the aquifer prevents a complete degradation to ethene. The type of carbon source, the cell number of D. mccartyi or the reductive dehalogenase genes, however, did not unequivocally explain the observed differences in degradation rates or the extent of dechlorination. Neither first-order, Michaelis-Menten nor Monod kinetics could perfectly simulate the dechlorination reactions in TCE spiked microcosms. A sensitivity analysis indicated that the inclusion of donor limitation would not significantly enhance the simulations without a clear process understanding. Results point to the role of the supporting microbial community but it remains to be verified how the complexity of the microbial (inter)actions should be represented in a model framework. (C) 2013 Elsevier B.V. All rights reserved.

Abstract: Identifying priority areas is an essential step in developing management strategies to reduce pesticide loads in surface water. A spatially explicit model-based approach was developed to detect priority areas for diffuse pesticide pollution at catchment scale. The method uses available datasets and considers different pesticide pathways in the environment post-application. The approach was applied in a catchment area in SE Flanders (Belgium) as a case study. Calculated risk areas were obtained using detailed landscape data and combining pesticide emissions and hydrological connectivity. The risk areas obtained were further compared with an alternative observation-based method, developed specifically for this study site that includes long-term field observations and local expert knowledge. Both methods equally classified 50% of the areas. The impact of crop rotation on the calculated risk was analysed. High-risk areas were identified and added to a cumulative map over all five years to evaluate temporal variations. The model-based approach was used for the initial identification of risk areas at the study site. The tool helps to prioritise zones and detect particular fields to target landscape mitigation measures to reduce diffuse pesticide pollution reaching surface water bodies.

Abstract: In this study, the aging behavior of microscale zerovalent iron (mZVI) particles was investigated by quantifying the hydrogen gas generated by anaerobic mZVI corrosion in batch degradation experiments. Granular iron and nanoscale zerovalent iron (nZVI) particles were included in this study as controls. Firstly, experiments in liquid medium (without aquifer material) were performed and revealed that mZV1 particles have approximately a 10-30 times lower corrosion rate than nZVI particles. A good correlation was found between surface area normalized corrosion rate (R-SA) and reaction rate constants (K-SA) of PCE, TCE, cDCE and 1,1,1-TCA. Generally, particles with higher degradation rates also have faster corrosion rates, but exceptions do exists. In a second phase, the hydrogen evolution was also monitored during batch tests in the presence of aquifer material and real groundwater. A 4-9 times higher corrosion rate of mZV1 particles was observed under the natural environment in comparison with the aquifer free artificial condition, which can be attributed to the low pH of the aquifer and its buffer capacity. A corrosion model was calibrated on the batch experiments to take into account the inhibitory effects of the corrosion products (dissolved iron, hydrogen and OH-) on the iron corrosion rate. (C) 2014 Elsevier B.V. All rights reserved.

Abstract: When applying hydrological models, different sources of uncertainty are present, and evaluations of model performances should take these into account to assess model outcomes correctly. Furthermore, uncertainty in the discharge observations complicates the model identification, both in terms of model structure and parameterization. In this paper, the authors compare two different lumped model structures (PDM and NAM) considering uncertainty coming from the rating curve. Limits of acceptability for the model simulations were determined based on derived uncertainty bounds of the discharge observations. The authors applied the DYNamic Identifiability Approach (DYNIA) to identify structural failure of both models and to evaluate the configuration of their structures. In general, similar model performances are observed. However, the model structures tend to behave differently in the course of time, as revealed by the DYNIA approach. Based on the analyses performed, the probability based soil storage representation of the PDM model outperforms the NAM structure. The incorporation of the observation error did not prevent the DYNIA analysis to identify potential model structural deficiencies that are limiting the representation of the seasonal variation, primarily indicated by shifting regions of parameter identifiability. As such, the proposed approach is able to indicate where deficiencies are found and model improvement is needed.

Abstract: Field water flow processes can be precisely delineated with proper sets of soil hydraulic properties derived from in situ and/or laboratory experiments. In this study we analyzed and compared soil hydraulic properties obtained by traditional laboratory experiments and inverse optimization tension infiltrometer data along the vertical direction within two typical Podzol profiles with sand texture in a potato field. The main goal was to identify proper sets of hydraulic parameters and to evaluate their relevance on hydrological model performance for irrigation management purposes. Tension disc infiltration experiments were carried out at four and five different depths for both profiles at consecutive negative pressure heads of 12, 6, 3 and 0.1 cm. At the same locations and depths undisturbed samples were taken to determine Mualem-van Genuchten (MVG) hydraulic parameters (theta(r), residual water content, theta(s), saturated water content, alpha and n, shape parameters and K-ls, saturated hydraulic conductivity) in the laboratory. Results demonstrated horizontal differences and vertical variability of hydraulic properties. The tension disc infiltration data fitted well in inverse modeling using Hydrus 2D/3D in combination with final water content at the end of the experiment, theta(f). Four MVG parameters (theta(s), alpha, n and field saturated hydraulic conductivity K-fs) were estimated (theta(r) set to zero), with estimated K-ls and alpha values being relatively similar to values from Wooding's solution which used as initial value and estimated theta(s) corresponded to (effective) field saturated water content, theta(f). The laboratory measurement of K-ls yielded 2-30 times higher values than the field method K-fs from top to subsoil layers, while there was a significant correlation between both K-s values (r = 0.75). We found significant differences of MVG parameters theta(s), n and alpha values between laboratory and field measurements, but again a significant correlation was observed between laboratory and field MVG parameters namely K-s, n, theta(s) (r >= 0.59). Assessment of the parameter relevance in 1-D model simulations, illustrated that the model over predicted and under predicted top soil-water content using laboratory and field experiments data sets respectively. The field MVG parameter data set resulted in better agreement to observed soil-water content as compared to the laboratory data set at nodes 10 and 20 cm. However, better simulation results were achieved using the laboratory data set at 30-60 cm depths. Results of our study do not confirm whether laboratory or field experiments data sets are most appropriate to predict soil water fluctuations in a complete soil profile, while field experiments are preferred in many studies. (C) 2016 Elsevier B.V. All rights reserved.